Perché questa trasformata di Fourier inversa non fornisce i risultati corretti?

Question

Voglio invertire la trasformata di Fourier di un'immagine in MATLAB, ma il risultato non è l'immagine originale (come dovrebbe essere). Vi sono ovviamente alcuni dettagli di implementazione che non so che causino il problema. Ecco il codice:

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

Answer 1

Dal FFT2 e IFFT2 entrambi ingressi solo di supporto di tipo double e single, il vostro image data (che rischia di tipo uint8) viene convertito in digitare double prima prima di essere elaborati da FFT2. Sarà quindi necessario convertire l'immagine in uscita inv di nuovo ad un unsigned 8-bit integer utilizzando la funzione di UINT8 per recuperare l'immagine originale:

>> 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: Come ulteriore suggerimento, vorrei evitare di nominare le variabili fft e inv poiché le funzioni con questi nomi esistono già in MATLAB.

Answer 2

Anche se si sta tentando di eseguire un'immagine FFT a colori (24 bit), tenere presente che imread() restituirà un array M x N x 3. Quindi dovresti eseguire FFT su ciascun canale R/G/B separatamente.

See this per dettagli.