En el procesamiento de imágenes digitales y la visión artificial , la segmentación de imágenes es el proceso de dividir una imagen digital en múltiples segmentos de imagen , también conocidos como regiones de imagen u objetos de imagen ( conjuntos de píxeles ). El objetivo de la segmentación es simplificar o cambiar la representación de una imagen para que sea más significativa y más fácil de analizar. [1] [2] La segmentación de imágenes se utiliza normalmente para localizar objetos y límites (líneas, curvas, etc.) en imágenes. Más precisamente, la segmentación de imágenes es el proceso de asignar una etiqueta a cada píxel de una imagen de forma que los píxeles con la misma etiqueta compartan determinadas características.
El resultado de la segmentación de imágenes es un conjunto de segmentos que cubren colectivamente toda la imagen, o un conjunto de contornos extraídos de la imagen (ver detección de bordes ). Cada uno de los píxeles de una región es similar con respecto a alguna característica o propiedad calculada, [3] como el color , la intensidad o la textura . Las regiones adyacentes son significativamente diferentes con respecto a las mismas características. [1] Cuando se aplica a una pila de imágenes, algo típico en las imágenes médicas , los contornos resultantes después de la segmentación de imágenes se pueden utilizar para crear reconstrucciones 3D con la ayuda de algoritmos de reconstrucción geométrica como los cubos de marcha . [4]
Algunas de las aplicaciones prácticas de la segmentación de imágenes son:
Se han desarrollado varios algoritmos y técnicas de uso general para la segmentación de imágenes. Para que sean útiles, estas técnicas normalmente deben combinarse con el conocimiento específico de un dominio para resolver eficazmente los problemas de segmentación del mismo.
Hay dos clases de técnicas de segmentación.
El método más simple de segmentación de imágenes se denomina método de umbralización . Este método se basa en un nivel de recorte (o un valor de umbral) para convertir una imagen en escala de grises en una imagen binaria.
La clave de este método es seleccionar el valor umbral (o valores cuando se seleccionan varios niveles). En la industria se utilizan varios métodos populares, incluidos el método de entropía máxima, el umbral de histograma equilibrado , el método de Otsu (varianza máxima) y el agrupamiento de k-medias .
Recientemente, se han desarrollado métodos para determinar el umbral de las imágenes de tomografía computarizada (TC). La idea clave es que, a diferencia del método de Otsu, los umbrales se derivan de las radiografías en lugar de la imagen (reconstruida). [21] [22]
Los nuevos métodos sugirieron el uso de umbrales no lineales basados en reglas difusas multidimensionales. En estos trabajos, la decisión sobre la pertenencia de cada píxel a un segmento se basa en reglas multidimensionales derivadas de la lógica difusa y algoritmos evolutivos basados en el entorno de iluminación de la imagen y la aplicación. [23]
El algoritmo K-means es una técnica iterativa que se utiliza para particionar una imagen en K grupos. [24] El algoritmo básico es
En este caso, la distancia es la diferencia al cuadrado o absoluta entre un píxel y el centro de un grupo. La diferencia se basa normalmente en el color , la intensidad , la textura y la ubicación del píxel, o en una combinación ponderada de estos factores. K se puede seleccionar de forma manual, aleatoria o mediante una heurística . Se garantiza que este algoritmo convergerá, pero es posible que no devuelva la solución óptima . La calidad de la solución depende del conjunto inicial de grupos y del valor de K.
El algoritmo Mean Shift es una técnica que se utiliza para dividir una imagen en un número de grupos desconocido a priori . Tiene la ventaja de no tener que empezar con una estimación inicial de dicho parámetro, lo que lo convierte en una mejor solución general para casos más diversos.
La segmentación basada en el movimiento es una técnica que se basa en el movimiento de la imagen para realizar la segmentación.
La idea es sencilla: observar las diferencias entre un par de imágenes. Suponiendo que el objeto de interés se esté moviendo, la diferencia será exactamente ese objeto.
Mejorando esta idea, Kenney et al. propusieron la segmentación interactiva [2]. Utilizan un robot para tocar objetos con el fin de generar la señal de movimiento necesaria para la segmentación basada en el movimiento.
La segmentación interactiva sigue el marco de percepción interactiva propuesto por Dov Katz [3] y Oliver Brock [4].
Otra técnica que se basa en el movimiento es la segmentación de movimiento rígido .
Los métodos basados en la compresión postulan que la segmentación óptima es la que minimiza, sobre todas las segmentaciones posibles, la longitud de codificación de los datos. [25] [26] La conexión entre estos dos conceptos es que la segmentación intenta encontrar patrones en una imagen y cualquier regularidad en la imagen puede usarse para comprimirla. El método describe cada segmento por su textura y forma de borde. Cada uno de estos componentes se modela mediante una función de distribución de probabilidad y su longitud de codificación se calcula de la siguiente manera:
Para cualquier segmentación dada de una imagen, este esquema produce el número de bits necesarios para codificar esa imagen en función de la segmentación dada. Por lo tanto, entre todas las segmentaciones posibles de una imagen, el objetivo es encontrar la segmentación que produzca la longitud de codificación más corta. Esto se puede lograr mediante un método de agrupamiento aglomerativo simple. La distorsión en la compresión con pérdida determina la tosquedad de la segmentación y su valor óptimo puede diferir para cada imagen. Este parámetro se puede estimar heurísticamente a partir del contraste de texturas en una imagen. Por ejemplo, cuando las texturas en una imagen son similares, como en las imágenes de camuflaje, se requiere una sensibilidad más fuerte y, por lo tanto, una cuantificación menor.
Los métodos basados en histogramas son muy eficientes en comparación con otros métodos de segmentación de imágenes porque normalmente requieren solo una pasada a través de los píxeles . En esta técnica, se calcula un histograma a partir de todos los píxeles de la imagen, y los picos y valles del histograma se utilizan para ubicar los grupos en la imagen. [1] Se puede utilizar el color o la intensidad como medida.
Una mejora de esta técnica consiste en aplicar de forma recursiva el método de búsqueda de histogramas a los grupos de la imagen para dividirlos en grupos más pequeños. Esta operación se repite con grupos cada vez más pequeños hasta que no se formen más grupos. [1] [27]
Una desventaja del método de búsqueda de histograma es que puede resultar difícil identificar picos y valles significativos en la imagen.
Los enfoques basados en histogramas también se pueden adaptar rápidamente para aplicarlos a múltiples fotogramas, manteniendo al mismo tiempo su eficiencia de pasada única. El histograma se puede realizar de múltiples maneras cuando se consideran múltiples fotogramas. El mismo enfoque que se adopta con un fotograma se puede aplicar a varios y, una vez que se fusionan los resultados, es más probable que se puedan distinguir los picos y valles que antes eran difíciles de identificar. El histograma también se puede aplicar por píxel, donde la información resultante se utiliza para determinar el color más frecuente para la ubicación del píxel. Este enfoque segmenta en función de los objetos activos y un entorno estático, lo que da como resultado un tipo diferente de segmentación útil en el seguimiento de video .
La detección de bordes es un campo bien desarrollado en el procesamiento de imágenes. Los límites de las regiones y los bordes están estrechamente relacionados, ya que a menudo hay un ajuste brusco de la intensidad en los límites de las regiones. Por lo tanto, las técnicas de detección de bordes se han utilizado como base de otra técnica de segmentación.
Los bordes identificados mediante la detección de bordes suelen estar desconectados. Sin embargo, para segmentar un objeto de una imagen, se necesitan límites de región cerrados. Los bordes deseados son los límites entre dichos objetos o taxones espaciales. [28] [29]
Los taxones espaciales [30] son gránulos de información [31] que consisten en una región de píxeles nítidos, ubicados en niveles de abstracción dentro de una arquitectura de escena anidada jerárquica. Son similares a la designación psicológica de la Gestalt de figura-fondo, pero se extienden para incluir el primer plano, los grupos de objetos, los objetos y las partes salientes de los objetos. Los métodos de detección de bordes se pueden aplicar a la región del taxón espacial, de la misma manera que se aplicarían a una silueta. Este método es particularmente útil cuando el borde desconectado es parte de un contorno ilusorio [32] [33]
Los métodos de segmentación también se pueden aplicar a los bordes obtenidos a partir de detectores de bordes. Lindeberg y Li [34] desarrollaron un método integrado que segmenta los bordes en segmentos de borde rectos y curvos para el reconocimiento de objetos basado en partes, con base en un criterio de longitud de descripción mínima (MDL ) que se optimizó mediante un método similar a dividir y fusionar con puntos de interrupción candidatos obtenidos a partir de señales de unión complementarias para obtener puntos más probables en los que considerar particiones en diferentes segmentos.
La detección de puntos aislados en una imagen es una parte fundamental de la segmentación de imágenes. Este proceso depende principalmente de la segunda derivada, lo que indica el uso del operador laplaciano. El laplaciano de una función viene dado por:
El operador laplaciano se emplea de modo que las derivadas parciales se deriven de una ecuación específica. La segunda derivada parcial de con respecto a y se da por:
Estas derivadas parciales se utilizan luego para calcular el laplaciano como:
Esta expresión matemática se puede implementar mediante una convolución con una máscara adecuada. Si ampliamos esta ecuación a tres dimensiones (x, y, z), la intensidad en cada ubicación de píxel alrededor de un píxel central en (x, y, z) se reemplaza por sus valores correspondientes. Esta ecuación se vuelve particularmente útil cuando asumimos que todos los píxeles tienen un espaciado unitario a lo largo de cada eje.
Se ha desarrollado una máscara esférica para su uso con conjuntos de datos tridimensionales. La máscara esférica está diseñada para utilizar únicamente aritmética de números enteros durante los cálculos, eliminando así la necesidad de hardware o software de punto flotante.
Al aplicar estos conceptos a imágenes reales representadas como matrices de números, debemos tener en cuenta lo que sucede cuando llegamos a un borde o región límite. La función se define como:
La ecuación anterior se utiliza para determinar si un punto de la imagen es un punto aislado en función de la magnitud de respuesta y un valor umbral . Si la magnitud de respuesta es mayor o igual que el umbral, la función devuelve 1, lo que indica la presencia de un punto aislado; de lo contrario, devuelve 0. Esto ayuda a detectar y segmentar eficazmente los puntos aislados de la imagen. [35]
La detección de puntos aislados tiene importantes aplicaciones en diversos campos, incluido el procesamiento de imágenes de rayos X. Por ejemplo, una imagen de rayos X original de un álabe de turbina se puede examinar píxel por píxel para detectar porosidad en el cuadrante superior derecho del álabe. El resultado de aplicar la respuesta de un detector de bordes a esta imagen de rayos X se puede aproximar. Esto demuestra la segmentación de puntos aislados en una imagen con la ayuda de sondas de un solo píxel. [36]
Este método es una combinación de tres características de la imagen: la partición de la imagen basada en el análisis del histograma se verifica por la alta compacidad de los cúmulos (objetos) y los altos gradientes de sus bordes. Para ello se deben introducir dos espacios: un espacio es el histograma unidimensional de brillo H = H ( B ); el segundo espacio es el espacio dual tridimensional de la propia imagen original B = B ( x , y ). El primer espacio permite medir la compactibilidad de la distribución del brillo de la imagen calculando un agrupamiento mínimo kmin. El umbral de brillo T correspondiente a kmin define la imagen binaria (en blanco y negro) – mapa de bits b = φ ( x , y ), donde φ ( x , y ) = 0, si B ( x , y ) < T , y φ ( x , y ) = 1, si B ( x , y ) ≥ T . El mapa de bits b es un objeto en el espacio dual. En ese mapa de bits se debe definir una medida que refleje la distribución compacta de los píxeles negros (o blancos). Por lo tanto, el objetivo es encontrar objetos con buenos bordes. Para todo T se debe calcular la medida M DC = G /( k × L ) (donde k es la diferencia de brillo entre el objeto y el fondo, L es la longitud de todos los bordes y G es el gradiente medio en los bordes). El máximo de MDC define la segmentación. [37]
Los métodos de crecimiento de regiones se basan principalmente en el supuesto de que los píxeles vecinos dentro de una región tienen valores similares. El procedimiento habitual consiste en comparar un píxel con sus vecinos. Si se cumple un criterio de similitud, el píxel puede configurarse para que pertenezca al mismo grupo que uno o más de sus vecinos. La selección del criterio de similitud es significativa y los resultados se ven influenciados por el ruido en todos los casos.
El método de Fusión Estadística de Regiones [38] (SRM) comienza construyendo el gráfico de píxeles utilizando 4-conectividad con bordes ponderados por el valor absoluto de la diferencia de intensidad. Inicialmente, cada píxel forma una única región de píxeles. Luego, SRM ordena esos bordes en una cola de prioridad y decide si fusionar o no las regiones actuales que pertenecen a los píxeles de borde utilizando un predicado estadístico.
Un método de crecimiento de regiones es el método de crecimiento de regiones con semillas. Este método toma un conjunto de semillas como entrada junto con la imagen. Las semillas marcan cada uno de los objetos que se van a segmentar. Las regiones se hacen crecer iterativamente mediante la comparación de todos los píxeles vecinos no asignados a las regiones. La diferencia entre el valor de intensidad de un píxel y la media de la región, , se utiliza como una medida de similitud . El píxel con la diferencia más pequeña medida de esta manera se asigna a la región respectiva. Este proceso continúa hasta que todos los píxeles se asignan a una región. Debido a que el crecimiento de regiones con semillas requiere semillas como entrada adicional, los resultados de la segmentación dependen de la elección de las semillas, y el ruido en la imagen puede hacer que las semillas se coloquen de manera incorrecta.
Otro método de crecimiento de regiones es el método de crecimiento de regiones sin semillas. Es un algoritmo modificado que no requiere semillas explícitas. Comienza con una sola región (el píxel elegido aquí no influye notablemente en la segmentación final). En cada iteración, considera los píxeles vecinos de la misma manera que el crecimiento de regiones con semillas. Se diferencia del crecimiento de regiones con semillas en que, si el mínimo es menor que un umbral predefinido , se agrega a la región respectiva . Si no lo es, el píxel se considera diferente de todas las regiones actuales y se crea una nueva región con este píxel.
Una variante de esta técnica, propuesta por Haralick y Shapiro (1985), [1] se basa en las intensidades de los píxeles . La media y la dispersión de la región y la intensidad del píxel candidato se utilizan para calcular una estadística de prueba. Si la estadística de prueba es suficientemente pequeña, el píxel se añade a la región y se vuelven a calcular la media y la dispersión de la región. De lo contrario, el píxel se rechaza y se utiliza para formar una nueva región.
Un método especial de crecimiento de regiones se denomina segmentación conectada (véase también conectividad lambda ). Se basa en las intensidades de los píxeles y en las rutas de enlace de vecindad. Se calcula un grado de conectividad (conectividad) en función de una ruta formada por píxeles. Para un cierto valor de , se dice que dos píxeles están conectados si existe una ruta que une esos dos píxeles y la conectividad de esta ruta es al menos . La conectividad es una relación de equivalencia. [39]
La segmentación por división y fusión se basa en una partición de árbol cuádruple de una imagen. A veces se la denomina segmentación de árbol cuádruple.
Este método comienza en la raíz del árbol que representa la imagen completa. Si se encuentra que no es uniforme (no es homogéneo), entonces se divide en cuatro cuadrados secundarios (el proceso de división), y así sucesivamente. Si, por el contrario, los cuatro cuadrados secundarios son homogéneos, se fusionan como varios componentes conectados (el proceso de fusión). El nodo en el árbol es un nodo segmentado. Este proceso continúa recursivamente hasta que no sean posibles más divisiones o fusiones. [40] [41] Cuando una estructura de datos especial está involucrada en la implementación del algoritmo del método, su complejidad temporal puede alcanzar , un algoritmo óptimo del método. [42]
Utilizando un método basado en ecuaciones diferenciales parciales (EDP) y resolviendo la ecuación de EDP mediante un esquema numérico, se puede segmentar la imagen. [43] La propagación de curvas es una técnica popular en esta categoría, con numerosas aplicaciones para la extracción de objetos, seguimiento de objetos, reconstrucción estéreo, etc. La idea central es desarrollar una curva inicial hacia el potencial más bajo de una función de costo, donde su definición refleja la tarea a abordar. Como en la mayoría de los problemas inversos , la minimización de la función de costo no es trivial e impone ciertas restricciones de suavidad en la solución, que en el presente caso se pueden expresar como restricciones geométricas en la curva en evolución.
Las técnicas lagrangianas se basan en parametrizar el contorno según alguna estrategia de muestreo y luego evolucionar cada elemento según la imagen y los términos internos. Estas técnicas son rápidas y eficientes, sin embargo, la formulación original "puramente paramétrica" (debida a Kass, Witkin y Terzopoulos en 1987 y conocida como " serpientes "), generalmente es criticada por sus limitaciones con respecto a la elección de la estrategia de muestreo, las propiedades geométricas internas de la curva, los cambios de topología (división y fusión de curvas), el abordaje de problemas en dimensiones superiores, etc. Hoy en día, se han desarrollado formulaciones "discretizadas" eficientes para abordar estas limitaciones manteniendo una alta eficiencia. En ambos casos, la minimización de la energía generalmente se lleva a cabo utilizando un descenso de gradiente más pronunciado, mediante el cual las derivadas se calculan utilizando, por ejemplo, diferencias finitas.
El método de conjunto de niveles fue propuesto inicialmente para rastrear interfaces móviles por Dervieux y Thomasset [44] [45] en 1979 y 1981 y luego fue reinventado por Osher y Sethian en 1988. [46] Esto se ha extendido a través de varios dominios de imágenes a fines de la década de 1990. Se puede utilizar para abordar de manera eficiente el problema de la propagación de curvas/superficies/etc. de manera implícita. La idea central es representar el contorno en evolución utilizando una función con signo cuyo cero corresponde al contorno real. Luego, de acuerdo con la ecuación de movimiento del contorno, se puede derivar fácilmente un flujo similar para la superficie implícita que, cuando se aplica al nivel cero, reflejará la propagación del contorno. El método de conjunto de niveles ofrece numerosas ventajas: es implícito, no tiene parámetros, proporciona una forma directa de estimar las propiedades geométricas de la estructura en evolución, permite el cambio de topología y es intrínseco. Se puede utilizar para definir un marco de optimización, como propusieron Zhao, Merriman y Osher en 1996. Se puede concluir que es un marco muy conveniente para abordar numerosas aplicaciones de visión por computadora y análisis de imágenes médicas. [47] La investigación en varias estructuras de datos de conjuntos de niveles ha llevado a implementaciones muy eficientes de este método.
El método de marcha rápida se ha utilizado en la segmentación de imágenes, [48] y este modelo se ha mejorado (permitiendo velocidades de propagación tanto positivas como negativas) en un enfoque llamado método de marcha rápida generalizado. [49]
El objetivo de los métodos variacionales es encontrar una segmentación que sea óptima con respecto a un funcional de energía específico. Los funcionales consisten en un término de ajuste de datos y términos de regularización. Un representante clásico es el modelo de Potts definido para una imagen por
Un minimizador es una imagen constante por partes que tiene un equilibrio óptimo entre la distancia L2 al cuadrado a la imagen dada y la longitud total de su conjunto de saltos. El conjunto de saltos de define una segmentación. El peso relativo de las energías se ajusta mediante el parámetro . La variante binaria del modelo de Potts, es decir, si el rango de está restringido a dos valores, a menudo se denomina modelo de Chan - Vese . [50] Una generalización importante es el modelo de Mumford-Shah [51] dado por
El valor funcional es la suma de la longitud total de la curva de segmentación , la suavidad de la aproximación y su distancia a la imagen original . El peso de la penalización por suavidad se ajusta mediante . El modelo de Potts a menudo se denomina modelo de Mumford-Shah constante por partes, ya que puede considerarse el caso degenerado . Se sabe que los problemas de optimización son NP-hard en general, pero las estrategias de minimización cercana funcionan bien en la práctica. Los algoritmos clásicos son la no convexidad graduada y la aproximación de Ambrosio-Tortorelli .
Los métodos de partición de grafos son herramientas eficaces para la segmentación de imágenes, ya que modelan el impacto de los vecindarios de píxeles en un grupo determinado de píxeles o píxeles, bajo el supuesto de homogeneidad en las imágenes. En estos métodos, la imagen se modela como un grafo ponderado, no dirigido . Por lo general, un píxel o un grupo de píxeles se asocian con nodos y los pesos de los bordes definen la (des)similitud entre los píxeles del vecindario. Luego, el grafo (imagen) se particiona de acuerdo con un criterio diseñado para modelar grupos "buenos". Cada partición de los nodos (píxeles) resultantes de estos algoritmos se considera un segmento de objeto en la imagen; consulte Categorización de objetos basada en segmentación . Algunos algoritmos populares de esta categoría son cortes normalizados, [52] caminante aleatorio , [53] corte mínimo, [54] partición isoperimétrica, [55] segmentación basada en árbol de expansión mínimo , [56] y categorización de objetos basada en segmentación .
La aplicación de campos aleatorios de Markov (MRF) para imágenes fue sugerida a principios de 1984 por Geman y Geman. [57] Su sólida base matemática y su capacidad para proporcionar un óptimo global incluso cuando se define en características locales demostraron ser la base para una nueva investigación en el dominio del análisis de imágenes, la eliminación de ruido y la segmentación. Los MRF se caracterizan completamente por sus distribuciones de probabilidad previa, distribuciones de probabilidad marginal, camarillas , restricción de suavizado, así como por un criterio para actualizar valores. El criterio para la segmentación de imágenes utilizando MRF se replantea como encontrar el esquema de etiquetado que tiene la máxima probabilidad para un conjunto dado de características. Las amplias categorías de segmentación de imágenes utilizando MRF son la segmentación supervisada y no supervisada.
En términos de segmentación de imágenes, la función que los MRF buscan maximizar es la probabilidad de identificar un esquema de etiquetado dado que se detecta un conjunto particular de características en la imagen. Esta es una reformulación del método de estimación máxima a posteriori .
El algoritmo genérico para la segmentación de imágenes utilizando MAP se detalla a continuación:
Cada algoritmo de optimización es una adaptación de modelos de distintos campos y se distinguen por sus funciones de costo únicas. El rasgo común de las funciones de costo es penalizar el cambio en el valor del píxel, así como la diferencia en la etiqueta del píxel en comparación con las etiquetas de los píxeles vecinos.
El algoritmo de modos condicionales iterados (ICM) intenta reconstruir el esquema de etiquetado ideal cambiando los valores de cada píxel en cada iteración y evaluando la energía del nuevo esquema de etiquetado utilizando la función de costo que se muestra a continuación.
donde α es la penalización por el cambio en la etiqueta del píxel y β es la penalización por la diferencia en la etiqueta entre los píxeles vecinos y el píxel elegido. Aquí se encuentra la vecindad del píxel i y δ es la función delta de Kronecker. Un problema importante con el ICM es que, de manera similar al descenso de gradiente, tiene una tendencia a permanecer sobre máximos locales y, por lo tanto, no obtener un esquema de etiquetado globalmente óptimo.
El recocido simulado (SA) , que se deriva como un análogo del recocido en metalurgia, utiliza el cambio en la etiqueta de píxel a lo largo de las iteraciones y estima la diferencia en energía de cada gráfico recién formado con respecto a los datos iniciales. Si el gráfico recién formado es más rentable, en términos de bajo costo de energía, se obtiene:
El algoritmo selecciona el gráfico recién formado. El recocido simulado requiere la entrada de programas de temperatura que afectan directamente la velocidad de convergencia del sistema, así como el umbral de energía para que se produzca la minimización.
Existe una variedad de otros métodos para resolver MRF simples y de orden superior. Incluyen la maximización del margen posterior, la estimación de MAP multiescala, [58] la segmentación de resolución múltiple [59] y más. Además de las estimaciones de verosimilitud, existen métodos de corte de gráficos con flujo máximo [60] y otros métodos basados en gráficos altamente restringidos [61] [62] para resolver MRF.
El algoritmo de expectativa-maximización se utiliza para estimar iterativamente las probabilidades y distribuciones a posteriori de etiquetado cuando no se dispone de datos de entrenamiento y no se puede formar una estimación del modelo de segmentación. Un enfoque general consiste en utilizar histogramas para representar las características de una imagen y proceder como se describe brevemente en este algoritmo de tres pasos:
1. Se utiliza una estimación aleatoria de los parámetros del modelo.
2. Paso E: Estimar las estadísticas de clase en función del modelo de segmentación aleatoria definido. Con ellas, calcular la probabilidad condicional de pertenecer a una etiqueta dado el conjunto de características, calculada mediante el teorema de Bayes ingenuo .
Aquí , el conjunto de todas las etiquetas posibles.
3. Paso M: La relevancia establecida de un conjunto de características dado para un esquema de etiquetado se utiliza ahora para calcular la estimación a priori de una etiqueta dada en la segunda parte del algoritmo. Dado que se desconoce la cantidad real de etiquetas totales (de un conjunto de datos de entrenamiento), en los cálculos se utiliza una estimación oculta de la cantidad de etiquetas proporcionada por el usuario.
donde es el conjunto de todas las características posibles.
La transformación de la cuenca hidrográfica considera la magnitud del gradiente de una imagen como una superficie topográfica. Los píxeles que tienen las intensidades de magnitud del gradiente (GMI) más altas corresponden a las líneas divisorias de aguas, que representan los límites de la región. El agua colocada en cualquier píxel encerrado por una línea divisoria de aguas común fluye cuesta abajo hasta un mínimo de intensidad local (LIM) común. Los píxeles que drenan hasta un mínimo común forman una cuenca de captación, que representa un segmento.
El supuesto central de los enfoques basados en modelos es que las estructuras de interés tienen una tendencia hacia una forma particular. Por lo tanto, se puede buscar un modelo probabilístico que caracterice la forma y su variación. Al segmentar una imagen, se pueden imponer restricciones utilizando este modelo como una prioridad. [63] Tal tarea puede implicar (i) el registro de los ejemplos de entrenamiento en una pose común, (ii) la representación probabilística de la variación de las muestras registradas y (iii) la inferencia estadística entre el modelo y la imagen. Otros métodos importantes en la literatura para la segmentación basada en modelos incluyen modelos de forma activa y modelos de apariencia activa .
Las segmentaciones de imágenes se calculan en múltiples escalas en el espacio de escala y, a veces, se propagan desde escalas gruesas a escalas finas; consulte segmentación en el espacio de escala .
Los criterios de segmentación pueden ser arbitrariamente complejos y pueden tener en cuenta criterios globales y locales. Un requisito común es que cada región debe estar conectada de algún modo.
El trabajo seminal de Witkin [64] [65] en el espacio de escala incluyó la noción de que una señal unidimensional podría segmentarse de manera inequívoca en regiones, con un parámetro de escala que controla la escala de segmentación.
Una observación clave es que los cruces por cero de las derivadas segundas (mínimos y máximos de la derivada primera o pendiente) de versiones suavizadas en múltiples escalas de una señal forman un árbol de anidamiento, que define relaciones jerárquicas entre segmentos a diferentes escalas. Específicamente, los extremos de pendiente a escalas gruesas se pueden rastrear hasta las características correspondientes a escalas finas. Cuando un máximo de pendiente y un mínimo de pendiente se aniquilan entre sí a una escala mayor, los tres segmentos que separaron se fusionan en un solo segmento, definiendo así la jerarquía de segmentos.
Se han realizado numerosos trabajos de investigación en este campo, de los cuales algunos han llegado a un punto en el que se pueden aplicar ya sea con intervención manual interactiva (normalmente en aplicaciones de imágenes médicas) o de forma totalmente automática. A continuación se presenta una breve descripción de algunas de las principales ideas de investigación en las que se basan los enfoques actuales.
Sin embargo, la estructura de anidamiento que describió Witkin es específica para señales unidimensionales y no se transfiere fácilmente a imágenes de dimensiones superiores. No obstante, esta idea general ha inspirado a varios otros autores a investigar esquemas de segmentación de imágenes de grueso a fino. Koenderink [66] propuso estudiar cómo evolucionan los contornos de isointensidad a lo largo de las escalas y este enfoque fue investigado con más detalle por Lifshitz y Pizer [67] . Desafortunadamente, sin embargo, la intensidad de las características de la imagen cambia a lo largo de las escalas, lo que implica que es difícil rastrear características de imágenes de escala gruesa a escalas más finas utilizando información de isointensidad.
Lindeberg [68] [69] estudió el problema de vincular los extremos locales y los puntos de silla sobre escalas, y propuso una representación de imagen llamada boceto primario del espacio de escala que hace explícitas las relaciones entre las estructuras a diferentes escalas, y también hace explícitas qué características de la imagen son estables en grandes rangos de escala, incluidas las escalas localmente apropiadas para ellas. Bergholm propuso detectar bordes en escalas gruesas en el espacio de escala y luego rastrearlos hasta escalas más finas con la elección manual tanto de la escala de detección gruesa como de la escala de localización fina.
Gauch y Pizer [70] estudiaron el problema complementario de crestas y valles en múltiples escalas y desarrollaron una herramienta para la segmentación interactiva de imágenes basada en cuencas hidrográficas multiescala. Olsen y Nielsen [71] también investigaron el uso de cuencas hidrográficas multiescala con aplicación al mapa de gradiente y Dam [ 72] lo trasladó al uso clínico. Vincken et al. [73] propusieron una hiperpila para definir relaciones probabilísticas entre estructuras de imágenes a diferentes escalas. Ahuja [74] [75] y sus colaboradores avanzaron en el uso de estructuras de imágenes estables en distintas escalas hasta convertirlas en un sistema completamente automatizado. Undeman y Lindeberg [76] presentaron un algoritmo de segmentación cerebral completamente automático basado en ideas estrechamente relacionadas de cuencas hidrográficas multiescala y lo probaron ampliamente en bases de datos cerebrales.
Florack y Kuijper también han retomado estas ideas para la segmentación de imágenes en múltiples escalas mediante la vinculación de estructuras de imágenes a lo largo de las escalas. [77] Bijaoui y Rué [78] asocian las estructuras detectadas en el espacio de escala por encima de un umbral de ruido mínimo en un árbol de objetos que abarca múltiples escalas y corresponde a un tipo de característica en la señal original. Las características extraídas se reconstruyen con precisión utilizando un método iterativo de matriz de gradiente conjugado.
En un tipo de segmentación, el usuario delinea la región de interés con los clics del mouse y se aplican algoritmos para que se muestre la ruta que mejor se ajusta al borde de la imagen.
En este tipo de segmentación se utilizan técnicas como SIOX , Livewire , Intelligent Scissors o IT-SNAPS. En un tipo alternativo de segmentación semiautomática, los algoritmos devuelven un taxón espacial (es decir, primer plano, grupo de objetos, objeto o parte de objeto) seleccionado por el usuario o designado a través de probabilidades previas. [79] [80]
La mayoría de los métodos de segmentación mencionados anteriormente se basan únicamente en la información de color de los píxeles de la imagen. Los seres humanos utilizan mucho más conocimiento al segmentar imágenes, pero implementar este conocimiento costaría una considerable cantidad de ingeniería humana y tiempo computacional, y requeriría una enorme base de datos de conocimiento del dominio que actualmente no existe. Los métodos de segmentación entrenables, como la segmentación de redes neuronales , superan estos problemas al modelar el conocimiento del dominio a partir de un conjunto de datos de píxeles etiquetados.
Una red neuronal de segmentación de imágenes puede procesar áreas pequeñas de una imagen para extraer características simples como los bordes. [81] Otra red neuronal, o cualquier mecanismo de toma de decisiones, puede combinar estas características para etiquetar las áreas de una imagen en consecuencia. Un tipo de red diseñada de esta manera es el mapa de Kohonen .
Las redes neuronales acopladas a pulsos (PCNN) son modelos neuronales propuestos mediante el modelado de la corteza visual de un gato y desarrollados para el procesamiento de imágenes biomiméticas de alto rendimiento . En 1989, Reinhard Eckhorn presentó un modelo neuronal para emular el mecanismo de la corteza visual de un gato. El modelo de Eckhorn proporcionó una herramienta simple y efectiva para estudiar la corteza visual de pequeños mamíferos, y pronto se reconoció que tenía un potencial de aplicación significativo en el procesamiento de imágenes. En 1994, el modelo de Eckhorn fue adaptado para ser un algoritmo de procesamiento de imágenes por John L. Johnson, quien denominó a este algoritmo Red neuronal acoplada a pulsos. [82] Durante la última década, las PCNN se han utilizado para una variedad de aplicaciones de procesamiento de imágenes, que incluyen: segmentación de imágenes, generación de características, extracción de rostros, detección de movimiento, crecimiento de regiones, reducción de ruido, etc. Una PCNN es una red neuronal bidimensional. Cada neurona en la red corresponde a un píxel en una imagen de entrada, recibiendo la información de color de su píxel correspondiente (por ejemplo, intensidad) como estímulo externo. Cada neurona también se conecta con sus neuronas vecinas, recibiendo estímulos locales de ellas. Los estímulos externos y locales se combinan en un sistema de activación interna, que acumula los estímulos hasta que supera un umbral dinámico, lo que da como resultado una salida de pulso. A través de un cálculo iterativo, las neuronas PCNN producen series temporales de salidas de pulso. Las series temporales de salidas de pulso contienen información de imágenes de entrada y se pueden utilizar para varias aplicaciones de procesamiento de imágenes, como la segmentación de imágenes y la generación de características. En comparación con los medios de procesamiento de imágenes convencionales, las PCNN tienen varias ventajas significativas, incluida la robustez frente al ruido, la independencia de las variaciones geométricas en los patrones de entrada, la capacidad de superar pequeñas variaciones de intensidad en los patrones de entrada, etc.
U-Net es una red neuronal convolucional que toma como entrada una imagen y genera una etiqueta para cada píxel. [83] U-Net se desarrolló inicialmente para detectar los límites de las células en imágenes biomédicas. U-Net sigue la arquitectura clásica de autocodificador , por lo que contiene dos subestructuras. La estructura del codificador sigue la pila tradicional de capas convolucionales y de agrupamiento máximo para aumentar el campo receptivo a medida que pasa por las capas. Se utiliza para capturar el contexto en la imagen. La estructura del decodificador utiliza capas de convolución transpuestas para el sobremuestreo de modo que las dimensiones finales sean cercanas a las de la imagen de entrada. Las conexiones de salto se colocan entre las capas de convolución y convolución transpuesta de la misma forma para preservar los detalles que se habrían perdido de otra manera.
Además de las tareas de segmentación semántica a nivel de píxel que asignan una categoría determinada a cada píxel, las aplicaciones de segmentación modernas incluyen tareas de segmentación semántica a nivel de instancia en las que cada individuo en una categoría determinada debe ser identificado de forma única, así como tareas de segmentación panóptica que combinan estas dos tareas para proporcionar una segmentación de escena más completa. [20]
Las imágenes relacionadas, como un álbum de fotos o una secuencia de fotogramas de vídeo, a menudo contienen objetos y escenas semánticamente similares, por lo que suele ser beneficioso explotar dichas correlaciones. [84] La tarea de segmentar simultáneamente escenas a partir de imágenes o fotogramas de vídeo relacionados se denomina cosegmentación , [16] que normalmente se utiliza en la localización de acciones humanas . A diferencia de la detección de objetos basada en cuadros delimitadores convencionales , los métodos de localización de acciones humanas proporcionan resultados más granulares, normalmente máscaras de segmentación por imagen que delimitan el objeto humano de interés y su categoría de acción (por ejemplo, Segment-Tube [17] ). A menudo se emplean técnicas como redes dinámicas de Markov , CNN y LSTM para explotar las correlaciones entre fotogramas.
Existen muchos otros métodos de segmentación, como la segmentación multiespectral o la segmentación basada en conectividad basada en imágenes DTI . [85] [86]
{{cite journal}}
: Requiere citar revista |journal=
( ayuda )