图像处理实验一实验二

1、实验一 空域图像增强一、实验目的（1）掌握基本的空域图像增强方法，观察图像增强的效果，加深理解；（2）了解空域平滑模板的特性及其对不同噪声的影响；（3）了解空域锐化模板的特性及其对边缘的影响。二、实验内容 （1）直方图处理：直方图均衡 （2）空域平滑：均值滤波、中值滤波；三、实验要求 （1）用matlab语言进行仿真实验； （2）递交实验报告，要求给出实验原理、源程序、实验结果及分析。四、实验图像 pollen.png lena.bmp bridge.gif五、具体实验内容及结果 1. 直方图均衡实验内容 （1）读入原图像pollen.png并显示原图像以及直方图 （2）对原图像进行直方图均衡

2、处理 （3）显示均衡后图像以及直方图。 要求：分析直方图能够改善视觉效果的原因实验代码：I=imread('D:实验图像pollen.png');subplot(4,5,1);imshow(I);title('原图');Gr=rgb2gray(I);subplot(4,5,2);imshow(Gr);title('gray');subplot(4,5,3);imhist(Gr);title('直方图');Grc=histeq(Gr);subplot(4,5,4);imshow(Grc);title('均衡化');s

3、ubplot(4,5,5);imhist(Grc);title('均衡化直方图');R=I(:,:,1);A=histeq(R);subplot(4,5,7);imshow(R);title('Red');subplot(4,5,8);imhist(R);title('直方图');subplot(4,5,9);imshow(A);title('均衡化');subplot(4,5,10);imhist(A);title('均衡化直方图');G=I(:,:,2);N=histeq(G);subplot(4,5,12);

4、imshow(G);title('Green');subplot(4,5,13);imhist(G);title('直方图');subplot(4,5,14);imshow(N);title('均衡化');subplot(4,5,15);imhist(N);title('均衡化直方图');B=I(:,:,3);C=histeq(B);subplot(4,5,17);imshow(B);title('Blue');subplot(4,5,18);imhist(B);title('直方图');subpl

5、ot(4,5,19);imshow(C);title('均衡化');subplot(4,5,20);imhist(C);title('均衡化直方图');结果：直方图均衡化能使原始图像的灰度直方图从比较集中的某个灰度区间变成在全部灰度范围内的均匀分布。故能够改善过暗或过亮的图像的视觉效果。2. 空域平滑实验内容 （1）读入原图像lena.bmp并显示； （2）对原图像分别添加高斯噪声和椒盐噪声，并显示加噪图像； （3）采用均值滤波进行去噪处理，并显示去噪图像； （4）采用中值滤波进行去噪处理，并显示去噪图像。 要求： （1）观察模板大小（窗口尺寸）对图像去噪效果的

6、影响（选用3×3、5×5、7×7）。 （2）观察均值滤波和中值滤波对不同噪声的滤除效果 实验代码：3*3窗口：I=imread('D:实验图像lena.bmp');G=imnoise(I,'gaussian',0,0.01);subplot(3,3,2);imshow(G);J=imnoise(I,'salt & pepper',0.05);subplot(3,3,3);imshow(J);h=(1/9)*ones(3);subplot(3,3,4);imshow(I);B1=conv2(G,h,'s

7、ame');subplot(3,3,5);imshow(B1,);B2=conv2(J,h,'same');subplot(3,3,6);imshow(B2,);subplot(3,3,7);imshow(I);C1=medfilt2(G,3 3);subplot(3,3,8);imshow(C1,);C2=medfilt2(J,3 3);subplot(3,3,9);imshow(C2,);5*5窗口：I=imread('D:实验图像lena.bmp');G=imnoise(I,'gaussian',0,0.01);subplot(3,

8、3,2);imshow(G);J=imnoise(I,'salt & pepper',0.05);subplot(3,3,3);imshow(J);h=(1/9)*ones(5);subplot(3,3,4);imshow(I);B1=conv2(G,h,'same');subplot(3,3,5);imshow(B1,);B2=conv2(J,h,'same');subplot(3,3,6);imshow(B2,);subplot(3,3,7);imshow(I);C1=medfilt2(G,5 5);subplot(3,3,8);im

9、show(C1,);C2=medfilt2(J,5 5);subplot(3,3,9);imshow(C2,);7*7窗口：I=imread('D:实验图像lena.bmp');G=imnoise(I,'gaussian',0,0.01);subplot(3,3,2);imshow(G);J=imnoise(I,'salt & pepper',0.05);subplot(3,3,3);imshow(J);h=(1/9)*ones(7);subplot(3,3,4);imshow(I);B1=conv2(G,h,'same'

10、);subplot(3,3,5);imshow(B1,);B2=conv2(J,h,'same');subplot(3,3,6);imshow(B2,);subplot(3,3,7);imshow(I);C1=medfilt2(G,7 7);subplot(3,3,8);imshow(C1,);C2=medfilt2(J,7 7);subplot(3,3,9);imshow(C2,);实验结果：3*3：5*5：7*7：通过比较以上三幅图可知中值滤波更适合滤除椒盐噪声均值滤波更适合滤除高斯噪声而窗口越大越能更好地滤除噪声但图像会越模糊。3. 空域锐化实验内容 （1）读入原图像br

11、idge.gif并显示； （2）采用sobel算子对图像进行处理，并显示结果； （3）尝试采用其他锐化模板进行处理。 实验代码：I=imread('D:实验图像brige.gif');subplot(1,3,1);imshow(I);H=fspecial('sobel');G=filter2(H,I);subplot(1,3,2);imshow(abs(G),);H1=fspecial('prewitt');J=filter2(H1,I);subplot(1,3,3);imshow(abs(J),);实验结果：实验二 频域图像增强一、实验目的 （

12、1）掌握基本的频域图像增强方法，观察图像增强的效果，加深理解； （2）了解不同类型频域滤波器的特性； （3）了解滤波器的截止频率对滤波效果的影响。 二、实验内容 （1）频率域平滑 （2）频率域锐化三、实验要求 （1）用matlab语言进行仿真实验； （2）递交实验报告，要求给出实验原理、源程序、实验结果及分析。四、实验图像 五、具体实验内容及结果 1. 频域平滑实验内容 （1）读入原图像test.tif并显示； （2）采用理想低通滤波器对图像进行处理，并显示结果； （3）采用二阶巴特沃斯低通滤波器对图像进行处理，并显示结果。 要求： （1）观察截止频率分别为15、30、80时对滤波结果的影响。

13、 （2）比较理想低通滤波器和巴特沃斯低通滤波器的滤波效果。 实验代码：主程序：理想低通：I=imread('D:实验图像test.tif');subplot(2,3,1);imshow(I,);G1=mylowpass(15);subplot(2,3,4);imshow(G1,);G2=mylowpass(30);subplot(2,3,5);imshow(G2,);G3=mylowpass(80);subplot(2,3,6);imshow(G3,);巴特沃斯：I=imread('D:实验图像test.tif');subplot(2,3,1);imshow(I

14、,);G1=mybutterw(15);subplot(2,3,4);imshow(G1);G2=mybutterw(30);subplot(2,3,5);imshow(G2);G3=mybutterw(80);subplot(2,3,6);imshow(G3);理想低通滤波器函数：function G=mylowpass(d)I=imread('D:实验图像test.tif');f=fft2(I);H=zeros(500);for i=1:500 for j=1:500 if(f(i,j)<=d) H(i,j)=1; end endendF=f.*H;G = ifft2(F);巴特沃斯滤波器函数：function G=mybutterw(d)I=imread('D:实验图像test.tif');f=fft2(I);H=ze