Registro de imagen invariable a escala basado en DFT Log-Polar

Estoy tratando de hacer el registro de imágenes usando la correlación de fase como se describe en el documento de Reddy Chatterji . En mi caso, las imágenes pueden ser escaladas y traducidas entre sí.

El algoritmo para encontrar la escala relativa, según tengo entendido, es (ver: el diagrama de flujo del documento ):

F1 = DFT(I1)
F2 = DFT(I2)
H1 = Highpass(F1)
H2 = Highpass(F2)
L1 = LogPolar(Magnitude(H1))
L2 = LogPolar(Magnitude(H2))
PC = PhaseCorrelate(L1,L2)
PM = norm(PC)
R = IDFT(PhaseCorr/PM)
P = Peak(R)
Scale = LogBase^P[1]

La escala me da valores aparentemente sin sentido (muy diferentes de una imagen a otra y nunca correctos).

Pero ignorando la escala, el mismo enfoque de correlación de fase funciona bien para la traducción; y entonces sospecho que tengo un problema con mi transformación log-polar. Aquí hay un ejemplo, donde resolví la traducción: la imagen de la izquierda es la original y la derecha se ha recortado y traducido, la solución se muestra en la parte superior del original:

Para el log-polar transformar, I primera transformo en el espacio dondees la imagen original,es el radio de la imagen (medio ancho) yes el número de muestras en ladirección. I a continuación muestra a partir de esta transformar en espacio de registro , dondecomo se describe en

\hat{yo} (ρ, θ) = yo (r + ρ \cos (\frac{2 π θ}{{norte}_{θ}}), r - ρ pecado (\frac{2 π θ}{{norte}_{θ}}))

$\hat{I}(\rho,\theta) = I\left(r+\rho\cos\left(\frac{2\pi\theta}{N_{\theta}}\right),r-\rho\sin\left(\frac{2\pi\theta}{N_{\theta}}\right) \right)$

I

$I$

r

$r$

N_{θ}

$N_{\theta}$

θ

$\theta$

{\hat{yo}}_{l o sol} (ρ, θ) = \hat{yo} ({Iniciar sesión}_{si} (ρ), θ)

$\hat{I}_{log}(\rho,\theta) = \hat{I}\left(\log_{b}(\rho),\theta\right)$

b = (2 r)^{- N_{ρ}}

$b = (2r)^{-N_{\rho}}$ 1 para que abarque todo el espacio polar.

$\rho=\theta=256$

Por último, esto muestra la transformación real por la que pasan las imágenes antes del paso de correlación de fase (la parte superior es el filtro de paso alto de magnitud DFT, la parte inferior es la del espacio polar logarítmico):

Estoy usando OpenCV, que tiene los métodos LogPolar y PhaseCorrelate. Mientras que PhaseCorrelate, como mi implementación manual, me da la respuesta correcta para la traducción, es incorrecta en escala. Dado que el uso de OpenCV LogPolar o el mío no afecta la corrección, me falta algo.

Cualquier ayuda sería apreciada.

image-processing dft opencv Drew Cummins
fuente

¿Te diste cuenta de lo que estaba mal?

Mr.WorshipMe

@ Mr.WorshipMe Lamentablemente no.

Drew Cummins el

@Drew Cummins, supongo que se debió a la imagen de prueba que usaste, porque hubo una transición nítida desde el fondo. ¿Qué tal otras imágenes de prueba? Además, a partir de la última figura, hubo diferencias aparentes entre las dos magnitudes, por lo que es mejor realizar un preprocesamiento de ventanas adecuado antes del DFT.

lxg

Un par de días antes encontré ese documento y he estado tratando de implementar el algoritmo sin éxito. Me preguntaba si podrías compartir tu implementación con un principiante :)

Alexis España

Registro de imagen invariable a escala basado en DFT Log-Polar

Respuestas: