Why image is shifted when using ifft2

Question

Khang Truong 2019-10-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/483287-why-image-is-shifted-when-using-ifft2

评论： Yulin WANG 2020-11-15

Hi, I'm implementing 2-D convolution by using FFT. Here is my code:

img = im2single(imread('dog.bmp'));
filter = fspecial('gaussian', 53, 3);
F = fft2(img);
mask = fft2(filter, size(img, 1), size(img, 2));
filtered_img = ifft2(F .* mask);
imshow(real(filtered_img));

Here is the original image:

And the result:

Why this happens? How can I fix it? Please help me.

Many thanks.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Yulin WANG 2020-11-15

Well, you are applying a Gaussian Filter on an image in frequency domain and convert it back to the spatial domain, and that's why the image becomes a bit blurry than the orginal one.

I did not see anything wrong with this.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Cris LaPierre 2019-10-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/483287-why-image-is-shifted-when-using-ifft2#answer_394569

编辑：Cris LaPierre 2019-10-3

在 MATLAB Online 中打开

The shift is related to your hsize value in fspecial (shifted ~0.5*hsize in both X and Y).

I'm not sure I can do any better explaining than what you can find googling, but you need to center your frequency domain in the center of the image. You can do this using fftshift and ifftshift. See here, and here.

Also note that, in the comment of the second link, there is a recommendation to use psf2otf to take the FFT of a point spread function.

After some playing around with a built in image, I came up with this code:

img = im2single(imread('autumn.tif'));
imshow(img)
filter = fspecial('gaussian', 53, 2);
F = fftshift(fft2(img));
mask = fftshift(psf2otf(filter,[size(img, 1), size(img, 2)]));
filtered_img = ifft2(ifftshift(F .* mask));
imshow(abs(filtered_img));