光信数字图像处理图像平滑实验l.docVIP

实验2 图像平滑实验 一、实验目的 1．通过实验掌握图像去噪的基本方法； 2．学会根据情况选用不同方法。 二、实验的硬件、软件平台 硬件：计算机 软件：操作系统：WINDOWS 7 应用软件：MATLAB 三、实验内容及要求 1．实验内容 请在如下面方法中选择多个，完成图像去噪操作，并进行分析、比较。 （1）对静态场景的多幅图片取平均； （2）空间域模板卷积（不同模板、不同尺寸）； （3）频域低通滤波器（不同滤波器模型、不同截止频率）； （4）中值滤波方法。 2．实验要求 （1）图片可根据需要选取； （2）对不同方法和同一方法的不同参数的实验结果进行分析和比较，如空间域卷积模板可有高斯型模板、矩形模板、三角形模板和自己根据需求设计的模板等；模板大小可以是3×3，5×5，7×7或更大。频域滤波可采用矩形或巴特沃斯等低通滤波器模型，截止频率也是可选的。 （3）分析比较不同方法的结果。 四、思考题 1．不同空间域卷积器模板的滤波效果有何不同？ 2．空间域卷积器模板的大小的滤波效果有何影响？ 3．用多幅图像代数平均的方法去噪对图像有何要求？ 4．不同频域滤波器的效果有何不同？ 五、实验报告要求 1．列出程序清单并进行功能注释； 2．说明不同方法去噪效果； 3．对去噪方法进行详细分析对比。 ******************对静态场景的多幅图片取平均************************ 实验代码如下 clc clear all a1=imread(lenna_noise1.bmp);%读入图片 a2=imread(lenna_noise2.bmp); a3=imread(lenna_noise3.bmp); a4=imread(lenna_noise4.bmp); a5=imread(lenna_noise5.bmp); a6=imread(lenna_noise6.bmp); a7=imread(lenna_noise7.bmp); a8=imread(lenna_noise8.bmp); b1=double(a1)/255;%变换图像数据类型 b2=double(a2)/255; b3=double(a3)/255; b4=double(a4)/255; b5=double(a5)/255; b6=double(a6)/255; b7=double(a7)/255; b8=double(a8)/255; c=(b1+b2+b3+b4+b5+b6+b7+b8)/8;%八幅图取平均 subplot(331),imshow(a1);%三行三列，显示图像1-8 subplot(332),imshow(a2); subplot(333),imshow(a3); subplot(334),imshow(a4); subplot(335),imshow(a5); subplot(336),imshow(a6); subplot(337),imshow(a7); subplot(338),imshow(a8); subplot(339),imshow(c) %%显示平均后的图像 ******************空间域模板卷积（不同模板、不同尺寸）*************** %%邻域模板 clc clear all i=imread(Panda_jiaoyan.bmp); %读入图像 i=double(i)/255;%变换数据类型 subplot(221); imshow(i);%显示图像 title(原图像); h=1/5*[0 1 0;1 0 1;0 1 0];%定义4邻域平均模板 a=filter2(h,i); %进行滤波 subplot(222); imshow(a);%显示图像 title(4邻域平均模板); subplot(223); h=1/12*[0 1 1 0;1 1 1 1;1 1 1 1;0 1 1 0]; %定义8邻域平均模板 a=filter2(h,i); imshow(a); title(8邻域平均模板); subplot(224); h=1/12*[0 1 1 1 0;1 1 1 1 1;1 1 1 1 1;0 1 1 1 0]; %定义12邻域平均模板 a=filter2(h,i); imshow(a); title(12邻域平均模板) %%高斯模板 clc clear all %高斯模板 i=imread(Panda_jiaoyan.bmp);%读入图像 i=double(i)/255; %数值转换 subplot(231); imshow(i); %显示图像 %选用3*3的高斯模板 title(原图像(高斯模板)) h=fspeci