¿Por qué esta transformación de Fourier inversa no da los resultados correctos?

Question

Quiero invertir la transformada de Fourier de una imagen en MATLAB, pero el resultado no es la imagen original (como debería ser). Obviamente hay algunos detalles de implementación que desconozco que están causando el problema. Aquí está el código:

img = imread('img.jpg'); 
fft = fft2(img); 
inv = ifft2(fft); 
imshow(inv); 

Answer 1

Desde FFT2 y IFFT2 dos únicas entradas de apoyo de tipo double y single, su image data (que es probable de tipo uint8) se convierte al tipo double primero antes de ser procesada por FFT2. Por lo tanto, tendrá que convertir la imagen de salida inv de nuevo a un entero sin signo de 8 bits utilizando la función UINT8 para recuperar la imagen original:

>> img = imread('peppers.png'); %# Load a sample image 
>> fft = fft2(img); %# Get the Fourier transform 
>> inv = ifft2(fft); %# Get the inverse Fourier transform 
>> inv = uint8(inv); %# Convert to uint8 
>> imshow(inv);  %# Show the image 
>> isequal(img,inv) %# Test if inv matches the original image img 

ans = 

    1    %# It does! 

NOTA: Como una propina adicional, me gustaría evitar nombrar las variables fft y inv ya que las funciones con esos nombres ya existen en MATLAB.

Answer 2

Además, si está tratando de hacer FFT en una imagen de color (24 bits), tenga en cuenta que imread() devolverá la matriz M x N x 3. Por lo tanto, debe realizar FFT en cada canal R/G/B por separado.

See this para más detalles.