MATLAB 图像处理：傅立叶变换、频谱分析与图像重建

% 读取图像文件 I = imread('cameraman.tif');

% 傅立叶变换 F = fft2(I);

% 零频率分量移位至图形中心 F_shift = fftshift(F);

% 傅立叶幅度谱 F_mag = abs(F_shift);

% 相位谱 F_phase = angle(F_shift);

% 显示图像和傅立叶变换结果 figure; subplot(2,2,1); imshow(I); title('Original Image'); subplot(2,2,2); imshow(log(1+F_mag),[]); title('Magnitude Spectrum'); subplot(2,2,3); imshow(F_phase,[]); title('Phase Spectrum'); subplot(2,2,4); imshow(ifft2(ifftshift(F_shift)),[]); title('Reconstructed Image');