图像处理MATLAB 实验三.pptVIP

* 图像处理的MATLAB实现 * 实验二 频域图像增强 一、实验目的 （1）掌握基本的频域图像增强方法，观察图像增 强的效果，加深理解； （2）了解不同类型频域滤波器的特性； （3）了解滤波器的截止频率对滤波效果的影响。 * 实验二 频域图像增强 二、实验内容 （1）频率域平滑 （2）频率域锐化 三、实验要求 （1）用matlab语言进行仿真实验； （2）递交实验报告，要求给出实验原理、源程序、 实验结果及分析。 * 实验二 频域图像增强 四、实验图像 test.tif * 实验二 频域图像增强 五、具体实验内容及要求 1. 频域平滑实验内容 （1）读入原图像test.tif并显示； （2）采用理想低通滤波器对图像进行处理，并显示 结果； （3）采用二阶巴特沃斯低通滤波器对图像进行处理， 并显示结果。 * 实验二 频域图像增强 要求： （1）观察截止频率分别为15、30、80时对滤波结果的 影响。 （2）比较理想低通滤波器和巴特沃斯低通滤波器的滤 波效果。 * 频率域平滑—低通滤波法 例：利用二阶巴特沃斯低通滤波器和理想低通滤波 器分别对灰度图像test.tif 进行平滑处理。 分析： 理想低通滤波器的传递函数为 n阶巴特沃斯低通滤波器的传递函数为 * 频率域平滑—低通滤波法 I=imread(test.tif); subplot(1,3,1),imshow(I) title(原始图像) f=double(I); g=fft2(f); % fft2函数用于计算图像的快速傅立叶变换 g=fftshift(g); % fftshift可以用于调整fft2的输出结果，将 变换后的图像频谱中心从矩阵的原点移到 矩阵的中心。 [M,N]=size(g); n=2; d0=30; % 截止频率设为30 m0=fix(M/2); n0=fix(N/2); * 频率域平滑—低通滤波法 for i=1:M for j=1:N d=sqrt((i-m0)^2+(j-n0)^2); h=1/(1+(d/d0)^(2*n)); % 计算二阶巴特沃斯低通 滤波器的传递函数 result1(i,j)=h*g(i,j); % 进行巴特沃斯低通滤波 if d30 % 进行理想低通滤波 result2(i,j)=0; else result2(i,j)=g(i,j); end end end * 频率域平滑—低通滤波法 result1=ifftshift(result1); J1=ifft2(result1); X1=uint8(real(J1)); subplot(1,3,2),imshow(X1) title(二阶巴特沃斯低通滤波图像) result2=ifftshift(result2); J2=ifft2(result2); X2=uint8(real(J2)); subplot(1,3,3),imshow(X2) title(理想低通滤波图像)