数字图像处理上机作业.docx

数字图像处理上机实验题 一、产生右图所示图像 f1(m,n)，其中图像大小为256×256，中间亮条为128×32，暗处=0，亮处=100。对其进行 FFT： 屏显示原图 f1(m,n)和FFT(f1)的幅度谱图； 2、令f2(m,n)=（-1）^（m+n）*f1(m,n)，重复以上过程，比较二者幅度谱的异同，简述理由； 3、若将f2(m,n)顺时针旋转 90 度得到f3(m,n)，试显示 FFT(f3)的幅度谱，并与 FFT(f2) 的幅度谱进行比较； 4、若将f1(m,n) 顺时针旋转 90 度得到f4(m,n)，令f5(m,n)=f1(m,n)+f4(m,n)，试显示 FFT(f5) 的幅度谱，并指出其与 FFT(f1)和FFT(f4)的关系； 5、若令f6(m,n)=f2(m,n)+f3(m,n)，试显示 FFT(f6)的幅度谱，并指出其与 FFT(f 2)和FFT(f3) 的关系，比较 FFT(f6)和FFT(f5)的幅度谱。 代码 f1=zeros(256,256); for i =64:1:191 for j = 112:1:143 f1(i,j) = 100; end end f2 = fft2(f1); %f2(m,n) = f3 f3 = ((-1)^(i+j))*f1; f4 = fft2(f3); %f3(m,n) = f5 f5 = imrotate(f3,90,bilinear); f6 = fft2(f5); %f4(m,n) = f7 f7 = imrotate(f1,90,bilinear); f8 = fft2(f7); %f5(m,n) = f8 f9 = f1 + f7; f10 = fft2 (f9); %f6(m,n) = f2(m,n)+f3(m,n) f11 = f3 + f5; f12 = fft2(f11); figure(1) subplot(1,2,1); imshow(abs(f1)); title(原图f1); subplot(1,2,2); imshow(abs(f2)); title(幅度谱fft2(f1)); figure(2) subplot(2,2,1) imshow(abs(f1)); title(原图f1) subplot(2,2,2) imshow(abs(f2)); title(幅度谱fft2(f1)); subplot(2,2,3); imshow(abs(f3)) title(变换谱f2); subplot(2,2,4); imshow(abs(f4)); title(幅度谱fft2(f2)); figure(3) subplot(2,2,1) imshow(abs(f3)) title(变换谱f2); subplot(2,2,2); imshow(abs(f4)); title(幅度谱fft2(f2)); subplot(2,2,3); imshow(abs(f5)) title(变换谱f3); subplot(2,2,4); imshow(abs(f6)); title(幅度谱fft2(f3)); figure(4) subplot(3,2,1); imshow(f7); title(f1旋转图f4); subplot(3,2,2); imshow(abs(f8)); title(幅度谱fft2(f4)); subplot(3,2,3); imshow(f9); title(f5(m,n)=f1+f4); subplot(3,2,4); imshow(abs(f10)); title(幅度谱fft2(f5)); subplot(3,2,5) imshow(abs(f1)); title(原图f1); subplot(3,2,6); imshow(abs(f2)); title(幅度谱fft2(f1)); figure(5) subplot(3,2,1) imshow(abs(f3)) title(变换谱f2); subplot(3,2,2); imshow(abs(f4)); title(幅度谱fft2(f2)); subplot(3,2,3); imshow(abs(f5)) title(变换谱f3); subplot(3,2,4); imshow(abs(f6)); title(幅度谱fft2(f3)); subplot(3,2,5) imshow(abs(f11)) title(变换谱f6=f2+f3); subplot(3,2,6); imshow(abs(f12)); title(幅度谱fft2(f6)); f