实验三 图像的空间域滤波.docVIP

1 A=imread(E:\pic\1.jpg); I=rgb2gray(A); subplot(1,3,1); imshow(I);title(原图); J=imnoise(I,salt pepper,0.05); subplot(1,3,2); imshow(J); title(加入椒盐噪声图象); K=imnoise(I,gaussian,0.01,0.02); subplot(1,3,3); imshow(K);title(加入高斯噪声图象); 2 A=imread(E:\pic\1.jpg); I=rgb2gray(A); Subplot(2,2,1); Imshow(I);title(原图); H=fspecial(motion,20,45); MotionBlur=imfilter(I,H,replicate); Subplot(2,2,2); Imshow(MotionBlur);title(MotionBlur image); H=fspecial(disk,10); blurred=imfilter(I,H,replicate); Subplot(2,2,3); Imshow(blurred);title(Blurred image); H=fspecial(unsharp,0.5); Sharpened=imfilter(I,H,replicate); Subplot(2,2,4); Imshow(Sharpened);title(sharpened image); 3 A=imread(E:\pic\1.jpg); I=rgb2gray(A); J=imnoise(I,salt pepper,0.05); Subplot(2,2,1);imshow(J);title(加入椒盐噪声图象); H=fspecial(motion,20,45); MotionBlur=imfilter(J,H,replicate); Subplot(2,2,2); Imshow(MotionBlur);title(replicate); MotionBlur=imfilter(J,H,symmetric); Subplot(2,2,3); Imshow(MotionBlur);title(symmetric); MotionBlur=imfilter(J,H,circular); Subplot(2,2,4); Imshow(MotionBlur);title(circular); 4 A=imread(E:\pic\1.jpg); I=rgb2gray(A); J=imnoise(I,salt pepper,0.05); Aver=J; H=fspecial(average); for i =1:10 Aver=imfilter(Aver,H,replicate); end Subplot(1,3,1);imshow(J);title(加有椒盐噪声的图象); Subplot(1,3,2);imshow(Aver);title(10次均值滤波图象¨); for i = 1:10 Aver=imfilter(Aver,H,replicate); end Subplot(1,3,3);imshow(Aver);title(20次均值滤波图象); 5 A=imread(E:\pic\1.jpg); I=rgb2gray(A); J=imnoise(I,salt pepper,0.05); Subplot(1,3,1);imshow(J);title(加有椒盐噪声的图象); H=fspecial(average); Aver=imfilter(J,H,replicate); Subplot(1,3,2);imshow(Aver);title(均值滤波后的图象); Med=medfilt2(J); Subplot(1,3,3);imshow(Med);title(中值滤波后的图象); 6 A=imread(E:\pic\1.jpg); f=rgb2gray(A); Subplot(1,2,1);imshow(f);title(原图); w=(1/16)*[1 2 1;2 4 2;1 2 1]; f=im2double(f); g=imfilter(f,w,replicate); Subplot(1,2,2);imshow(g);title(平滑滤波后的图象); 1 A=imread(E:\pic\1.jpg); f=rgb2gray(A); w8=[1 1 1;1 -8 1;1 1 1]; f=im2double(f); g8=f-imfilter(f,w8,replicate); fig