matlab - 图像处理——使用FFT放大图像(Matlab代码)

Question

我正在尝试编写一个简单的 matlab 代码，它使用 fft 放大图像。我尝试了已知的图像扩展算法，该算法计算图像的傅立叶变换，用零填充它并计算填充图像的逆傅立叶。 然而，逆傅里叶变换返回包含复数的图像。 因此，当我尝试使用 imshow 显示结果时，出现以下错误:

Warning: Displaying real part of complex input.

你知道我做错了什么吗？

我的代码:

im = imread('fruit.jpg');
imFFT = fft2(im);
bigger = padarray(imFFT,[10,10]);
imEnlarged = ifft2(bigger);

谢谢!

Answer 1

这是因为 FFT 返回对应于从 0 到 Fs 的离散(空间)频率的值，其中 Fs 是(空间)采样率。您需要在高频处插入零，它们位于返回的 FFT 的中心，而不是末端。

您可以使用 fftshift 将高频移动到末尾，用零填充，然后使用 ifftshift 移回(感谢@Shai 的更正):

bigger = ifftshift(padarray(fftshift(imFFT),[10,10]));

另请注意，用零填充会降低放大图像中的值。您可以使用合适的放大系数 amp 更正它，在本例中等于 (1+2*10/length(im))^2:

bigger = ifftshift(padarray(fftshift(amp*imFFT),[10,10]));