图像频谱的相位谱谱恢复图像

要通过相位谱恢复一张图像，可以使用 Matlab 中的傅里叶反变换函数。具体步骤如下：

1. 读入图片并转换为灰度图像。

img = imread('image.png');
gray_img = rgb2gray(img);

2. 对灰度图像进行二维傅里叶变换，并将频谱移到中心。

f = fft2(double(gray_img));
fshift = fftshift(f);

3. 计算频谱的幅度谱和相位谱。

magnitude_spectrum = log(abs(fshift));
phase_spectrum = angle(fshift);

4. 将相位谱恢复为频谱。

fshift_new = abs(fshift) .* exp(1i * phase_spectrum);

5. 对新的频谱进行傅里叶反变换。

f_new = ifftshift(fshift_new);
img_new = real(ifft2(f_new));

6. 显示结果。

subplot(1,2,1), imshow(gray_img), title('原图');
subplot(1,2,2), imshow(uint8(img_new)), title('恢复图像');

上述代码中的 `subplot` 函数用于将两张图片显示在同一个窗口中。其中，第一个参数表示行数，第二个参数表示列数，第三个参数表示当前图片的位置。`imshow` 函数用于显示图像，`title` 函数用于设置标题。通过这些步骤，就可以通过相位谱恢复一张图像。