数字图像处理作业合编.doc

数字图像处理上机实验题 基础部分 1. 产生右图所示图像 f1(m,n)，其中图像大小为256×256，中间亮条为128×32，暗处=0，亮处=100。对其进行FFT： 同屏显示原图f1(m,n)和FFT(f1)的幅度谱图； 若令f2(m,n)=（-1）m+n f1(m,n)，重复以上过程，比较二者幅度谱的异同，简述理由； 若将f2(m,n)顺时针旋转90度得到f3(m,n)，试显示FFT(f3)的幅度谱，并与FFT(f2)的幅度谱进行比较； 若将f1(m,n) 顺时针旋转90度得到f4(m,n)，令f5(m,n)=f1(m,n)+f4(m,n)，试显示FFT(f5)的幅度谱，并指出其与 FFT(f1)和FFT(f4)的关系； 若令f6(m,n)=f2(m,n)+f3(m,n)，试显示FFT(f6)的幅度谱，并指出其与 FFT(f2)和FFT(f3)的关系，比较FFT(f6)和FFT(f5)的幅度谱。 程序如下： m=256; n=256; f1=zeros(m,n); f1(65:192,113:144)=100; FFTf1=fft2(f1); ff1=log(1+abs(FFTf1)); figure(1),subplot(131),imshow(f1),title(原图像f1); subplot(132),imshow(abs(FFTf1),[]),title(原图像的傅里叶变换); subplot(133),imshow(abs(ff1),[]),title(以对数坐标显示); for x=1:m for y=1:n f2(x,y)=(-1)^(x+y)*f1(x,y); end end FFTf2=fft2(f2); ff2=log(1+abs(FFTf2)); figure(2),subplot(131),imshow(f2),title(原图像乘-1的m+n次方得到f2); subplot(132),imshow(abs(FFTf2),[]),title(f2的傅里叶变换); subplot(133),imshow(abs(ff2),[]),title(以对数坐标显示); f3=imrotate(f2,-90); FFTf3=fft2(f3); ff3=log(1+abs(FFTf3)); figure(3),subplot(131),imshow(f3),title(f2顺时针旋转90度得到f3); subplot(132),imshow(abs(FFTf3),[]),title(f3的傅里叶变换); subplot(133),imshow(abs(ff3),[]),title(以对数坐标显示); f4=imrotate(f1,-90); f5=f1+f4; FFTf4=fft2(f4); FFTf5=fft2(f5); ff4=log(1+abs(FFTf4)); ff5=log(1+abs(FFTf5)); figure(4),subplot(131),imshow(f4),title(f1顺时针旋转90度得到f4); subplot(132),imshow(abs(FFTf3),[]),title(f4的傅里叶变换); subplot(133),imshow(abs(ff3),[]),title(以对数坐标显示); figure(5),subplot(131),imshow(f5),title(f5=f1+f4); subplot(132),imshow(abs(FFTf5),[]),title(f5的傅里叶变换); subplot(133),imshow(abs(ff5),[]),title(以对数坐标显示); f6=f2+f3; FFTf6=fft2(f6); ff6=log(1+abs(FFTf6)); figure(6),subplot(131),imshow(f6),title(f6=f2+f3); subplot(132),imshow(abs(FFTf6),[]),title(f6的傅里叶变换); subplot(133),imshow(abs(ff6),[]),title(以对数坐标显示); 结果分析： ① 同屏显示原图f1(m,n)和FFT(f1)的幅度谱图； f1的幅度谱 ② 若令f2(m,n)=(-1)^(m+n)*f1(m,n)，重复以上过程，比较二者幅度谱的异同，简述理由； f2的幅度谱 f1(m,n)的幅度谱中，能量集中在四个角上，f2(m,n)的幅度谱能量集中在中心位置附近。 根据傅里叶变换的频域移位特性，f1(m,n)乘(-1）^(m+n),在频域里面对应其幅度谱平移M/2和N/2，即进行了频谱中心化。 ③ 若将f2(m,n)顺