实验四 图像的空域平滑滤波.doc

实验四 图像的空域平滑滤波 实验目的 理解图像空域平滑滤波的基本定义及目的； 掌握图像空域平滑滤波的基本原理及方法； 掌握用MATLAB语言实现图像的空域平滑滤波的方法。 实验内容 程序 I=imread(girl.bmp); figure; subplot(1,2,1); imshow(I); title(原图); J=imnoise(I,gaussian); subplot(1,2,2); imshow(J); title(加高斯噪声); 图像： 程序 I=imread(girl.bmp); title(原图); J=imnoise(I,gaussian); h1=fspecial(average,3); h2=fspecial(average,4); h3=fspecial(average,5); h4=fspecial(average,6); I1=filter2(h1,J); I2=filter2(h2,J); I3=filter2(h3,J); I4=filter2(h4,J); subplot(2,3,1), imshow(I), title(原图像); subplot(2,3,2), imshow(J), title(加噪声图像); subplot(2,3,3), imshow(I1, [ ]), title(3*3); subplot(2,3,4), imshow(I2, [ ]), title(4*4); subplot(2,3,5), imshow(I3, [ ]), title(5*5); subplot(2,3,6), imshow(I4, [ ]), title(6*6); 图像 程序： I=imread(girl.bmp); J=imnoise(I,gaussian); I1=medfilt2(J,[3 3]); I2=medfilt2(J,[4 4]); I3=medfilt2(J,[5 5]); I4=medfilt2(J,[6 6]); subplot(2,3,1), imshow(I), title(原图像); subplot(2,3,2), imshow(J), title(加噪声图像); subplot(2,3,3), imshow(I1), title(3*3); subplot(2,3,4), imshow(I2), title(4*4); subplot(2,3,5), imshow(I3), title(5*5); subplot(2,3,6), imshow(I4), title(6*6); 图像： 加高斯噪声程序： I=imread(girl.bmp); J=imnoise(I,gaussian); h=fspecial(average,3); I1=filter2(h,J); I2=medfilt2(J,[3 3]); subplot(2,2,1), imshow(I), title(原图像); subplot(2,2,2), imshow(J), title(加高斯噪声图像); subplot(2,2,3), imshow(I1,[ ]), title(均值滤波); subplot(2,2,4), imshow(I2), title(中值滤波); 图像： 加椒盐噪声程序： I=imread(girl.bmp); J=imnoise(I,salt pepper); h=fspecial(average,3); I1=filter2(h,J); I2=medfilt2(J,[3 3]); subplot(2,2,1), imshow(I), title(原图像); subplot(2,2,2), imshow(J), title(加椒盐噪声图像); subplot(2,2,3), imshow(I1,[ ]), title(均值滤波); subplot(2,2,4), imshow(I2), title(中值滤波); 图像： 均值滤波： I=imread(girl.bmp); J=imnoise(I,gaussian); J=double(J); for a=2:127 for b=2:127 I1(a,b)=1/9*(J(a-1,b-1)+J(a-1,b)+J(a-1,b+1)+J(a,b-1)+J(a,b)+J(a,b+1)+J(a+1,b-1)+J(a,b)+J(a,b+1)); end end subplot(1,3,1), imshow(I), title(原图像); subplot(1,3,2), imshow(J,[ ]), title(加高斯噪声图像); subplot(1,3,3), imsh