clear, clc
I = imread('...');
F = fftshift(fft2(I)); % 对图像进⾏⼆维 DFT(fft2)，并移⾄中⼼位置magn = log(abs(F)); % 加 log 是便于显⽰，缩⼩值域phase = log(angle(F)*180/pi); % 转换为度数subplot(2, 2, 1), imshow(I, []), title('原图像')
subplot(2, 2, 2), imshow(log(F), []), title('图像频谱')subplot(2, 2, 3), imshow(magn, []), title('图像幅度谱')subplot(2, 2, 4), imshow(phase, []), title('图像相位谱')
abs(F) == sqrt(real(F).^2 + imag(F).^2)
angle(F) == log(imag(F)) == atan2(imag(F), real(F))
因篇幅问题不能全部显示，请点此查看更多更全内容