MATLAB中脑部MRI图像傅里叶变换（Fourier transform）结果异常求助

兄弟，我看了你的代码和问题描述，马上就发现几个关键问题，这就是为啥你看不到预期的中心白十字，逆变换还出噪声的原因：

问题1：直接显示复数形式的傅里叶变换结果 #

fft2()输出的是复数矩阵，里面包含幅度和相位信息。你直接用imshow()显示复数矩阵时，MATLAB只会取实部来显示，而且傅里叶变换的幅度值动态范围极大——低频成分（对应图像的整体结构，也就是你要的中心十字）的幅度远远高于高频细节，直接显示的话，大部分像素都会因为数值太小被压成黑色，少数极大值变成纯白，看起来就像噪声。

问题2：没对幅度谱做对数压缩 #

就算你取了幅度，直接显示还是看不到清晰的十字，因为低频和高频的幅度差了好几个数量级。必须用对数变换来压缩这个动态范围，让不同幅度的频率成分都能被肉眼分辨。

问题3：逆变换的步骤错误 #

你直接对移位后的brain_DFT_shifted做ifft2()，相当于对已经移到中心的频谱做逆变换，得到的图像是移位后的版本，而且因为复数的虚部（数值计算误差），显示出来也会异常。正确的做法是先把移位的频谱移回去，再做逆变换。

修正后的完整代码 #

% 读取并预处理图像
brain = imread('brain.jpeg');
brain_bw = rgb2gray(brain);

% 计算傅里叶变换并移位
brain_DFT = fft2(brain_bw);
brain_DFT_shifted = fftshift(brain_DFT);

% 计算幅度谱并做对数压缩
brain_amp = abs(brain_DFT_shifted);
brain_amp_log = log(1 + brain_amp); % 加1是为了避免log(0)的问题

% 显示原始图像和处理后的傅里叶幅度谱
figure;
subplot(1,2,1);
imshow(brain);
title('Original Image');
subplot(1,2,2);
imshow(brain_amp_log, []); % 加[]让MATLAB自动调整显示范围
title('Fourier Transform (Log Amplitude)');

% 正确的逆变换步骤
brain_DFT_unshifted = ifftshift(brain_DFT_shifted); % 先移回原始频谱位置
inverse = ifft2(brain_DFT_unshifted);
inverse_real = real(inverse); % 取实部消除数值误差带来的虚部

% 显示逆变换结果
figure();
imshow(inverse_real, []);
title('Inverse Fourier Transform');

补充说明 #

• 用abs()取幅度谱是关键，这才是反映图像频率分布的直观指标；
• 对数变换log(1 + x)是处理傅里叶幅度谱的常规操作，能把大范围的数值压缩到肉眼可辨的区间；
• 逆变换时必须用ifftshift()把移位后的频谱还原到原始位置，再做ifft2，最后取实部是为了消除浮点运算带来的微小虚部，保证图像正常显示。

备注：内容来源于stack exchange，提问作者Shika93