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=fft
2、2(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);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)title
3、(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=fft2(f6)subplot(122)imshow(abs(F4
4、)title(f6幅度谱图 )2 题作 256*256 的图像如下从中值滤波和均值滤波后的图像可以看出中值滤波把小黑色正方形快的角抹去了，而均值滤波则把图像变模糊了尤其是边缘模糊的比较厉害。且可以看出中值滤波后图像的直方图没有变化，而均值滤波后的直方图变化了。附 2 题代码：A=zeros(256,256);for i=(1:64:256);for j=(1:64:256);A(j:j+32,i:i+32)=1;endendfor m=(32:64:256);for n=(32:64:256);A(n:n+32,m:m+32)=1;endendfigureimshow(A)title(原 图像
5、)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(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(均 值滤波后的直方图)3 题整体而言中值滤波的结果比均值滤
6、波的结果要好，均值滤波模糊了轮廓边缘，但是中值滤波对椒盐噪声的滤波破坏了轮廓边缘。中值滤波对高斯噪声的水平滤波不是很好。附 4 题程序：I=zeros(256,256);for i=(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(高斯白噪声)B=imnoise(I,salt subplot(122)imshow(B)title(椒 盐噪声)h=1 1 1;1 1 1;1 1 1;%领域平均法消除噪声模
7、板为1/91 1 1;1 1 1;1 1 1h=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)imshow(Ka,)title(对 高斯白噪声的中值滤波)Kb=medfilt2(B,3,3);subplot(122)imshow(Kb,)title(对 椒盐噪声的中值波)4 题roberts 自动选择的阈值为：0.1305sob
8、el 自动选择的阈值为：0.1162prewitt 自动选择的阈值为：0.1138附 4 题程序：I=imread(rice.png);figureimshow(I)title(原始 图像)BW1,thresh1=edge(I,roberts);disp(roberts自动选择 的阈值为 ：)disp(thresh1)figuresubplot(131)imshow(BW1)title(roberts算子)BW2,thresh2=edge(I,sobel);disp(sobel自动选择的阈值为：)disp(thresh2)subplot(132)imshow(BW2)title(sobel算子
9、)BW3,thresh3=edge(I,prewitt);disp(prewitt自动选择 的阈值为：)disp(thresh3)subplot(133)imshow(BW3)title(prwitte算子)BW4=filter2(fspecial(Prewitt),I);BW5=filter2(fspecial(Prewitt),I);BW6=filter2(fspecial(Sobel),I);figuresubplot(2,2,1);imshow(I);title(灰度 图像);subplot(2,2,2);imshow(BW4); title(Roberts锐化结果); subplot
10、(2,2,3);imshow(BW5);title(Prewitt锐化 结果);subplot(2,2,4);imshow(BW6);title(Sobel锐化 结果);5 题迭代后的阈值：131附 5 题程序：clear all;I=imread(rice.png);ZMAX=max(max(I);ZMIN=min(min(I);TK=(ZMAX+ZMIN)/2;bcal=1;ISIZE=size(I);while(bcal)ifground=0;ibground=0;FgroundS=0;BgroundS=0;for i=1:ISIZE(1)for j=1:ISIZE(2)tmp=I(i,
11、j);if(tmp=TK)ifground=ifground+1;FgroundS= FgroundS+double(tmp);elseibground=ibground+1;BgroundS= BgroundS+double(tmp);endendendZO=FgroundS/ifground;ZB=BgroundS/ibground;TKTmp=uint8(ZO+ZB)/2);if(TKTmp=TK)bcal=0;elseTK=TKTmp;endenddisp(strcat(迭代后的阈值：,num2str(TK);newI=im2bw(I,double(TK)/255);subplot(121)imshow(I)title(原始 图像)subplot(122)imshow(newI)title(迭代法分割得到的二值图像)
