图像二维离散傅里叶变换、幅度谱、相位谱

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('图像相位谱')