La figura 1.(c) muestra la imagen de prueba reconstruida únicamente a partir del espectro de MAGNITUD. Podemos decir que los valores de intensidad de los píxeles de BAJA frecuencia son comparativamente más altos que los píxeles de ALTA frecuencia.
La figura 1.(d) muestra la imagen de prueba reconstruida únicamente a partir 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 altos que los píxeles de BAJA frecuencia.
¿Por qué está presente esta contradicción mágica de cambio de intensidad (o intercambio) entre la imagen de prueba reconstruida solo a partir del espectro de MAGNITUD y la imagen de prueba reconstruida solo a partir del espectro de FASE, que cuando se combinan forman la imagen de prueba original?
La naturaleza de las imágenes C y D se explica por la imagen de potencia espectral, B.
Puede ver en la imagen de potencia espectral que la mayor parte de la energía se concentra en bajas frecuencias (la parte central más brillante de la imagen). Las frecuencias más altas (que se alejan del centro) tienen menos energía porque hay menos brillo en esta parte de la imagen.
Este tipo de espectro de potencia es típico en imágenes naturales . La caída de energía con la frecuencia tiende a ajustarse bien a una función exponencial.
Dado esto, la imagen C (abajo a la izquierda) tiene sentido. Esperarías ver formas de manchas que son predominantemente de baja frecuencia, porque ahí es donde está el poder en la imagen. Su simetría se debe a que las fases están todas puestas a cero.
Ha proporcionado la función utilizada para reconstruir la función de solo fase. Si elimina las funciones de suma y base, se queda con:
F_phase_only(u,v) = F(u,v) / |F(u,v)|
F(u,v) es un número complejo con fase y magnitud. Está dividiendo por la magnitud para eliminarlo, para que sea 1 para todos los componentes. Entonces, F_phase_only(u,v) solo tiene las fases. La magnitud de F_phase_only(u,v) = 1, para todo u,v.
Ahora, mire la imagen espectral nuevamente (arriba a la derecha). Esto es realmente solo un gráfico de |F(u,v)| – es una imagen de magnitud única.
Podemos ver que para u,v pequeños, las bajas frecuencias que están en el medio de la imagen, |F(u,v)| es grande, porque el medio es brillante.
Para u,v más grandes, las frecuencias más altas que están lejos del medio, |F(u,v)| es relativamente pequeño, porque estas partes de la imagen son más oscuras.
Entonces, volviendo a F_phase_only(u,v), lo reescribiré muy sutilmente.
F_phase_only(u,v) = F(u,v) . BOOST(u,v).
¿Qué es impulsar?
BOOST(u,v) = 1 / |F(u,v)|
¿Qué podemos determinar sobre la naturaleza de BOOST?
Para pequeñas u,v las frecuencias bajas , BOOST(u,v) debe ser un número pequeño. Debe ser pequeño, porque BOOST es el recíproco de la magnitud. Sabemos por el espectro que la magnitud es grande para bajas frecuencias. Por lo tanto, BOOST(u,v) debe ser pequeño y así disminuir las bajas frecuencias .
Para u,v más grandes las frecuencias altas , BOOST(u,v) debe ser un número grande. Debe ser grande, porque BOOST es el recíproco de la magnitud. Sabemos por el espectro que la magnitud es pequeña para frecuencias altas . Por lo tanto, BOOST(u,v) debe ser grande y, por lo tanto, mejora las frecuencias altas .
Entonces, la imagen de "solo fase" da a todos los componentes del dominio de frecuencia la misma magnitud. Para un componente individual tomado en relación con los demás, esto enfatiza enormemente los componentes de frecuencia más alta y disminuye en gran medida los componentes de baja frecuencia .
Entonces, si bien lo ha llamado una "imagen de solo fase", en realidad es una imagen en la que los términos de alta frecuencia se mejoran enormemente a expensas de los términos de baja frecuencia.
Por esta razón, la imagen resultante muestra las características de borde que existen en las frecuencias más altas ampliadas.
Teniendo en cuenta el papel que mencioné anteriormente , podría preparar una imagen de fase más justa.
F_fairer_phase(u,v) = Expected(u,v) . F(u,v) / |F(u,v)|
¿Dónde Expected(u,v)
está la potencia espectral estadísticamente esperada en u,v? Se deja como ejercicio :-)
Alternativamente, podría derivar Expected(u,v)
promediando todo |F(w,x)|
donde sqrt(w.w + x.x) = sqrt(u.u, v.v)
.
benjohn