Articulo de referencia

Diferencia de gaussianas

En el campo de la imagen , la diferencia de gaussianas ( DoG ) es un algoritmo de mejora de características que consiste en la resta de una versión desenfocada gaussiana de una ...

En el campo de la imagen , la diferencia de gaussianas ( DoG ) es un algoritmo de mejora de características que consiste en la resta de una versión desenfocada gaussiana de una imagen original de otra versión menos desenfocada de la misma. En el caso simple de imágenes en escala de grises , las imágenes desenfocadas se obtienen mediante la convolución de las imágenes originales en escala de grises con núcleos gaussianos de diferente ancho (desviación estándar). El desenfoque de una imagen mediante un núcleo gaussiano suprime únicamente la información espacial de alta frecuencia . La resta de una imagen de la otra preserva la información espacial que se encuentra dentro del rango de frecuencias preservadas en las dos imágenes desenfocadas. Por lo tanto, la DoG es un filtro de paso de banda espacial que atenúa las frecuencias en la imagen original en escala de grises que están lejos del centro de la banda. [ 1 ]

Formulación

Comparación entre la diferencia de gaussianas y la ondícula de sombrero mexicano.

DejarΦt:RnorteR{\displaystyle \Phi _{t}:\mathbb {R} ^{n}\rightarrow \mathbb {R} }denotamos la función gaussiana radialΦt(incógnita)=norte(incógnita|0,t){\displaystyle \Phi _{t}(x)={\mathcal {N}}(x|0,t)}con significa0{\displaystyle 0}y varianzat{\displaystyle t}, es decir, la función gaussiana multivariadaΦt(incógnita)=norte(incógnita|0,tI){\displaystyle \Phi _ {t}(x)={\mathcal {N}}(x|0,tI)}con significa0{\displaystyle 0}y covarianzatI{\displaystyle tI}. Más explícitamente, tenemos

Φt(incógnita)=1(2πt)norte/2miincógnita22t.{\displaystyle \Phi _{t}(x)={\frac {1}{(2\pi {}t)^{n/2}}}e^{-{\frac {\|x\|^{2}}{2t}}}.}

La diferencia de gaussianas con varianzast1<t2{\displaystyle t_{1}<t_{2}}es la función kernel

Kt1,t2=Φt1Φt2{\displaystyle K_{t_{1},t_{2}}=\Phi _{t_{1}}-\Phi _{t_{2}}}

obtenido al restar la gaussiana de mayor varianza de la gaussiana de menor varianza. El operador de diferencia de gaussianas es el operador convolucional asociado con esta función kernel. Por lo tanto, dada una imagen en escala de grises n -dimensionalI:RnorteR{\displaystyle I:\mathbb {R} ^{n}\rightarrow \mathbb {R} }, la diferencia de gaussianas de la imagenI{\displaystyle I}es la imagen n -dimensional

IKt1,t2=I(Φt1Φt2)=IΦt1IΦt2.{\displaystyle I*K_{t_{1},t_{2}}=I*(\Phi _{t_{1}}-\Phi _{t_{2}})=I*\Phi _{t_{1}}-I*\Phi _{t_{2}}.}

Dado que la convolución es bilineal, convolucionar contra la diferencia de gaussianas equivale a aplicar dos desenfoques gaussianos diferentes y luego calcular la diferencia. En la práctica, esto es más rápido porque el desenfoque gaussiano es un filtro separable .

Aproximación

La diferencia de gaussianas puede pensarse como una aproximación de la función núcleo de sombrero mexicano utilizada para el laplaciano del operador gaussiano. La observación clave es que la familia de gaussianasΦt{\displaystyle \Phi _{t}}es la solución fundamental de la ecuación del calor

tΦt(incógnita)=12ΔΦt(incógnita).{\displaystyle \partial _{t}\Phi _{t}(x)={\frac {1}{2}}\Delta \Phi _{t}(x).}

El lado izquierdo se puede aproximar mediante el cociente de diferencias.

Φt+δt(incógnita)Φt(incógnita)δt=1δtKt+δt,t(incógnita).{\displaystyle {\frac {\Phi _{t+\delta {t}}(x)-\Phi _{t}(x)}{\delta {t}}}={\frac {1}{\delta {t}}}K_{t+\delta {t},t}(x).}

Mientras tanto, el lado derecho es precisamente el laplaciano de la función gaussiana. Nótese que el laplaciano de la gaussiana se puede utilizar como filtro para producir un desenfoque gaussiano del laplaciano de la imagen porqueIΔΦt=ΔIΦt{\displaystyle I*\Delta \Phi _{t}=\Delta {I}*\Phi _{t}}por propiedades estándar de convolución. La relación entre el operador de diferencia de gaussianas y el laplaciano del operador gaussiano se explica con más detalle en el Apéndice A de Lindeberg (2015). [ 2 ]

Detalles y solicitudes

Ejemplo antes de la diferencia de gaussianas
Después del filtrado de diferencia de gaussianas en blanco y negro

Como algoritmo de mejora de características , la diferencia de gaussianas se puede utilizar para aumentar la visibilidad de los bordes y otros detalles presentes en una imagen digital. Una amplia variedad de filtros alternativos de nitidez de bordes funcionan mejorando los detalles de alta frecuencia, pero debido a que el ruido aleatorio también tiene una alta frecuencia espacial, muchos de estos filtros tienden a aumentar el ruido, lo que puede ser un artefacto indeseable. El algoritmo de diferencia de gaussianas elimina los detalles de alta frecuencia que a menudo incluyen ruido aleatorio, lo que hace que este enfoque sea uno de los más adecuados para procesar imágenes con un alto grado de ruido. Una desventaja importante de la aplicación del algoritmo es una reducción inherente en el contraste general de la imagen producida por la operación. [ 1 ]

Cuando se utiliza para mejorar imágenes, el algoritmo de diferencia de gaussianas se aplica normalmente cuando la relación de tamaño del núcleo (2) al núcleo (1) es de 4:1 o 5:1. En las imágenes de ejemplo, los tamaños de los núcleos gaussianos empleados para suavizar la imagen de muestra fueron de 10 píxeles y 5 píxeles.

El algoritmo también puede utilizarse para obtener una aproximación del Laplaciano de Gauss cuando la relación entre el tamaño 2 y el tamaño 1 es aproximadamente igual a 1,6. [ 3 ] El Laplaciano de Gauss es útil para detectar bordes que aparecen en diversas escalas de imagen o grados de enfoque de la imagen. Los valores exactos de los tamaños de los dos núcleos que se utilizan para aproximar el Laplaciano de Gauss determinarán la escala de la imagen de diferencia, que puede aparecer borrosa como resultado.

Las diferencias de gaussianas también se han utilizado para la detección de blobs en la transformación de características invariantes a la escala (SIFT). De hecho, la DoG, como diferencia de dos distribuciones normales multivariadas, siempre tiene una suma nula total y su convolución con una señal uniforme no genera respuesta. Se aproxima bien a la segunda derivada de una gaussiana ( laplaciano de gaussiana ) con K~1,6 y a los campos receptivos de las células ganglionares de la retina con K~5. Puede utilizarse fácilmente en esquemas recursivos y se emplea como operador en algoritmos en tiempo real para la detección de blobs y la selección automática de escala.

Más información

In its operation, the difference of Gaussians algorithm is believed to mimic how neural processing in the retina of the eye extracts details from images destined for transmission to the brain.[4][5][6]

See also

References

  1. 12"Molecular Expressions Microscopy Primer: Digital Image Processing – Difference of Gaussians Edge Enhancement Algorithm", Olympus America Inc., and Florida State University Michael W. Davidson, Mortimer Abramowitz
  2. Lindeberg, Tony (2015). "Image Matching Using Generalized Scale-Space Interest Points". Journal of Mathematical Imaging and Vision. 52: 3–36. doi:10.1007/s10851-014-0541-0. S2CID 254657377.
  3. D. Marr; E. Hildreth (29 February 1980). "Theory of Edge Detection". Proceedings of the Royal Society of London. Series B, Biological Sciences. 207 (1167): 215–217. Bibcode:1980RSPSB.207..187M. doi:10.1098/rspb.1980.0020. JSTOR 35407. PMID 6102765. S2CID 2150419. — A difference of Gaussians of any scale is an approximation to the laplacian of the Gaussian (see the entry for difference of Gaussians under Blob detection). However, Marr and Hildreth recommend the ratio of 1.6 because of design considerations balancing bandwidth and sensitivity. The url for this reference may only make the first page and abstract of the article available depending on if you are connecting through an academic institution or not.
  4. Christina Enroth-Cugell; J. G. Robson (1966). "The Contrast Sensitivity of Retinal Ganglion Cells of the Cat". Journal of Physiology. 187 (3): 517–552. doi:10.1113/jphysiol.1966.sp008107. PMC 1395960. PMID 16783910.
  5. Matthew J. McMahon; Orin S. Packer; Dennis M. Dacey (April 14, 2004). "The Classical Receptive Field Surround of Primate Parasol Ganglion Cells Is Mediated Primarily by a Non-GABAergic Pathway"(PDF). Journal of Neuroscience. 24 (15): 3736–3745. doi:10.1523/JNEUROSCI.5252-03.2004. PMC 6729348. PMID 15084653.
  6. Young, Richard (1987). "The Gaussian derivative model for spatial vision: I. Retinal mechanisms". Spatial Vision. 2 (4): 273–293(21). doi:10.1163/156856887X00222. PMID 3154952.

Further reading

  • Notes by Melisa Durmuş on Edge Detection and Gaussian related mathematics from the University of Edinburgh.