Reconstrucción de imagen: fase vs. magnitud

La figura 1. (c) muestra la imagen de prueba reconstruida solo a partir del espectro MAGNITUD. Podemos decir que los valores de intensidad de los píxeles de frecuencia BAJA son comparativamente más que los píxeles de frecuencia ALTA.

La figura 1. (d) muestra la imagen de prueba reconstruida solo del espectro de FASE. Podemos decir que los valores de intensidad de los píxeles de ALTA frecuencia (bordes, líneas) son comparativamente más que los píxeles de BAJA frecuencia.

¿Por qué esta contradicción mágica de cambio de intensidad (o intercambio) está presente entre la imagen de prueba reconstruida solo del espectro MAGNITUD y la imagen de prueba reconstruida solo del espectro de FASE, que cuando se combinan forman la imagen de prueba original?

clc;
clear all;
close all;
i1=imread('C:\Users\Admin\Desktop\rough\Capture1.png');
i1=rgb2gray(i1);

f1=fftn(i1);
mag1=abs(f1);
s=log(1+fftshift(f1));
phase1=angle(f1);

r1=ifftshift(ifftn(mag1));
r2=ifftn(exp(1i*phase1));
figure,imshow(i1);
figure,imshow(s,[]);
figure,imshow(uint8(r1));
figure,imshow(r2,[]);
r2=histeq(r2);
r3=histeq(uint8(r2));     
figure,imshow(r2);
figure,imshow(r3);

La figura 1. (c) muestra la imagen de prueba reconstruida solo a partir del espectro MAGNITUD. Podemos decir que los valores de intensidad de los píxeles de frecuencia BAJA son comparativamente más que los píxeles de frecuencia ALTA.

En realidad, esto no es correcto. Los valores de fase determinan el cambio en los componentes sinusoides de la imagen. Con la fase cero, todas las sinusoides se centran en la misma ubicación y se obtiene una imagen simétrica cuya estructura no tiene correlación real con la imagen original. Estar centrado en la misma ubicación significa que las sinusoides son máximas en esa ubicación, y es por eso que hay un gran parche blanco en el centro de la Figura 1.c.

La reconstrucción solo de fase conserva las características debido al principio de congruencia de fase . En la ubicación de bordes y líneas, la mayoría de los componentes sinusoides tienen la misma fase. Consulte http://homepages.inf.ed.ac.uk/rbf/CVonline/LOCAL_COPIES/OWENS/LECT7/node2.html. Esto solo puede usarse para detectar líneas y bordes, http: //www.csse.uwa. edu.au/~pk/research/pkpapers/phasecorners.pdf , sin importar su magnitud. Entonces puede ver que la información de la fase es más importante.

Cambiar la magnitud de los diversos sinusoides componentes cambia la forma de la entidad. Cuando realiza una reconstrucción solo de fase, establece todas las magnitudes en una, lo que cambia la forma de las características, pero no su ubicación. En muchas imágenes, los componentes de baja frecuencia tienen una magnitud mayor que los componentes de alta frecuencia, por lo que la reconstrucción de fase única parece un filtro de paso alto.

En resumen, la fase contiene la información sobre la ubicación de las características.

No puede agregar las imágenes solo de fase y solo de magnitud para obtener el original. Puede multiplicarlos en el dominio de Fourier y transformarlos de nuevo para obtener el original.

En su línea mag1=abs(f1); , está dejando la intensidad total de la imagen sin cambios (pruebe esto resumiendo las intensidades en todos los píxeles). Rechazar la información de fase en el espacio de Fourier solo conduce a una redistribución espacial de la intensidad en el espacio real, de modo que r1 tendrá la misma insistencia total que i1.

En su línea phase1=angle(f1); está normalizando las amplitudes de cada píxel (en el espacio de Fourier) a 1, por lo que la intensidad total de la imagen cambiará. Como la fase lleva una gran parte de la información espacial de la imagen, se conservan las principales características de la imagen.