matlab傅里叶变换相位谱 相位重建

在Matlab中，可以使用fft2函数对图像进行二维离散傅里叶变换。变换后得到的结果是包含幅度谱和相位谱的复数矩阵。相位谱可以通过angle函数获取，单位为弧度。要进行相位重建，可以将相位谱与幅度谱结合起来，然后使用ifft2函数进行傅里叶逆变换，即可得到重建后的图像。

以下是一个示例代码片段，展示了如何实现相位重建：

I = imread('exp2.tif'); % 读入原图像
I1 = I(:,:,1:3); % 四通道转为三通道去掉透明层
I2 = rgb2gray(I1); % 对图像进行灰度处理
fftl = fft2(I2); % 二维离散傅里叶变换

% 获取幅度谱和相位谱
A = abs(fftl); % 幅度谱
T = angle(fftl); % 相位谱

% 进行相位重建
reconstructed = ifft2(A .* exp(1i*T));

% 显示重建结果
imshow(reconstructed);
title('相位重建图像');

请注意，相位谱的单位为弧度，要将其转换为角度可以使用angle函数乘以180再除以pi。