Warum liefert diese inverse Fourier-Transformation nicht die korrekten Ergebnisse?

Question

Ich möchte die Fourier-Transformation eines Bildes in MATLAB invertieren, aber das Ergebnis ist nicht das Originalbild (wie es sein sollte). Es gibt offensichtlich einige Implementierungsdetails, von denen ich nicht weiß, dass dies das Problem verursacht. Hier ist der Code:

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

Answer 1

Seit FFT2 und IFFT2 beide nur Unterstützung Eingänge des Typs double und single, Ihre image data (die uint8 wahrscheinlich vom Typ ist) umgewandelt wird double bevor eingeben, indem FFT2 verarbeitet werden. Sie werden also Ihre Ausgangsbild inv zurück zu einem unsigned 8-Bit-Integer mit der Funktion konvertieren müssen UINT8 das Originalbild wiederherzustellen:

>> 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! 

HINWEIS: Als zusätzliche Spitze, würde ich vermeiden, dass Ihre Benennung von Variablen fft und inv, da Funktionen mit diesen Namen bereits in MATLAB vorhanden sind.

Answer 2

Auch wenn Sie versuchen, FFT auf Farbe (24-Bit) Bild - beachten Sie, dass imread() wird M x N x 3-Array zurückgeben. Sie sollten also FFT für jeden R/G/B-Kanal separat ausführen.

See this für Details.