理想低通滤波及高通滤波实现高频增强实验报告.pdf

实验作业 5： 1、 用理想低通滤波器在频率域实现低通滤波 程序代码： I=imread(d:/3.jpg); for j=1:b I=rgb2gray(I); figure(1),imshow(I); distance=sqrt((i-a0)^2+(j-b0)^2); title( 原图像 ); if distance=d h=1; s=fftshift(fft2(I)); else h=0; figure(2); end; imshow(abs(s),[]); s(i,j)=h*s(i,j); title( 图像傅里叶变换所得频谱 ); end; [a,b]=size(s); end; a0=round(a/2); s=uint8(real(ifft2(ifftshift(s)))); b0=round(b/2); figure(3); d=10; imshow(s); for i=1:a title( 低通滤波所得图像 ); 滤波 d=10 2、用理想高通滤波器在频率域实现高频 增强 I=imread(d:/3.jpg); for i=1:a I=rgb2gray(I); for j=1:b figure(1),imshow(I); title( 原图像 ); distance=sqrt((i-a0)^2+(j-b0)^2); s=fftshift(fft2(I)); if distance=d h=0; figure(2); else h=1; imshow(abs(s),[]); end; title( 图像傅里叶变换所得频谱 ); s(i,j)=(p+q*h)*s(i,j); figure(3); end; imshow(log(abs(s)),[]); end; title( 图像傅里叶变换取对数所得频谱 s=uint8(real(ifft2(ifftshift(s)))); ); figure(4); [a,b]=size(s); imshow(s);title( 高通滤波所得图像 ); a0=round(a/2); figure(5); b0=round(b/2); imshow(s+I);title( 高通 滤波所得高频 d=10; 增强图像 ); p=0.2;q=0.5; 心得体会： 这次按照低通滤波器和高