Servicios Personalizados
Revista
Articulo
Indicadores
- Citado por SciELO
- Accesos
Links relacionados
- Similares en SciELO
Compartir
Ingeniería, investigación y tecnología
versión On-line ISSN 2594-0732versión impresa ISSN 1405-7743
Ing. invest. y tecnol. vol.12 no.1 Ciudad de México ene./mar. 2011
Registro de imágenes videoendoscópicas para reconstruir cartografías de órganos huecos
VideoEndoscopic Image Registration for Cartography Reconstruction of Hollow Organs
MirandaLuna R.1, PosadaGómez R.2, AlorHernández G.3, MartínezSibaja A4. y CortésRobles G.5
1 Universidad Tecnológica de la Mixteca, Oaxaca. Email: rmiranda@yahoo.com
2 Departamento de Posgrado e Investigación Instituto Tecnológico de Orizaba, Veracruz. Email: rposada@itorizaba.edu.mx
3 Departamento de Posgrado e Investigación Instituto Tecnológico de Orizaba, Veracruz. Email: galor@itorizaba.edu.mx
4 Departamento de Posgrado e Investigación Instituto Tecnológico de Orizaba, Veracruz. Email: asibaja@itorizaba.edu.mx
5 Departamento de Posgrado e Investigación Instituto Tecnológico de Orizaba, Veracruz. Email: grobles@itorizaba.edu.mx
Información del artículo: Recibido: septiembre de 2006.
Reevaluado: noviembre de 2007.
Aceptado: septiembre de 2010.
Resumen
Los sistemas videoendoscópicos se usan para examinar órganos huecos permitiendo observar las capas que recubren las cavidades internas de los órganos (superficies internas). Comúnmente, un especialista debe analizar la secuencia de imágenes adquiridas durante el examen endoscópico, ya sea para confirmar el diagnóstico efectuado durante la endoscopia o para valorar la evolución de una patología. En una secuencia de imágenes videoendoscópicas (en donde cada imagen muestra una pequeña área de la superficie tisular examinada), puede resultar difícil localizar el área de interés debido a: 1) la gran cantidad de imágenes que componen la secuencia, 2) la naturaleza secuencial del archivo (que obliga al médico a repasar toda la secuencia para poder encontrar la zona de interés) y 3) la posibilidad de que la zona de interés se encuentre o esté contenida en varios tramos no contiguos de la secuencia. En este trabajo, se implementa una metodología que permite la reconstrucción automática de una cartografía de la región de interés, permitiendo su visualización panorámica y eliminando información redundante. Además, se presenta un método basado en un algoritmo automático de registro de imágenes para la reconstrucción de cartografías. Finalmente, se presenta un protocolo de evaluación y los resultados obtenidos.
Descriptores: cartografías de órganos internos, videoendoscopia, registro de imágenes, corrección de la distorsión radial.
Abstract
In examinations of hollow organs, videoendoscopic systems are used. These systems allow to analyze layers recovering the internal organ cavities (internal surfaces). In this kind of analysis, an expert must review the acquired image sequence during an endoscopic analysis, either to confirm the diagnosis issue of endoscope, or to evaluate the pathology's evolution. For these reasons, it is difficult to find the surface of interest in a videoendoscopic image sequence (where each image shows a little area of examined tissue surface) due to: (1) the big number of images into the sequence, (2) the sequential nature of file (forcing doctor to review all sequence to find the interest region), and (3) frequently the interest region that could be physically contiguous is shown in frames faraway from each other. In this work, a methodology for automatic cartography reconstruction of the interest region is proposed. This methodology provides a panoramic view and eliminates redundant information. Furthermore, an automatic image registration algorithm for cartography construction is presented. Finally, an evaluation protocol and the obtained results are also presented.
Keywords: internal organs cartography, videoendoscopy, image registration, radial distortion correction.
Introducción
Actualmente, los sistemas videoendoscópicos utilizados en aplicaciones médicas generan secuencias de imágenes en los que cada cuadro contiene un fragmento de la superficie total escaneada. Una cartografía (por la analogía con el trazado de mapas geográficos) del área escaneada, constituye así una representación fácil de manejar, la cual es útil en un análisis postoperatorio. En una secuencia de video, las imágenes adyacentes comparten una superficie común (información redundante) que permiten construir el panorama de la superficie observada utilizando las imágenes de la secuencia. Para construir cartografías a partir de imágenes sucesivas, se propone efectuar un registro de estas imágenes y combinarlas para constituir la imagen panorámica.
El registro de imágenes endoscópicas se ve afectado por tres características inherentes a la naturaleza de las imágenes adquiridas (figura 1):
1) La distorsión geométrica no lineal de tipo radial causada por un sistema óptico diseñado para brindar un amplio campo de visión distorsiona las imágenes haciendo que un mismo punto de la escena ocupe posiciones relativas diferentes en dos imágenes sucesivas a causa de esta distorsión (Tian et al., 2002), (Smith et al, 1992);
2) El patrón fijo de alta frecuencia introducido por la fibra óptica de los fibroscopios (endoscopios que utilizan fibras ópticas para recuperar la imagen) y que corresponde al arreglo de las fibras, se superpone a la imagen generada y;
3) Un patrón de baja frecuencia generado por una iluminación no uniforme al interior del órgano (normalmente más importante en el centro que en los extremos de la imagen).
En este trabajo se presenta un método para reconstruir una cartografía de un órgano explorado a través de un registro monomodal.
La estructura de este trabajo es la siguiente: primero, se describe la implementación de un método de registro que corrige las distorsiones propias al sistema de adquisición y efectúa la cartografía de una región de interés. Después, se describe la metodología empleada en el desarrollo de la cartografía explicando cada etapa involucrada en este proceso. Posteriormente, se describen las pruebas efectuadas al sistema y los resultados obtenidos. Finalmente, se presenta la conclusión donde se analizan los resultados y se discuten las posibilidades del sistema propuesto.
Metodología para la creación de la cartografía
La metodología para la realización de la cartografía se divide en tres etapas principales:
1. El acondicionamiento de las imágenes que corrige la distorsión radial y reduce la presencia de los patrones de alta y baja frecuencia;
2. El registro de las imágenes previamente corregidas usando un algoritmo automático basado en la información mutua de las imágenes y;
3. La creación de la imagen panorámica conteniendo todas las imágenes registradas en un mismo referencial.
Acondicionamiento de las imágenes endoscópicas
Para reducir los problemas que afectan el registro de las imágenes se procede a acondicionar las imágenes efectuando primeramente la corrección de la distorsión radial (mediante un algoritmo de registro de imágenes) y filtrando posteriormente las imágenes resultantes para eliminar los patrones de alta y baja frecuencia.
A continuación, se describen las técnicas empleadas en la metodología propuesta para corregir los principales tipos de distorsión generada por el endoscopio (la distorsión radial y la corrección de patrones de alta y baja frecuencia).
Corrección de la distorsión radial
La transformación geométrica que relaciona las imágenes adquiridas a través del endoscopio, se afecta tanto por la posición del endoscopio desde la cual se adquieren como por la óptica del endoscopio. La distorsión radial que afecta las imágenes endoscópicas (figura 1a) se modela con las coordenadas de la proyección del centro óptico en la imagen (Cx, Cy) y la relación polinomial r' expresada en la siguiente fórmula:
donde
r es la distancia radial entre un punto de coordenadas (x, y) (Cx, Cy),
ki son los coeficientes de corrección usados para obtener una imagen sin distorsión,
M es el grado del polinomio.
El número de parámetros ki que deben calcularse depende de la importancia de la distorsión, así como de la precisión requerida. En este trabajo, se seleccionó M = 5 que representa un buen compromiso entre exactitud y tiempo de procesamiento.
El problema de la corrección radial se observa como un problema de calibración del endoscopio. Así, sea U una imagen ideal libre de deformaciones, se considera que la imagen V es la imagen U adquirida por medio del endoscopio a calibrar y se afectó por una transformación T, a su vez conformada por una transformación Ti (parámetros intrínsecos del endoscopio) dada por la distorsión radial y una transformación Te (parámetros extrínsecos) causada por la posición desde la cual se efectuó la adquisición de la imagen U, de manera que:
De esta manera, efectuando un registro entre las imágenes U y V se calculan los valores de la matriz T. Los parámetros de la matriz proyectiva (traslaciones: tx, ty, rotación: θ, escala: α) y del polinomio (Cx, Cy y ki) son independientes; de esta forma, los parámetros de corrección radial se aíslan después del registro, con el fin de aplicar la corrección radial de manera sistemática a todas las imágenes adquiridas con el endoscopio, siempre y cuando no se modifique la óptica del mismo. Una explicación detallada del procedimiento de corrección radial o calibración del endoscopio se encuentra en (Miranda et al, 2004).
Una vez corregida la distorsión radial en toda la secuencia de imágenes, se procede a reducir la presencia de los artefactos de alta y baja frecuencia generados por el arreglo de fibras ópticas y la iluminación no uniforme, respectivamente.
Reducción de la presencia de patrones de alta y baja frecuencia
La reducción de la presencia de los patrones de alta y baja frecuencia se obtiene filtrando la imagen. Para eliminar las altas frecuencias añadidas a las imágenes por el arreglo de fibras ópticas (figura 1b), se emplea un filtro con una función gaussiana, cuya función de transferencia en el dominio de la frecuencia (González, 2002) se expresa como:
donde
D(u, v) es la distancia espacial desde la componente de frecuencia (u, v) hasta el centro del espectro (u = 0, v = 0) y;
D0 es la desviación estándar del filtro.
Una vez analizado el espectro de Fourier correspondiente al patrón de fibras, se fija D0 = 4, obteniendo así un filtro para bajas. Para corregir el problema de la iluminación no uniforme, primero se obtiene una versión filtrada de la imagen endoscópica conteniendo las componentes de muy baja frecuencia, producidas esencialmente por la distribución no uniforme de la iluminación, mediante el filtro gaussiano con D0 = 16. La imagen obtenida se resta de la imagen previamente filtrada con D0 = 4. De esta manera, se tiene una imagen corregida que finalmente se utiliza en el registro de imágenes.
Registro de imágenes
El registro de imágenes se efectúa con pares de imágenes en niveles de gris. Dadas dos imágenes U y V, su registro se expresa como:
En la fórmula (4), representa la transformación buscada en el espacio de transformaciones T, es decir, la transformación que optimiza la medida de similitud S calculada con las estructuras homologas extraídas de las imágenes U y V con ayuda de algoritmos de segmentación representados por ft y fm. El algoritmo de registro emplea los niveles de grises como estructuras homologas, evitando la necesidad de efectuar una segmentación de las imágenes. En este trabajo, se implementó como medida de similitud la información mutua (Viola, 1995). Sea W el resultado de la transformación T aplicada sobre la imagen V, la información Mutua I(U, W) se expresa por la siguiente fórmula:
donde
H(U) y H(W) representan la entropía de las imágenes U y W, respectivamente dadas por:
donde
xn es el nesimo valor del espacio Ωx (Ωx es el intervalo de niveles de gris Ε [0255]).
H(U, W) representa la entropía conjunta de las imágenes y se define como:
donde:
xn y yn son los nesimos valores de Ωx y Ωy (niveles de gris: 0 255) para las variables x y y.
P(xn, yn) es la probabilidad de concurrencia de los eventos xn, yn.
La información mutua es una medida de similitud robusta que permite conocer la información que comparten las partes superpuestas de dos imágenes caracterizadas únicamente por la distribución de sus niveles de gris. Esta medida de similitud es poco sensible a los problemas de áreas con poco recubrimiento espacial común al utilizar, tanto la entropía de cada imagen como la entropía conjunta de ambas.
El modelo de transformación usado considera translaciones, rotaciones, factor de escala y perspectiva (que representan los cambios que pueden darse entre dos imágenes consecutivas por la manipulación del endoscopio). Para la optimización de los parámetros de la matriz de transformación que permiten optimizar la información mutua, se implementó el algoritmo de gradiente descendente estocástico, debido a la existencia de mínimos locales en el espacio de búsqueda generados esencialmente por la interpolación de los niveles de gris. Este algoritmo es un método de optimización iterativo bien adaptado para superar los pequeños óptimos locales capaces de impedir la convergencia global. Los valores iniciales de la matriz de transformación son los de la matriz identidad y se actualizan en cada iteración, de manera que:
donde λ es el parámetro de rapidez de convergencia, el cual se reduce progresivamente a medida que se alcanza el máximo de la información mutua.
El registro se efectúa entre pares de imágenes consecutivos en la secuencia generada por el de videoendoscopio. Una vez hecho esto, es necesario que todas las imágenes transformadas después del registro se sitúen en un referencial común para construir la cartografía.
Creación del mosaico
El proceso de registro entre pares de imágenes consecutivas genera una matriz de transformación Tr que optimiza la información mutua de dos imágenes adyacentes i e i 1. Posteriormente, se calcula para cada imagen la transformación global Tg que permita tener todas las imágenes en un mismo sistema de coordenadas, mediante la ecuación:
Esto permite situar todas las imágenes en el referencial de la primera imagen. Después de aplicar la transformación global a todas las imágenes, sólo resta combinarlas para formar una sola imagen que constituye el panorama de toda la superficie explorada (cartografía 2D).
Pruebas efectuadas
Con el fin de validar el método de creación de cartografías y las diferentes etapas que lo componen, se efectuaron pruebas para determinar la eficiencia en la corrección de la distorsión radial. Así también, se efectuó la reconstrucción de una cartografía a partir de imágenes adquiridas empleando un sistema de posicionamiento que permite conocer los desplazamientos reales del endoscopio y una imagen panorámica que simula la superficie de la cavidad interna a explorar. Al tratarse de una imagen bidimensional se efectúa una comparación precisa con los resultados de la cartografía reconstruida.
Sistema de posicionamiento
Con el fin de conocer los desplazamientos efectuados por el endoscopio al momento de la adquisición de la secuencia de imágenes, se construyó un sistema de posicionamiento (figura 2) constituido por una plataforma en la que se montaron 3 tablas micrométricas (NEWPORT modelo MUMR8.51). Las tablas se ensamblaron de manera tal que permiten obtener desplazamientos lineales en los 3 ejes (,,). El equipo de posicionamiento permite controlar de manera precisa los desplazamientos manualmente efectuados del endoscopio dentro de un volumen de 50 × 50 × 40 mm3.
Para la adquisición de las imágenes, se empleó una cámara Leritier 760 ULL y un fibroscopio (FUJINONFG100FP), el cual se fijó al sistema de posicionamiento con ayuda de una guía (figura 2b). De esta manera, la punta del endoscopio se alinea con el borde de la guía. El sistema permite colocar la punta del endoscopio a cualquier distancia entre 10 mm y 60 mm con respecto a la escena. Además, el sistema permite girar el endoscopio para simular adquisiciones perspectivas. En la parte inferior del sistema se sitúa la imagen a adquirir (fotografía de la región a explorar o patrón para evaluar cuantitativamente la calibración). La imagen a adquirir puede rotarse alrededor del eje a pasos de 2°, debido a un pivote que se encuentra en el centro de la base de adquisición y a un patrón con marcas colocado debajo de la imagen.
El sistema descrito permite efectuar, tanto las pruebas de la calibración como de la reconstrucción de la cartografía.
Calibración del endoscopio (corrección de la distorsión radial)
Para evaluar el método de corrección de la distorsión radial, se utilizó una imagen patrón (figura 3a). La imagen patrón, se adquirió situando el endoscopio perpendicular a ésta. El posicionamiento del endoscopio se efectuó empleando el sistema de posicionamiento antes descrito.
Para la corrección de distorsión radial, se inicializan los parámetros del polinomio ki = 0 y Cx y Cy toman las coordenadas del centro de la imagen. La matriz proyectiva, se inicializa con tx = ty = 1, θ = 0 y α = 0.85.
El error en la corrección de la distorsión radial se evalúa calculando la distancia euclidiana entre puntos de referencia identificados, tanto en la imagen patrón sin distorsión como en la imagen endoscópica binarizada. Estos puntos se definen por los centros de los cuadros blancos.
Con el fin de evaluar la influencia del ángulo de perspectiva en el método de calibración, se realizaron (con ayuda del sistema de posicionamiento) adquisiciones de la imagen patrón a diferentes ángulos (0°, 5° y 10°), para generar las perspectivas con respecto a la perpendicular al plano del patrón.
Debido a que el campo de visión del endoscopio se calculó en 449 píxeles, finalmente se normalizan los errores promedio (eprom) y máximos (emáx) que representan el comportamiento típico del sistema y el peor escenario, respectivamente (Shahidi et al, 2002).
Adquisición y reconstrucción de una secuencia
Empleando el sistema de posicionamiento se adquirió una secuencia de imágenes, simulando una exploración de endoscopia. Primeramente, se colocó sobre la plataforma la imagen a adquirir. Durante la adquisición, el endoscopio se situó a una altura de 40 mm con respecto a la plataforma. Posteriormente, se adquirieron 20 imágenes, desplazando progresivamente la punta del endoscopio de izquierda a derecha (dirección positiva del eje horizontal) (figura 4a). Cada imagen se adquirió después de un desplazamiento de 2.5 mm. Este valor representa el desplazamiento máximo de la punta del endoscopio durante un examen de endoscopia clínica (se determinó de manera empírica tomando en cuenta la frecuencia de adquisición de 25 cuadros por segundo). Esta secuencia permite probar la respuesta del algoritmo, ante un desplazamiento puro.
A partir de la última posición (imagen I20) se efectúa una segunda etapa de 10 adquisiciones, en las que la punta del endoscopio se desplaza en dirección diagonal sobre el plano (, ) con incrementos de 2.5 y 2.5 mm, respectivamente, a lo largo de cada eje. Estos desplazamientos se acompañaron de rotaciones de la imagen alrededor del eje de +2° entre cada adquisición. Finalmente, se efectuaron 10 adquisiciones desplazando el endoscopio de forma diagonal sobre el plano (, ) con incrementos de 2.5 y +2.5 mm, respectivamente a lo largo de cada eje, al mismo tiempo se efectuaron rotaciones alrededor de los ejes y de 2° y 2°, respectivamente, produciéndose así un efecto de perspectiva en las imágenes.
Para evaluar el método propuesto, se efectuó primeramente la adquisición de una secuencia utilizando una fotografía digital de la superficie interna de la vejiga de un puerco (figura 4a) y posteriormente un diseño constituido de cuadros alineados en forma de rejilla con una separación de 10 mm entre cada línea (figura 4b). Las 40 imágenes de la primera adquisición (fotografía de la vejiga) se registraron con el fin de obtener los parámetros de transformación. Los parámetros de transformación obtenidos con la primera prueba se aplicaron a la segunda secuencia (patrón de puntos), a fin de reconstruir la cartografía de dicho patrón. La calidad de la cartografía de la primera adquisición (vejiga) se evalúa así, comparando el patrón de puntos reconstruido contra el patrón de puntos original. El error se calcula como la suma de las distancias euclidianas entre los centros de puntos correspondientes del patrón reconstruido y del patrón original, esto se expresa mediante la siguiente ecuación:
donde
N representa los puntos a comparar,
xi y yi son los valores absolutos de los errores en las direcciones vertical y horizontal entre el centro del iésimo punto en el patrón original y su correspondiente en el mosaico reconstruido.
Resultados
Corrección de la distorsión radial
La tabla 1 muestra los valores de los parámetros calculados por el algoritmo de calibración para la imagen patrón adquirida desde 3 ángulos diferentes (figura 3). Los valores de k4 y k5 son similares independientemente del ángulo de adquisición, lo cual confirma la irrelevancia del ángulo de perspectiva. Se realizó una comparación de los errores normalizados obtenidos para todas las imágenes con el fin de evaluar con mayor exactitud la influencia del ángulo de inclinación obteniendo un error promedio eprom de 0.21 a 0.34%. Los valores del error máximo emáx son de 2 a 4 veces mayores que eprom. Sin embargo, los errores obtenidos son comparables con otros resultados (Asari et al., 1999).
De esta manera, el método propuesto sólo requiere que el endoscopio se posicione de manera perpendicular al patrón de la figura 3a para efectuar la calibración, el error final de la reconstrucción de la cartografía no varía, aún cuando se haya adquirido la imagen en un rango de 10°.
Eliminación de los patrones de alta y baja frecuencia
Los resultados del filtro aplicado a las imágenes adquiridas para eliminar las componentes correspondientes al patrón de fibras y la iluminación no homogénea se muestran en la figura 5. Aunque la calidad visual de la imagen después del filtrado no resulta apropiada para un análisis visual, la información que contiene es adecuada para su uso en el procedimiento de registro, por lo que una vez calculado los parámetros de la transformación pueden aplicarse a las imágenes sin filtrar y permitir una mejor visualización de la cartografía.
Reconstrucción de la cartografía (Mosaico)
El mosaico mostrado en la figura 6 se construyó a partir de una secuencia de 40 imágenes en niveles de gris de tamaño 320 × 385 píxeles adquiridas, situando la punta del endoscopio a 40 mm de altura. La figura 7 muestra la reconstrucción del patrón de prueba utilizado para la evaluación cuantitativa del método de registro. El error promedio obtenido entre imágenes adyacentes fue de 0.09 píxel para la translación en el eje ,0.33 píxel para la translación en el eje , 0.001 para el factor de escala y 0.073° para el ángulo de rotación. El error global de reconstrucción del mosaico (error entre el patrón original y el patrón reconstruido) fue de 6.37 pixeles para el error promedio, con un error mínimo y máximo de 1.1 y 30.1053 pixeles, respectivamente. Para efectuar el registro entre 2 imágenes se requirió aproximadamente 250 iteraciones (aproximadamente 3 minutos utilizando una Pentium III a 850 MHz y 64 Mb de RAM). Cabe mencionar que se si el grado de la transformación en la ecuación (4) se incrementa, se requerirá un mayor número de iteraciones para encontrar la transformación optima , o incluso puede impedir que ésta pueda ser encontrada (proceso de optimización atrapado por un óptimo local).
Discusión
Los errores obtenidos en el proceso de corrección radial para diferentes ángulos de inclinación, indican que este procedimiento puede efectuarse sujetando el endoscopio con la mano (en posición aproximadamente perpendicular con respecto al patrón), con un error final comparable con los reportados por otros métodos que necesitan de un patrón de calibración específico y de una plataforma de posicionamiento, resultando versátil y sencillo.
Es importante observar que aún cuando el error del registro de imágenes adyacentes es relativamente pequeño, este error es acumulativo, por lo que la cartografía final puede requerir de un ajuste final como podría ser un registro en las secuencias de imágenes de manera inversa a la original con lo que se esperaría que los errores fuesen promediados y distribuidos de manera uniforme en toda la secuencia. Adicionalmente, deberá en futuros trabajos tomarse en cuenta que la superficie reconstruida no se encuentra en un plano uniforme. De la misma manera, es importante mencionar que puede reducirse significativamente el tiempo para la reconstrucción al no emplear toda la secuencia de imágenes en el video original.
Conclusiones
Los resultados obtenidos muestran que la metodología propuesta permite la reconstrucción 2D de la superficie de órganos internos. El preprocesamiento empleado permite mejorar las imágenes endoscópicas para su uso en el proceso de registro de imágenes. Como puede verse en los resultados obtenidos, la información lograda después de aplicar el filtro y la corrección de la distribución no homogénea de la iluminación permiten obtener los parámetros de la transformación global para aplicarlos posteriormente en las imágenes originales y efectuar la cartografía correspondiente.
Referencias
Asari K.V., Kumar S., Radhakrishnan D. A New Approach for Nonlinear Distortion Correction in Endoscopic Images Based on Least Squares Method. IEEE Trans. Med. Imaging, 18:345354. 1999 [ Links ]
González R.C. Digital Image Processing. Prentice Hall. 2002. [ Links ]
Miranda Luna R., Blondel W.C.P.M., Daul, Hernandez Y, Posada R., Wolf D. A Simplified Method of Endoscopic Image Distortion, Correction Based on Grey Level Registration. International Conference on Image Processing (ICIP 04). Singapore, Octubre 2004. [ Links ]
Shahidi R., Bax M.R., Maurer C.R., Jhonson J.A., Wilkinson E.P., Wang B., West J.P., Citardi M.J., Manwaring K.H., Khadem R. Implementation, Calibration and Accuracy Testing of ImageEnhanced Endoscopy System. IEEE Trans. Med. Imaging, 11:15241535. 2002. [ Links ]
Smith W.E., Vakil N., Maislin A. Correction of Distortion in Endoscope Images. IEEE Trans. Med. Imaging, 11:117122. 1992. [ Links ]
Tian H., Srikanthan T., Asari K.V., Lam S.K. Study on the Effect of Object to Camera Distance on Polynomial Expansion Coefficients in Barrel Distortion Correction. Fifth IEEE Southwest Symposium on Image Analysis and Interpretation (SSIAI'02). New Mexico. Abril 2002. [ Links ]
Viola P.A. Alignment by Maximization of Mutual Information. Tesis (Doctorado en Ingeniería). Boston, MA, USA. Massachusetts Institute of Technology. 1995. [ Links ]
Semblanza de los autores
Rosebet MirandaLuna. Obtuvo el grado de doctor en automatización y procesamiento de señales por el Institut National Politechnique de Lorraine (Nancy, Francia). Es profesor en el Instituto Tecnológico de la Mixteca. Sus temas de interés son el procesamiento de señales e imágenes con aplicación médica, así como la instrumentación biomédica.
Rubén PosadaGómez. Es doctor en automatización y procesamiento de señales por el Institut National Politechnique de Lorraine (Nancy, Francia). Actualmente es profesor investigador del Instituto Tecnológico de Orizaba. Sus áreas de interés son el procesamiento de imágenes médicas, la instrumentación, el tratamiento de señales y la bioelectrónica
Giner AlorHernández. Es doctor en ciencias en ingeniería eléctrica con especialidad en computación por el Centro de Investigación y de Estudios Avanzados del IPN. Pertenece al Sistema Nacional de Investigadores. Actualmente es profesor investigador del Instituto Tecnológico de Orizaba.
Albino MartínezSibaja. Es maestro en ciencias en ingeniería electrónica con especialidad en electrónica de potencia por el CENIDET, sus líneas de investigación son las aplicaciones de la electrónica en sistemas digitales.
Guillermo CortésRobles. Es doctor en ciencias en sistemas industriales por el Institut National Polytechnique de Toulouse. Pertenece al Sistema Nacional de Investigadores. Actualmente es profesor investigador del Instituto Tecnológico de Orizaba.