<?xml version="1.0" encoding="ISO-8859-1"?><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
<front>
<journal-meta>
<journal-id>1405-7743</journal-id>
<journal-title><![CDATA[Ingeniería, investigación y tecnología]]></journal-title>
<abbrev-journal-title><![CDATA[Ing. invest. y tecnol.]]></abbrev-journal-title>
<issn>1405-7743</issn>
<publisher>
<publisher-name><![CDATA[Universidad Nacional Autónoma de México, Facultad de Ingeniería]]></publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id>S1405-77432011000100009</article-id>
<title-group>
<article-title xml:lang="es"><![CDATA[Registro de imágenes video-endoscópicas para reconstruir cartografías de órganos huecos]]></article-title>
<article-title xml:lang="en"><![CDATA[Video-Endoscopic Image Registration for Cartography Reconstruction of Hollow Organs]]></article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Miranda-Luna]]></surname>
<given-names><![CDATA[R]]></given-names>
</name>
<xref ref-type="aff" rid="A01"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Posada-Gómez]]></surname>
<given-names><![CDATA[R]]></given-names>
</name>
<xref ref-type="aff" rid="A02"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Alor-Hernández]]></surname>
<given-names><![CDATA[G]]></given-names>
</name>
<xref ref-type="aff" rid="A03"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Martínez-Sibaja]]></surname>
<given-names><![CDATA[A]]></given-names>
</name>
<xref ref-type="aff" rid="A04"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Cortés-Robles]]></surname>
<given-names><![CDATA[G]]></given-names>
</name>
<xref ref-type="aff" rid="A05"/>
</contrib>
</contrib-group>
<aff id="A01">
<institution><![CDATA[,Universidad Tecnológica de la Mixteca  ]]></institution>
<addr-line><![CDATA[ Oaxaca]]></addr-line>
</aff>
<aff id="A02">
<institution><![CDATA[,Instituto Tecnológico de Orizaba Departamento de Posgrado e Investigación ]]></institution>
<addr-line><![CDATA[ Veracruz]]></addr-line>
</aff>
<aff id="A03">
<institution><![CDATA[,Instituto Tecnológico de Orizaba Departamento de Posgrado e Investigación ]]></institution>
<addr-line><![CDATA[ Veracruz]]></addr-line>
</aff>
<aff id="A04">
<institution><![CDATA[,Instituto Tecnológico de Orizaba Departamento de Posgrado e Investigación ]]></institution>
<addr-line><![CDATA[ Veracruz]]></addr-line>
</aff>
<aff id="A05">
<institution><![CDATA[,Instituto Tecnológico de Orizaba Departamento de Posgrado e Investigación ]]></institution>
<addr-line><![CDATA[ Veracruz]]></addr-line>
</aff>
<pub-date pub-type="pub">
<day>00</day>
<month>03</month>
<year>2011</year>
</pub-date>
<pub-date pub-type="epub">
<day>00</day>
<month>03</month>
<year>2011</year>
</pub-date>
<volume>12</volume>
<numero>1</numero>
<fpage>83</fpage>
<lpage>91</lpage>
<copyright-statement/>
<copyright-year/>
<self-uri xlink:href="http://www.scielo.org.mx/scielo.php?script=sci_arttext&amp;pid=S1405-77432011000100009&amp;lng=en&amp;nrm=iso"></self-uri><self-uri xlink:href="http://www.scielo.org.mx/scielo.php?script=sci_abstract&amp;pid=S1405-77432011000100009&amp;lng=en&amp;nrm=iso"></self-uri><self-uri xlink:href="http://www.scielo.org.mx/scielo.php?script=sci_pdf&amp;pid=S1405-77432011000100009&amp;lng=en&amp;nrm=iso"></self-uri><abstract abstract-type="short" xml:lang="es"><p><![CDATA[Los sistemas video-endoscó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 video-endoscó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.]]></p></abstract>
<abstract abstract-type="short" xml:lang="en"><p><![CDATA[In examinations of hollow organs, video-endoscopic 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 video-endoscopic 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 far-away 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.]]></p></abstract>
<kwd-group>
<kwd lng="es"><![CDATA[cartografías de órganos internos]]></kwd>
<kwd lng="es"><![CDATA[video-endoscopia]]></kwd>
<kwd lng="es"><![CDATA[registro de imágenes]]></kwd>
<kwd lng="es"><![CDATA[corrección de la distorsión radial]]></kwd>
<kwd lng="en"><![CDATA[internal organs cartography]]></kwd>
<kwd lng="en"><![CDATA[video-endoscopy]]></kwd>
<kwd lng="en"><![CDATA[image registration]]></kwd>
<kwd lng="en"><![CDATA[radial distortion correction]]></kwd>
</kwd-group>
</article-meta>
</front><body><![CDATA[ <p align="center"><font face="verdana" size="4"><b>Registro de im&aacute;genes video&#150;endosc&oacute;picas para reconstruir cartograf&iacute;as de &oacute;rganos huecos</b></font></p>     <p align="center"><font face="verdana" size="2">&nbsp;</font></p>     <p align="center"><font size="3" face="verdana"><b>Video&#150;Endoscopic Image Registration for Cartography Reconstruction of Hollow Organs</b></font></p>     <p align="center"><font face="verdana" size="2">&nbsp;</font></p>     <p align="center"><font face="verdana" size="2"><b>Miranda&#150;Luna R.<sup>1</sup>, Posada&#150;G&oacute;mez R.<sup>2</sup>, Alor&#150;Hern&aacute;ndez G.<sup>3</sup>, Mart&iacute;nez&#150;Sibaja A<sup>4</sup>. y Cort&eacute;s&#150;Robles G.<sup>5</sup></b></font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     <p align="justify"><font size="2" face="verdana"><i><sup>1</sup> Universidad Tecnol&oacute;gica de la Mixteca, Oaxaca. </i>E&#150;mail: <a href="mailto:rmiranda@yahoo.com">rmiranda@yahoo.com</a></font></p>     <p align="justify"><font face="verdana" size="2"><i><sup>2</sup> Departamento de Posgrado e Investigaci&oacute;n Instituto Tecnol&oacute;gico de Orizaba, Veracruz. </i>E&#150;mail: <a href="mailto:rposada@itorizaba.edu.mx">rposada@itorizaba.edu.mx</a></font></p>     <p align="justify"><font face="verdana" size="2"><i><sup>3</sup> Departamento de Posgrado e Investigaci&oacute;n Instituto Tecnol&oacute;gico de Orizaba, Veracruz. </i>E&#150;mail: <a href="mailto:galor@itorizaba.edu.mx">galor@itorizaba.edu.mx</a></font></p>     <p align="justify"><font face="verdana" size="2"><i><sup>4</sup> Departamento de Posgrado e Investigaci&oacute;n Instituto Tecnol&oacute;gico de Orizaba, Veracruz. </i>E&#150;mail: <a href="mailto:asibaja@itorizaba.edu.mx">asibaja@itorizaba.edu.mx</a></font></p>     ]]></body>
<body><![CDATA[<p align="justify"><font face="verdana" size="2"><i><sup>5</sup> Departamento de Posgrado e Investigaci&oacute;n Instituto Tecnol&oacute;gico de Orizaba, Veracruz. </i>E&#150;mail: <a href="mailto:grobles@itorizaba.edu.mx">grobles@itorizaba.edu.mx</a></font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     <p align="justify"><font face="verdana" size="2">Informaci&oacute;n del art&iacute;culo: Recibido: septiembre de 2006.     <br> Reevaluado: noviembre de 2007.    <br> Aceptado: septiembre de 2010.</font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     <p align="justify"><font size="2" face="verdana"><b>Resumen</b></font></p>     <p align="justify"><font size="2" face="verdana">Los sistemas video&#150;endosc&oacute;picos se usan para examinar &oacute;rganos huecos permitiendo observar las capas que recubren las cavidades internas de los &oacute;rganos (superficies internas). Com&uacute;nmente, un especialista debe analizar la secuencia de im&aacute;genes adquiridas durante el examen endosc&oacute;pico, ya sea para confirmar el diagn&oacute;stico efectuado durante la endoscopia o para valorar la evoluci&oacute;n de una patolog&iacute;a. En una secuencia de im&aacute;genes video&#150;endosc&oacute;picas (en donde cada imagen muestra una peque&ntilde;a &aacute;rea de la superficie tisular examinada), puede resultar dif&iacute;cil localizar el &aacute;rea de inter&eacute;s debido a: 1) la gran cantidad de im&aacute;genes que componen la secuencia, 2) la naturaleza secuencial del archivo (que obliga al m&eacute;dico a repasar toda la secuencia para poder encontrar la zona de inter&eacute;s) y 3) la posibilidad de que la zona de inter&eacute;s se encuentre o est&eacute; contenida en varios tramos no contiguos de la secuencia. En este trabajo, se implementa una metodolog&iacute;a que permite la reconstrucci&oacute;n autom&aacute;tica de una cartograf&iacute;a de la regi&oacute;n de inter&eacute;s, permitiendo su visualizaci&oacute;n panor&aacute;mica y eliminando informaci&oacute;n redundante. Adem&aacute;s, se presenta un m&eacute;todo basado en un algoritmo autom&aacute;tico de registro de im&aacute;genes para la reconstrucci&oacute;n de cartograf&iacute;as. Finalmente, se presenta un protocolo de evaluaci&oacute;n y los resultados obtenidos.</font></p>     <p align="justify"><font face="verdana" size="2"><b>Descriptores: </b>cartograf&iacute;as de &oacute;rganos internos, video&#150;endoscopia, registro de im&aacute;genes, correcci&oacute;n de la distorsi&oacute;n radial.</font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     ]]></body>
<body><![CDATA[<p align="justify"><font face="verdana" size="2"><b>Abstract</b></font></p>     <p align="justify"><font face="verdana" size="2"><i>In examinations of hollow organs, video&#150;endoscopic 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 video&#150;endoscopic 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 far&#150;away 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.</i></font></p>     <p align="justify"><font face="verdana" size="2"><b><i>Keywords: </i></b><i>internal organs cartography, video&#150;endoscopy, image registration, radial distortion correction.</i></font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     <p align="justify"><font size="2" face="verdana"><b>Introducci&oacute;n</b></font></p>     <p align="justify"><font face="verdana" size="2">Actualmente, los sistemas video&#150;endosc&oacute;picos utilizados en aplicaciones m&eacute;dicas generan secuencias de im&aacute;genes en los que cada cuadro contiene un fragmento de la superficie total escaneada. Una cartograf&iacute;a (por la analog&iacute;a con el trazado de mapas geogr&aacute;ficos) del &aacute;rea escaneada, constituye as&iacute; una representaci&oacute;n f&aacute;cil de manejar, la cual es &uacute;til en un an&aacute;lisis post&#150;operatorio. En una secuencia de video, las im&aacute;genes adyacentes comparten una superficie com&uacute;n (informaci&oacute;n redundante) que permiten construir el panorama de la superficie observada utilizando las im&aacute;genes de la secuencia. Para construir cartograf&iacute;as a partir de im&aacute;genes sucesivas, se propone efectuar un registro de estas im&aacute;genes y combinarlas para constituir la imagen panor&aacute;mica.</font></p>     <p align="justify"><font face="verdana" size="2">El registro de im&aacute;genes endosc&oacute;picas se ve afectado por tres caracter&iacute;sticas inherentes a la naturaleza de las im&aacute;genes adquiridas (figura 1):</font></p>     <blockquote>       <p align="justify"><font face="verdana" size="2">1) La distorsi&oacute;n geom&eacute;trica no lineal de tipo radial causada por un sistema &oacute;ptico dise&ntilde;ado para brindar un amplio campo de visi&oacute;n distorsiona las im&aacute;genes haciendo que un mismo punto de la escena ocupe posiciones relativas diferentes en dos im&aacute;genes sucesivas a causa de esta distorsi&oacute;n (Tian <i>et al.,</i> 2002), (Smith <i>et al,</i> 1992);</font></p>       <p align="justify"><font face="verdana" size="2">2) El patr&oacute;n fijo de alta frecuencia introducido por la fibra &oacute;ptica de los fibroscopios (endoscopios que utilizan fibras &oacute;pticas para recuperar la imagen) y que corresponde al arreglo de las fibras, se superpone a la imagen generada y;</font></p>       ]]></body>
<body><![CDATA[<p align="justify"><font face="verdana" size="2">3) Un patr&oacute;n de baja frecuencia generado por una iluminaci&oacute;n no uniforme al interior del &oacute;rgano (normalmente m&aacute;s importante en el centro que en los extremos de la imagen).</font></p> </blockquote>     <p align="justify"><font face="verdana" size="2">En este trabajo se presenta un m&eacute;todo para reconstruir una cartograf&iacute;a de un &oacute;rgano explorado a trav&eacute;s de un registro monomodal.</font></p>     <p align="justify"><font face="verdana" size="2">La estructura de este trabajo es la siguiente: primero, se describe la implementaci&oacute;n de un m&eacute;todo de registro que corrige las distorsiones propias al sistema de adquisici&oacute;n y efect&uacute;a la cartograf&iacute;a de una regi&oacute;n de inter&eacute;s. Despu&eacute;s, se describe la metodolog&iacute;a empleada en el desarrollo de la cartograf&iacute;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&oacute;n donde se analizan los resultados y se discuten las posibilidades del sistema propuesto.</font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     <p align="justify"><font size="2" face="verdana"><b>Metodolog&iacute;a para la creaci&oacute;n de la cartograf&iacute;a</b></font></p>     <p align="justify"><font face="verdana" size="2">La metodolog&iacute;a para la realizaci&oacute;n de la cartograf&iacute;a se divide en tres etapas principales:</font></p>     <blockquote>       <p align="justify"><font face="verdana" size="2">1. El acondicionamiento de las im&aacute;genes que corrige la distorsi&oacute;n radial y reduce la presencia de los patrones de alta y baja frecuencia;</font></p>       <p align="justify"><font face="verdana" size="2">2. El registro de las im&aacute;genes previamente corregidas usando un algoritmo autom&aacute;tico basado en la informaci&oacute;n mutua de las im&aacute;genes y;</font></p>       <p align="justify"><font face="verdana" size="2">3. La creaci&oacute;n de la imagen panor&aacute;mica conteniendo todas las im&aacute;genes registradas en un mismo referencial.</font></p> </blockquote>     ]]></body>
<body><![CDATA[<p align="justify"><font face="verdana" size="2">Acondicionamiento de las im&aacute;genes endosc&oacute;picas</font></p>     <p align="justify"><font face="verdana" size="2">Para reducir los problemas que afectan el registro de las im&aacute;genes se procede a acondicionar las im&aacute;genes efectuando primeramente la correcci&oacute;n de la distorsi&oacute;n radial (mediante un algoritmo de registro de im&aacute;genes) y filtrando posteriormente las im&aacute;genes resultantes para eliminar los patrones de alta y baja frecuencia.</font></p>     <p align="justify"><font face="verdana" size="2">A continuaci&oacute;n, se describen las t&eacute;cnicas empleadas en la metodolog&iacute;a propuesta para corregir los principales tipos de distorsi&oacute;n generada por el endoscopio (la distorsi&oacute;n radial y la correcci&oacute;n de patrones de alta y baja frecuencia).</font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     <p align="justify"><font size="2" face="verdana"><b>Correcci&oacute;n de la distorsi&oacute;n radial</b></font></p>     <p align="justify"><font face="verdana" size="2">La transformaci&oacute;n geom&eacute;trica que relaciona las im&aacute;genes adquiridas a trav&eacute;s del endoscopio, se afecta tanto por la posici&oacute;n del endoscopio desde la cual se adquieren como por la &oacute;ptica del endoscopio. La distorsi&oacute;n radial que afecta las im&aacute;genes endosc&oacute;picas (<a href="#f1">figura 1a</a>) se modela con las coordenadas de la proyecci&oacute;n del centro &oacute;ptico en la imagen (<i>Cx, Cy</i>) y la relaci&oacute;n polinomial <i>r'</i> expresada en la siguiente f&oacute;rmula:</font></p>     <p align="center"><font size="2" face="verdana"><img src="/img/revistas/iit/v12n1/a9s1.jpg"></font></p>     <p align="justify"><font face="verdana" size="2">donde</font></p>     <p align="justify"><font face="verdana" size="2"><i>r </i>es la distancia radial entre un punto de coordenadas <i>(x, y) (Cx, Cy),</i></font></p>     <p align="justify"><font face="verdana" size="2"><i>k<sub>i</sub></i> son los coeficientes de correcci&oacute;n usados para obtener una imagen sin distorsi&oacute;n,</font></p>     ]]></body>
<body><![CDATA[<p align="justify"><font face="verdana" size="2"><i>M</i> es el grado del polinomio.</font></p>     <p align="center"><font size="2" face="verdana"><a name="f1"></a></font></p>     <p align="center"><font size="2" face="verdana"><img src="/img/revistas/iit/v12n1/a9f1.jpg"></font></p>     <p align="justify"><font face="verdana" size="2">El n&uacute;mero de par&aacute;metros <i>k<sub>i</sub></i> que deben calcularse depende de la importancia de la distorsi&oacute;n, as&iacute; como de la precisi&oacute;n requerida. En este trabajo, se seleccion&oacute; <i>M =</i> 5 que representa un buen compromiso entre exactitud y tiempo de procesamiento.</font></p>     <p align="justify"><font face="verdana" size="2">El problema de la correcci&oacute;n radial se observa como un problema de calibraci&oacute;n del endoscopio. As&iacute;, sea <i>U</i> una imagen ideal libre de deformaciones, se considera que la imagen <i>V</i> es la imagen <i>U</i> adquirida por medio del endoscopio a calibrar y se afect&oacute; por una transformaci&oacute;n <i>T,</i> a su vez conformada por una transformaci&oacute;n <i>T<sub>i</sub></i> (par&aacute;metros intr&iacute;nsecos del endoscopio) dada por la distorsi&oacute;n radial y una transformaci&oacute;n <i>T<sub>e</sub></i> (par&aacute;metros extr&iacute;nsecos) causada por la posici&oacute;n desde la cual se efectu&oacute; la adquisici&oacute;n de la imagen U, de manera que:</font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9s2.jpg"></font></p>     <p align="justify"><font size="2" face="verdana">De esta manera, efectuando un registro entre las im&aacute;genes <i>U</i> y <i>V</i> se calculan los valores de la matriz <i>T.</i> Los par&aacute;metros de la matriz proyectiva (traslaciones: <i>tx, ty,</i> rotaci&oacute;n: &#952;, escala: &#945;) y del polinomio (<i>Cx</i>, <i>Cy</i> y <i>k<sub>i</sub></i>) son independientes; de esta forma, los par&aacute;metros de correcci&oacute;n radial se a&iacute;slan despu&eacute;s del registro, con el fin de aplicar la correcci&oacute;n radial de manera sistem&aacute;tica a todas las im&aacute;genes adquiridas con el endoscopio, siempre y cuando no se modifique la &oacute;ptica del mismo. Una explicaci&oacute;n detallada del procedimiento de correcci&oacute;n radial o calibraci&oacute;n del endoscopio se encuentra en (Miranda <i>et al,</i> 2004).</font></p>     <p align="justify"><font face="verdana" size="2">Una vez corregida la distorsi&oacute;n radial en toda la secuencia de im&aacute;genes, se procede a reducir la presencia de los artefactos de alta y baja frecuencia generados por el arreglo de fibras &oacute;pticas y la iluminaci&oacute;n no uniforme, respectivamente.</font></p>     <p align="justify"><font face="verdana" size="2">Reducci&oacute;n de la presencia de patrones de alta y baja frecuencia</font></p>     <p align="justify"><font face="verdana" size="2">La reducci&oacute;n de la presencia de los patrones de alta y baja frecuencia se obtiene filtrando la imagen. Para eliminar las altas frecuencias a&ntilde;adidas a las im&aacute;genes por el arreglo de fibras &oacute;pticas (<a href="#f1">figura 1b</a>), se emplea un filtro con una funci&oacute;n gaussiana, cuya funci&oacute;n de transferencia en el dominio de la frecuencia (Gonz&aacute;lez, 2002) se expresa como:</font></p>     ]]></body>
<body><![CDATA[<p align="center"><font size="2" face="verdana"><img src="/img/revistas/iit/v12n1/a9s3.jpg"></font></p>     <p align="justify"><font face="verdana" size="2">donde</font></p>     <p align="justify"><font face="verdana" size="2">D(<i>u, v</i>) es la distancia espacial desde la componente de frecuencia (<i>u, v</i>) hasta el centro del espectro (<i>u =</i> 0, <i>v =</i> 0) y;</font></p>     <p align="justify"><font face="verdana" size="2">D<sub><i>0</i></sub> es la desviaci&oacute;n est&aacute;ndar del filtro.</font></p>     <p align="justify"><font face="verdana" size="2">Una vez analizado el espectro de Fourier correspondiente al patr&oacute;n de fibras, se fija <i>D<sub>0</sub></i> = 4, obteniendo as&iacute; un filtro para bajas. Para corregir el problema de la iluminaci&oacute;n no uniforme, primero se obtiene una versi&oacute;n filtrada de la imagen endosc&oacute;pica conteniendo las componentes de muy baja frecuencia, producidas esencialmente por la distribuci&oacute;n no uniforme de la iluminaci&oacute;n, mediante el filtro gaussiano con <i>D<sub>0</sub></i> = 16. La imagen obtenida se resta de la imagen previamente filtrada con <i>D<sub>0</sub></i> = 4. De esta manera, se tiene una imagen corregida que finalmente se utiliza en el registro de im&aacute;genes.</font></p>     <p align="justify"><font face="verdana" size="2">Registro de im&aacute;genes</font></p>     <p align="justify"><font face="verdana" size="2">El registro de im&aacute;genes se efect&uacute;a con pares de im&aacute;genes en niveles de gris. Dadas dos im&aacute;genes <i>U</i> y <i>V,</i> su registro se expresa como:</font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9s4.jpg"></font></p>     <p align="justify"><font size="2" face="verdana">En la f&oacute;rmula (4), <img src="/img/revistas/iit/v12n1/a9s11.jpg"> representa la transformaci&oacute;n buscada en el espacio de transformaciones <i>T,</i> es decir, la transformaci&oacute;n que optimiza la medida de similitud <i>S</i> calculada con las estructuras homologas extra&iacute;das de las im&aacute;genes <i>U</i> y <i>V</i> con ayuda de algoritmos de segmentaci&oacute;n representados por <i>f<sub>t</sub></i> y <i>f<sub>m</sub>.</i> El algoritmo de registro emplea los niveles de grises como estructuras homologas, evitando la necesidad de efectuar una segmentaci&oacute;n de las im&aacute;genes. En este trabajo, se implement&oacute; como medida de similitud la informaci&oacute;n mutua (Viola, 1995). Sea <i>W</i> el resultado de la transformaci&oacute;n <i>T</i> aplicada sobre la imagen <i>V,</i> la informaci&oacute;n Mutua <i>I</i>(<i>U, W</i>) se expresa por la siguiente f&oacute;rmula:</font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9s5.jpg"></font></p>     ]]></body>
<body><![CDATA[<p align="justify"><font size="2" face="verdana">donde</font></p>     <p align="justify"><font face="verdana" size="2"><i>H(U)</i> y <i>H(W)</i> representan la entrop&iacute;a de las im&aacute;genes <i>U</i> y <i>W,</i> respectivamente dadas por:</font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9s6.jpg"></font></p>     <p align="justify"><font size="2" face="verdana">donde</font></p>     <p align="justify"><font size="2" face="verdana"><i>x<sub>n</sub></i> es el <i>n<sup>esimo</sup></i> valor del espacio &#937;<i>x</i> (&#937;<i>x</i> es el intervalo de niveles de gris &#917; &#91;0&#150;255&#93;).</font></p>     <p align="justify"><font face="verdana" size="2"><i>H(U, W)</i> representa la entrop&iacute;a conjunta de las im&aacute;genes y se define como:</font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9s7.jpg"></font></p>     <p align="justify"><font size="2" face="verdana">donde:</font></p>     <p align="justify"><font size="2" face="verdana"><i>x<sub>n</sub></i> y <i>y<sub>n</sub></i> son los <i>n<sup>esimos</sup></i> valores de &#937;<i>x</i> y &#937;<i>y</i> (niveles de gris: 0 &#150; 255) para las variables <i>x</i> y y.</font></p>     <p align="justify"><font face="verdana" size="2"><i>P</i>(<i>x<sub>n</sub>, y<sub>n</sub></i>) es la probabilidad de concurrencia de los eventos <i>x<sub>n</sub>, y<sub>n</sub>.</i></font></p>     ]]></body>
<body><![CDATA[<p align="justify"><font face="verdana" size="2">La informaci&oacute;n mutua es una medida de similitud robusta que permite conocer la informaci&oacute;n que comparten las partes superpuestas de dos im&aacute;genes caracterizadas &uacute;nicamente por la distribuci&oacute;n de sus niveles de gris. Esta medida de similitud es poco sensible a los problemas de &aacute;reas con poco recubrimiento espacial com&uacute;n al utilizar, tanto la entrop&iacute;a de cada imagen como la entrop&iacute;a conjunta de ambas.</font></p>     <p align="justify"><font face="verdana" size="2">El modelo de transformaci&oacute;n usado considera translaciones, rotaciones, factor de escala y perspectiva (que representan los cambios que pueden darse entre dos im&aacute;genes consecutivas por la manipulaci&oacute;n del endoscopio). Para la optimizaci&oacute;n de los par&aacute;metros de la matriz de transformaci&oacute;n que permiten optimizar la informaci&oacute;n mutua, se implement&oacute; el algoritmo de gradiente descendente estoc&aacute;stico, debido a la existencia de m&iacute;nimos locales en el espacio de b&uacute;squeda generados esencialmente por la interpolaci&oacute;n de los niveles de gris. Este algoritmo es un m&eacute;todo de optimizaci&oacute;n iterativo bien adaptado para superar los peque&ntilde;os &oacute;ptimos locales capaces de impedir la convergencia global. Los valores iniciales de la matriz de transformaci&oacute;n son los de la matriz identidad y se actualizan en cada iteraci&oacute;n, de manera que:</font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9s8.jpg"></font></p>     <p align="justify"><font size="2" face="verdana">donde <i>&#955;</i> es el par&aacute;metro de rapidez de convergencia, el cual se reduce progresivamente a medida que se alcanza el m&aacute;ximo de la informaci&oacute;n mutua.</font></p>     <p align="justify"><font face="verdana" size="2">El registro se efect&uacute;a entre pares de im&aacute;genes consecutivos en la secuencia generada por el de video&#150;endoscopio. Una vez hecho esto, es necesario que todas las im&aacute;genes transformadas despu&eacute;s del registro se sit&uacute;en en un referencial com&uacute;n para construir la cartograf&iacute;a.</font></p>     <p align="justify"><font face="verdana" size="2">Creaci&oacute;n del mosaico</font></p>     <p align="justify"><font face="verdana" size="2">El proceso de registro entre pares de im&aacute;genes consecutivas genera una matriz de transformaci&oacute;n <i>Tr</i> que optimiza la informaci&oacute;n mutua de dos im&aacute;genes adyacentes <i>i</i> e <i>i</i> &#150; 1. Posteriormente, se calcula para cada imagen la transformaci&oacute;n global <i>Tg</i> que permita tener todas las im&aacute;genes en un mismo sistema de coordenadas, mediante la ecuaci&oacute;n:</font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9s9.jpg"></font></p>     <p align="justify"><font size="2" face="verdana">Esto permite situar todas las im&aacute;genes en el referencial de la primera imagen. Despu&eacute;s de aplicar la transformaci&oacute;n global a todas las im&aacute;genes, s&oacute;lo resta combinarlas para formar una sola imagen que constituye el panorama de toda la superficie explorada (cartograf&iacute;a 2D).</font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     ]]></body>
<body><![CDATA[<p align="justify"><font size="2" face="verdana"><b>Pruebas efectuadas</b></font></p>     <p align="justify"><font size="2" face="verdana">Con el fin de validar el m&eacute;todo de creaci&oacute;n de cartograf&iacute;as y las diferentes etapas que lo componen, se efectuaron pruebas para determinar la eficiencia en la correcci&oacute;n de la distorsi&oacute;n radial. As&iacute; tambi&eacute;n, se efectu&oacute; la reconstrucci&oacute;n de una cartograf&iacute;a a partir de im&aacute;genes adquiridas empleando un sistema de posicionamiento que permite conocer los desplazamientos reales del endoscopio y una imagen panor&aacute;mica que simula la superficie de la cavidad interna a explorar. Al tratarse de una imagen bidimensional se efect&uacute;a una comparaci&oacute;n precisa con los resultados de la cartograf&iacute;a reconstruida.</font></p>     <p align="justify"><font face="verdana" size="2">Sistema de posicionamiento</font></p>     <p align="justify"><font face="verdana" size="2">Con el fin de conocer los desplazamientos efectuados por el endoscopio al momento de la adquisici&oacute;n de la secuencia de im&aacute;genes, se construy&oacute; un sistema de posicionamiento (<a href="#f2">figura 2</a>) constituido por una plataforma en la que se montaron 3 tablas microm&eacute;tricas (NEWPORT modelo MUMR8.51). Las tablas se ensamblaron de manera tal que permiten obtener desplazamientos lineales en los 3 ejes (<img src="/img/revistas/iit/v12n1/a9s12.jpg">,<img src="/img/revistas/iit/v12n1/a9s13.jpg">,<img src="/img/revistas/iit/v12n1/a9s14.jpg">). El equipo de posicionamiento permite controlar de manera precisa los desplazamientos manualmente efectuados del endoscopio dentro de un volumen de 50 &times; 50 &times; 40 mm<sup>3</sup>.</font></p>     <p align="center"><font size="2" face="verdana"><a name="f2"></a></font></p>     <p align="center"><font size="2" face="verdana"><img src="/img/revistas/iit/v12n1/a9f2.jpg"></font></p>     <p align="justify"><font face="verdana" size="2">Para la adquisici&oacute;n de las im&aacute;genes, se emple&oacute; una c&aacute;mara Leritier 760 ULL y un fibroscopio (FUJINONFG&#150;100FP), el cual se fij&oacute; al sistema de posicionamiento con ayuda de una gu&iacute;a (<a href="#f2">figura 2b</a>). De esta manera, la punta del endoscopio se alinea con el borde de la gu&iacute;a. El sistema permite colocar la punta del endoscopio a cualquier distancia entre 10 mm y 60 mm con respecto a la escena. Adem&aacute;s, el sistema permite girar el endoscopio para simular adquisiciones perspectivas. En la parte inferior del sistema se sit&uacute;a la imagen a adquirir (fotograf&iacute;a de la regi&oacute;n a explorar o patr&oacute;n para evaluar cuantitativamente la calibraci&oacute;n). La imagen a adquirir puede rotarse alrededor del eje <img src="/img/revistas/iit/v12n1/a9s14.jpg"> a pasos de 2&deg;, debido a un pivote que se encuentra en el centro de la base de adquisici&oacute;n y a un patr&oacute;n con marcas colocado debajo de la imagen.</font></p>     <p align="justify"><font face="verdana" size="2">El sistema descrito permite efectuar, tanto las pruebas de la calibraci&oacute;n como de la reconstrucci&oacute;n de la cartograf&iacute;a.</font></p>     <p align="justify"><font face="verdana" size="2">Calibraci&oacute;n del endoscopio (correcci&oacute;n de la distorsi&oacute;n radial)</font></p>     <p align="justify"><font face="verdana" size="2">Para evaluar el m&eacute;todo de correcci&oacute;n de la distorsi&oacute;n radial, se utiliz&oacute; una imagen patr&oacute;n (<a href="#f3">figura 3a</a>). La imagen patr&oacute;n, se adquiri&oacute; situando el endoscopio perpendicular a &eacute;sta. El posicionamiento del endoscopio se efectu&oacute; empleando el sistema de posicionamiento antes descrito.</font></p>     ]]></body>
<body><![CDATA[<p align="center"><font size="2" face="verdana"><a name="f3"></a></font></p>     <p align="center"><font size="2" face="verdana"><img src="/img/revistas/iit/v12n1/a9f3.jpg"></font></p>     <p align="justify"><font face="verdana" size="2">Para la correcci&oacute;n de distorsi&oacute;n radial, se inicializan los par&aacute;metros del polinomio <i>k<sub>i</sub></i> = 0 y <i>Cx</i> y <i>Cy</i> toman las coordenadas del centro de la imagen. La matriz proyectiva, se inicializa con <i>tx</i> = <i>ty</i> = 1, &#952; = 0 y &#945; = 0.85.</font></p>      <p align="justify"><font face="verdana" size="2">El error en la correcci&oacute;n de la distorsi&oacute;n radial se eval&uacute;a calculando la distancia euclidiana entre puntos de referencia identificados, tanto en la imagen patr&oacute;n sin distorsi&oacute;n como en la imagen endosc&oacute;pica binarizada. Estos puntos se definen por los centros de los cuadros blancos.</font></p>     <p align="justify"><font face="verdana" size="2">Con el fin de evaluar la influencia del &aacute;ngulo de perspectiva en el m&eacute;todo de calibraci&oacute;n, se realizaron (con ayuda del sistema de posicionamiento) adquisiciones de la imagen patr&oacute;n a diferentes &aacute;ngulos (0&deg;, 5&deg; y 10&deg;), para generar las perspectivas con respecto a la perpendicular al plano del patr&oacute;n.</font></p>     <p align="justify"><font face="verdana" size="2">Debido a que el campo de visi&oacute;n del endoscopio se calcul&oacute; en 449 p&iacute;xeles, finalmente se normalizan los errores promedio (<i>e<sub>prom</sub></i>) y m&aacute;ximos (<i>e<sub>m&aacute;x</sub></i>) que representan el comportamiento t&iacute;pico del sistema y el peor escenario, respectivamente (Shahidi <i>et al,</i> 2002).</font></p>     <p align="justify"><font face="verdana" size="2">Adquisici&oacute;n y reconstrucci&oacute;n de una secuencia</font></p>     <p align="justify"><font face="verdana" size="2">Empleando el sistema de posicionamiento se adquiri&oacute; una secuencia de im&aacute;genes, simulando una exploraci&oacute;n de endoscopia. Primeramente, se coloc&oacute; sobre la plataforma la imagen a adquirir. Durante la adquisici&oacute;n, el endoscopio se situ&oacute; a una altura de 40 mm con respecto a la plataforma. Posteriormente, se adquirieron 20 im&aacute;genes, desplazando progresivamente la punta del endoscopio de izquierda a derecha (direcci&oacute;n positiva del eje horizontal) (<a href="#f4a">figura 4a</a>). Cada imagen se adquiri&oacute; despu&eacute;s de un desplazamiento de 2.5 mm. Este valor representa el desplazamiento m&aacute;ximo de la punta del endoscopio durante un examen de endoscopia cl&iacute;nica (se determin&oacute; de manera emp&iacute;rica tomando en cuenta la frecuencia de adquisici&oacute;n de 25 cuadros por segundo). Esta secuencia permite probar la respuesta del algoritmo, ante un desplazamiento puro.</font></p>     <p align="center"><font size="2" face="verdana"><a name="f4a" id="f4a"></a></font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9f4a.jpg"></font></p>     ]]></body>
<body><![CDATA[<p align="center"><font face="verdana" size="2"><a name="f4b"></a></font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9f4b.jpg"></font></p>     <p align="justify"><font face="verdana" size="2">A partir de la &uacute;ltima posici&oacute;n (imagen <i>I<sub>20</sub></i>) se efect&uacute;a una segunda etapa de 10 adquisiciones, en las que la punta del endoscopio se desplaza en direcci&oacute;n diagonal sobre el plano (<img src="/img/revistas/iit/v12n1/a9s12.jpg">, <img src="/img/revistas/iit/v12n1/a9s13.jpg">) con incrementos de &#150;2.5 y &#150;2.5 mm, respectivamente, a lo largo de cada eje. Estos desplazamientos se acompa&ntilde;aron de rotaciones de la imagen alrededor del eje <img src="/img/revistas/iit/v12n1/a9s14.jpg"> de +2&deg; entre cada adquisici&oacute;n. Finalmente, se efectuaron 10 adquisiciones desplazando el endoscopio de forma diagonal sobre el plano (<img src="/img/revistas/iit/v12n1/a9s12.jpg">, <img src="/img/revistas/iit/v12n1/a9s13.jpg">) con incrementos de &#150;2.5 y +2.5 mm, respectivamente a lo largo de cada eje, al mismo tiempo se efectuaron rotaciones alrededor de los ejes <img src="/img/revistas/iit/v12n1/a9s13.jpg"> y <img src="/img/revistas/iit/v12n1/a9s14.jpg"> de 2&deg; y &#150;2&deg;, respectivamente, produci&eacute;ndose as&iacute; un efecto de perspectiva en las im&aacute;genes.</font></p>     <p align="justify"><font face="verdana" size="2">Para evaluar el m&eacute;todo propuesto, se efectu&oacute; primeramente la adquisici&oacute;n de una secuencia utilizando una fotograf&iacute;a digital de la superficie interna de la vejiga de un puerco (<a href="#f4a">figura 4a</a>) y posteriormente un dise&ntilde;o constituido de cuadros alineados en forma de rejilla con una separaci&oacute;n de 10 mm entre cada l&iacute;nea (<a href="#f4b">figura 4b</a>). Las 40 im&aacute;genes de la primera adquisici&oacute;n (fotograf&iacute;a de la vejiga) se registraron con el fin de obtener los par&aacute;metros de transformaci&oacute;n. Los par&aacute;metros de transformaci&oacute;n obtenidos con la primera prueba se aplicaron a la segunda secuencia (patr&oacute;n de puntos), a fin de reconstruir la cartograf&iacute;a de dicho patr&oacute;n. La calidad de la cartograf&iacute;a de la primera adquisici&oacute;n (vejiga) se eval&uacute;a as&iacute;, comparando el patr&oacute;n de puntos reconstruido contra el patr&oacute;n de puntos original. El error se calcula como la suma de las distancias euclidianas entre los centros de puntos correspondientes del patr&oacute;n reconstruido y del patr&oacute;n original, esto se expresa mediante la siguiente ecuaci&oacute;n:</font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9s10.jpg"></font></p>     <p align="justify"><font size="2" face="verdana">donde</font></p>     <p align="justify"><font face="verdana" size="2"><i>N</i> representa los puntos a comparar, </font></p>     <p align="justify"><font face="verdana" size="2"><i>x<sub>i</sub></i> y <i>y<sub>i</sub></i> son los valores absolutos de los errores en las direcciones vertical y horizontal entre el centro del i&eacute;simo punto en el patr&oacute;n original y su correspondiente en el mosaico reconstruido.</font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     <p align="justify"><font size="2" face="verdana"><b>Resultados</b></font></p>     ]]></body>
<body><![CDATA[<p align="justify"><font size="2" face="verdana">Correcci&oacute;n de la distorsi&oacute;n radial</font></p>     <p align="justify"><font face="verdana" size="2">La <a href="#t1">tabla 1</a> muestra los valores de los par&aacute;metros calculados por el algoritmo de calibraci&oacute;n para la imagen patr&oacute;n adquirida desde 3 &aacute;ngulos diferentes (<a href="#f3">figura 3</a>). Los valores de <i>k<sub>4</sub></i> y <i>k<sub>5</sub></i> son similares independientemente del &aacute;ngulo de adquisici&oacute;n, lo cual confirma la irrelevancia del &aacute;ngulo de perspectiva. Se realiz&oacute; una comparaci&oacute;n de los errores normalizados obtenidos para todas las im&aacute;genes con el fin de evaluar con mayor exactitud la influencia del &aacute;ngulo de inclinaci&oacute;n obteniendo un error promedio <i>e<sub>prom</sub></i> de 0.21 a 0.34%. Los valores del error m&aacute;ximo <i>e<sub>m&aacute;x</sub></i> son de 2 a 4 veces mayores que <i>e<sub>prom</sub>.</i> Sin embargo, los errores obtenidos son comparables con otros resultados (Asari <i>et al.,</i> 1999).</font></p>     <p align="center"><font size="2" face="verdana"><a name="t1"></a></font></p>     <p align="center"><font size="2" face="verdana"><img src="/img/revistas/iit/v12n1/a9t1.jpg"></font></p>     <p align="justify"><font face="verdana" size="2">De esta manera, el m&eacute;todo propuesto s&oacute;lo requiere que el endoscopio se posicione de manera perpendicular al patr&oacute;n de la <a href="#f3">figura 3a</a> para efectuar la calibraci&oacute;n, el error final de la reconstrucci&oacute;n de la cartograf&iacute;a no var&iacute;a, a&uacute;n cuando se haya adquirido la imagen en un rango de 10&deg;.</font></p>     <p align="justify"><font face="verdana" size="2">Eliminaci&oacute;n de los patrones de alta y baja frecuencia</font></p>     <p align="justify"><font face="verdana" size="2">Los resultados del filtro aplicado a las im&aacute;genes adquiridas para eliminar las componentes correspondientes al patr&oacute;n de fibras y la iluminaci&oacute;n no homog&eacute;nea se muestran en la <a href="#f5">figura 5</a>. Aunque la calidad visual de la imagen despu&eacute;s del filtrado no resulta apropiada para un an&aacute;lisis visual, la informaci&oacute;n que contiene es adecuada para su uso en el procedimiento de registro, por lo que una vez calculado los par&aacute;metros de la transformaci&oacute;n pueden aplicarse a las im&aacute;genes sin filtrar y permitir una mejor visualizaci&oacute;n de la cartograf&iacute;a.</font></p>     <p align="center"><font size="2" face="verdana"><a name="f5"></a></font></p>     <p align="center"><font size="2" face="verdana"><img src="/img/revistas/iit/v12n1/a9f5.jpg"></font></p>     <p align="justify"><font face="verdana" size="2">Reconstrucci&oacute;n de la cartograf&iacute;a (Mosaico)</font></p>     ]]></body>
<body><![CDATA[<p align="justify"><font face="verdana" size="2">El mosaico mostrado en la <a href="#f6">figura 6</a> se construy&oacute; a partir de una secuencia de 40 im&aacute;genes en niveles de gris de tama&ntilde;o 320 &#215; 385 p&iacute;xeles adquiridas, situando la punta del endoscopio a 40 mm de altura. La <a href="#f7">figura 7</a> muestra la reconstrucci&oacute;n del patr&oacute;n de prueba utilizado para la evaluaci&oacute;n cuantitativa del m&eacute;todo de registro. El error promedio obtenido entre im&aacute;genes adyacentes fue de 0.09 p&iacute;xel para la translaci&oacute;n en el eje <img src="/img/revistas/iit/v12n1/a9s12.jpg">,&#150;0.33 p&iacute;xel para la translaci&oacute;n en el eje <img src="/img/revistas/iit/v12n1/a9s13.jpg">, 0.001 para el factor de escala y 0.073&deg; para el &aacute;ngulo de rotaci&oacute;n. El error global de reconstrucci&oacute;n del mosaico (error entre el patr&oacute;n original y el patr&oacute;n reconstruido) fue de 6.37 pixeles para el error promedio, con un error m&iacute;nimo y m&aacute;ximo de 1.1 y 30.1053 pixeles, respectivamente. Para efectuar el registro entre 2 im&aacute;genes se requiri&oacute; 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&oacute;n en la ecuaci&oacute;n (4) se incrementa, se requerir&aacute; un mayor n&uacute;mero de iteraciones para encontrar la transformaci&oacute;n optima <img src="/img/revistas/iit/v12n1/a9s11.jpg">, o incluso puede impedir que &eacute;sta pueda ser encontrada (proceso de optimizaci&oacute;n atrapado por un &oacute;ptimo local).</font></p>     <p align="center"><font size="2" face="verdana"><a name="f6"></a></font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9f6.jpg"></font></p>     <p align="center"><font face="verdana" size="2"><a name="f7"></a></font></p>     <p align="center"><font face="verdana" size="2"><img src="/img/revistas/iit/v12n1/a9f7.jpg"></font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     <p align="justify"><font size="2" face="verdana"><b>Discusi&oacute;n</b></font></p>     <p align="justify"><font face="verdana" size="2">Los errores obtenidos en el proceso de correcci&oacute;n radial para diferentes &aacute;ngulos de inclinaci&oacute;n, indican que este procedimiento puede efectuarse sujetando el endoscopio con la mano (en posici&oacute;n aproximadamente perpendicular con respecto al patr&oacute;n), con un error final comparable con los reportados por otros m&eacute;todos que necesitan de un patr&oacute;n de calibraci&oacute;n espec&iacute;fico y de una plataforma de posicionamiento, resultando vers&aacute;til y sencillo.</font></p>     <p align="justify"><font face="verdana" size="2">Es importante observar que a&uacute;n cuando el error del registro de im&aacute;genes adyacentes es relativamente peque&ntilde;o, este error es acumulativo, por lo que la cartograf&iacute;a final puede requerir de un ajuste final como podr&iacute;a ser un registro en las secuencias de im&aacute;genes de manera inversa a la original con lo que se esperar&iacute;a que los errores fuesen promediados y distribuidos de manera uniforme en toda la secuencia. Adicionalmente, deber&aacute; 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&oacute;n al no emplear toda la secuencia de im&aacute;genes en el video original.</font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     ]]></body>
<body><![CDATA[<p align="justify"><font size="2" face="verdana"><b>Conclusiones</b></font></p>     <p align="justify"><font size="2" face="verdana">Los resultados obtenidos muestran que la metodolog&iacute;a propuesta permite la reconstrucci&oacute;n 2D de la superficie de &oacute;rganos internos. El pre&#150;procesamiento empleado permite mejorar las im&aacute;genes endosc&oacute;picas para su uso en el proceso de registro de im&aacute;genes. Como puede verse en los resultados obtenidos, la informaci&oacute;n lograda despu&eacute;s de aplicar el filtro y la correcci&oacute;n de la distribuci&oacute;n no homog&eacute;nea de la iluminaci&oacute;n permiten obtener los par&aacute;metros de la transformaci&oacute;n global para aplicarlos posteriormente en las im&aacute;genes originales y efectuar la cartograf&iacute;a correspondiente.</font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     <p align="justify"><font size="2" face="verdana"><b>Referencias</b></font></p>     <!-- ref --><p align="justify"><font size="2" face="verdana">Asari K.V., Kumar S., Radhakrishnan D. A New Approach for Nonlinear Distortion Correction in Endoscopic Images Based on Least Squares Method. <i>IEEE Trans. Med. Imaging,</i> 18:345&#150;354. 1999</font>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=4283005&pid=S1405-7743201100010000900001&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --><!-- ref --><p align="justify"><font face="verdana" size="2">Gonz&aacute;lez R.C. <i>Digital Image Processing.</i> Prentice Hall. 2002.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=4283006&pid=S1405-7743201100010000900002&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></font></p>     <!-- ref --><p align="justify"><font face="verdana" size="2">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.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=4283008&pid=S1405-7743201100010000900003&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></font></p>     <!-- ref --><p align="justify"><font face="verdana" size="2">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 Image&#150;Enhanced Endoscopy System. <i>IEEE Trans. Med. Imaging, </i>11:1524&#150;1535. 2002.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=4283010&pid=S1405-7743201100010000900004&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></font></p>     <!-- ref --><p align="justify"><font face="verdana" size="2">Smith W.E., Vakil N., Maislin A. Correction of Distortion in Endoscope Images. <i>IEEE Trans. Med. Imaging,</i> 11:117&#150;122. 1992.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=4283012&pid=S1405-7743201100010000900005&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></font></p>     <!-- ref --><p align="justify"><font face="verdana" size="2">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.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=4283014&pid=S1405-7743201100010000900006&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></font></p>     <!-- ref --><p align="justify"><font face="verdana" size="2">Viola P.A. Alignment by Maximization of Mutual Information. Tesis (Doctorado en Ingenier&iacute;a). Boston, MA, USA. Massachusetts Institute of Technology. 1995.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=4283016&pid=S1405-7743201100010000900007&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></font></p>     <p align="justify"><font face="verdana" size="2">&nbsp;</font></p>     <p align="justify"><font size="2" face="verdana"><b>Semblanza de los autores</b></font></p>     <p align="justify"><font size="2" face="verdana"><i>Rosebet Miranda&#150;Luna.</i> Obtuvo el grado de doctor en automatizaci&oacute;n y procesamiento de se&ntilde;ales por el Institut National Politechnique de Lorraine (Nancy, Francia). Es profesor en el Instituto Tecnol&oacute;gico de la Mixteca. Sus temas de inter&eacute;s son el procesamiento de se&ntilde;ales e im&aacute;genes con aplicaci&oacute;n m&eacute;dica, as&iacute; como la instrumentaci&oacute;n biom&eacute;dica.</font></p>     ]]></body>
<body><![CDATA[<p align="justify"><font face="verdana" size="2"><i>Rub&eacute;n Posada&#150;G&oacute;mez.</i> Es doctor en automatizaci&oacute;n y procesamiento de se&ntilde;ales por el Institut National Politechnique de Lorraine (Nancy, Francia). Actualmente es profesor investigador del Instituto Tecnol&oacute;gico de Orizaba. Sus &aacute;reas de inter&eacute;s son el procesamiento de im&aacute;genes m&eacute;dicas, la instrumentaci&oacute;n, el tratamiento de se&ntilde;ales y la bioelectr&oacute;nica</font></p>     <p align="justify"><font face="verdana" size="2"><i>Giner Alor&#150;Hern&aacute;ndez.</i> Es doctor en ciencias en ingenier&iacute;a el&eacute;ctrica con especialidad en computaci&oacute;n por el Centro de Investigaci&oacute;n y de Estudios Avanzados del IPN. Pertenece al Sistema Nacional de Investigadores. Actualmente es profesor investigador del Instituto Tecnol&oacute;gico de Orizaba.</font></p>     <p align="justify"><font face="verdana" size="2"><i>Albino Mart&iacute;nez&#150;Sibaja.</i> Es maestro en ciencias en ingenier&iacute;a electr&oacute;nica con especialidad en electr&oacute;nica de potencia por el CENIDET, sus l&iacute;neas de investigaci&oacute;n son las aplicaciones de la electr&oacute;nica en sistemas digitales.</font></p>     <p align="justify"><font face="verdana" size="2"><i>Guillermo Cort&eacute;s&#150;Robles.</i> 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&oacute;gico de Orizaba.</font></p>      ]]></body><back>
<ref-list>
<ref id="B1">
<nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Asari]]></surname>
<given-names><![CDATA[K.V.]]></given-names>
</name>
<name>
<surname><![CDATA[Kumar]]></surname>
<given-names><![CDATA[S]]></given-names>
</name>
<name>
<surname><![CDATA[Radhakrishnan]]></surname>
<given-names><![CDATA[D. A]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA[New Approach for Nonlinear Distortion Correction in Endoscopic Images Based on Least Squares Method]]></article-title>
<source><![CDATA[IEEE Trans. Med. Imaging]]></source>
<year>1999</year>
<volume>18</volume>
<page-range>345-354</page-range></nlm-citation>
</ref>
<ref id="B2">
<nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[González]]></surname>
<given-names><![CDATA[R.C.]]></given-names>
</name>
</person-group>
<source><![CDATA[Digital Image Processing]]></source>
<year>2002</year>
<publisher-name><![CDATA[Prentice Hall]]></publisher-name>
</nlm-citation>
</ref>
<ref id="B3">
<nlm-citation citation-type="confpro">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Miranda Luna]]></surname>
<given-names><![CDATA[R]]></given-names>
</name>
<name>
<surname><![CDATA[Blondel]]></surname>
<given-names><![CDATA[W.C.P.M.]]></given-names>
</name>
<name>
<surname><![CDATA[Daul, Hernandez]]></surname>
<given-names><![CDATA[Y]]></given-names>
</name>
<name>
<surname><![CDATA[Posada]]></surname>
<given-names><![CDATA[R]]></given-names>
</name>
<name>
<surname><![CDATA[Wolf]]></surname>
<given-names><![CDATA[D]]></given-names>
</name>
</person-group>
<source><![CDATA[A Simplified Method of Endoscopic Image Distortion, Correction Based on Grey Level Registration]]></source>
<year></year>
<conf-name><![CDATA[ International Conference on Image Processing]]></conf-name>
<conf-date>2004</conf-date>
<conf-loc>Singapore </conf-loc>
</nlm-citation>
</ref>
<ref id="B4">
<nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Shahidi]]></surname>
<given-names><![CDATA[R]]></given-names>
</name>
<name>
<surname><![CDATA[Bax]]></surname>
<given-names><![CDATA[M.R.]]></given-names>
</name>
<name>
<surname><![CDATA[Maurer]]></surname>
<given-names><![CDATA[C.R.]]></given-names>
</name>
<name>
<surname><![CDATA[Jhonson]]></surname>
<given-names><![CDATA[J.A.]]></given-names>
</name>
<name>
<surname><![CDATA[Wilkinson]]></surname>
<given-names><![CDATA[E.P.]]></given-names>
</name>
<name>
<surname><![CDATA[Wang]]></surname>
<given-names><![CDATA[B]]></given-names>
</name>
<name>
<surname><![CDATA[West]]></surname>
<given-names><![CDATA[J.P.]]></given-names>
</name>
<name>
<surname><![CDATA[Citardi]]></surname>
<given-names><![CDATA[M.J.]]></given-names>
</name>
<name>
<surname><![CDATA[Manwaring]]></surname>
<given-names><![CDATA[K.H.]]></given-names>
</name>
<name>
<surname><![CDATA[Khadem]]></surname>
<given-names><![CDATA[R]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA[Implementation, Calibration and Accuracy Testing of Image-Enhanced Endoscopy System]]></article-title>
<source><![CDATA[IEEE Trans. Med. Imaging]]></source>
<year>2002</year>
<volume>11</volume>
<page-range>1524-1535</page-range></nlm-citation>
</ref>
<ref id="B5">
<nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Smith]]></surname>
<given-names><![CDATA[W.E.]]></given-names>
</name>
<name>
<surname><![CDATA[Vakil]]></surname>
<given-names><![CDATA[N]]></given-names>
</name>
<name>
<surname><![CDATA[Maislin]]></surname>
<given-names><![CDATA[A]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA[Correction of Distortion in Endoscope Images]]></article-title>
<source><![CDATA[IEEE Trans. Med. Imaging]]></source>
<year>1992</year>
<volume>11</volume>
<page-range>117-122</page-range></nlm-citation>
</ref>
<ref id="B6">
<nlm-citation citation-type="confpro">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Tian]]></surname>
<given-names><![CDATA[H]]></given-names>
</name>
<name>
<surname><![CDATA[Srikanthan]]></surname>
<given-names><![CDATA[T]]></given-names>
</name>
<name>
<surname><![CDATA[Asari]]></surname>
<given-names><![CDATA[K.V.]]></given-names>
</name>
<name>
<surname><![CDATA[Lam]]></surname>
<given-names><![CDATA[S.K.]]></given-names>
</name>
</person-group>
<source><![CDATA[Study on the Effect of Object to Camera Distance on Polynomial Expansion Coefficients in Barrel Distortion Correction]]></source>
<year></year>
<conf-name><![CDATA[Fifth IEEE Southwest Symposium on Image Analysis and Interpretation]]></conf-name>
<conf-date>2002</conf-date>
<conf-loc>New Mexico </conf-loc>
</nlm-citation>
</ref>
<ref id="B7">
<nlm-citation citation-type="">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Viola]]></surname>
<given-names><![CDATA[P.A.]]></given-names>
</name>
</person-group>
<source><![CDATA[Alignment by Maximization of Mutual Information]]></source>
<year></year>
</nlm-citation>
</ref>
</ref-list>
</back>
</article>
