¿Por qué parte real de FFT convierte la imagen en rotación + original?

16

He leído esta imagen:

ingrese la descripción de la imagen aquí

tomó su FFT (2D) y luego FFT inversa para recuperar exactamente la imagen. El código se proporciona para referencia:

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

imshow(im); %Output is same image

Pero cuando cambio el fourier y tomo solo la parte real,

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

Obtengo una imagen como esta ( tenga en cuenta que el cambio de tamaño es irrelevante y solo debido a guardarlo desde el controlador de figura de Matlab ):

ingrese la descripción de la imagen aquí

¿Alguien puede explicarme la teoría (matemáticas) detrás de esto? Gracias

Científico fallido
fuente

Respuestas:

16

Digamos que tu imagen está dada por . Entonces su transformación de Fourier viene dada por I f ( ω x , ω y ) = x y I ( x , y ) e j ω x x e j ω y y d x d yyo(X,y)

yoF(ωX,ωy)=Xyyo(X,y)mijωXXmijωyyreXrey

Ahora tomas la parte real y realizas la inversa:

yometro(α,β)=ωXωy{yoF(ωX,ωy)}mijωXαmijωyβreωXreωy=ωXωy{Xyyo(X,y)mijωXXmijωyyreXrey}mijωXαmijωyβreωXreωy=Xyyo(X,y)ωXωy{mijωXXmijωyy}mijωXαmijωyβreωXreωyreXrey

Puede ver claramente que la integral interna es la transformación 2D de Fourier de que es

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

Sustituyendo el resultado por obtiene yometro

yometro(X,y)=12[yo(X,y)+yo(-X,-y)]

Por supuesto, en su caso , sin embargo, la transformada discreta de Fourier supone que su señal es -periódica y obtiene donde son las dimensiones de su imagen. Creo que ahora puedes ver por qué obtuviste ese resultado.X,y>0 0norte

yometro(X,y)=12[yo(X,y)+yo(norte-X,METRO-y)]
norte,METRO
ThP
fuente
¡Buena respuesta! +1
Peter K.
3
I think you can see now why got that result.Si. Sin embargo, dado que esta pregunta llegó a la lista HNQ, tal vez consideraría agregar el paso final para aquellos que vienen de sitios inclinados menos matemáticos.
Mástil
9

z(X,y)z(-X,-y)Sobre el origen. Tenga en cuenta que el origen aquí sería el centro del espacio de Fourier. Esto puede reformularse, por supuesto, si el componente DC no está en el centro de su implementación de FFT. Y esto es lo que ves en tu imagen: una versión de punto reflejado está superponiendo la imagen verdadera, porque forzaste a un espacio a ser realmente valorado.

En realidad, esta propiedad se está utilizando para acelerar la resonancia magnética (MRI) en algunos casos: la MRI adquiere los datos directamente en el espacio de Fourier. Dado que una imagen de RM ideal puede describirse solo con valores reales (todos los vectores de magnetización excitados tienen fase 0), solo tiene que adquirir la mitad del espacio de datos, lo que le ahorra la mitad del tiempo de imagen. Por supuesto, las imágenes de resonancia magnética no están completamente valoradas debido a las limitaciones de la realidad ... pero con algunos trucos aún puede utilizar esta técnica de manera ventajosa.

M529
fuente
2
Me gustó la forma simple de indicar la misma respuesta que proporcionó ThP. Y gracias por la información sobre MRI. No sabía sobre eso.
Falló el científico