图像处理MATLAB实验三

《图像处理MATLAB实验三》由会员分享，可在线阅读，更多相关《图像处理MATLAB实验三（10页珍藏版）》请在金锄头文库上搜索。

1、1 图像处理的图像处理的MATLAB实现实现骄任肝说以抓猛装示膝特碑艳季躺狰筑素迫擞戒蝉彤抽蓑泛摆玩何芍放纳图像处理MATLAB-实验三图像处理MATLAB-实验三2 实验二 频域图像增强一、实验目的一、实验目的 （1 1）掌掌握握基基本本的的频频域域图图像像增增强强方方法法，观观察察图图像像增增 强的效果，加深理解；强的效果，加深理解； （2 2）了解不同类型频域滤波器的特性；）了解不同类型频域滤波器的特性； （3 3）了解滤波器的截止频率对滤波效果的影响。）了解滤波器的截止频率对滤波效果的影响。 圈残锨煞晒獭暖棱且瑞恨月果姬汐菌杠会仓缀郭社枣斧锄宗梧符典喉痴乱图像处理MATLAB-实验三图

2、像处理MATLAB-实验三3 实验二 频域图像增强二、实验内容二、实验内容 （1 1）频率域平滑）频率域平滑 （2 2）频率域锐化）频率域锐化三、实验要求三、实验要求 （1 1）用）用matlabmatlab语言进行仿真实验；语言进行仿真实验； （2 2）递交实验报告，要求给出实验原理、源程序、递交实验报告，要求给出实验原理、源程序、 实验结果及分析。实验结果及分析。 予恤忆取爬染氨滦痛钻淬芬驴怕急羔矛噎湛孤烬慈柴整鸿豺欺啦呛猎豆搪图像处理MATLAB-实验三图像处理MATLAB-实验三4 实验二 频域图像增强四、实验图像四、实验图像 test.tif调氧伎迫周沫情沦控丛蠢瘩踢妻鱼阀昏划饱渗忿

3、曾春琼蓄蹲动轨宵份西墓图像处理MATLAB-实验三图像处理MATLAB-实验三5 实验二 频域图像增强 五、具体实验内容及要求五、具体实验内容及要求 1. 1. 频域平滑实验内容频域平滑实验内容 （1 1）读入原图像）读入原图像test.tiftest.tif并显示；并显示； （2 2）采用理想低通滤波器对图像进行处理，并显示）采用理想低通滤波器对图像进行处理，并显示 结果；结果； （3 3）采用二阶巴特沃斯低通滤波器对图像进行处理，）采用二阶巴特沃斯低通滤波器对图像进行处理， 并显示结果。并显示结果。 潦绎镁评瞻舅莱褪崇穴炮袁奈溅隐逗递冕黎较侯绳想抠圭啤彩田佩件盅铭图像处理MATLAB-实验

4、三图像处理MATLAB-实验三6 实验二 频域图像增强 要求：要求： （1 1）观察截止频率分别为观察截止频率分别为15、30、80时对滤波结果的时对滤波结果的 影响。影响。 （2 2）比较理想低通滤波器和巴特沃斯低通滤波器的滤比较理想低通滤波器和巴特沃斯低通滤波器的滤 波效果。波效果。砖奸胎翱缆袖嘱藉虫蝎涟羔辙碘晾腻解般纶硬头卿叁役雷佣堂念虞薛钥萤图像处理MATLAB-实验三图像处理MATLAB-实验三7频率域平滑低通滤波法例：例：利用二阶巴特沃斯低通滤波器和理想低通滤波利用二阶巴特沃斯低通滤波器和理想低通滤波器分别对灰度图像器分别对灰度图像test.tif test.tif 进行平滑处理。

5、进行平滑处理。 分析：分析： 理想低通滤波器的传递函数为理想低通滤波器的传递函数为 n阶巴特沃斯低通滤波器的传递函数为阶巴特沃斯低通滤波器的传递函数为雪炕咱合眼坎俏缠弦或良秉侗蓄勿黄拘佑蓉剃记卡卓呵砌号雹室睬存惺汀图像处理MATLAB-实验三图像处理MATLAB-实验三8频率域平滑低通滤波法 I=imread(test.tif); I=imread(test.tif); subplot(1,3,1),imshow(I) subplot(1,3,1),imshow(I) title( title(原始图像原始图像) f=double(I); f=double(I); g=fft2(f); g=f

6、ft2(f); % fft2% fft2函数用于计算图像的快函数用于计算图像的快速傅立叶变换速傅立叶变换 g=fftshift(g); g=fftshift(g); % fftshift% fftshift可以用于调整可以用于调整fft2fft2的的输出结果，将输出结果，将 变换后的变换后的图像频谱中心从矩阵的原点移到图像频谱中心从矩阵的原点移到 矩阵的中矩阵的中心。心。 M,N=size(g); M,N=size(g); n=2; n=2; d0=30; d0=30; % % 截止频率设为截止频率设为3030 m0=fix(M/2); m0=fix(M/2); n0=fix(N/2); n0

7、=fix(N/2); 烬楔且汲醒牛赠纂郧吼酋玖崖麓谬煮四汀谋啸颧茄论全篮衅肤潦忍苟勃填图像处理MATLAB-实验三图像处理MATLAB-实验三9频率域平滑低通滤波法 for i=1:M for i=1:M for j=1:N for j=1:N d=sqrt(i-m0)2+(j-n0)2); d=sqrt(i-m0)2+(j-n0)2); h=1/(1+(d/d0)(2*n); h=1/(1+(d/d0)(2*n); % % 计计算二阶巴特沃斯低通算二阶巴特沃斯低通 滤波器的传递函数滤波器的传递函数 result1(i,j)=h*g(i,j); result1(i,j)=h*g(i,j); %

8、 % 进行巴特沃斯低通滤波进行巴特沃斯低通滤波 if d30 if d30 % % 进行理想低通滤波进行理想低通滤波 result2(i,j)=0; result2(i,j)=0; else else result2(i,j)=g(i,j); result2(i,j)=g(i,j); end end end end end end 胚居羹肩谴剐捡拥吐洗蹲穗旨喀兆淄界阻骇匀支宵慎寸恰勉佳墅终泅枪息图像处理MATLAB-实验三图像处理MATLAB-实验三10频率域平滑低通滤波法 result1=ifftshift(result1); result1=ifftshift(result1); J1=i

9、fft2(result1); J1=ifft2(result1); X1=uint8(real(J1); X1=uint8(real(J1); subplot(1,3,2),imshow(X1) subplot(1,3,2),imshow(X1) title( title(二阶巴特沃斯低通滤波图像二阶巴特沃斯低通滤波图像) result2=ifftshift(result2); result2=ifftshift(result2); J2=ifft2(result2); J2=ifft2(result2); X2=uint8(real(J2); X2=uint8(real(J2); subplot(1,3,3),imshow(X2) subplot(1,3,3),imshow(X2) title( title(理想低通滤波图像理想低通滤波图像) 圭吨缨锗撩躬箱叮遁蛮郴凤尝雁谩砚洪汪握友岔葫赫犯勒旺像怔瓶民狗肥图像处理MATLAB-实验三图像处理MATLAB-实验三