Pourquoi une partie réelle de la FFT convertit l'image en rotation + original?

16

J'ai lu cette image:

entrez la description de l'image ici

pris sa FFT (2D) puis FFT inverse pour récupérer exactement l'image. Le code est fourni pour référence:

imfft = fft2(photographer);
im = uint8(ifft2(imfft));

imshow(im); %Output is same image

Mais quand je change le fourier et ne prends que la vraie partie,

imfft = real(fft2(photographer));
im = uint8(ifft2(imfft));
imshow(im);

J'obtiens une image comme celle-ci ( notez que le changement de taille n'est pas pertinent et uniquement en raison de son enregistrement à partir du gestionnaire de figures Matlab ):

entrez la description de l'image ici

Quelqu'un peut-il m'expliquer la théorie (mathématiques) derrière cela? Merci

Scientifique en échec
la source

Réponses:

16

Disons que votre image est donnée par je(X,y)

jeF(ωX,ωy)=Xyje(X,y)ejωXXejωyyXy

Maintenant, vous prenez la partie réelle et effectuez l'inverse:

jem(α,β)=ωXωy{jeF(ωX,ωy)}ejωXαejωyβωXωy=ωXωy{Xyje(X,y)ejωXXejωyyXy}ejωXαejωyβωXωy=Xyje(X,y)ωXωy{ejωXXejωyy}ejωXαejωyβωXωyXy

cos(ωXX)cos(ωyy)+péché(ωXX)péché(ωyy)
12[δ(X-α)δ(y-β)+δ(X+α)δ(y+β)]

jem

jem(X,y)=12[je(X,y)+je(-X,-y)]

X,y>0N

jem(X,y)=12[je(X,y)+je(N-X,M-y)]
N,M
ThP
la source
Bonne réponse! +1
Peter K.
3
I think you can see now why got that result.Oui. Cependant, puisque cette question a frappé la liste HNQ, vous pourriez peut-être envisager d'ajouter l'étape finale pour ceux qui viennent de sites inclinés moins mathématiques.
Mât
9

z(X,y)z(-X,-y)sur l'origine. Notez que l'origine ici serait le centre de l'espace de Fourier. Cela peut être reformulé, bien sûr, si le composant DC n'est pas au centre de votre implémentation FFT. Et voici ce que vous voyez dans votre image: une version à réflexion ponctuelle recouvre la véritable image - parce que vous avez forcé un espace à être réellement valorisé.

Cette propriété est en fait utilisée pour accélérer l'imagerie par résonance magnétique (IRM) dans certains cas: l'IRM acquiert les données directement dans l'espace de Fourier. Puisqu'une image MR idéale ne peut être décrite que par des valeurs réelles (tous les vecteurs de magnétisation excités ont la phase 0), vous n'avez qu'à acquérir la moitié de l'espace de données, ce qui vous fait gagner la moitié du temps d'imagerie. Bien sûr, les images IRM ne sont pas complètement réelles en raison des limites de la réalité ... mais avec quelques astuces, vous pouvez toujours utiliser cette technique de manière avantageuse.

M529
la source
2
J'ai aimé la manière simple de dire la même réponse que ThP a fournie. Et merci pour les informations sur l'IRM. Je n'en savais rien.
Échec scientifique