数字图像实验2解析.docVIP

实验2 图像平滑实验 一、实验目的 ?? ?1．通过实验掌握图像去噪的基本方法； ? ?2．学会根据情况选用不同方法。 ? 二、实验内容及要求 ? 1．实验内容 ?请在如下面方法中选择多个，完成图像去噪操作，并进行分析、比较。? （1）空间域模板卷积（不同模板、不同尺寸）； ? （2）频域低通滤波器（不同滤波器模型、不同截止频率）；? （3）中值滤波方法。? 2．实验要求 ? （1）图片可根据需要选取； ? （2）对不同方法和同一方法的不同参数的实验结果进行分析和比较，如空间域卷积模板可有矩形模板、三角形模板和自己根据需求设计的模板等；模板大小可以是3×3，5×5，7×7或更大。频域滤波可采用矩形或巴特沃斯等低通滤波器模型，截止频率也是可选的。 ? （3）分析比较不同方法的结果。 ?? 三、思考题 ? 1．不同空间域卷积器模板的滤波效果有何不同？? 2．空间域卷积器模板的大小的滤波效果有何影响？? 3．不同频域滤波器的效果有何不同？ ??四、实验报告要求 ? 1．列出程序清单并进行功能注释，递交.m文件； 2．说明不同方法去噪效果； ? 3．对去噪方法进行详细分析对比。 I=imread(Q:\4.tiff); subplot(3,3,1);%三行三列第一幅 imshow(I);%显示图片 title(原图) B=medfilt2(I,[3,3]);%中值滤波函数 subplot(3,3,2); imshow(B); title(中值滤波) C=ordfilt2(I,6,ones(3,3));%二维顺序统计滤波函数 subplot(3,3,3); imshow(C); title(二维顺序统计滤波) h=fspecial(average,[3,3]);%均值滤波器 D=filter2(h,I)/255;%滤波 subplot(3,3,4); imshow(D)%显示图片 title(均值滤波) h=fspecial(sobel);%边缘增强算子 D=filter2(h,I)/255;%滤波 subplot(3,3,5); imshow(D);%显示图片 title(sobel算子) h=fspecial(prewitt);%边缘增强算子 D=filter2(h,I)/255;%滤波 subplot(3,3,6); imshow(D)%显示图片 title(Prewitt算子) figure(2) Q=medfilt2(I,[3,3]);%中值率波函数 subplot(2,2,1); imshow(Q); title(模版大小为[3,3]的中值滤波函数) W=medfilt2(I,[5,5]);%中值率波函数 subplot(2,2,2); imshow(W); title(模版大小为[5,5]的中值滤波函数) E=medfilt2(I,[7,7]);%中值率波函数 subplot(2,2,3); imshow(E); title(模版大小为[7,7]的中值滤波函数) R=medfilt2(I,[9,9]);%中值率波函数 subplot(2,2,4); imshow(R); title(模版大小为[9,9]的中值滤波函数) 换图片只需改编号 (3)巴斯低通滤波 [I,map]=imread(Q:\8.jpg); [M,N]=size(I); F=fft2(I);%进行傅里叶变换 fftshift(F); %将灰度图像的二维不连续Frourier变换的零频率成分移到频谱的中心 Dcut=100; D0=30; %定义D0 for u=1:M for v=1:N D(u,v)=sqrt(u^2+v^2); BUTTERH(u,v)=1/(1+(sqrt(2)-1)*(D0/Dcut)^2); %巴特罗斯滤波器的函数 end end BUTTERG=BUTTERH.*F; BUTTERfiltered=ifft2(BUTTERG); %进行二维反离散的Fourier变换后，取复数的实部转化为无符号8位整数 subplot(121) imshow(I); title(源图) subplot(122) imshow(BUTTERfiltered,map) title(巴特罗斯低通滤波后的图像) 思考题 1．不同空间域卷积器模板的滤波效果有何不同？? 2．空间域卷积器模板的大小的滤波效果有何影响？?3．不同频域滤波器的效果有何不同？