1、 数 字 图 像 处 理 上机作业 教 师 ： 刘清华 学 院 ： 电子工程 学院 班 级 ： 021141 专 业 ： 探测制导与控制技术 姓 名 ： 苏二鹏 学 号 ： 02114030(1) (2) f1 与 f2 的幅度谱相同，因为 f1 与 f2 振幅的绝对值是相同的。 （ 3） f3 的幅度谱与 f2 的幅度谱按顺时针旋转 90度后相同 （ 4） f4 的幅度谱与 f1 的幅度谱按顺时针旋转 90度后 相同 （ 5） f5 的幅度谱是 f1 与 f4 的幅度谱相加， f6 的幅度谱是 f2 与 f3 的幅度谱相加 程序： f1=zeros(256,256); f1(65:192,1
2、13:144)=100; figure subplot(121) imshow(f1) title(原图 f1) F=fft2(f1); subplot(122) imshow(abs(F) title(f1幅度谱图 ) for m=(1:256); for n=(1:256); f2(m,n)=power(-1,m+n)*f1(m,n); end end figure subplot(121) imshow(f2) title(f2图 ) F1=fft2(f2); subplot(122) imshow(abs(F1) title(f2幅度谱图 ) f3=imrotate(f2,-90);
3、figure subplot(121) imshow(f3) title(f3图 ) F2=fft2(f3); subplot(122) imshow(abs(F2) title(f3幅度谱图 ) f4=imrotate(f1,-90); figure subplot(121) imshow(f4) title(f4图 ) F5=fft2(f4); subplot(122) imshow(abs(F5) title(f4幅度谱图 ) f5=f1+f4; figure subplot(121) imshow(f5) title(f5图 ) F3=fft2(f5) subplot(122) ims
4、how(abs(F3) title(f5幅度谱图 ) f6=f2+f3; figure subplot(121) imshow(f6) title(f6图 ) F4=fft2(f6) subplot(122) imshow(abs(F4) title(f6幅度谱图 ) 256*256图像如下： 从中值滤波和均值滤波后的图像可以看出，中值滤波把黑色方快完全完全隔离开，而白色方块则通过边角的两个多出来的像素点相连，而均值滤波则把图像变模糊了，边缘模糊的比较厉害。中值滤波后图像的直方图没有变化，而均值滤波后的直方图变化了。 程序： A=zeros(256,256); for i=(1:64:256)
5、; for j=(1:64:256); A(j:j+32,i:i+32)=1; end end for m=(32:64:256); for n=(32:64:256); A(n:n+32,m:m+32)=1; end end figure imshow(A) title(原图像 ) figure K=medfilt2(A,3,3) subplot(121) imshow(K) title(3*3中值滤波 ) M=filter2(fspecial(average,3),A) subplot(122) imshow(M) title(3*3均值值滤波 ) figure subplot(131)
6、imhist(A) axis(0 2 0 2) title(原图直方图 ) subplot(132) imhist(K) axis(0 2 0 2) title(中值滤波后的直方图 ) subplot(133) imhist(M) axis(0 2 0 2) title(均值滤波后的直方图 ) 整体而言中值滤波的结果比均值滤波的结果要好，均值滤波模糊了轮廓边缘，但是中值滤波对椒盐噪声的滤波破坏了轮廓边缘。中值滤波对高斯噪声的水平滤波不是很好。 程序： I=zeros(256,256); for i=(32:24:224); I(23:233,i:i+7)=1; end figure imsho
7、w(I) title(原图像 ) A=imnoise(I,gaussian,0,0.025); figure subplot(121) imshow(A) title(高斯白噪声 ) B=imnoise(I,salt subplot(122) imshow(B) title(椒盐噪声 ) h=1 1 1;1 1 1;1 1 1;%领域平均法消除噪声模板为 1/91 1 1;1 1 1;1 1 1 h=9./h; Ja=conv2(A,h); figure subplot(121) imshow(Ja,) title(对高斯白噪声的均值滤波 ) Jb=conv2(B,h); subplot(122) imshow(Jb,) title(对椒盐噪声的均值滤波 ) Ka=medfilt2(A,3,3) figure subplot(121) imshow(Ka,) title(对高斯白噪声的中值滤波 ) Kb=medfilt2(B,3,3); subplot(122) imshow(Kb,) title(对椒盐噪声的中值 波 )
