数字图像处理大作业

《数字图像处理大作业》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数字图像处理上机实验学院：电子工程学院学号：1202121443姓名：王美妙1题1题11题2f1与f2的幅度谱相同，因为f1与f2的振幅绝对只相同。1题3f3的幅度谱与f2的幅度谱顺时钟旋转90度后相同1题4f4的幅度谱与f1的幅度谱顺时钟旋转90度后相同1题5f5的幅度谱是f1与f4的幅度谱相加，f6的幅度谱是f2与f3的幅度谱相加附1题代码：f1=zeros(256,256);f1(65:192,113:144)=100;figuresubplot(121)imshow(f1)title('原图f1')F=fft2(f1);subplot(122)imsh

2、ow(abs(F))title('f1幅度谱图')form=(1:256);forn=(1:256);f2(m,n)=power(-1,m+n)*f1(m,n);endendfiguresubplot(121)imshow(f2)title('f2图')F1=fft2(f2);subplot(122)imshow(abs(F1))title('f2幅度谱图')f3=imrotate(f2,-90);figuresubplot(121)imshow(f3)title('f3图')F2=fft2(f3);subplot(122)imshow(abs(F2))tit

3、le('f3幅度谱图')f4=imrotate(f1,-90);figuresubplot(121)imshow(f4)title('f4图')F5=fft2(f4);subplot(122)imshow(abs(F5))title('f4幅度谱图')f5=f1+f4;figuresubplot(121)imshow(f5)title('f5图')F3=fft2(f5)subplot(122)imshow(abs(F3))title('f5幅度谱图')f6=f2+f3;figuresubplot(121)imshow(f6)title('f6图')F4=fft

4、2(f6)subplot(122)imshow(abs(F4))title('f6幅度谱图')2题作256*256的图像如下从中值滤波和均值滤波后的图像可以看出中值滤波把小黑色正方形快的角抹去了，而均值滤波则把图像变模糊了尤其是边缘模糊的比较厉害。且可以看出中值滤波后图像的直方图没有变化，而均值滤波后的直方图变化了。附2题代码：A=zeros(256,256);fori=(1:64:256);forj=(1:64:256);A(j:j+32,i:i+32)=1;endendform=(32:64:256);forn=(32:64:256);A(n:n+32,m

5、:m+32)=1;endendfigureimshow(A)title('原图像')figureK=medfilt2(A,[3,3])subplot(121)imshow(K)title('3*3中值滤波')M=filter2(fspecial('average',3),A)subplot(122)imshow(M)title('3*3均值值滤波')figuresubplot(131)imhist(A)axis([0202])title('原图直方图')subplot(132)imhist(K)axis([0202])title('中值滤波后的直方图')sub

6、plot(133)imhist(M)axis([0202])title('均值滤波后的直方图')3题整体而言中值滤波的结果比均值滤波的结果要好，均值滤波模糊了轮廓边缘，但是中值滤波对椒盐噪声的滤波破坏了轮廓边缘。中值滤波对高斯噪声的水平滤波不是很好。附4题程序：I=zeros(256,256);fori=(32:24:224);I(23:233,i:i+7)=1;endfigureimshow(I)title('原图像')A=imnoise(I,'gaussian',0,0.025);figuresubplot(121)imshow(A)title('高斯白噪

7、声')B=imnoise(I,'salt&pepper',0.025);subplot(122)imshow(B)title('椒盐噪声')h=[111;111;111];%%领域平均法消除噪声模板为1/9[111;111;111]h=9./h;Ja=conv2(A,h);figuresubplot(121)imshow(Ja,[])title('对高斯白噪声的均值滤波')Jb=conv2(B,h);subplot(122)imshow(Jb,[])title('对椒盐噪声的均值滤波')Ka=medfilt2(A,[3,3])figuresubplot(121

8、)imshow(Ka,[])title