实验七 图像增强—平滑和滤波.docxVIP

实验七 图像增强—平滑与滤波****************共5小题***************付蕊一、实验目的及要求在熟悉图像平滑的基本原理和方法的基础上，在理论指导下，能在MATLAB环境下对图像进行平滑处理。本实验要求用线性平滑滤波、中值平滑滤波、小波降噪和频域低通滤波的方法进行程序设计。通过平滑处理，对结果图像加以比较，得出自己的实验结论。学生应当完成对于给定图像+噪声，使用平均滤波器、中值滤波器对不同强度的高斯噪声和椒盐噪声，进行滤波处理；能够正确地评价处理的结果；能够从理论上作出合理的解释。二、实验设备与软件1．PC计算机2．MATLAB软件/语言包括图像处理工具箱(Image Processing Toolbox) 3．实验所需要的图片三．实验内容1．对给定的图像添加均值为0，方差分别为200、400的高斯噪声，以及概率分别为0.1、0.2的椒盐噪声，显示添加噪声后的图像。加高斯 I=imread(rice.jpg); J=imnoise(I,gaussian,0,200); K=imnoise(I,gaussian,0,400); subplot(2,2,1);imshow(I);title(原图); subplot(2,2,2);imshow(J);title(200高斯噪声); subplot(2,2,3);imshow(I);title(原图); subplot(2,2,4);imshow(K);title(400高斯噪声);加椒盐 I=imread(rice.jpg); J=imnoise(I,salt pepper,0.2); K=imnoise(I,salt pepper,0.1); subplot(2,2,1);imshow(I);title(原图); subplot(2,2,2);imshow(J);title(0.2椒盐噪声); subplot(2,2,3);imshow(I);title(原图); subplot(2,2,4);imshow(K);title(0.1椒盐噪声);2．对加噪图像分别进行3×3、5×5、9×9的均值滤波，显示并比较滤波结果。 I=imread(rice.jpg); J=imnoise(I,salt pepper,0.2); ave1=fspecial(average,3); ave2=fspecial(average,5); ave3=fspecial(average,9); K = filter2(ave1,J)/255; L = filter2(ave2,J)/255; M = filter2(ave3,J)/255; subplot(2,2,1);imshow(J);title(椒盐噪声图); subplot(2,2,2);imshow(K);title(3*3均值滤波); subplot(2,2,3);imshow(L);title(5*5均值滤波); subplot(2,2,4);imshow(M);title(9*9均值滤波);3．对加噪图像分别进行3×3、5×5、9×9的中值滤波，显示并比较滤波结果。 K=medfilt2(J,[3 3]); L=medfilt2(J,[5 5]); M=medfilt2(J,[9 9]); subplot(2,2,1);imshow(J);title(椒盐噪声图); subplot(2,2,2);imshow(K);title(3*3中值滤波); subplot(2,2,3);imshow(L);title(5*5中值滤波); subplot(2,2,4);imshow(M);title(9*9中值滤波);4．对加噪图像进行截至频率D0分别为10、25的高斯低通滤波，显示并比较滤波结果。I=imread(rice.jpg);J=imnoise(I,salt pepper,0.2);f=double(J);ff=fft2(f);g=fftshift(ff);[n1,n2]=size(g);d0=25;u0=round(n1/2);v0=round(n2/2);for i=1:n1 for j=1:n2 d=sqrt((i-u0)^2+(j-v0)^2); h=exp(-d^2/(2*d0^2)); y(i,j)=g(i,j)*h; end;end;y1=ifftshift(y);e1=ifft2(y1);e2=uint8(real(e1));for i=1:n1 for j=1:n2 d=sqrt((i-u0)^2+(j-v0)^2); h=exp(-d^2/(2*d1^2)); y(i,j)=g(i,