image processing - enlarging an image using FFT (Matlab code)

Question

I'm trying to write a simple matlab code which enlarges an image using fft. I tried the known algorithm for image expansion, which computes the Fourier transform of the image, pads it with zeros and computes the inverse Fourier of the padded image. However, the inverse Fourier transform returns an image which contains complex numbers. Therefore, when I'm trying to show the result using imshow, I'm getting the following error:

Warning: Displaying real part of complex input.

Do you have an idea what am I doing wrong?

my code:

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

Thanks!

Luis Mendo · Accepted Answer · 2014-04-10T21:47:03.793

6

That's because the FFT returns values corresponding to the discrete (spatial) frequencies from 0 through Fs, where Fs is the (spatial) sampling rate. You need to insert zeros at high frequencies, which are located at the center of the returned FFT, not in its end.

You can use fftshift to shift the high frequencies to the end, pad with zeros, and then shift back with ifftshift (thanks to @Shai for the correction):

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

Also, note that padding with zeros decreases the values in the enlarged image. You can correct that using a suitable amplification factor amp, which in this case would be equal to (1+2*10/length(im))^2:

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

edited Apr 10 '14 at 21:47

answered Nov 05 '13 at 11:23

Luis Mendo

110,752
13
76
147

1

nice answer. However, shouldn't you be using `ifftshift` after padding? – Shai Nov 05 '13 at 11:40
@Shai You're right. `fftshift` and `ifftshift` are the same only for even size – Luis Mendo Nov 05 '13 at 12:50
Nice. I've added this question as a favorite just to bookmark your answer. +1 – chappjc Nov 08 '13 at 19:57

score 3 · Answer 2 · edited May 23 '17 at 11:52

3

You can pad at the higher frequencies directly (without fftshift suggested by Luis Mendo)

>> BIG = padarray( amp*imFFT, [20 20], 0, 'post' );
>> big = ifft2( BIG );

edited May 23 '17 at 11:52

Community

1
1

answered Nov 05 '13 at 11:39

Shai

111,146
38
238
371

Thank you both for your answers, however, i'm still getting complex results. is that ok? or should I just show the absolute value of the given image? – ibezito Nov 05 '13 at 11:53
@user2955997 if the ansewr is no good - do not "accept" it. you can use comments to ask for calrifications. It might be the case that you have a very small imagenary component due to numerical inaccurecies. look at the imagenary part: `>> figure;imagesc(imag(big));axis image;impixelinfo` and see how large is the error. – Shai Nov 05 '13 at 12:15

hotpaw2 · Answer 3 · 2013-11-05T16:32:44.837

If you want a strictly real result, then before you do the IFFT you need to make sure the zero-padded array is exactly conjugate symmetric. Adding the zeros off-center could prevent this required symmetry.

Due to finite numerical precision, you may still end up with a complex IFFT result, but the imaginary components will all be tiny values that are essentially equivalent to zero.

Your FFT library may contain a half-to-real (quarter-size input for 2D) version that enforces symmetry and throws away the almost-zero numerical noise for you.

image processing - enlarging an image using FFT (Matlab code)

3 Answers3