Articulo de referencia

Momento de la imagen

En el procesamiento de imágenes , la visión por computadora y campos relacionados, un momento de imagen es un promedio ponderado particular ( momento ) de las intensidades de lo...

En el procesamiento de imágenes , la visión por computadora y campos relacionados, un momento de imagen es un promedio ponderado particular ( momento ) de las intensidades de los píxeles de la imagen, o una función de dichos momentos, generalmente elegido para tener alguna propiedad o interpretación atractiva.

Los momentos de imagen son útiles para describir objetos después de la segmentación . Las propiedades simples de la imagen que se encuentran a través de los momentos de imagen incluyen el área (o intensidad total), su centroide e información sobre su orientación .

Momentos crudos

Para una función continua bidimensional f ( x , y ), el momento (a veces llamado "momento bruto") de orden ( p + q ) se define como

METROpagq=incógnitapagyqF(incógnita,y)dincógnitady{\displaystyle M_{pq}=\int \limits _{-\infty }^{\infty }\int \limits _{-\infty }^{\infty }x^{p}y^{q}f(x,y)\,dx\,dy}

para p , q = 0,1,2,... Adaptando esto a una imagen escalar ( escala de grises ) con intensidades de píxeles I ( x , y ), los momentos de imagen sin procesar M ij se calculan mediante

METROij=incógnitayincógnitaiyjI(incógnita,y){\displaystyle M_{ij}=\sum _{x}\sum _{y}x^{i}y^{j}I(x,y)\,\!}

En algunos casos, esto se puede calcular considerando la imagen como una función de densidad de probabilidad , es decir , dividiendo lo anterior por

incógnitayI(incógnita,y){\displaystyle \sum _{x}\sum _{y}I(x,y)\,\!}

Un teorema de unicidad [ 1 ] establece que si f ( x , y ) es continua a trozos y tiene valores distintos de cero solo en una parte finita del plano xy , existen momentos de todos los órdenes, y la secuencia de momentos ( M pq ) está determinada unívocamente por f ( x , y ). [ 2 ] Recíprocamente, ( M pq ) determina unívocamente f ( x , y ). En la práctica, la imagen se resume con funciones de algunos momentos de orden inferior.

Ejemplos

Las propiedades simples de la imagen derivadas a través de momentos brutos incluyen:

  • Área (para imágenes binarias) o suma de niveles de gris (para imágenes en escala de grises):METRO00{\displaystyle M_{00}}
  • Centroide:{incógnita¯, y¯}={METRO10METRO00,METRO01METRO00}{\displaystyle \{{\bar {x}},\ {\bar {y}}\}=\left\{{\frac {M_{10}}{M_{00}}},{\frac {M_{01}}{M_{00}}}\right\}}

Momentos centrales

Los momentos centrales se definen como

μpagq=(incógnitaincógnita¯)pag(yy¯)qF(incógnita,y)dincógnitady{\displaystyle \mu _{pq}=\int \limits _{-\infty }^{\infty }\int \limits _{-\infty }^{\infty }(x-{\bar {x}})^{p}(y-{\bar {y}})^{q}f(x,y)\,dx\,dy}

dóndeincógnita¯=METRO10METRO00{\displaystyle {\bar {x}}={\frac {M_{10}}{M_{00}}}}yy¯=METRO01METRO00{\displaystyle {\bar {y}}={\frac {M_{01}}{M_{00}}}}son los componentes del centroide .

Si ƒ ( x , y ) es una imagen digital, entonces la ecuación anterior se convierte en: 

μpagq=incógnitay(incógnitaincógnita¯)pag(yy¯)qF(incógnita,y){\displaystyle \mu _{pq}=\sum _{x}\sum _{y}(x-{\bar {x}})^{p}(y-{\bar {y}})^{q}f(x,y)}

Los momentos centrales de orden hasta 3 son:

μ00=METRO00,μ01=0,μ10=0,μ11=METRO11incógnita¯METRO01=METRO11y¯METRO10,μ20=METRO20incógnita¯METRO10,μ02=METRO02y¯METRO01,μ21=METRO212incógnita¯METRO11y¯METRO20+2incógnita¯2METRO01,μ12=METRO122y¯METRO11incógnita¯METRO02+2y¯2METRO10,μ30=METRO303incógnita¯METRO20+2incógnita¯2METRO10,μ03=METRO033y¯METRO02+2y¯2METRO01.{\displaystyle {\begin{aligned}\mu _{00}&=M_{00},&\mu _{01}&=0,\\\mu _{10}&=0,&\mu _{11}&=M_{11}-{\bar {x}}M_{01}=M_{11}-{\bar {y}}M_{10},\\\mu _{20}&=M_{20}-{\bar {x}}M_{10},&\mu _{02}&=M_{02}-{\bar {y}}M_{01},\\\mu _{21}&=M_{21}-2{\bar {x}}M_{11}-{\bar {y}}M_{20}+2{\bar {x}}^{2}M_{01},&\mu _{12}&=M_{12}-2{\bar {y}}M_{11}-{\bar {x}}M_{02}+2{\bar {y}}^{2}M_{10},\\\mu _{30}&=M_{30}-3{\bar {x}}M_{20}+2{\bar {x}}^{2}M_{10},&\mu _{03}&=M_{03}-3{\bar {y}}M_{02}+2{\bar {y}}^{2}M_{01}.\end{aligned}}}

Se puede demostrar que:

μpagq=metropagnorteq(pagmetro)(qnorte)(incógnita¯)(pagmetro)(y¯)(qnorte)METROmetronorte{\displaystyle \mu _{pq}=\sum _{m}^{p}\sum _{n}^{q}{p \choose m}{q \choose n}(-{\bar {x}})^{(pm)}(-{\bar {y}})^{(qn)}M_{mn}}

Los momentos centrales son invariantes traslacionales .

Ejemplos

La información sobre la orientación de la imagen se puede obtener utilizando primero los momentos centrales de segundo orden para construir una matriz de covarianza .

μ20=μ20/μ00=METRO20/METRO00incógnita¯2μ02=μ02/μ00=METRO02/METRO00y¯2μ11=μ11/μ00=METRO11/METRO00incógnita¯y¯{\displaystyle {\begin{aligned}\mu '_{20}&=\mu _{20}/\mu _{00}=M_{20}/M_{00}-{\bar {x}}^{2}\\\mu '_{02}&=\mu _{02}/\mu _{00}=M_{02}/M_{00}-{\bar {y}}^{2}\\\mu '_{11}&=\mu _{11}/\mu _{00}=M_{11}/M_{00}-{\bar {x}}{\bar {y}}\end{aligned}}}

La matriz de covarianza de la imagenI(incógnita,y){\displaystyle I(x,y)}ahora es

cobertura[I(incógnita,y)]=[μ20μ11μ11μ02].{\displaystyle \operatorname {cov} [I(x,y)]={\begin{bmatrix}\mu '_{20}&\mu '_{11}\\\mu '_{11}&\mu '_{02}\end{bmatrix}}.}

Los autovectores de esta matriz corresponden a los ejes mayor y menor de la intensidad de la imagen, por lo que la orientación se puede extraer del ángulo del autovector asociado al autovalor más grande hacia el eje más cercano a dicho autovector. Se puede demostrar que este ángulo Θ viene dado por la siguiente fórmula:

Θ=12arctan(2μ11μ20μ02){\displaystyle \Theta ={\frac {1}{2}}\arctan \left({\frac {2\mu '_{11}}{\mu '_{20}-\mu '_{02}}}\right)}

La fórmula anterior es válida siempre que:

μ20μ020{\displaystyle \mu '_{20}-\mu '_{02}\neq 0}

Se puede demostrar fácilmente que los valores propios de la matriz de covarianza son:

λi=μ20+μ022±4μ112+(μ20μ02)22,{\displaystyle \lambda _{i}={\frac {\mu '_{20}+\mu '_{02}}{2}}\pm {\frac {\sqrt {4{\mu '}_{11}^{2}+({\mu '}_{20}-{\mu '}_{02})^{2}}}{2}},}

y son proporcionales al cuadrado de la longitud de los ejes del vector propio. La diferencia relativa en magnitud de los valores propios es, por lo tanto, una indicación de la excentricidad de la imagen, o cuán alargada está. La excentricidad es

1λ2λ1.{\displaystyle {\sqrt {1-{\frac {\lambda _{2}}{\lambda _{1}}}}}.}

invariantes de momento

Los momentos son bien conocidos por su aplicación en el análisis de imágenes, ya que pueden utilizarse para derivar invariantes con respecto a clases de transformación específicas.

El término momentos invariantes se usa con frecuencia de forma incorrecta en este contexto. Sin embargo, si bien los invariantes de momento son invariantes que se forman a partir de momentos, los únicos momentos que son invariantes en sí mismos son los momentos centrales.

Cabe destacar que los invariantes que se detallan a continuación son exactamente invariantes solo en el dominio continuo. En un dominio discreto, ni el escalado ni la rotación están bien definidos: una imagen discreta transformada de esta manera suele ser una aproximación, y la transformación no es reversible. Por lo tanto, estos invariantes son solo aproximadamente invariantes al describir una forma en una imagen discreta.

invariantes de traslación

Los momentos centrales μ i j de cualquier orden son, por construcción, invariantes con respecto a las traslaciones .

invariantes de escala

Los invariantes η i j con respecto a la traslación y la escala se pueden construir a partir de los momentos centrales dividiendo por un momento central de orden cero debidamente escalado:

ηij=μijμ00(1+i+j2){\displaystyle \eta _{ij}={\frac {\mu _{ij}}{\mu _{00}^{\left(1+{\frac {i+j}{2}}\right)}}}\,\!}

donde i + j ≥ 2. Nótese que la invariancia traslacional se deduce directamente utilizando únicamente momentos centrales.

Invariantes de rotación

Como se muestra en el trabajo de Hu, [ 3 ] [ 4 ] se pueden construir invariantes con respecto a la traslación , la escala y la rotación :

I1=η20+η02{\displaystyle I_{1}=\eta _{20}+\eta _{02}}

I2=(η20η02)2+4η112{\displaystyle I_{2}=(\eta _{20}-\eta _{02})^{2}+4\eta _{11}^{2}}

I3=(η303η12)2+(3η21η03)2{\displaystyle I_{3}=(\eta _{30}-3\eta _{12})^{2}+(3\eta _{21}-\eta _{03})^{2}}

I4=(η30+η12)2+(η21+η03)2{\displaystyle I_{4}=(\eta _{30}+\eta _{12})^{2}+(\eta _{21}+\eta _{03})^{2}}

I5=(η303η12)(η30+η12)[(η30+η12)23(η21+η03)2]+(3η21η03)(η21+η03)[3(η30+η12)2(η21+η03)2]{\displaystyle I_{5}=(\eta _{30}-3\eta _{12})(\eta _{30}+\eta _{12})[(\eta _{30}+\eta _{12})^{2}-3(\eta _{21}+\eta _{03})^{2}]+(3\eta _{21}-\eta _{03})(\eta _{21}+\eta _{03})[3(\eta _{30}+\eta _{12})^{2}-(\eta _{21}+\eta _{03})^{2}]}

I6=(η20η02)[(η30+η12)2(η21+η03)2]+4η11(η30+η12)(η21+η03){\displaystyle I_{6}=(\eta _{20}-\eta _{02})[(\eta _{30}+\eta _{12})^{2}-(\eta _{21}+\eta _{03})^{2}]+4\eta _{11}(\eta _{30}+\eta _{12})(\eta _{21}+\eta _{03})}

I7=(3η21η03)(η30+η12)[(η30+η12)23(η21+η03)2](η303η12)(η21+η03)[3(η30+η12)2(η21+η03)2].{\displaystyle I_{7}=(3\eta _{21}-\eta _{03})(\eta _{30}+\eta _{12})[(\eta _{30}+\eta _{12})^{2}-3(\eta _{21}+\eta _{03})^{2}]-(\eta _{30}-3\eta _{12})(\eta _{21}+\eta _{03})[3(\eta _{30}+\eta _{12})^{2}-(\eta _{21}+\eta _{03})^{2}].}

Estos son bien conocidos como invariantes de momento de Hu .

La primera, I 1 , es análoga al momento de inercia alrededor del centroide de la imagen, donde las intensidades de los píxeles son análogas a la densidad física. Las seis primeras, I 1 ... I 6 , son simétricas respecto a la reflexión, es decir, no cambian si la imagen se convierte en una imagen especular. La última, I 7 , es antisimétrica respecto a la reflexión (cambia de signo al reflejarse), lo que le permite distinguir imágenes especulares de imágenes que, de otro modo, serían idénticas.

J. Flusser propuso una teoría general sobre la derivación de conjuntos completos e independientes de invariantes de momento de rotación. [ 5 ] Demostró que el conjunto tradicional de invariantes de momento de Hu no es ni independiente ni completo. I 3 no es muy útil ya que depende de los demás (I3=(I52+I72)/I43{\displaystyle I_{3}=(I_{5}^{2}+I_{7}^{2})/I_{4}^{3}}). En el conjunto original de Hu falta un invariante de momento independiente de tercer orden:

I8=η11[(η30+η12)2(η03+η21)2](η20η02)(η30+η12)(η03+η21){\displaystyle I_{8}=\eta _{11}[(\eta _{30}+\eta _{12})^{2}-(\eta _{03}+\eta _{21})^{2}]-(\eta _{20}-\eta _{02})(\eta _{30}+\eta _{12})(\eta _{03}+\eta _{21})}

Al igual que I 7 , I 8 también es antisimétrico de reflexión.

Más tarde, J. Flusser y T. Suk [ 6 ] especializaron la teoría para el caso de formas N-rotacionalmente simétricas.

Aplicaciones

Zhang et al. aplicaron invariantes de momento de Hu para resolver el problema de detección patológica del cerebro (PBD). [ 7 ] Doerr y Florence utilizaron información de la orientación del objeto relacionada con los momentos centrales de segundo orden para extraer eficazmente secciones transversales de objetos invariantes a la traslación y la rotación a partir de datos de imágenes de tomografía de micro-rayos X. [ 8 ]

DA Hoeltzel y Wei-Hua Chieng utilizaron el invariante de momento de Hu para realizar un análisis en un mecanismo de cuatro barras parametrizado dimensionalmente, lo que produjo 15 grupos distintos de curvas de acoplamiento (patrones) de un total de 356 curvas de acoplamiento generadas. [ 9 ]

  • Análisis de imágenes binarias , Universidad de Edimburgo
  • Momentos estadísticos , Universidad de Edimburgo
  • Momentos variantes , página de Percepción Artificial y Visión por Computadora (código fuente en Matlab y Python)
  • Vídeo de presentación de Hu Moments en YouTube
  • Implementación básica de esta página, Jupyter y Python.

Referencias

  1. Hu, Ming-Kuei (1962). "Reconocimiento de patrones visuales mediante invariantes de momento". IRE Transactions on Information Theory . 8 (2): 179– 187. doi : 10.1109/TIT.1962.1057692 . ISSN 2168-2712 . 
  2. González, Rafael C.; Woods, Richard E. (2001). Procesamiento de imágenes digitales . Prentice Hall. pág. 672. ISBN  0-201-18075-8.
  3. MK Hu, "Reconocimiento de patrones visuales mediante invariantes de momento", IRE Trans. Info. Theory, vol. IT-8, pp.179 187, 1962
  4. http://docs.opencv.org/modules/imgproc/doc/structural_analysis_and_shape_descriptors.html?highlight=cvmatchshapes#humoments Método OpenCV de Hu Moments
  5. J. Flusser: " Sobre la independencia de los invariantes del momento de rotación ", Pattern Recognition, vol. 33, pp. 1405 1410, 2000.
  6. J. Flusser y T. Suk, " Invariantes del momento de rotación para el reconocimiento de objetos simétricos ", IEEE Trans. Image Proc., vol. 15, pp. 3784 3790, 2006.
  7. Zhang, Y. (2015). "Detección patológica del cerebro basada en la entropía wavelet y los invariantes del momento de Hu" . Bio-Medical Materials and Engineering . 26 : 1283–1290 . doi : 10.3233/BME-151426 . PMID 26405888 . 
  8. Doerr, Frederik; Florence, Alastair (2020). "Una metodología de análisis de imágenes micro-XRT y aprendizaje automático para la caracterización de formulaciones de cápsulas multiparticuladas" . International Journal of Pharmaceutics: X. 2 100041. doi : 10.1016 /j.ijpx.2020.100041 . PMC 6997304. PMID 32025658 .  
  9. Hoeltzel, DA; Chieng, Wei-Hua (1990). "Síntesis de coincidencia de patrones como un enfoque automatizado para el diseño de mecanismos" . Journal of Mechanical Design . 112 (2): 190-199. doi : 10.1115/1.2912592 .