test_img11_3.m

来自「对于普通图象的一些频域的分析小程序,可以对图象处理参考作用.」· M 代码 · 共 36 行

clear;

I=imread('img11.bmp');
I2=rgb2gray(I);
figure;
subplot(3,3,1)
imshow(I2);
title('原图像')
I3=fft2(I2);
subplot(3,3,4)
imshow(I3);
title('FFT变换结果');
subplot(3,3,7)
I4=fftshift(I3);
W=abs(I3);
imshow(I4);
title('零点平移FFT');
subplot(3,3,2)
imshow(log(abs(I3)),[]);
title('幅度谱')
subplot(3,3,6)
imshow(log(abs(I4)),[]);
title('零点平移的幅度谱');
FUF=ifft2(I3);          %幅度的IDFT
subplot(3,3,3)
imshow(log(abs(FUF)),[]);
title('幅频IDFT图');
AN=angle(I3);
subplot(3,3,8)
imshow(AN);
title('相位谱');
FUX=ifft2(AN);
subplot(3,3,9)
imshow(log(abs(FUX)),[]);
title('相频IDFT图');