(15 intermediate revisions by 2 users not shown) | |||
Line 25: | Line 25: | ||
mencionar la medición del interior del cilindro maestro que forma parte del sistema de frenado de los autos, el escaneo 3D del | mencionar la medición del interior del cilindro maestro que forma parte del sistema de frenado de los autos, el escaneo 3D del | ||
interior del cuerpo humano, la medición en el interior de tuberías, la medición del interior de botellas hechas de plástico, entre | interior del cuerpo humano, la medición en el interior de tuberías, la medición del interior de botellas hechas de plástico, entre | ||
− | otros (Figura 1). El presente trabajo propone el uso de un sistema de medición capaz de obtener la posición de la luz activa | + | otros ([[#img-1|Figura 1]]). El presente trabajo propone el uso de un sistema de medición capaz de obtener la posición de la luz activa |
para poder dimensionar el interior del objeto, este sistema consta de dos módulos: un sistema de dimensionamiento 2D de visión | para poder dimensionar el interior del objeto, este sistema consta de dos módulos: un sistema de dimensionamiento 2D de visión | ||
activa, y un sistema de posicionamiento que calcula la profundidad del sistema de dimensionamiento 2D de visión activa. La | activa, y un sistema de posicionamiento que calcula la profundidad del sistema de dimensionamiento 2D de visión activa. La | ||
Line 31: | Line 31: | ||
obtener una medición de mayor precisión y exactitud. | obtener una medición de mayor precisión y exactitud. | ||
+ | <div id='img-1'></div> | ||
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 65%;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 65%;" | ||
|- | |- | ||
Line 42: | Line 43: | ||
consiste en la sustitución del sistema de iluminación original del fibroscopio por un sistema que proyecta una línea de luz. Esta | consiste en la sustitución del sistema de iluminación original del fibroscopio por un sistema que proyecta una línea de luz. Esta | ||
línea de luz revela una sección transversal del interior del objeto bajo estudio. La línea de luz se obtiene mediante un LED de | línea de luz revela una sección transversal del interior del objeto bajo estudio. La línea de luz se obtiene mediante un LED de | ||
− | potencia encerrado entre 2 placas y unido por una extensión a la punta del fibroscopio ( | + | potencia encerrado entre 2 placas y unido por una extensión a la punta del fibroscopio ([[#img-2|Figuras 2]](a) y (b)). Las imágenes de |
diferentes secciones transversales del objeto son capturadas por una cámara conectada al ocular del fibroscopio por medio de un | diferentes secciones transversales del objeto son capturadas por una cámara conectada al ocular del fibroscopio por medio de un | ||
− | video acoplador (Figura 2 (c)). | + | video acoplador ([[#img-2|Figura 2]](c)). |
− | + | <div id='img-2'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 70%;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 70%;" | ||
|- | |- | ||
Line 123: | Line 124: | ||
posicionamiento. | posicionamiento. | ||
− | La Tabla 1 muestra las principales características que presenta nuestro sistema para escanear el interior de piezas comparado con | + | La [[#tab-1|Tabla 1]] muestra las principales características que presenta nuestro sistema para escanear el interior de piezas comparado con |
otros sistemas. Los artículos presentados en esta tabla utilizan luz estructurada con uno o varios planos, o proyectan luz con | otros sistemas. Los artículos presentados en esta tabla utilizan luz estructurada con uno o varios planos, o proyectan luz con | ||
formas conocidas. Los trabajos que no usan luz estructurada, [23], detectan características en la imagen que son seguidas mientras | formas conocidas. Los trabajos que no usan luz estructurada, [23], detectan características en la imagen que son seguidas mientras | ||
Line 133: | Line 134: | ||
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;font-size: 75%;"> | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;font-size: 75%;"> | ||
'''Tabla 1'''. Principales características de nuestro sistema de medición de interiores comparado con otras propuestas</div> | '''Tabla 1'''. Principales características de nuestro sistema de medición de interiores comparado con otras propuestas</div> | ||
− | + | <div id='tab-1'></div> | |
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:50%;" | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:50%;" | ||
Line 317: | Line 318: | ||
El sistema cuenta con una cámara CCD monocromática (C1) con una resolución de 2448<math>\times></math>2046 pixeles y un sensor de 1/1.45" que se | El sistema cuenta con una cámara CCD monocromática (C1) con una resolución de 2448<math>\times></math>2046 pixeles y un sensor de 1/1.45" que se | ||
− | conecta al ocular del fibroscopio utilizando un video acoplador óptico (VA) | + | conecta al ocular del fibroscopio utilizando un video acoplador óptico (VA) ([[#img-2|Figura 2]]). Un fibroscopio de 7 mm de diámetro (FIBRO) |
con una longitud efectiva de 1320 mm, dispone de una distancia mínima de trabajo de 30.48 mm, y un campo de visión de | con una longitud efectiva de 1320 mm, dispone de una distancia mínima de trabajo de 30.48 mm, y un campo de visión de | ||
aproximadamente 50 grados alrededor del eje vertical y 40 grados alrededor del eje horizontal. La configuración del sistema del | aproximadamente 50 grados alrededor del eje vertical y 40 grados alrededor del eje horizontal. La configuración del sistema del | ||
− | sistema de medición propuesto se muestra en la Figura 3. | + | sistema de medición propuesto se muestra en la [[#img-3|Figura 3]]. |
− | + | <div id='img-3'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: 80%;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: 80%;" | ||
|- | |- | ||
Line 330: | Line 331: | ||
− | La Figura 4 muestra la cámara endoscópica (mC2), que utiliza tecnología optoelectrónica para inspeccionar áreas de difícil acceso. | + | La [[#img-4|Figura 4]] muestra la cámara endoscópica (mC2), que utiliza tecnología optoelectrónica para inspeccionar áreas de difícil acceso. |
La cámara endoscópica (mC2) tiene un diámetro de 7 mm, las dimensiones son de 510.54 mm <math>\times</math> 7.62 mm <math>\times</math> 7.62 mm, la resolución es de | La cámara endoscópica (mC2) tiene un diámetro de 7 mm, las dimensiones son de 510.54 mm <math>\times</math> 7.62 mm <math>\times</math> 7.62 mm, la resolución es de | ||
640x480 pixeles (JPEG), adquiere a 30 cuadros por segundo (FPS) y usa una luz LED ajustable en la punta de la cámara para iluminar | 640x480 pixeles (JPEG), adquiere a 30 cuadros por segundo (FPS) y usa una luz LED ajustable en la punta de la cámara para iluminar | ||
Line 341: | Line 342: | ||
cámara 2 (mC2) localizan el plano de luz al interior del objeto. | cámara 2 (mC2) localizan el plano de luz al interior del objeto. | ||
+ | <div id='img-4'></div> | ||
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;" | ||
|- | |- | ||
Line 358: | Line 360: | ||
===3.1 Visión Activa=== | ===3.1 Visión Activa=== | ||
− | El diagrama de flujo de la Figura 5 muestra la metodología para realizar el dimensionamiento o escaneo 2D. En él se muestran tres | + | El diagrama de flujo de la [[#img-5|Figura 5]] muestra la metodología para realizar el dimensionamiento o escaneo 2D. En él se muestran tres |
− | grandes bloques de operaciones delimitados con línea punteada. Procedimiento de calibración, trazabilidad, posición y orientación | + | grandes bloques de operaciones delimitados con línea punteada. Procedimiento de calibración, trazabilidad, posición y orientación (POSE). |
− | (POSE). | + | <div id='img-5'></div> |
− | + | ||
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;" | ||
|- | |- | ||
Line 377: | Line 378: | ||
El último bloque del diagrama describe como el contorno de las secciones transversales 2D definidas en el plano imagen, se | El último bloque del diagrama describe como el contorno de las secciones transversales 2D definidas en el plano imagen, se | ||
transforma en una nube de puntos (coordenadas 2D en el mundo) utilizando los parámetros intrínsecos de la cámara y la POSE del plano | transforma en una nube de puntos (coordenadas 2D en el mundo) utilizando los parámetros intrínsecos de la cámara y la POSE del plano | ||
− | de luz. Cada uno de los bloques mostrados en la Figura 5 se encuentra descritos en [3]. | + | de luz. Cada uno de los bloques mostrados en la [[#img-5|Figura 5]] se encuentra descritos en [3]. |
===3.2 Determinación de la profundidad=== | ===3.2 Determinación de la profundidad=== | ||
Para calcular la profundidad se usó un sistema de visión, un odómetro y un sensor magnético. Primeramente determinamos el modelo del | Para calcular la profundidad se usó un sistema de visión, un odómetro y un sensor magnético. Primeramente determinamos el modelo del | ||
− | sensor, después hicimos una corrección, la caracterización del sensor (Figura 6) y finalmente la fusión de los datos obtenidos | + | sensor, después hicimos una corrección, la caracterización del sensor ([[#img-6|Figura 6]]) y finalmente la fusión de los datos obtenidos |
por los sensores. | por los sensores. | ||
− | + | <div id='img-6'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;" | ||
|- | |- | ||
Line 405: | Line 406: | ||
El odómetro es un sistema que está compuesto por un conjunto de tarjetas electrónicas STM 32F4 discovery que se encargan de adquirir | El odómetro es un sistema que está compuesto por un conjunto de tarjetas electrónicas STM 32F4 discovery que se encargan de adquirir | ||
− | y procesar la señal de los pulsos generados (Figura 7). Entre las partes que componen el sistema se encuentra un módulo de | + | y procesar la señal de los pulsos generados ([[#img-7|Figura 7]]). Entre las partes que componen el sistema se encuentra un módulo de |
medición lineal. El módulo lineal usa un encoder opto-electrónico con un carrete acoplado al mismo, se jala un cable por un sistema | medición lineal. El módulo lineal usa un encoder opto-electrónico con un carrete acoplado al mismo, se jala un cable por un sistema | ||
dentado con una cadena en la que cada vuelta el encoder hace un conteo de los pulsos generados. El conteo de los pulsos generados | dentado con una cadena en la que cada vuelta el encoder hace un conteo de los pulsos generados. El conteo de los pulsos generados | ||
Line 412: | Line 413: | ||
Coordenadas (CMM) y la metodología presentada en [52] para una sola dimensión. La caracterización del sensor se hizo de forma | Coordenadas (CMM) y la metodología presentada en [52] para una sola dimensión. La caracterización del sensor se hizo de forma | ||
estadística aproximando la medición de los sensores a una distribución paramétrica normal. | estadística aproximando la medición de los sensores a una distribución paramétrica normal. | ||
− | + | <div id='img-7'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | ||
|- | |- | ||
Line 427: | Line 428: | ||
− | El sensor magnético ( | + | El sensor magnético ([[#img-8|Figura 8]]) está formado por un imán de posición sobre un eje de movimiento del dispositivo. El sistema de |
medición propone el uso de una escala magnética que funciona como un potenciómetro con una resistencia de 10K cuyo desplazamiento | medición propone el uso de una escala magnética que funciona como un potenciómetro con una resistencia de 10K cuyo desplazamiento | ||
depende de las posiciones del imán, para calcular el desplazamiento lineal al hacer presión sobre la capa superior del sensor, un | depende de las posiciones del imán, para calcular el desplazamiento lineal al hacer presión sobre la capa superior del sensor, un | ||
Line 434: | Line 435: | ||
Medición de Coordenadas (CMM) como en [52]. La caracterización del sensor de desplazamiento magnético se hizo de forma estadística | Medición de Coordenadas (CMM) como en [52]. La caracterización del sensor de desplazamiento magnético se hizo de forma estadística | ||
considerando una distribución normal. | considerando una distribución normal. | ||
− | + | <div id='img-8'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;" | ||
|- | |- | ||
Line 520: | Line 521: | ||
La alineación de los sensores es importante ya que la posición inicial del sistema se estima con esta referencia. Los datos | La alineación de los sensores es importante ya que la posición inicial del sistema se estima con esta referencia. Los datos | ||
− | obtenidos de los sensores son fusionados para mejorar la estimación de la posición del generador de luz activa (Figura 9). El | + | obtenidos de los sensores son fusionados para mejorar la estimación de la posición del generador de luz activa ([[#img-9|Figura 9]]). El |
odómetro y magneto adquieren datos los cuales son corregidos usando una curva polinómica previamente calculada. Las imágenes de la | odómetro y magneto adquieren datos los cuales son corregidos usando una curva polinómica previamente calculada. Las imágenes de la | ||
cámara son corregidas en cuanto la distorsión radial usando los coeficientes de distorsión. | cámara son corregidas en cuanto la distorsión radial usando los coeficientes de distorsión. | ||
+ | <div id='img-9'></div> | ||
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;" | ||
|- | |- | ||
Line 595: | Line 597: | ||
===3.3.1 Filtro de Kalman === | ===3.3.1 Filtro de Kalman === | ||
+ | |||
El Filtro de Kalman(FK) se utiliza para modelar sistemas lineales en el que las mediciones del sistema se ven afectadas por ruido | El Filtro de Kalman(FK) se utiliza para modelar sistemas lineales en el que las mediciones del sistema se ven afectadas por ruido | ||
aleatorio provenientes del comportamiento del sistema. Para hacer una estimación de los valores de la posición, el filtro requiere | aleatorio provenientes del comportamiento del sistema. Para hacer una estimación de los valores de la posición, el filtro requiere | ||
Line 685: | Line 688: | ||
El fibroscopio genera estructuras en forma de panal de abeja que se manifiestan en las imágenes capturadas. Estas estructuras deben | El fibroscopio genera estructuras en forma de panal de abeja que se manifiestan en las imágenes capturadas. Estas estructuras deben | ||
ser filtradas de lo contrario el procedimiento de calibración no puede realizarse a menos que se desarrolle un extractor de | ser filtradas de lo contrario el procedimiento de calibración no puede realizarse a menos que se desarrolle un extractor de | ||
− | características inmune a este tipo de ruido. La Figura 10 muestra la imagen de un patrón de calibración capturado a través del | + | características inmune a este tipo de ruido. La [[#img-10|Figura 10]] muestra la imagen de un patrón de calibración capturado a través del |
fibroscopio. En este trabajo se utilizaron dos filtros, uno propuesto por [56] y el otro propuesto por [57]. Después de algunas | fibroscopio. En este trabajo se utilizaron dos filtros, uno propuesto por [56] y el otro propuesto por [57]. Después de algunas | ||
pruebas, se estableció que la aplicación de ambos filtros en secuencia, primero Winter y luego Xu, da excelentes resultados. La | pruebas, se estableció que la aplicación de ambos filtros en secuencia, primero Winter y luego Xu, da excelentes resultados. La | ||
− | secuencia opuesta no funciona pues los resultados son similares con solo aplicar el filtro de Winter, Figura 10 (b). La Figura 10 | + | secuencia opuesta no funciona pues los resultados son similares con solo aplicar el filtro de Winter, [[#img-10|Figura 10]] (b). La [[#img-10|Figura 10]] |
muestra los resultados individuales de ambos filtros y los resultados cuando se aplican en secuencia. | muestra los resultados individuales de ambos filtros y los resultados cuando se aplican en secuencia. | ||
+ | <div id='img-10'></div> | ||
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;" | ||
|- | |- | ||
Line 703: | Line 707: | ||
utilizar como patrón de referencia. | utilizar como patrón de referencia. | ||
− | La Figura 11 representa un escaneo 3D de interior del cilindro de frenado de autos. La Figura 11 (a) muestra el escaneo de varias | + | La [[#img-11|Figura 11]] representa un escaneo 3D de interior del cilindro de frenado de autos. La [[#img-11|Figura 11]](a) muestra el escaneo de varias |
− | secciones transversales en el interior del cilindro maestro. La Figura 11 (b) muestra el error del diámetro al comparar el diámetro | + | secciones transversales en el interior del cilindro maestro. La [[#img-11|Figura 11]](b) muestra el error del diámetro al comparar el diámetro |
del cilindro ajustado con todas las secciones transversales medidas por el sistema propuesto y el diámetro medido por una CMM. El | del cilindro ajustado con todas las secciones transversales medidas por el sistema propuesto y el diámetro medido por una CMM. El | ||
− | error es de 0.028 mm con una incertidumbre de 0.050 mm la cual cubre dicho error. En la Figura 11 (b) también se aprecia la | + | error es de 0.028 mm con una incertidumbre de 0.050 mm la cual cubre dicho error. En la [[#img-11|Figura 11]](b) también se aprecia la |
cilindricidad de 0.186 mm. El error de diámetro es una prueba de la exactitud del sistema propuesto en el dimensionamiento 2D. Los | cilindricidad de 0.186 mm. El error de diámetro es una prueba de la exactitud del sistema propuesto en el dimensionamiento 2D. Los | ||
− | resultados de la Figura 11 (b) se obtuvieron con el software de medición certificado GOM Inspect V8. | + | resultados de la [[#img-11|Figura 11]](b) se obtuvieron con el software de medición certificado GOM Inspect V8. |
− | + | <div id='img-11'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | ||
|- | |- | ||
|style="padding:10px;"| [[File:Draft_Gonzalez-Barbosa_236306777_2789_cilindro1.jpg|600px]] | |style="padding:10px;"| [[File:Draft_Gonzalez-Barbosa_236306777_2789_cilindro1.jpg|600px]] | ||
|- | |- | ||
− | | style="text-align: center;font-size: 75%;"|(a) Escaneo 3D del interior de un cilindro de frenado de autos como el mostrado en la Figura | + | | style="text-align: center;font-size: 75%;"|(a) Escaneo 3D del interior de un cilindro de frenado de autos como el mostrado en la [[#img-1|Figura 1]](b) |
|- | |- | ||
|style="padding:10px;"| [[File:Draft_Gonzalez-Barbosa_236306777_9073_cilindro2.jpg|600px]] | |style="padding:10px;"| [[File:Draft_Gonzalez-Barbosa_236306777_9073_cilindro2.jpg|600px]] | ||
Line 744: | Line 748: | ||
probabilidad para las lecturas de los sensores. | probabilidad para las lecturas de los sensores. | ||
− | La Figura 12 muestra los resultados de las posiciones de la luz activa provenientes de datos de diferentes sensores, 15 posiciones | + | La [[#img-12|Figura 12]] muestra los resultados de las posiciones de la luz activa provenientes de datos de diferentes sensores, 15 posiciones |
diferentes se miden. La descripción del modelo del sistema está relacionado con estas mediciones, con los resultados obtenidos se | diferentes se miden. La descripción del modelo del sistema está relacionado con estas mediciones, con los resultados obtenidos se | ||
observa el comportamiento de un sistema dinámico lineal de tres variables independientes con la misma tendencia. A partir de conocer | observa el comportamiento de un sistema dinámico lineal de tres variables independientes con la misma tendencia. A partir de conocer | ||
Line 750: | Line 754: | ||
odométrico, se demuestra que a partir de este intervalo se tienen un valor de 0.1 mm de error acumulado con referencia al valor | odométrico, se demuestra que a partir de este intervalo se tienen un valor de 0.1 mm de error acumulado con referencia al valor | ||
medido para cada posición. | medido para cada posición. | ||
− | + | <div id='img-12'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | ||
|- | |- | ||
Line 760: | Line 764: | ||
Cuando se hacen varias mediciones se puede obtener un valor estadístico del error, el error asociado a las mediciones de la posición | Cuando se hacen varias mediciones se puede obtener un valor estadístico del error, el error asociado a las mediciones de la posición | ||
− | se presenta en la Figura 13, con estos valores sabemos que la medición de la posición es proporcional a la desviación que presenta, | + | se presenta en la [[#img-13|Figura 13]], con estos valores sabemos que la medición de la posición es proporcional a la desviación que presenta, |
y nos da la exactitud de éstas medidas, es decir, a mayor repetibilidad mayor exactitud. Al tomar en cuenta estas medidas nos dirá | y nos da la exactitud de éstas medidas, es decir, a mayor repetibilidad mayor exactitud. Al tomar en cuenta estas medidas nos dirá | ||
− | que tan preciso es el sensor con el que estamos midiendo. En la Figura 13 podemos apreciar que las incertidumbre en los datos | + | que tan preciso es el sensor con el que estamos midiendo. En la [[#img-13|Figura 13]] podemos apreciar que las incertidumbre en los datos |
presentados están contenidas en el rango de 0 hasta 3 mm . El intervalo en el que el sensor odométrico mostró un mayor error | presentados están contenidas en el rango de 0 hasta 3 mm . El intervalo en el que el sensor odométrico mostró un mayor error | ||
sistemático es de 50 mm a 100 mm que presenta incertidumbres cercanas a 3 mm, mientras que el sensor magnético en el intervalo de 20 | sistemático es de 50 mm a 100 mm que presenta incertidumbres cercanas a 3 mm, mientras que el sensor magnético en el intervalo de 20 | ||
Line 769: | Line 773: | ||
mayor capacidad de adquirir los datos, esta información nos es útil analizarla ya que posteriormente al hacer la fusión de los datos | mayor capacidad de adquirir los datos, esta información nos es útil analizarla ya que posteriormente al hacer la fusión de los datos | ||
esperamos que los sensores con menor variabilidad incluyan en la estimación final de la posición. | esperamos que los sensores con menor variabilidad incluyan en la estimación final de la posición. | ||
− | + | <div id='img-13'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | ||
|- | |- | ||
Line 786: | Line 790: | ||
variación se ajusta a una distribución normal que representamos con PDF's centrado en la media que es el mejor valor de la posición | variación se ajusta a una distribución normal que representamos con PDF's centrado en la media que es el mejor valor de la posición | ||
medido en un intervalo de medición de 10 mm y con una dispersión en proporción del error de esta medición, conocer la distribución | medido en un intervalo de medición de 10 mm y con una dispersión en proporción del error de esta medición, conocer la distribución | ||
− | del conjunto de datos para cada sensor nos proporciona información de la posición. En la Figura 14 se muestran tres PDF's, la | + | del conjunto de datos para cada sensor nos proporciona información de la posición. En la [[#img-14|Figura 14]] se muestran tres PDF's, la |
PDF1(x) representa al sensor óptico en rojo, en verde la PDF2(x) corresponde al sensor odométrico y en azul PDF3(x) al sensor magnético. Nuestro trabajo está enfocado en fusionar los resultados de estas distribuciones para dar una mejor estimación de la posición. | PDF1(x) representa al sensor óptico en rojo, en verde la PDF2(x) corresponde al sensor odométrico y en azul PDF3(x) al sensor magnético. Nuestro trabajo está enfocado en fusionar los resultados de estas distribuciones para dar una mejor estimación de la posición. | ||
− | + | <div id='img-14'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | ||
|- | |- | ||
Line 798: | Line 802: | ||
===4.3 Interpolación de PDF's=== | ===4.3 Interpolación de PDF's=== | ||
− | En la Figura 14 se muestra el conjunto de datos medidos en posiciones discretas definidas con distribuciones de probabilidad (PDF's), sin embargo, en la práctica cuando se realizan las mediciones hay datos en los que el sistema se encuentra en alguna posición en la que no se guarda ese valor de medición, ya que la presencia de información faltante lleva consigo perdida de precisión en las estimaciones se debe tomar en cuenta alguna metodología para estimar estos datos faltantes. Para encontrar una solución a esto, proponemos el uso del método de la interpolación lineal con ella llegamos a la obtención del modelado de la distribución de la posición cuando toma un valor fuera de un intervalo discreto donde no hay mediciones disponibles. | + | En la [[#img-14|Figura 14]] se muestra el conjunto de datos medidos en posiciones discretas definidas con distribuciones de probabilidad (PDF's), sin embargo, en la práctica cuando se realizan las mediciones hay datos en los que el sistema se encuentra en alguna posición en la que no se guarda ese valor de medición, ya que la presencia de información faltante lleva consigo perdida de precisión en las estimaciones se debe tomar en cuenta alguna metodología para estimar estos datos faltantes. Para encontrar una solución a esto, proponemos el uso del método de la interpolación lineal con ella llegamos a la obtención del modelado de la distribución de la posición cuando toma un valor fuera de un intervalo discreto donde no hay mediciones disponibles. |
La interpolación consiste en hallar un dato dentro de un intervalo en el que conocemos los valores extremos. Los puntos <math>X_1(N,\mu ,\sigma )</math> y | La interpolación consiste en hallar un dato dentro de un intervalo en el que conocemos los valores extremos. Los puntos <math>X_1(N,\mu ,\sigma )</math> y | ||
Line 809: | Line 813: | ||
sobre el resultado de la medición cada distribución aporta información complementaria que puede incrementar la exactitud de la | sobre el resultado de la medición cada distribución aporta información complementaria que puede incrementar la exactitud de la | ||
estimación de la posición del LED. Los resultados de las distribuciones que se fusionan mediante el producto de las PDF individuales | estimación de la posición del LED. Los resultados de las distribuciones que se fusionan mediante el producto de las PDF individuales | ||
− | de cada uno de los valores proporcionados por los diferentes sensores se muestran en la Figura 15. El producto final muestra una | + | de cada uno de los valores proporcionados por los diferentes sensores se muestran en la [[#img-15|Figura 15]]. El producto final muestra una |
incertidumbre menor que cuando el error es proveniente de cada uno de los sensores de manera individual, ya que es posible que el | incertidumbre menor que cuando el error es proveniente de cada uno de los sensores de manera individual, ya que es posible que el | ||
sensor que presente una menor incertidumbre tendrá mayor impacto sobre el resultado de la medición, esta distribución total está | sensor que presente una menor incertidumbre tendrá mayor impacto sobre el resultado de la medición, esta distribución total está | ||
asociada a un nuevo valor de posición. | asociada a un nuevo valor de posición. | ||
− | + | <div id='img-15'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | ||
|- | |- | ||
Line 825: | Line 829: | ||
Una vez obtenida la distribución total se calcula la media y la desviación estándar estimados asociados a esta distribución, el | Una vez obtenida la distribución total se calcula la media y la desviación estándar estimados asociados a esta distribución, el | ||
estrechamiento de la distribución indica que la incertidumbre implícita se reduce, y la media se centra con mayor precisión con | estrechamiento de la distribución indica que la incertidumbre implícita se reduce, y la media se centra con mayor precisión con | ||
− | respecto a la posición objetivo real. En la Tabla 2 se recopilan los resultados. Con este método modelamos el error de la posición | + | respecto a la posición objetivo real. En la [[#tab-2|Tabla 2]] se recopilan los resultados. Con este método modelamos el error de la posición |
para la localización del LED de potencia, hemos demostrado que este método es capaz de encontrar una estimación de la posición del | para la localización del LED de potencia, hemos demostrado que este método es capaz de encontrar una estimación de la posición del | ||
LED y podemos concluir que este método es capaz de realizar una estimación de la posición aun cuando los valores de los sensores | LED y podemos concluir que este método es capaz de realizar una estimación de la posición aun cuando los valores de los sensores | ||
Line 846: | Line 850: | ||
tiempo. | tiempo. | ||
− | La Figura 16 representa los valores de medición de los datos iniciales que son procesados por ambos filtros. Los filtros procesan la | + | La [[#img-16|Figura 16]] representa los valores de medición de los datos iniciales que son procesados por ambos filtros. Los filtros procesan la |
medición de 15 muestras dando como resultado la señal estimada X. | medición de 15 muestras dando como resultado la señal estimada X. | ||
− | La Figura 17 muestra la estimación de las posiciones después de que los datos son procesados por los filtros, junto con los valores | + | La [[#img-17|Figura 17]] muestra la estimación de las posiciones después de que los datos son procesados por los filtros, junto con los valores |
del movimiento real del objeto (rojo) con sus respectivos valores de posición estimados con ambos filtros, el filtro de kalman | del movimiento real del objeto (rojo) con sus respectivos valores de posición estimados con ambos filtros, el filtro de kalman | ||
(azul) y el filtro de partículas(negro). Los resultados muestran la estimación de la posición en cada punto de medición, se observa | (azul) y el filtro de partículas(negro). Los resultados muestran la estimación de la posición en cada punto de medición, se observa | ||
que los valores de posición son estimados son muy similares a la real. | que los valores de posición son estimados son muy similares a la real. | ||
− | + | <div id='img-16'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | ||
|- | |- | ||
Line 860: | Line 864: | ||
|} | |} | ||
− | + | <div id='img-17'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | ||
|- | |- | ||
Line 869: | Line 873: | ||
− | En la Figura 18, se muestra la incertidumbre de la posición con los métodos implementados, la línea roja representa la incertidumbre | + | En la [[#img-18|Figura 18]], se muestra la incertidumbre de la posición con los métodos implementados, la línea roja representa la incertidumbre |
generada al fusionar las incertidumbres procedentes de los sensores utilizados observamos que existe una diferencia considerable con | generada al fusionar las incertidumbres procedentes de los sensores utilizados observamos que existe una diferencia considerable con | ||
respecto a la incertidumbre obtenida de manera individual de cada sensor, los resultados con este método oscilan entre 1 mm a 3 mm, | respecto a la incertidumbre obtenida de manera individual de cada sensor, los resultados con este método oscilan entre 1 mm a 3 mm, | ||
Line 876: | Line 880: | ||
respectivamente, los valores obtenidos nos indican qué tan fiable es la medición que estamos obteniendo con nuestro sistema de | respectivamente, los valores obtenidos nos indican qué tan fiable es la medición que estamos obteniendo con nuestro sistema de | ||
posicionamiento, al tomar una decisión del método a utilizar para mejorar la posición del sistema. | posicionamiento, al tomar una decisión del método a utilizar para mejorar la posición del sistema. | ||
− | + | <div id='img-18'></div> | |
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | ||
|- | |- | ||
Line 885: | Line 889: | ||
− | Los resultados de los algoritmos propuestos se muestran en la Tabla 2, donde se observan las posiciones y el error de medición, | + | Los resultados de los algoritmos propuestos se muestran en la [[#tab-2|Tabla 2]], donde se observan las posiciones y el error de medición, |
tanto para el algoritmo de mezcla de gaussianas con EM, FK y FP, los resultados muestran que los métodos implementados tienden a | tanto para el algoritmo de mezcla de gaussianas con EM, FK y FP, los resultados muestran que los métodos implementados tienden a | ||
aproximarse a la posición real del sistema y ambos filtros pueden ser utilizados para la estimación de la posición, sin embargo, el | aproximarse a la posición real del sistema y ambos filtros pueden ser utilizados para la estimación de la posición, sin embargo, el | ||
filtro de partículas es el algoritmo que realiza una mejor estimación de la posición en comparación con los otros métodos, además, | filtro de partículas es el algoritmo que realiza una mejor estimación de la posición en comparación con los otros métodos, además, | ||
− | también se puede ver en los resultados de la Tabla 2 como la desviación estándar es menor en el algoritmo filtro de partículas lo | + | también se puede ver en los resultados de la [[#tab-2|Tabla 2]] como la desviación estándar es menor en el algoritmo filtro de partículas lo |
que indica que los resultados son más estables manteniéndose el LED localizado a lo largo de todo el recorrido del escaneo. Se puede | que indica que los resultados son más estables manteniéndose el LED localizado a lo largo de todo el recorrido del escaneo. Se puede | ||
observar como todos los valores de posiciones y el error comprendido entre los valores obtenidos se reducen drásticamente con la | observar como todos los valores de posiciones y el error comprendido entre los valores obtenidos se reducen drásticamente con la | ||
Line 895: | Line 899: | ||
otros métodos. | otros métodos. | ||
− | Para llevar a cabo la evaluación del proceso de medición de nuestro sensor, se calcularon los valores de incertidumbre de la posición en cada punto de medición, siguiendo los cálculos establecidos en [58], mostrados en la Tabla 2 y en la Figura 19. Finalmente, la | + | Para llevar a cabo la evaluación del proceso de medición de nuestro sensor, se calcularon los valores de incertidumbre de la posición en cada punto de medición, siguiendo los cálculos establecidos en [58], mostrados en la [[#tab-2|Tabla 2]] y en la [[#img-19|Figura 19]]. Finalmente, la [[#tab-3|Tabla 3]] muestra el error promedio de los valores de incertidumbre de las mediciones obtenidas en la [[#tab-2|Tabla 2]], analizando los resultados obtenidos se puede ver que estos valores representan el 1% del error en la medición calculado a partir de la dispersión de los datos con respecto al valor promedio de las medidas obtenidas por el sensor. |
− | + | ||
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;font-size: 75%;"> | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;font-size: 75%;"> | ||
− | '''Tabla 2'''. Datos obtenidos de todas las estimaciones (mm)</div> | + | '''Tabla 2'''. Datos obtenidos de todas las estimaciones (mm)</div><br /> |
− | + | <center> | |
− | {| class="wikitable" align="center" style="margin: | + | <div id='tab-2'></div> |
− | |- | + | {| class="wikitable" align="center" style="margin:0;text-align: center;font-size:85%;" |
+ | |-style="border-top:2px solid; border-bottom:2px solid;" | ||
| '''Distancia de referencia''' || '''Mezcla de gaussianas''' || '''Incertidumbre'''|| '''FK''' || '''Incertidumbre''' || '''FP''' || '''Incertidumbre''' | | '''Distancia de referencia''' || '''Mezcla de gaussianas''' || '''Incertidumbre'''|| '''FK''' || '''Incertidumbre''' || '''FP''' || '''Incertidumbre''' | ||
|- | |- | ||
Line 932: | Line 936: | ||
|- | |- | ||
|130 || 131.388 || 1.1 || 130.786 || 0.6 || 130.326 || 0.2 | |130 || 131.388 || 1.1 || 130.786 || 0.6 || 130.326 || 0.2 | ||
− | |- | + | |-style="border-bottom:2px solid;" |
|140 || 142.619 || 2.1 || 141.854 || 1.5 || 141.056 || 0.8 | |140 || 142.619 || 2.1 || 141.854 || 1.5 || 141.056 || 0.8 | ||
|} | |} | ||
+ | </center> | ||
− | + | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;font-size: 75%;"> | |
− | <div class=" | + | '''Tabla 3'''. Incertidumbres obtenidas de estimaciónes de posición discretas (mm) |
− | '''Tabla 3''' | + | |
</div> | </div> | ||
− | {| class="wikitable" style="align: center;" | + | <div id='tab-3'></div> |
− | |- | + | <center> |
+ | {| class="wikitable" align="center" style="margin:0;text-align: center;font-size:85%;" | ||
+ | |-style="border-top:2px solid; border-bottom:2px solid;" | ||
! incertidumbre !! Mezcla de gausssianas !! FK !! FP | ! incertidumbre !! Mezcla de gausssianas !! FK !! FP | ||
− | |- | + | |-style="border-bottom:2px solid;" |
| || 1.8 || 1.64 || 0.73 | | || 1.8 || 1.64 || 0.73 | ||
|} | |} | ||
+ | </center> | ||
− | |||
+ | Los resultados experimentales muestran que el error de medición del sensor propuesto está dentro del 1%. Este resultado se compara a continuación con el obtenido en otros trabajos, en los que se presentan sistemas que proporcionan información útil para aplicaciones de navegación e inspección de tuberías y modelado 3D. Autores como [59] presentan un sistema de visión capaz de obtener información de la profundidad con un error de estimación menor al 2%. Un sistema de visión para el escaneo de tuberías 3D como el propuesto en [60], en el método de escaneo consiste en proyectar un haz de luz circular sobre la superficie del tubo, que permite calcular la posición de escaneo con una resolución de menos del 1%. En [61] demuestran cómo se puede realizar la medición del perfil interno de piezas cilíndricas con una resolución del 2%. Mientras que en [62] proponen un sistema de medición para la reconstrucción 3D de tuberías, el método consiste en estimar el movimiento de la cámara con referencia las coordenadas 3D de la sección de luz, este sistema puede estimar el movimiento de la cámara con un error del 3%. En [63] proponen un método de localización simultánea para estimar con precisión la posición de un instrumento utilizado con aplicaciones quirúrgicas, con una precisión menor del 2%. En [64] proponen un sistema para la inspección de tuberías que consiste en un robot móvil que se desplaza en el interior del ducto, la posición del móvil es medida a través de un distanciómetro láser con una precisión del 3%. En [65] presentan una técnica de estimación de la posición de la cámara dentro de la tubería con una precisión del 1%. De los resultados obtenidos con el sistema de medición propuesto en nuestro trabajo, se deduce que la precisión de nuestro sensor puede considerarse del mismo nivel que los trabajados mencionados, lo que demuestra que este sensor es adecuado para utilizarse en procesos donde es posible adaptar mediciones ópticas donde se espera que el sistema sea capaz de evaluar el perfil de la superficie interna de la pieza a medir. | ||
+ | <div id='img-19'></div> | ||
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | ||
|- | |- | ||
Line 968: | Line 976: | ||
comportamientos similares, se encontró que existe un error sistemático de 1.64 mm y 0.73 mm respectivamente logrando de esta forma | comportamientos similares, se encontró que existe un error sistemático de 1.64 mm y 0.73 mm respectivamente logrando de esta forma | ||
obtener resultados que permiten mejorar la medición de la posición del LED. Los resultados del posicionamiento obtenidos con el | obtener resultados que permiten mejorar la medición de la posición del LED. Los resultados del posicionamiento obtenidos con el | ||
− | sistema de medición propuesto se comparan con los resultados de trabajos existentes y se presentan en la Figura 19, la solución que | + | sistema de medición propuesto se comparan con los resultados de trabajos existentes y se presentan en la [[#img-19|Figura 19]], la solución que |
proporcionan estas técnicas son cercanos a los presentados por nuestro sistema, lo que indica que el sistema de escaneo se puede | proporcionan estas técnicas son cercanos a los presentados por nuestro sistema, lo que indica que el sistema de escaneo se puede | ||
utilizar para diferentes aplicaciones | utilizar para diferentes aplicaciones | ||
Line 977: | Line 985: | ||
<span id='_Ref449083719'></span> | <span id='_Ref449083719'></span> | ||
<div id="1"></div> | <div id="1"></div> | ||
− | [[#cite-1|[1]]] Z. | + | [[#cite-1|[1]]] Zhang Z. Flexible camera calibration by viewing a plane from unknown orientations. Proceedings of the Seventh IEEE International Conference on Computer Vision, pp. 666-673, 1999. |
<div id="2"></div> | <div id="2"></div> | ||
− | [[#cite-2|[2]]] S. | + | [[#cite-2|[2]]] Rupp S., Winter C., Wittenberg T. Camera calibration from fiberscopic views with accuracy evaluation. Bildverarbeitung für die Medizin, pp. 424-428, 2006. |
<div id="3"></div> | <div id="3"></div> | ||
− | [[#cite-3|[3]]] | + | [[#cite-3|[3]]] Icasio-Hernández O., Gonzalez-Barbosa J.-J., Hurtado-Ramos J., Viliesid-Alonso M. 3d reconstruction of hollow parts analyzing images acquired by a fiberscope. Measurement Science and Technology 25 (7), 2014. |
<div id="4"></div> | <div id="4"></div> | ||
− | [[#cite-4|[4]]] L. | + | [[#cite-4|[4]]] Chen L., Bavigadda V., Kofidis T., Howe R. Fiber optic projection-imaging system for shape measurement in confined space. The ScientificWorld Journal, 28:206569, 2014. |
<div id="5"></div> | <div id="5"></div> | ||
− | [[#cite-5|[5]]] T. | + | [[#cite-5|[5]]] Stehle T., Truhn D., Aach T., Trautwein C., Tischendorf J. Camera calibration for fish-eye lenses in endoscopy with an application to 3D reconstruction. In Proceedings IEEE International Symposium on Biomedical Imaging (ISBI), IEEE, Washington D.C., pp. 1176-1179, 2007. |
<div id="6"></div> | <div id="6"></div> | ||
− | [[#cite-6|[6]]] J. | + | [[#cite-6|[6]]] Barreto J., Santos J. Ray-based calibration of rigid medical endoscopes. The 8th Workshop Omnidirectional Vision Camera Networks and Non Classical Cameras (OMNIVIS), Marseille, France, 2008. |
<div id="7"></div> | <div id="7"></div> | ||
− | [[#cite-7|[7]]] J. P. | + | [[#cite-7|[7]]] Barreto J.P., Roquette J., Sturm P., Fonseca F. Automatic camera calibration applied to medical endoscopy. Procedings of the British Machine Vision Conference, 52.1-52.10, 2009. |
<div id="8"></div> | <div id="8"></div> | ||
− | [[#cite-8|[8]]] R. | + | [[#cite-8|[8]]] Melo R., Barreto J., Falcão G. A new solution for camera calibration and real-time image distortion correction in medical endoscopy-initial technical evaluation. IEEE Transactions on Bio-Medical Engineering, 59(3):634-44, 2012. |
<div id="9"></div> | <div id="9"></div> | ||
− | [[#cite-9|[9]]] C. | + | [[#cite-9|[9]]] Schmalz C., Forster F., Schick A., Angelopoulou E. An endoscopic 3d scanner based on structured light. Medical Image Analysis, 16(5):1063-1072, 2012. |
<div id="10"></div> | <div id="10"></div> | ||
− | [[#cite-10|[10]]] M. | + | [[#cite-10|[10]]] Hayashibe M., Suzuki N., Nakamura Y. Laser-scan endoscope system for intraoperative geometry acquisition and surgical robot safety management. Medical Image Analysis, 10(4):509-19, 2006. |
<div id="11"></div> | <div id="11"></div> | ||
− | [[#cite-11|[11]]] C. | + | [[#cite-11|[11]]] Wu C.H., Sun Y.N., Chang C.C. Three-dimensional modeling from endoscopic video using geometric constraints via feature positioning. IEEE Transactions on Biomedical Engineering, 54(7):1199-1211, 2007. |
<div id="12"></div> | <div id="12"></div> | ||
− | [[#cite-12|[12]]] ] C. Wu, S. | + | [[#cite-12|[12]]] ] C. Wu, Narasimhan S., Jaramaz B. Shape-from-shading under near point lighting and partial views for orthopedic endoscopy. In Proceedings of the First Int. Workshop on Photometric Analysis For Computer Vision - PACV 2007, 8 pp., Rio de Janeiro, Brazil, 2007. |
<div id="13"></div> | <div id="13"></div> | ||
− | [[#cite-13|[13]]] G. | + | [[#cite-13|[13]]] Wang G.H., Han J.Q., Zhang X.M. Three-dimensional reconstruction of endoscope images by a fast shape from shading method. Measurement Science and Technology, 20 (12):125801, 2009. |
<div id="14"></div> | <div id="14"></div> | ||
− | [[#cite-14|[14]]] L. | + | [[#cite-14|[14]]] Criante L., Lucesoli A., Farabollini B., Bonifazi F., Rozzi T., Simoni F. Size measurement in endoscopic images by low coherence interferometry. Journal of Optics A: Pure and Applied Optics, 11(3):034007, 2009. |
<div id="15"></div> | <div id="15"></div> | ||
− | [[#cite-15|[15]]] J. | + | [[#cite-15|[15]]] Penne J., Höller K., Stürmer M., Schrauder T., Schneider A., Engelbrecht R., Feussner H., Schmauss B., Hornegger J. Time-of-flight 3-D endoscopy. International Conference on Medical Image Computing and Computer-Assisted Intervention (MICCAI 2009), 12 (Pt 1), 467-474, 2009. |
− | Computer-Assisted Intervention 12 (Pt 1) | + | |
<div id="16"></div> | <div id="16"></div> | ||
− | [[#cite-16|[16]]] C. | + | [[#cite-16|[16]]] Daul C., Blondel W. From 2D towards 3D cartography of hollow organs. 7th IEEE International Conference on ElectricalEngineering, Computing Science and Automatic Control, CCE 2010, CCE 201, pp. 285-293, Tuxtla Gutiérrez, Mexico, Sept. 2010. |
<div id="17"></div> | <div id="17"></div> | ||
− | [[#cite-17|[17]]] T. | + | [[#cite-17|[17]]] Seki T., Oka K., Naganawa A., Yamashita H., Kim K., ChibaT.T. Laser distance measurement using a newly developed composite-type optical fiberscope for fetoscopic laser surgery. Optics and Lasers in Engineering, 48(10):974-977, 2010. |
<div id="18"></div> | <div id="18"></div> | ||
− | [[#cite-18|[18]]] K. | + | [[#cite-18|[18]]] Arai K. Method for 3D object reconstruction using several portions of 2D images from the diferent aspects acquired with image scopes included in the fiber retractor. International Journal of Advanced Research in Artificial Intelligence, 1(9):14-19, 2012. |
<div id="19"></div> | <div id="19"></div> | ||
− | [[#cite-19|[19]]]C. | + | [[#cite-19|[19]]] Schmalz C., Forster F., Schick A., Angelopoulou E. An endoscopic 3d scanner based on structured light. Medical Image Analysis, 16(5):1063-1072, 2012. |
<div id="20"></div> | <div id="20"></div> | ||
− | [[#cite-20|[20]]] M. | + | [[#cite-20|[20]]] Hartrumpf M., Munser R. Optical three-dimensional measurements by radially symmetric structured light projection. Applied Optics, 36(13):2923-2928, 1997. |
<div id="21"></div> | <div id="21"></div> | ||
− | [[#cite-21|[21]]] D. | + | [[#cite-21|[21]]] Schaack D. Making accurate, three-dimensional measurements through a standard borescope. Proc. SPIE 3397, Nondestructive Evaluation of Aging Aircraft, Airports, and Aerospace Hardware II, 3397:264-276, 1998. |
<div id="22"></div> | <div id="22"></div> | ||
− | [[#cite-22|[22]]] K. | + | [[#cite-22|[22]]] Hasegawa K., Noda K., Sato Y. Electronic endoscope system for shape measurement. IEEE: Object Recognition Supported by User interaction for Service Robots, 1:792,795, Quebec, Canada, 2002. |
<div id="23"></div> | <div id="23"></div> | ||
− | [[#cite-23|[23]]] C. | + | [[#cite-23|[23]]] Winter C., Scholz I. Reconstruction of tubes from monocular fiberscopic images. Application and first results. In Vision, Modeling, and Visualization, pp. 57-64, 2005. |
<div id="24"></div> | <div id="24"></div> | ||
− | [[#cite-24|[24]]] G. | + | [[#cite-24|[24]]] Biegelbauer G., Vincze M. 3D vision-guided bore inspection system. In Fourth IEEE International Conference on Computer Vision Systems (ICVS'06), pp. 22- 22, 2006. |
<div id="25"></div> | <div id="25"></div> | ||
− | [[#cite-25|[25]]] H. D. | + | [[#cite-25|[25]]] Ford H.D., Tatam R.P. Passive OCT probe head for 3D duct inspection. Measurement Science and Technology, 24(9):094001, 2013. |
<div id="26"></div> | <div id="26"></div> | ||
− | [[#cite-26|[26]]] A. | + | [[#cite-26|[26]]] Müller A., Schubert M., Beleites E. Noncontact three-dimensional laser measuring device for tracheoscopy. Annals of Otology, Rhinology and Laryngology, 111(9):821-827, 2002. |
<div id="27"></div> | <div id="27"></div> | ||
− | [[#cite-27|[27]]] W. | + | [[#cite-27|[27]]] Dörfiel W., Sugano Y., Stalling D., Coconu L., Hentschel D., Linÿ G., Witt C. Laser-based endoscopic measurement of airway dimensions. Pneumologie, 57(9):503-509, 2003. |
<div id="28"></div> | <div id="28"></div> | ||
− | [[#cite-28|[28]]] J. | + | [[#cite-28|[28]]] Williamson J., James A., Phillips M., Sampson D., Hillman D., Eastwood P. Quantifying tracheobronchial tree dimensions: Methods, limitations and emerging techniques. European Respiratory Journal, 34(1):42-55, (2009). |
<div id="29"></div> | <div id="29"></div> | ||
− | [[#cite-29|[29]]] J. | + | [[#cite-29|[29]]] Geng J., Xie J. Review of 3-d endoscopic surface imaging techniques. IEEE Sensors Journal, 14(4):945-960, 2014. |
<div id="30"></div> | <div id="30"></div> | ||
− | [[#cite-30|[30]]] D. | + | [[#cite-30|[30]]] Yang D., Liu Z., Sun F., Zhang J., Liu H., Wang S. Recursive depth parametrization of monocular visual navigation: Observability analysis and performance evaluation. Information Sciences, 287:38-49, 2014. |
<div id="31"></div> | <div id="31"></div> | ||
− | [[#cite-31|[31]]] H. | + | [[#cite-31|[31]]] Firouzi H., Najjaran H. Real-time monocular vision-based object tracking with object distance and motion estimation. IEEE/ASME International Conference on Advanced Intelligent Mechatronics, (AIM 2010), pp. 987-992, 2010. |
<div id="32"></div> | <div id="32"></div> | ||
− | [[#cite-32|[32]]] Z. | + | [[#cite-32|[32]]] Zhang Z., Han Y., Zhou Y., Dai M. A novel absolute localization estimation of a target with monocular vision. Optik, 124(12):1218-1223, 2013. |
<div id="33"></div> | <div id="33"></div> | ||
− | [[#cite-33|[33]]]A. | + | [[#cite-33|[33]]] Assa A., Janabi-Sharifi F. Virtual visual servoing for multicamera pose estimation. IEEE/ASME Transactions on Mechatronics, 20(2):789-798, 2014. |
<div id="34"></div> | <div id="34"></div> | ||
− | [[#cite-34|[34]]] D. | + | [[#cite-34|[34]]] Burschka D., Mair E. Direct pose estimation with a monocular camera. Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics), 4931 LNCS, pp. 440-453, 2008. |
<div id="35"></div> | <div id="35"></div> | ||
− | [[#cite-35|[35]]] B. | + | [[#cite-35|[35]]] Williams B., Klein G., Reid I. Real-time slam relocalisation. IEEE 11th International Conference on Computer Vision, pp. 1-8, Rio de Janeiro, 2007. |
<div id="36"></div> | <div id="36"></div> | ||
− | [[#cite-36|[36]]] P. | + | [[#cite-36|[36]]] Owczarek P., Goslinski J., Rybarczyk D., Kubacki A. Control of an electro-hydraulic manipulator by vision system using central point of a marker estimated via kalman filter. ICA 2016. Advances in Intelligent Systems and Computing, Springer International Publishing, pp. 587-596, 2016. |
<div id="37"></div> | <div id="37"></div> | ||
− | [[#cite-37|[37]]] A. W. | + | [[#cite-37|[37]]] Stroupe A.W., Martin M.C., Balch T. Merging Gaussian distributions for object localization in multi-robot systems. In Experimental Robotics VII, Lecture Notes in Control and Information Sciences, Vol. 271, pp. 343-352, Springer, Berlin, Heidelberg, 2001. |
<div id="38"></div> | <div id="38"></div> | ||
− | [[#cite-38|[38]]] D. | + | [[#cite-38|[38]]] Fox D., Burgard W., Thrun S. Markov localization for mobile robots in dynamic environments. Journal of Artificial Intelligence Research, 11:391-427, 1999. |
<div id="39"></div> | <div id="39"></div> | ||
− | [[#cite-39|[39]]] | + | [[#cite-39|[39]]] O'Kane J. Global localization using odometry. Proceedings 2006 IEEE International Conference on Robotics and Automation (ICRA 2006), pp. 37-42, 2006. |
<div id="40"></div> | <div id="40"></div> | ||
− | [[#cite-40|[40]]] L. | + | [[#cite-40|[40]]] Delahoche L., Pégard C., Marhic B. Vasseur P. Navigation system based on an omnidirectional vision sensor. In Proc. of the Int. Conf. on Intelligent Robotics and Systems 1997 (IROS’97), Vol. 2, pp. 718–724, Grenoble, France, 1997. |
<div id="41"></div> | <div id="41"></div> | ||
− | [[#cite-41|[41]]] P. | + | [[#cite-41|[41]]] Eberhard P., Tang Q. Sensor data fusion for the localization and position control of one kind of omnidirectional mobile robots. Multibody System Dynamics, Robotics and Control, Springer, Chapter 4, pp. 45-73, 2013. |
<div id="6"></div> | <div id="6"></div> | ||
− | [[#cite-42|[42]]] M. | + | [[#cite-42|[42]]] Caruso M., Bratland T., Smith C., Schneider R. Anisotropic magnetoresistive sensors: Theory and applications. Sensors (Peterborough, NH), 16(3):18-26, 1999. |
<div id="43"></div> | <div id="43"></div> | ||
− | [[#cite-43|[43]]] C. | + | [[#cite-43|[43]]] Schott C., Racz R., Betschart F., Popovic R. A new two-axis magnetic position sensor. SENSORS, 2002 IEEE, Vol. 2, pp. 911-915, 2002. |
<div id="44"></div> | <div id="44"></div> | ||
− | [[#cite-44|[44]]] | + | [[#cite-44|[44]]] Bosch-Jorge M., Soriano A., Valera A., Sánchez Salmerón A.J. Detección y seguimiento de objetos móviles utilizando un escáner láser mediante filtros de kalman. Comité Español de Automática (CEA-IFAC), Valencia (Spain), 2014. |
<div id="45"></div> | <div id="45"></div> | ||
− | [[#cite-45|[45]]] J. | + | [[#cite-45|[45]]] Woo J., Kim Y., Lee J., Lim M. Localization of mobile robot using particle filter. In 2006 SICE-ICASE International Joint Conference, pp. 3031-3034, 2006. |
<div id="46"></div> | <div id="46"></div> | ||
− | [[#cite-46|[46]]] P. | + | [[#cite-46|[46]]] Vadakkepat P., Jing L. Improved particle filter in sensor fusion for tracking randomly moving object. IEEE Transactions on Instrumentation and Measurement, 55(5):1823-1832, 2006. |
<div id="47"></div> | <div id="47"></div> | ||
− | [[#cite-47|[47]]] C. | + | [[#cite-47|[47]]] Yuan C., Lianchuan M., Ma W. Mobile target tracking based on hybrid open-loop monocular vision motion control strategy. Discrete Dynamics in Nature and Society, 10 pp., 2015. |
<div id="48"></div> | <div id="48"></div> | ||
− | [[#cite-48|[48]]] Z. | + | [[#cite-48|[48]]] Boukhers Z., Shirahama K., Li F., Grzegorzek M. Object detection and depth estimation for 3d trajectory extraction. In Proceedings of International Workshop on Content-Based Multimedia Indexing, Vol. 2015, July 2015. |
<div id="49"></div> | <div id="49"></div> | ||
− | [[#cite-49|[49]]] | + | [[#cite-49|[49]]] Sakharkar S.S., Kamble S.D., Khobragade A.S. Object detection and tracking using particle filters. International Journal of Computer Trends and Technology(IJCTT), Vol. 22(1), April 2015. |
<div id="50"></div> | <div id="50"></div> | ||
− | [[#cite-50|[50]]] J.-Y. | + | [[#cite-50|[50]]] Bouguet J.-Y. Camera calibration toolbox for Matlab, 2015. |
<div id="51"></div> | <div id="51"></div> | ||
− | [[#cite-51|[51]]] J. Y. | + | [[#cite-51|[51]]] Bouguet J.-Y. Camera calibration toolbox for Matlab, 2008. |
<div id="52"></div> | <div id="52"></div> |
En este trabajo presentamos un sistema de posicionamiento de visión activa para el escaneo 3D del interior de piezas. El diseño del sistema propuesto consta de dos módulos: un sistema de dimensionamiento 2D de visión activa, y un sistema que posiciona el módulo de visión activa. El sistema de posicionamiento es capaz de determinar la profundidad del sistema de dimensionamiento 2D de visión activa en el interior del objeto a escanear usando varios sensores. Las principales contribuciones de este trabajo son la caracterización del sistema de dimensionamiento 2D, y el desarrollo de algoritmos de posicionamiento de la luz activa con énfasis en el modelado y fusión de sensores. El sistema puede utilizarse como un sistema de dimensionamiento en aplicaciones industriales como la industria metal mecánica, la aeronáutica, la medicina, en el control de calidad y en áreas de visión por computadora.
Palabras clave: Estimación de la posición, procesamiento de imágenes, fusión de datos, modelado sensorial
In this work we present an active positioning system for 3D scan of interior parts. The design of the proposed system consists of two modules: an active 2D dimensional system and positional system based on active vision. The active 2D dimensional system is able to determine the depth of the 2D dimensional system inside the object to be scanned using several sensors. The main contributions of this work are the characterization of the 2D dimensional system and the development of active light positioning algorithms with emphasis on the modeling and fusion of the sensors. The system can be used as a dimensional system in industrial applications such as the metal mechanical industry, aeronautics industry, medicine, quality control and computer vision.
Keywords: Pose estimation, image processing, data fusion, sensor models
En este trabajo presentamos una propuesta de un sistema de posicionamiento de luz activa para el escaneo 3D del interior de piezas industriales con gran utilidad en diferentes tipos de industria como la metal mecánica, la aeronáutica, la militar y la medicina que requieren constantemente la medición de dimensiones internas de la pieza bajo inspección. Como ejemplos de aplicación podemos mencionar la medición del interior del cilindro maestro que forma parte del sistema de frenado de los autos, el escaneo 3D del interior del cuerpo humano, la medición en el interior de tuberías, la medición del interior de botellas hechas de plástico, entre otros (Figura 1). El presente trabajo propone el uso de un sistema de medición capaz de obtener la posición de la luz activa para poder dimensionar el interior del objeto, este sistema consta de dos módulos: un sistema de dimensionamiento 2D de visión activa, y un sistema de posicionamiento que calcula la profundidad del sistema de dimensionamiento 2D de visión activa. La profundidad es calculada por la información proveniente de varios sensores de forma independiente los cuales son fusionados para obtener una medición de mayor precisión y exactitud.
El sistema de dimensionamiento del interior de piezas con difícil acceso está basado en un fibroscopio modificado. La modificación
consiste en la sustitución del sistema de iluminación original del fibroscopio por un sistema que proyecta una línea de luz. Esta
línea de luz revela una sección transversal del interior del objeto bajo estudio. La línea de luz se obtiene mediante un LED de
potencia encerrado entre 2 placas y unido por una extensión a la punta del fibroscopio (Figuras 2(a) y (b)). Las imágenes de
diferentes secciones transversales del objeto son capturadas por una cámara conectada al ocular del fibroscopio por medio de un
video acoplador (Figura 2(c)).
El sistema que determina la profundidad de la visión activa, está conformado por un sistema de visión, un sensor magnético y un
odómetro. El sistema de visión consta de una cámara y de un patrón de referencia pegado en las placas del sistema de iluminación,
unidos al sistema de guiado del fibroscopio que componen al sistema de visión activa. El patrón está constituido por una matriz de 8
columnas de círculos y 8 renglones de círculos de color blanco con fondo negro. La cámara es del tipo endoscopio digital flexible de
7 mm de diámetro y una distancia de trabajo de 150 mm y captura imágenes con una resolución de 640x480 pixeles. El sensor magnético
funciona como un potenciómetro cuyo valor de resistencia está relacionado con el valor de la posición, este puede adherirse a la
pieza bajo medición, la lectura del sensor magnético se realiza mediante un módulo Bluetooth. El odómetro se coloca en forma
paralela a la cámara y se desplaza en dirección al eje óptico de la cámara. Nuestra configuración del sistema está diseñado para una
distancia óptima de 150 mm de profundidad. La principal aportación de este trabajo está centrada en el módulo de posicionamiento del
sistema de dimensionamiento. La contribución de este módulo es la adaptación de un sistema de visión junto con la integración de
sensores de desplazamiento lineal de bajo costo y fácil adaptación; un sensor magnético y un odómetro que es desarrollado por
nosotros, para calcular el posicionamiento del sistema. En particular, el enfoque consiste en fusionar las mediciones provenientes
de los sensores con el fin de lograr estimaciones del posicionamiento del sistema más exactas, el cálculo de la posición es formulada utilizando técnicas de estimación como los filtros FK y FP que modelan el comportamiento lineal del sistema, ambos métodos
han sido implementados con el propósito de reducir el error de la estimación resultante. La novedad de este trabajo es el desarrollo
de un método de posicionamiento para el sistema de visión activa para la mejora del dimensionamiento de objetos.
El sistema de escaneo tridimensional presentado está compuesto por un módulo de dimensionamiento de visión activa y un módulo de posicionamiento. El módulo de visión activa consta de un fibroscopio y un generador de luz. El módulo de posicionamiento está formado por un sistema de visión, un sensor magnético y un odómetro. La aplicación del sistema de medición se basa en la determinación de la profundidad del sistema de dimensionamiento 2D de visión activa, para poder dimensionar el interior del objeto con este sistema, cada componente ó todo el sistema de visión debe ser calibrado. Son pocos los artículos que especifican métodos particulares para calibrar una cámara conectada a un fibroscopio utilizada como parte del sistema de escaneo y que proponen un sistema de posicionamiento de luz activa con varios sensores para el dimensionamiento de objetos, en esta sección presentamos algunos trabajos comparativos con el presente trabajo.
Los artículos que describen la calibración de la cámara conectada al ocular de un endoscopio los podemos dividir en dos grupos: los que se basan en la técnica propuesta por [1] como [2,3,4], y los que realizan la calibración considerando una lente de pescado (lente fish-eye) con un ángulo de visión superior a los 45 grados como los propuestos por [5,6,7,8]. El trabajo presentado por [9] consta de un sistema catadióptrico formado por un endoscopio y un espejo esférico. Los autores calibran el sistema de dos formas: usando el modelo pinhole presentado en [1] donde la distorsión radial modela la curvatura del espejo y proponen un modelo propio para el sistema.
Varias de las publicaciones para el escaneo 3D de piezas huecas utilizando endoscopios rígidos o flexibles se encuentran en el campo de la medicina, por ejemplo [10-19] utilizan piezas de referencia para validar los escaneos 3D en el interior del cuerpo humano, por lo que bien pudieran ser utilizadas en el sector industrial. Los trabajos que utilizan endoscopios rígidos o flexibles para escaneos 3D de utilidad en la industria son: [20,21,22,23,24,25]. En [20] se desarrolló un robot para inspeccionar tuberías con diámetros entre 20 cm y 50 cm, su sistema emite luz estructurada para escanear tridimensionalmente el interior de las tuberías. En [21] se menciona que la mayoría de los trabajos dedicados a escaneos 3D a través de endoscopios han sido enfocados a patentes, el método de calibración de la cámara descrito es similar a nuestro trabajo, incluso la manera de comprobar la exactitud del sistema, ya que sigue la ideología de nuestro trabajo usando piezas de referencia. En [22] se utiliza un videoscopio y un escáner láser para realizar los escaneos 3D, la manera de comprobar la exactitud es con piezas de referencia y su trabajo está enfocado al área médica.
De acuerdo con [23], existen pocas publicaciones que describen el escaneo 3D de piezas huecas a partir de imágenes adquiridas con fibroscopios monoculares. Los trabajos [26,27], al igual que nosotros, utilizan un fibroscopio y escanean en 2D secciones transversales que luego posicionan mediante otro sistema. Sin embargo, [26,27] se basan en un enfoque de iluminación activa utilizando un espejo cónico para proyectar un rayo láser en forma de anillo sobre la superficie interior de la tráquea con el fin de encontrar estenosis. En el trabajo de [10] al igual que nosotros, proyectan un haz de luz para resaltar la superficie a medir, para ello utilizan un galvanómetro con el fin de resaltar diferentes zonas en el interior del órgano. Este trabajo encuentra también la relación angular entre la cámara y el láser con el fin de calcular coordenadas 3D. Para lograr el escaneo 3D, [16], proyecta un láser de color verde con el fin de triangular haciendo uso del láser y la cámara. [19] desarrolló su propio escáner endoscópico, su trabajo al igual que el nuestro proyecta uno o más anillos en el interior a escanear, con dichos anillos y una lente catadióptrica conectada a la cámara puede capturar puntos en 3D. En [23] el escaneo 3D es a partir de la extracción de puntos característicos de las imágenes y al igual que en este trabajo filtran las estructuras de panal de abeja para mejorar la calibración de la cámara y los escaneos 3D. [24] desarrollaron un sistema para inspeccionar barrenos que utiliza un endoscopio rígido y un proyector láser que proyecta un patrón circular. El sistema es guiado por un brazo robótico e inspecciona diámetros de barrenos entre 4 mm y 50 mm con una profundidad de 100 mm. [25] presentan una técnica conocida como Tomografía de Coherencia Óptica (OCT) para el escaneo 3D, también utilizaron un espejo cónico para la iluminación. Esta técnica no incluye componentes electro-mecánicos en la punta de la sonda. Dos revisiones que tratan el escaneo 3D utilizando endoscopios y análisis de imágenes tanto en el área médica como la industrial son presentados en [28,29]. Ninguna de estas revisiones presenta el diseño de un sistema de posicionamiento similar a nuestro trabajo, por lo que se puede considerar a nuestro sistema como innovador tomando en cuenta el sistema de iluminación propuesto, las metodologías para calibrar el sistema, la integración de los sensores y la forma de calcular el posicionamiento de la luz activa y evaluar los errores de medición y su incertidumbre, donde la principal aportación es el desarrollo del sistema de posicionamiento.
La Tabla 1 muestra las principales características que presenta nuestro sistema para escanear el interior de piezas comparado con otros sistemas. Los artículos presentados en esta tabla utilizan luz estructurada con uno o varios planos, o proyectan luz con formas conocidas. Los trabajos que no usan luz estructurada, [23], detectan características en la imagen que son seguidas mientras el sistema de visión calibrado se desplaza de manera controlada, usando el método de factorización las características seguidas son escaneadas tridimensionalmente. La mayoría de los trabajos usan una referencia con el propósito de transformar las unidades de pixeles a unidades métricas. Pocos son los trabajos que hacen un análisis del error y de las incertidumbres del sistema de medición.
Ref. | A) | B) | C) | D) | E) | F) | G) |
[26] | x | x | x | ||||
[27] | x | x | |||||
[10] | x | x | x | ||||
[16] | x | x | x | ||||
[19] | x | x | |||||
[20] | x | x | x | ||||
[21] | x | x | x | x | x | ||
[22] | x | x | |||||
[23] | x | x | |||||
[24] | x | x | x | x | |||
[25] | x | x | x | x | |||
Nuestro trabajo | x | x | x | x | x | x | x |
A) Uso de luz estructurada, B) Uso endoscopio o fibroscopio, C) Uso de referencias, D) Aplicación Industrial, E) Cálculo POSE de luz, F) Cálculo del error e incertidumbre, G) Cálculo de profundidad.
Existen en la literatura trabajos relacionados con la estimación de la posición basados en la visión y la integración de sensores. Algunos de los trabajos que utilizan sistemas de visión para calcular la posición de los objetos utilizando información proveniente de las imágenes capturadas por la cámara son [30-32]. Autores como [33] realizan una estimación de la posición utilizando un sistema de múltiples cámaras, la configuración de estos sistemas combinan las observaciones de todas las cámaras para calcular los parámetros de rotación y traslación del sistema a partir de una secuencia de imágenes. En [34] presentan un método que permite calcular los valores de la posición y rotación de la navegación de una cámara monocular, calculan directamente los parámetros de movimiento a partir de una serie de imágenes y estiman la precisión de los parámetros. El trabajo propuesto por [35] aborda el problema de localización utilizando la aplicación de un algoritmo de mapeo probabilístico y localización simultáneos (SLAM) utilizando una sola cámara. El algoritmo SLAM permite generar mapas probabilísticos del entorno dados los puntos de las proyecciones de la imagen; con los mapas obtenidos es posible generar trayectorias de navegación.
Otros trabajos utilizan múltiples sensores como láser, giroscopios, GPS, etc. Para obtener información del posicionamiento de un objeto y mejorar la estimación de la posición de este, que se basan en métodos probabilísticos para fusionar los datos de los sensores como el filtro de Kalman [36], mezcla de gaussianas [37], localización de markov para entornos dinámicos como [38], etc. En el trabajo de [39] se describe la implementación de algoritmos de navegación para determinar la posición de un robot móvil a partir de la estimación del odómetro cuando este se mueve en un entorno conocido, donde muestran la integración de odómetros en aplicaciones de robótica móvil, el robot puede rotar y trasladarse a través del entorno y a su vez medir estos movimientos para su localización. Otro sistema de localización de un robot móvil en un ambiente estructurado es el desarrollado por [40] su sistema combina un sistema de visión omnidirecional junto con un odómetro, calculan una estimación fiable de la posición debido a que realiza un ajuste en las mediciones al combinar los datos de los sensores utilizando un método de filtrado estadístico (EKF). En [41], plantean el problema de la localización y control de posición de un robot móvil mediante el uso de sensores y la fusión de datos obtenidos de estos sensores, primero son corregidos los datos obtenidos por los odómetros del robot y después retroalimentan el EKF para tener mayor control de la posición. Otro enfoque para resolver el problema de localización de un objeto es el trabajo presentado por [42] donde se utiliza un sistema de detección de posición lineal, utilizando una serie de sensores magneto resistivos (MR) y señales de acondicionamiento electrónico. Los sensores se utilizan para determinar la posición de un imán que está unido a un objeto en movimiento, estos sensores miden la relación de la posición del imán y la función del campo magnético para determinar la posición del objeto [43], en este trabajo se describe un nuevo sensor de desplazamiento magnético que mide la posición de un imán a lo largo de dos ejes. El cálculo de la posición se consigue midiendo la dirección del campo magnético. Los sensores utilizados en este trabajo son dos circuitos de montaje superficial, a diferencia del sensor que utilizamos en nuestro trabajo que es un sensor magnético de película delgada. En [44] se expone la implementación de un algoritmo de detección y seguimiento de objetos mediante un escáner láser, su sistema funciona a partir de la segmentación de los datos del sensor que permite determinar el número de objetos presentes en el entorno y determinar sus movimientos para poder llevar a cabo una navegación de forma segura. Otros autores como [45] describen la auto localización de un robot utilizando las mediciones de un sensor ultrasónico basado en el filtro de partículas, los autores introducen varios métodos de localización para eliminar el ruido del proceso y de la medición, muestran resultados comparativos entre el filtro EKF y FP para la localización del robot con ambos métodos.
En [46] proponen un método de localización de un objeto en movimiento a través de una secuencia de imágenes tomadas por una cámara, utiliza el filtro de partículas como algoritmo de seguimiento basado en la información de múltiples sensores, integran un sistema de visión junto con un sistema de control para seguir de forma continua el objeto. Al igual que autores como [47] que proponen un método de seguimiento objetos basado en el control de movimiento de una cámara utilizan el filtro de partículas para predecir la posición del objeto en movimiento correspondiente a la información obtenida de una imagen. Autores como [48] analizan el problema entre la detección de un objeto y la estimación de la profundidad y realizan un seguimiento del objeto usando filtro de partículas, donde la posición de este objeto es modelado para la estimación de la profundidad. En [49] los autores presentan el seguimiento de un objeto en una secuencia de video, en sustitución de un modelo dinámico del objeto utilizan el método de substracción de fondo para predecir las ubicaciones del objeto y la aplicación del filtro de partículas para el seguimiento del objeto. En este trabajo presentamos el desarrollo de algoritmos para dar una solución al problema de la estimación de la posición, para mejorar el dimensionamiento 2D del interior del objeto a escanear.
En esta sección se presenta un sistema capaz de escanear el interior de objetos en dos etapas. La primer etapa realiza el dimensionamiento 2D de una sección transversal usando visión activa. La segunda etapa localiza el sistema de visión activa dentro del interior del objeto para asignarle al dimensionamiento 2D el valor de profundidad que corresponde a la posición de la luz activa.
El sistema cuenta con una cámara CCD monocromática (C1) con una resolución de 24482046 pixeles y un sensor de 1/1.45" que se conecta al ocular del fibroscopio utilizando un video acoplador óptico (VA) (Figura 2). Un fibroscopio de 7 mm de diámetro (FIBRO) con una longitud efectiva de 1320 mm, dispone de una distancia mínima de trabajo de 30.48 mm, y un campo de visión de aproximadamente 50 grados alrededor del eje vertical y 40 grados alrededor del eje horizontal. La configuración del sistema del sistema de medición propuesto se muestra en la Figura 3.
La Figura 4 muestra la cámara endoscópica (mC2), que utiliza tecnología optoelectrónica para inspeccionar áreas de difícil acceso.
La cámara endoscópica (mC2) tiene un diámetro de 7 mm, las dimensiones son de 510.54 mm 7.62 mm 7.62 mm, la resolución es de
640x480 pixeles (JPEG), adquiere a 30 cuadros por segundo (FPS) y usa una luz LED ajustable en la punta de la cámara para iluminar
la zona de inspección. En la punta del fibroscopio se adjuntó un dispositivo de iluminación, omitiendo el sistema de iluminación
proporcionado por el fibroscopio. Este dispositivo está formado por un LED de potencia de color rojo de 1 watt (LP) y dos placas
rectangulares opacas unidas a la circunferencia del LED. Las placas rectangulares se colocan de tal manera que el componente emisor
de luz del diodo está entre las placas (véase región del LP ampliada que muestra el LED entre placas sujetas al fibroscopio). De
esta manera, la luz sólo puede salir a través de la ranura formada por las placas, formando un plano de iluminación que toma la
forma de la sección transversal del interior donde se proyecte. Las placas tienen un patrón de círculos impresos, que junto con la
cámara 2 (mC2) localizan el plano de luz al interior del objeto.
![]() |
Figura 4. Cámara Endoscopica USB Super Eyes N005 |
Otro sistema de posicionamiento usado es un sistema llamado magneto-pot o escala magnética (ME), que es un potenciómetro con una
resistencia variable, su operación se fundamenta en la variación de la resistencia en función de la posición. La ME puede adherirse
a la pieza bajo medición. La lectura de posición de la ME y el control de intensidad del LP se realizan mediante un módulo Bluetooth
(BT). Una interfaz gráfica de usuario (PC) es utilizada para controlar y leer todos los dispositivos del sistema y procesar la
información necesaria para realizar el escaneo 3D. Un tercer método para calcular el posicionamiento de las secciones transversales
es utilizando un odómetro (ODO), el cual calcula el desplazamiento del LP uniendo la cuerda del odómetro a una de las placas del LP.
El diagrama de flujo de la Figura 5 muestra la metodología para realizar el dimensionamiento o escaneo 2D. En él se muestran tres grandes bloques de operaciones delimitados con línea punteada. Procedimiento de calibración, trazabilidad, posición y orientación (POSE).
![]() |
Figura 5. Diagrama de flujo para el escaneo 2D del interior de piezas huecas |
El procedimiento de calibración que se muestra en el bloque superior calcula los parámetros intrínsecos de la cámara conectada al
fibroscopio y los coeficientes de distorsión de la lente acoplada a la cámara. Los parámetros y los coeficientes son evaluados con
el toolbox de Bouguet et al. [50] modificado para detectar círculos. El bloque de trazabilidad permite encontrar la POSE del plano
de luz generado por el LED de potencia encerrado entre placas y la POSE de los anillos patrón, ambas POSEs respecto al sistema de
coordenadas de referencia ubicado en la cámara. Los anillos proporcionan la trazabilidad y corrección de la POSE del plano de luz.
El último bloque del diagrama describe como el contorno de las secciones transversales 2D definidas en el plano imagen, se
transforma en una nube de puntos (coordenadas 2D en el mundo) utilizando los parámetros intrínsecos de la cámara y la POSE del plano
de luz. Cada uno de los bloques mostrados en la Figura 5 se encuentra descritos en [3].
Para calcular la profundidad se usó un sistema de visión, un odómetro y un sensor magnético. Primeramente determinamos el modelo del sensor, después hicimos una corrección, la caracterización del sensor (Figura 6) y finalmente la fusión de los datos obtenidos por los sensores.
![]() |
Figura 6. Modelado, corrección y caracterización de los sensores para el cálculo de profundidad |
El sistema de visión utilizado consta de una cámara y un patrón de referencia pegado en una de las placas del LED de potencia. El
modelo de la cámara usado es el pinhole. Para calcular los parámetros de la cámara usamos un patrón de calibración con círculos
impresos. En la caracterización, el patrón se mueve respecto a la cámara en intervalos de 10 mm usando un sistema de
posicionamiento calibrado, por ejemplo una Máquina de Medición por Coordenadas (CMM), el sistema es capaz de medir hasta una
distancia de 150 mm. Se toman tres secuencias de medición para caracterizar el sistema de visión.
Con el propósito de dar trazabilidad a los parámetros del modelo, se utiliza un patrón de calibración compuesto por una matriz de 8 x 8 círculos. La posición de cada centro del círculo es calculado por una Máquina de Medición de Coordenadas Ópticas (OCMM). Con este procedimiento extra, eliminamos los errores que puede introducir la impresora que imprime el patrón de círculos. Se modificó el toolbox de Bouguet [51] para utilizar los círculos en el cálculo de los parámetros intrínsecos.
El odómetro es un sistema que está compuesto por un conjunto de tarjetas electrónicas STM 32F4 discovery que se encargan de adquirir y procesar la señal de los pulsos generados (Figura 7). Entre las partes que componen el sistema se encuentra un módulo de medición lineal. El módulo lineal usa un encoder opto-electrónico con un carrete acoplado al mismo, se jala un cable por un sistema dentado con una cadena en la que cada vuelta el encoder hace un conteo de los pulsos generados. El conteo de los pulsos generados por el encoder relaciona las vueltas y el desplazamiento lineal de la cuerda que está sujeta al generador de luz estructurada. El modelo utilizado para el odómetro es un polinomio. Para calcular la corrección del polinomio usamos una Máquina de Medición de Coordenadas (CMM) y la metodología presentada en [52] para una sola dimensión. La caracterización del sensor se hizo de forma estadística aproximando la medición de los sensores a una distribución paramétrica normal.
![]() | |
(a) Odómetro completo | |
![]() | |
(b) Tarjeta STM 32F4 discovery que permite la lectura del odómetro y la interfaz de los datos | |
Figura 7. Sensor odométrico |
El sensor magnético (Figura 8) está formado por un imán de posición sobre un eje de movimiento del dispositivo. El sistema de
medición propone el uso de una escala magnética que funciona como un potenciómetro con una resistencia de 10K cuyo desplazamiento
depende de las posiciones del imán, para calcular el desplazamiento lineal al hacer presión sobre la capa superior del sensor, un
segundo imán se coloca sobre las placas del generador de luz (LP), lo que permite conocer el desplazamiento del sistema de luz
activa. El modelo utilizado para el sensor magnético es un polinomio. Para la corrección del polinomio usamos una Máquina de
Medición de Coordenadas (CMM) como en [52]. La caracterización del sensor de desplazamiento magnético se hizo de forma estadística
considerando una distribución normal.
![]() |
Figura 8. Arriba: Módulo bluetooth. Abajo: Escala magnética |
Sea una variable que toma como posibles valores , , que representan donde se encuentra posicionada la luz activa en un instante de tiempo , donde a cada valor de la variable se le asocia un valor de probabilidad, entonces para el conjunto de valores de junto con sus probabilidades se define la distribución de probabilidad como:
|
(1) |
Si es la localización del sistema en el instante actual , dadas una secuencia de observaciones medidas por los tres sensores , entonces para calcular la posición del sistema en el instante actual a partir de cierta medida se define un valor de probabilidad condicionada a todos los datos de los sensores disponibles mediante :
|
(2) |
El problema general es estimar la posición del generador de luz activa con el sistema propuesto. Las mediciones del sistema son datos con incertidumbres que se pueden modelar mediante distribuciones de probabilidad (PDF), en este caso como una distribución gaussiana. La ecuación (3) describe la probabilidad de que el generador de luz activa se encuentre en la posición medida por el sensor [53]
|
(3) |
donde es el dato medido por el sensor (sistema de visión, odómetro y magneto), es la media y es la desviación estándar de la distribución para la medición del sensor. Las PDF's resultantes para cada se utilizan para representar cada una de las mediciones de los sensores con sus incertidumbres, cada medición produce una distribución gaussiana de la posición del generador de la luz activa. Para cada sensor la media y la desviación estándar han sido calculadas para posiciones discretas, el desplazamiento de referencia de la luz activa fue hecho con la Máquina de Medición de Coordenadas (CMM). Los datos son interpolados para posiciones donde no se conoce la media y desviación estándar de referencia. Existen diferentes métodos para la interpolación de los datos, en este trabajo utilizamos el algoritmo de maximización de la esperanza (expectation-maximization, EM). La interpolación lineal de los datos está definida por [54]:
|
(4) |
|
(5) |
donde , corresponde a la media y desviación estándar para las posiciones desplazadas por la Máquina de Medición de Coordenadas (CMM), y los datos estadísticos interpolados, es el parámetro de interpolación, indica la diferencia entre dos valores observados, dadas dos muestras aleatorias , la es la varianza de , variable aleatoria unidimensional .
La alineación de los sensores es importante ya que la posición inicial del sistema se estima con esta referencia. Los datos obtenidos de los sensores son fusionados para mejorar la estimación de la posición del generador de luz activa (Figura 9). El odómetro y magneto adquieren datos los cuales son corregidos usando una curva polinómica previamente calculada. Las imágenes de la cámara son corregidas en cuanto la distorsión radial usando los coeficientes de distorsión.
![]() |
Figura 9. Metodología utilizada para el cálculo de la profundidad |
Usando la CMM, el sistema de dimensionamiento activo se desplazó varias veces a diferentes posiciones y se caracterizaron los datos
de los sensores a diferentes posiciones de la CMM. Estas posiciones de la CMM son discretas, cuando una posición del sensor no ha
sido caracterizada se usa la interpolación descrita en la Sección 3.2.1.
El siguiente paso será el de fusionar todos los sensores. Cuando se tienen datos de diferentes sensores lo conveniente es realizar una fusión de datos con el fin de disminuir el error en la estimación de la posición, la idea de la fusión de dos o más datos es combinar las PDF's de estas mediciones para obtener un solo conjunto de datos. Las correspondientes distribuciones se fusionan mediante la multiplicación de las PDF's dada la ecuación 6, que representa la multiplicación de los datos y esta definida como la distribución de probabilidad conjunta
|
(6) |
De la resultante de la combinación de las PDF de cada sensor, extraemos los datos del valor esperado o media, , que representa la mejor estimación de la posición y la desviación estándar que es la incertidumbre asociada al valor esperado de la posición [55]. Entonces, la posición estimada de la luz activa es corregida mediante la fusión de estas distribuciones de probabilidad.
Los filtros de Kalman y Partículas son algoritmos que actualizan recursivamente una estimación del estado del sistema de visión activa dada una secuencia de observaciones. Utilizamos la aplicación de ambos filtros como estimadores del posicionamiento del sistema de dimensionamiento basado en luz activa para tener un mayor control de la medición de la posición. Ambos filtros comienzan con una representación del modelo dinámico del sistema, el filtro de Kalman estima la posición en base a retroalimentación del modelo lineal del sistema, mientras que el filtro de partículas lo hace en términos de la predicción de muestras aleatorias, la idea principal es representar la distribución de la posición del sistema. El modelo dinámico del sistema que describe el comportamiento del sistema es presentado por:
|
(7) |
con una medida:
|
(8) |
donde es el vector de estado que describe la posición del sistema, representa el error de la estimación del estado, es el error aleatorio asociado a las mediciones de los sensores y es la medición conocida dada por el sistema en la posición actual. y son constantes iguales a 1 debido a que nuestro sistema de posicionamiento se mueve sobre un solo eje. El objetivo es dar una solución al problema de la estimación de la posición de la luz activa aplicando los algoritmos propuestos el FK y FP para estimar el vector de estado dado un conjunto de mediciones provenientes del sistema de medición, así, los filtros procesan las mediciones sobre un conjunto de muestras, dando finalmente como resultado el valor estimado .
El Filtro de Kalman(FK) se utiliza para modelar sistemas lineales en el que las mediciones del sistema se ven afectadas por ruido aleatorio provenientes del comportamiento del sistema. Para hacer una estimación de los valores de la posición, el filtro requiere de los valores estadísticos del ruido del proceso y de las mediciones, las cuales son utilizados para la estimación de la posición de la luz activa para producir estimaciones que están mas cerca del valor real, a partir de las mediciones obtenidas de la mezcla de gaussianas se obtiene un vector de medición , estas mediciones son utilizadas para inicializar el filtro.
Las ecuaciones definidas en la sección 3.3 representan a como la ecuación del estado del sistema, como la ecuación de las mediciones, a como el error asociado a la medida y wn como el error asociado al sistema. es la covarianza del sistema relacionada con la diferencia entre el valor real del estado y la estimación, es la covarianza del error de medición que se estima a partir de la desviación estándar de la medida de los sensores. Experimentalmente el filtro toma en cuenta los siguientes parámetros: el estado inicial del sistema se definió como , la covarianza del sistema y la covarianza del error asociado a la medición fue . Posteriormente se actualizan la ecuación de observación usando los estados estimados y la covarianza mediante las ecuaciones:
|
(9) |
|
(10) |
El término minimiza la covarianza del error denominado como Ganancia de Kalman:
|
(11) |
La ecuación de actualización de las medidas por el filtro esta definida por:
|
(12) |
|
(13) |
Una vez descritas las ecuaciones matemáticas del filtro, se procedió a realizar la simulación correspondiente de las estimaciones. Las gráficas en la sección de resultados ilustran el comportamiento típico del filtro de Kalman, el objetivo de utilizar el filtro de Kalman es mostrar la estimación de la posición de la luz activa.
Presentamos la aplicación del filtro de partículas (FP) como un modelo de estimación similar al filtro de Kalman, en cada estimación usa como entrada el valor de posición en el estado anterior y la medición actual de los sensores. El algoritmo consiste en definir el posicionamiento de la luz activa a partir de una PDF inicial con un número de () partículas utilizando los datos provenientes de los sensores, hemos decidido inicializar el filtro como sigue y ; con un número total de 100 partículas, estas partículas representan muestras aleatorias de posibles valores cercanos a la posición del sistema, en un segundo instante de tiempo cuando el sistema se mueve a la siguiente posición las partículas se desplazan intentando predecir la posición del sistema en ese instante de tiempo, las partículas dispersadas alrededor de esa posición son evaluadas y se determina la probabilidad de cada partícula asociada a un peso , remuestrea anulando aquellas partículas con poco valor de probabilidad, después de estos pasos las partículas distribuidas con mayor peso tendrán mayor probabilidad de obtener una PDF en torno a una única solución para la estimación del posicionamiento de la luz activa, estimamos recursivamente esta distribución en cada momento en que una medida es tomada.
El sistema de visión activa propuesto es capaz de dimensionar el interior de objetos que van desde los 20 mm a 60 mm en dimensiones transversales. Cualquier perfil interior uniforme y no uniforme se puede escanear siempre y cuando la cámara conectada al fibroscopio tenga una visión clara de la línea de luz LED que alcanza las paredes. El cálculo de la profundidad del dimensionamiento produce una nube de puntos 3D de la pieza bajo medición.
El fibroscopio genera estructuras en forma de panal de abeja que se manifiestan en las imágenes capturadas. Estas estructuras deben ser filtradas de lo contrario el procedimiento de calibración no puede realizarse a menos que se desarrolle un extractor de características inmune a este tipo de ruido. La Figura 10 muestra la imagen de un patrón de calibración capturado a través del fibroscopio. En este trabajo se utilizaron dos filtros, uno propuesto por [56] y el otro propuesto por [57]. Después de algunas pruebas, se estableció que la aplicación de ambos filtros en secuencia, primero Winter y luego Xu, da excelentes resultados. La secuencia opuesta no funciona pues los resultados son similares con solo aplicar el filtro de Winter, Figura 10 (b). La Figura 10 muestra los resultados individuales de ambos filtros y los resultados cuando se aplican en secuencia.
Una comparación de las mediciones de una sección transversal del cilindro de frenado de autos se realizó entre los resultados del
sistema propuesto y los obtenidos por medio de una CMM. La CMM es un equipo calibrado trazable a patrones nacionales y se puede
utilizar como patrón de referencia.
La Figura 11 representa un escaneo 3D de interior del cilindro de frenado de autos. La Figura 11(a) muestra el escaneo de varias secciones transversales en el interior del cilindro maestro. La Figura 11(b) muestra el error del diámetro al comparar el diámetro del cilindro ajustado con todas las secciones transversales medidas por el sistema propuesto y el diámetro medido por una CMM. El error es de 0.028 mm con una incertidumbre de 0.050 mm la cual cubre dicho error. En la Figura 11(b) también se aprecia la cilindricidad de 0.186 mm. El error de diámetro es una prueba de la exactitud del sistema propuesto en el dimensionamiento 2D. Los resultados de la Figura 11(b) se obtuvieron con el software de medición certificado GOM Inspect V8.
![]() |
(a) Escaneo 3D del interior de un cilindro de frenado de autos como el mostrado en la Figura 1(b) |
![]() |
(b) Medición de las secciones transversales del cilindro |
Figura 11. Escanea 3D del Cilindro |
El LED de potencia se mueve a lo largo de una trayectoria dentro del objeto a escanear y se quiere medir su posición. Los tres sensores se mueven junto con el LED y toman mediciones, estas se agrupan en un vector de posición , donde representa el número de veces que son tomadas las mediciones, es el número de mediciones para una posición específica. Los datos se adquieren por medio de mediciones repetidas, para este experimento medimos tres secuencias de datos diferentes para cada sensor, con la repetición de las medidas observamos que se producen resultados iguales o similares debido a que estas están sujetas a errores asociados al sistema, una mejor manera de estimar estas mediciones es calculando la media de los datos para tener una aproximación al valor medido.
Una vez realizado el escaneo, se debe determinar el modelo del sistema de acuerdo con las lecturas recogidas por los sensores. Para realizar una estimación correcta de todas las posiciones con sus incertidumbres que ocurren durante la adquisición de datos, asumimos que nuestros datos representan un conjunto de variables aleatorias independientes, uniformemente distribuidas que siguen una distribución normal gaussiana de acuerdo a con una media centrada en el punto medio y una desviación estándar que modela la incertidumbre.
Por lo que debemos encontrar los valores de sus parámetros estadísticos de la media y la desviación estándar y a partir de estos podemos realizar una caracterización del modelo que pueda adecuarse a las posiciones que están más cerca de las posiciones reales. En secciones posteriores definimos formalmente el modelo del sensor definiendo cada una de las distribuciones de probabilidad para las lecturas de los sensores.
La Figura 12 muestra los resultados de las posiciones de la luz activa provenientes de datos de diferentes sensores, 15 posiciones diferentes se miden. La descripción del modelo del sistema está relacionado con estas mediciones, con los resultados obtenidos se observa el comportamiento de un sistema dinámico lineal de tres variables independientes con la misma tendencia. A partir de conocer el comportamiento del sistema vemos que en la zona de 80 a 150 mm es donde una variabilidad de medida está presente en el sensor odométrico, se demuestra que a partir de este intervalo se tienen un valor de 0.1 mm de error acumulado con referencia al valor medido para cada posición.
![]() |
Figura 12. Gráfica de la posición del LED de potencia cada 10 mm, para cada sensor |
Cuando se hacen varias mediciones se puede obtener un valor estadístico del error, el error asociado a las mediciones de la posición
se presenta en la Figura 13, con estos valores sabemos que la medición de la posición es proporcional a la desviación que presenta,
y nos da la exactitud de éstas medidas, es decir, a mayor repetibilidad mayor exactitud. Al tomar en cuenta estas medidas nos dirá
que tan preciso es el sensor con el que estamos midiendo. En la Figura 13 podemos apreciar que las incertidumbre en los datos
presentados están contenidas en el rango de 0 hasta 3 mm . El intervalo en el que el sensor odométrico mostró un mayor error
sistemático es de 50 mm a 100 mm que presenta incertidumbres cercanas a 3 mm, mientras que el sensor magnético en el intervalo de 20
a 50 mm se aproxima a los 0.1 mm, a su vez el sensor óptico después de los 110 mm presenta variaciones de 0.1 mm debido a que en
esas posiciones el campo visión de la cámara ya es limitado. Así pues podemos concluir que los sensores de menor error muestran
mayor capacidad de adquirir los datos, esta información nos es útil analizarla ya que posteriormente al hacer la fusión de los datos
esperamos que los sensores con menor variabilidad incluyan en la estimación final de la posición.
![]() |
Figura 13. Desviación estándar σ de los sensores a diferentes posiciones |
Al iniciar el análisis de una serie de datos un primer paso consiste en describir el comportamiento de estos. En nuestro caso estamos interesados en conocer un valor de posición después de la realización de una medición y la determinación de la distribución a la cual se refiere ese valor. Ya que nuestros resultados de medición son datos que miden una misma variable de manera repetida en intervalos de tiempo de 500 ms y que presentan cierta variación con respecto a la cantidad real de posición se asume que esta variación se ajusta a una distribución normal que representamos con PDF's centrado en la media que es el mejor valor de la posición medido en un intervalo de medición de 10 mm y con una dispersión en proporción del error de esta medición, conocer la distribución del conjunto de datos para cada sensor nos proporciona información de la posición. En la Figura 14 se muestran tres PDF's, la PDF1(x) representa al sensor óptico en rojo, en verde la PDF2(x) corresponde al sensor odométrico y en azul PDF3(x) al sensor magnético. Nuestro trabajo está enfocado en fusionar los resultados de estas distribuciones para dar una mejor estimación de la posición.
![]() |
Figura 14. La estimación de la posición se representa como una distribución de probabilidad |
En la Figura 14 se muestra el conjunto de datos medidos en posiciones discretas definidas con distribuciones de probabilidad (PDF's), sin embargo, en la práctica cuando se realizan las mediciones hay datos en los que el sistema se encuentra en alguna posición en la que no se guarda ese valor de medición, ya que la presencia de información faltante lleva consigo perdida de precisión en las estimaciones se debe tomar en cuenta alguna metodología para estimar estos datos faltantes. Para encontrar una solución a esto, proponemos el uso del método de la interpolación lineal con ella llegamos a la obtención del modelado de la distribución de la posición cuando toma un valor fuera de un intervalo discreto donde no hay mediciones disponibles.
La interpolación consiste en hallar un dato dentro de un intervalo en el que conocemos los valores extremos. Los puntos y , sean dos variables independientes asociadas a una distribución normal, entonces estimamos un valor de ubicado dentro del intervalo con sus estadísticos y (Sección 3.2.1) asociados a una distribución interpolada.
La combinación de las distribuciones de los sensores depende de la asociación de las incertidumbres de estos y del impacto que den sobre el resultado de la medición cada distribución aporta información complementaria que puede incrementar la exactitud de la estimación de la posición del LED. Los resultados de las distribuciones que se fusionan mediante el producto de las PDF individuales de cada uno de los valores proporcionados por los diferentes sensores se muestran en la Figura 15. El producto final muestra una incertidumbre menor que cuando el error es proveniente de cada uno de los sensores de manera individual, ya que es posible que el sensor que presente una menor incertidumbre tendrá mayor impacto sobre el resultado de la medición, esta distribución total está asociada a un nuevo valor de posición.
![]() |
Figura 15. Fusión de los sensores a diferentes posiciones |
Una vez obtenida la distribución total se calcula la media y la desviación estándar estimados asociados a esta distribución, el estrechamiento de la distribución indica que la incertidumbre implícita se reduce, y la media se centra con mayor precisión con respecto a la posición objetivo real. En la Tabla 2 se recopilan los resultados. Con este método modelamos el error de la posición para la localización del LED de potencia, hemos demostrado que este método es capaz de encontrar una estimación de la posición del LED y podemos concluir que este método es capaz de realizar una estimación de la posición aun cuando los valores de los sensores sean muy ruidosos y aun cuando las contribuciones individuales de los sensores no sean buenas.
Nuestro objetivo con la implementación de los filtros como estimadores de la posición es el de calcular de forma más precisa la posición del sistema. Primero fusionamos directamente las mediciones de los sensores en un solo vector de mediciones y después estos valores son utilizados por ambos filtros. En el caso del FK resuelve el problema de la estimación de la posición, tomando el vector de mediciones para retroalimentar con estos valores a las ecuaciones que definen el comportamiento del sistema descritas en la sección 3.3. La principal ventaja de los filtros de Kalman es su eficiencia computacional en la realización de las operaciones al utilizar datos con distribuciones unidimensionales. Aunque a diferencia de otros filtros, el filtro de Kalman no permite la posibilidad de añadir o quitar sensores ya que se basa en un modelo específico del comportamiento del sistema. La fiabilidad y respuesta del filtro está ligado a la formulación del modelo y una vez definido el sistema con las variables que retroalimentan al filtro ningún otro valor puede ser quitado o añadido debido a que a partir de estos valores se retroalimenta el filtro. En cambio, el filtro de partículas realiza una estimación de la posición en base al muestreo aleatorio de partículas, considerando una distribución inicial propuesta que se utiliza para realizar una estimación de las distribuciones combinando las observaciones obtenidas con el sistema a lo largo del tiempo.
La Figura 16 representa los valores de medición de los datos iniciales que son procesados por ambos filtros. Los filtros procesan la medición de 15 muestras dando como resultado la señal estimada X. La Figura 17 muestra la estimación de las posiciones después de que los datos son procesados por los filtros, junto con los valores del movimiento real del objeto (rojo) con sus respectivos valores de posición estimados con ambos filtros, el filtro de kalman (azul) y el filtro de partículas(negro). Los resultados muestran la estimación de la posición en cada punto de medición, se observa que los valores de posición son estimados son muy similares a la real.
![]() |
Figura 16. Datos iniciales |
![]() |
Figura 17. Estimación de las posiciones con FK y FP |
En la Figura 18, se muestra la incertidumbre de la posición con los métodos implementados, la línea roja representa la incertidumbre
generada al fusionar las incertidumbres procedentes de los sensores utilizados observamos que existe una diferencia considerable con
respecto a la incertidumbre obtenida de manera individual de cada sensor, los resultados con este método oscilan entre 1 mm a 3 mm,
adicionalmente, implementamos los algoritmos FK (azul) y FP (negra) con el propósito de disminuir la incertidumbre de la estimación
de la posición, observamos que con ambas implementaciones se produce una reducción de la incertidumbre de 1.64 mm y 0.73 mm
respectivamente, los valores obtenidos nos indican qué tan fiable es la medición que estamos obteniendo con nuestro sistema de
posicionamiento, al tomar una decisión del método a utilizar para mejorar la posición del sistema.
![]() |
Figura 18: Desviaciones de los filtros |
Los resultados de los algoritmos propuestos se muestran en la Tabla 2, donde se observan las posiciones y el error de medición,
tanto para el algoritmo de mezcla de gaussianas con EM, FK y FP, los resultados muestran que los métodos implementados tienden a
aproximarse a la posición real del sistema y ambos filtros pueden ser utilizados para la estimación de la posición, sin embargo, el
filtro de partículas es el algoritmo que realiza una mejor estimación de la posición en comparación con los otros métodos, además,
también se puede ver en los resultados de la Tabla 2 como la desviación estándar es menor en el algoritmo filtro de partículas lo
que indica que los resultados son más estables manteniéndose el LED localizado a lo largo de todo el recorrido del escaneo. Se puede
observar como todos los valores de posiciones y el error comprendido entre los valores obtenidos se reducen drásticamente con la
implementación del algoritmo FP, ya que el estimado que se obtiene sobre el error de su posición se mejora en comparación con los
otros métodos.
Para llevar a cabo la evaluación del proceso de medición de nuestro sensor, se calcularon los valores de incertidumbre de la posición en cada punto de medición, siguiendo los cálculos establecidos en [58], mostrados en la Tabla 2 y en la Figura 19. Finalmente, la Tabla 3 muestra el error promedio de los valores de incertidumbre de las mediciones obtenidas en la Tabla 2, analizando los resultados obtenidos se puede ver que estos valores representan el 1% del error en la medición calculado a partir de la dispersión de los datos con respecto al valor promedio de las medidas obtenidas por el sensor.
Distancia de referencia | Mezcla de gaussianas | Incertidumbre | FK | Incertidumbre | FP | Incertidumbre |
0 | 0 | 0 | 0 | 0 | 0 | 0 |
10 | 11.072 | 0.8 | 10.881 | 0.7 | 10.481 | 0.3 |
20 | 20.87 | 0.7 | 17.521 | 2.0 | 21.131 | 0.9 |
30 | 31.215 | 0.9 | 28.324 | 1.3 | 30.796 | 0.6 |
40 | 41.565 | 1.2 | 41.119 | 0.9 | 40.41 | 0.3 |
50 | 52.148 | 1.7 | 51.555 | 1.2 | 50.903 | 0.7 |
60 | 62.774 | 2.2 | 62.508 | 2.0 | 61.989 | 1.6 |
70 | 72.5820 | 2.1 | 72.363 | 1.9 | 69.76 | 0.1 |
80 | 83.915 | 3.1 | 83.677 | 3.0 | 80.884 | 0.7 |
90 | 93.688 | 3.0 | 93.401 | 2.7 | 92.349 | 1.9 |
100 | 103.133 | 2.5 | 102.634 | 2.1 | 99.963 | 0.3 |
110 | 112.371 | 1.9 | 112.03 | 1.6 | 111.415 | 1.1 |
120 | 122.1810 | 1.7 | 121.251 | 1.0 | 119.466 | 0.4 |
130 | 131.388 | 1.1 | 130.786 | 0.6 | 130.326 | 0.2 |
140 | 142.619 | 2.1 | 141.854 | 1.5 | 141.056 | 0.8 |
Tabla 3. Incertidumbres obtenidas de estimaciónes de posición discretas (mm)
incertidumbre | Mezcla de gausssianas | FK | FP |
---|---|---|---|
1.8 | 1.64 | 0.73 |
Los resultados experimentales muestran que el error de medición del sensor propuesto está dentro del 1%. Este resultado se compara a continuación con el obtenido en otros trabajos, en los que se presentan sistemas que proporcionan información útil para aplicaciones de navegación e inspección de tuberías y modelado 3D. Autores como [59] presentan un sistema de visión capaz de obtener información de la profundidad con un error de estimación menor al 2%. Un sistema de visión para el escaneo de tuberías 3D como el propuesto en [60], en el método de escaneo consiste en proyectar un haz de luz circular sobre la superficie del tubo, que permite calcular la posición de escaneo con una resolución de menos del 1%. En [61] demuestran cómo se puede realizar la medición del perfil interno de piezas cilíndricas con una resolución del 2%. Mientras que en [62] proponen un sistema de medición para la reconstrucción 3D de tuberías, el método consiste en estimar el movimiento de la cámara con referencia las coordenadas 3D de la sección de luz, este sistema puede estimar el movimiento de la cámara con un error del 3%. En [63] proponen un método de localización simultánea para estimar con precisión la posición de un instrumento utilizado con aplicaciones quirúrgicas, con una precisión menor del 2%. En [64] proponen un sistema para la inspección de tuberías que consiste en un robot móvil que se desplaza en el interior del ducto, la posición del móvil es medida a través de un distanciómetro láser con una precisión del 3%. En [65] presentan una técnica de estimación de la posición de la cámara dentro de la tubería con una precisión del 1%. De los resultados obtenidos con el sistema de medición propuesto en nuestro trabajo, se deduce que la precisión de nuestro sensor puede considerarse del mismo nivel que los trabajados mencionados, lo que demuestra que este sensor es adecuado para utilizarse en procesos donde es posible adaptar mediciones ópticas donde se espera que el sistema sea capaz de evaluar el perfil de la superficie interna de la pieza a medir.
![]() |
Figura 19. Comparación con otras técnicas en términos de precisión |
Se ha presentado un sistema para el escaneo de interiores en 3D, donde definimos un sistema para el dimensionamiento en 2D de las secciones transversales del interior de partes, nos enfocamos en realizar un sistema de medición para determinar de la posición del LED de potencia que fusiona sensores odométrico, magnético y óptico. Los resultados iniciales de las posiciones del LED llevaron a la observación de que es preciso que la referencia del sistema mantenga una posición de referencia fija, para que la asociación de datos no falle y las mediciones del sistema lleven una incorrecta estimación o error de la posición. En los resultados mostramos cómo las posiciones del LED pueden ser estimadas. El modelo lineal del sistema utilizado se basa en la asociación de datos, resuelto mediante mezcla de gaussianas con EM y después se desarrollaron y probaron los métodos de FK y FP con el propósito de estimar la posición del sistema de dimensionamiento, mostramos los resultados con la finalidad de destacar que ambos filtros tienen comportamientos similares, se encontró que existe un error sistemático de 1.64 mm y 0.73 mm respectivamente logrando de esta forma obtener resultados que permiten mejorar la medición de la posición del LED. Los resultados del posicionamiento obtenidos con el sistema de medición propuesto se comparan con los resultados de trabajos existentes y se presentan en la Figura 19, la solución que proporcionan estas técnicas son cercanos a los presentados por nuestro sistema, lo que indica que el sistema de escaneo se puede utilizar para diferentes aplicaciones
[1] Zhang Z. Flexible camera calibration by viewing a plane from unknown orientations. Proceedings of the Seventh IEEE International Conference on Computer Vision, pp. 666-673, 1999.
[2] Rupp S., Winter C., Wittenberg T. Camera calibration from fiberscopic views with accuracy evaluation. Bildverarbeitung für die Medizin, pp. 424-428, 2006.
[3] Icasio-Hernández O., Gonzalez-Barbosa J.-J., Hurtado-Ramos J., Viliesid-Alonso M. 3d reconstruction of hollow parts analyzing images acquired by a fiberscope. Measurement Science and Technology 25 (7), 2014.
[4] Chen L., Bavigadda V., Kofidis T., Howe R. Fiber optic projection-imaging system for shape measurement in confined space. The ScientificWorld Journal, 28:206569, 2014.
[5] Stehle T., Truhn D., Aach T., Trautwein C., Tischendorf J. Camera calibration for fish-eye lenses in endoscopy with an application to 3D reconstruction. In Proceedings IEEE International Symposium on Biomedical Imaging (ISBI), IEEE, Washington D.C., pp. 1176-1179, 2007.
[6] Barreto J., Santos J. Ray-based calibration of rigid medical endoscopes. The 8th Workshop Omnidirectional Vision Camera Networks and Non Classical Cameras (OMNIVIS), Marseille, France, 2008.
[7] Barreto J.P., Roquette J., Sturm P., Fonseca F. Automatic camera calibration applied to medical endoscopy. Procedings of the British Machine Vision Conference, 52.1-52.10, 2009.
[8] Melo R., Barreto J., Falcão G. A new solution for camera calibration and real-time image distortion correction in medical endoscopy-initial technical evaluation. IEEE Transactions on Bio-Medical Engineering, 59(3):634-44, 2012.
[9] Schmalz C., Forster F., Schick A., Angelopoulou E. An endoscopic 3d scanner based on structured light. Medical Image Analysis, 16(5):1063-1072, 2012.
[10] Hayashibe M., Suzuki N., Nakamura Y. Laser-scan endoscope system for intraoperative geometry acquisition and surgical robot safety management. Medical Image Analysis, 10(4):509-19, 2006.
[11] Wu C.H., Sun Y.N., Chang C.C. Three-dimensional modeling from endoscopic video using geometric constraints via feature positioning. IEEE Transactions on Biomedical Engineering, 54(7):1199-1211, 2007.
[12] ] C. Wu, Narasimhan S., Jaramaz B. Shape-from-shading under near point lighting and partial views for orthopedic endoscopy. In Proceedings of the First Int. Workshop on Photometric Analysis For Computer Vision - PACV 2007, 8 pp., Rio de Janeiro, Brazil, 2007.
[13] Wang G.H., Han J.Q., Zhang X.M. Three-dimensional reconstruction of endoscope images by a fast shape from shading method. Measurement Science and Technology, 20 (12):125801, 2009.
[14] Criante L., Lucesoli A., Farabollini B., Bonifazi F., Rozzi T., Simoni F. Size measurement in endoscopic images by low coherence interferometry. Journal of Optics A: Pure and Applied Optics, 11(3):034007, 2009.
[15] Penne J., Höller K., Stürmer M., Schrauder T., Schneider A., Engelbrecht R., Feussner H., Schmauss B., Hornegger J. Time-of-flight 3-D endoscopy. International Conference on Medical Image Computing and Computer-Assisted Intervention (MICCAI 2009), 12 (Pt 1), 467-474, 2009.
[16] Daul C., Blondel W. From 2D towards 3D cartography of hollow organs. 7th IEEE International Conference on ElectricalEngineering, Computing Science and Automatic Control, CCE 2010, CCE 201, pp. 285-293, Tuxtla Gutiérrez, Mexico, Sept. 2010.
[17] Seki T., Oka K., Naganawa A., Yamashita H., Kim K., ChibaT.T. Laser distance measurement using a newly developed composite-type optical fiberscope for fetoscopic laser surgery. Optics and Lasers in Engineering, 48(10):974-977, 2010.
[18] Arai K. Method for 3D object reconstruction using several portions of 2D images from the diferent aspects acquired with image scopes included in the fiber retractor. International Journal of Advanced Research in Artificial Intelligence, 1(9):14-19, 2012.
[19] Schmalz C., Forster F., Schick A., Angelopoulou E. An endoscopic 3d scanner based on structured light. Medical Image Analysis, 16(5):1063-1072, 2012.
[20] Hartrumpf M., Munser R. Optical three-dimensional measurements by radially symmetric structured light projection. Applied Optics, 36(13):2923-2928, 1997.
[21] Schaack D. Making accurate, three-dimensional measurements through a standard borescope. Proc. SPIE 3397, Nondestructive Evaluation of Aging Aircraft, Airports, and Aerospace Hardware II, 3397:264-276, 1998.
[22] Hasegawa K., Noda K., Sato Y. Electronic endoscope system for shape measurement. IEEE: Object Recognition Supported by User interaction for Service Robots, 1:792,795, Quebec, Canada, 2002.
[23] Winter C., Scholz I. Reconstruction of tubes from monocular fiberscopic images. Application and first results. In Vision, Modeling, and Visualization, pp. 57-64, 2005.
[24] Biegelbauer G., Vincze M. 3D vision-guided bore inspection system. In Fourth IEEE International Conference on Computer Vision Systems (ICVS'06), pp. 22- 22, 2006.
[25] Ford H.D., Tatam R.P. Passive OCT probe head for 3D duct inspection. Measurement Science and Technology, 24(9):094001, 2013.
[26] Müller A., Schubert M., Beleites E. Noncontact three-dimensional laser measuring device for tracheoscopy. Annals of Otology, Rhinology and Laryngology, 111(9):821-827, 2002.
[27] Dörfiel W., Sugano Y., Stalling D., Coconu L., Hentschel D., Linÿ G., Witt C. Laser-based endoscopic measurement of airway dimensions. Pneumologie, 57(9):503-509, 2003.
[28] Williamson J., James A., Phillips M., Sampson D., Hillman D., Eastwood P. Quantifying tracheobronchial tree dimensions: Methods, limitations and emerging techniques. European Respiratory Journal, 34(1):42-55, (2009).
[29] Geng J., Xie J. Review of 3-d endoscopic surface imaging techniques. IEEE Sensors Journal, 14(4):945-960, 2014.
[30] Yang D., Liu Z., Sun F., Zhang J., Liu H., Wang S. Recursive depth parametrization of monocular visual navigation: Observability analysis and performance evaluation. Information Sciences, 287:38-49, 2014.
[31] Firouzi H., Najjaran H. Real-time monocular vision-based object tracking with object distance and motion estimation. IEEE/ASME International Conference on Advanced Intelligent Mechatronics, (AIM 2010), pp. 987-992, 2010.
[32] Zhang Z., Han Y., Zhou Y., Dai M. A novel absolute localization estimation of a target with monocular vision. Optik, 124(12):1218-1223, 2013.
[33] Assa A., Janabi-Sharifi F. Virtual visual servoing for multicamera pose estimation. IEEE/ASME Transactions on Mechatronics, 20(2):789-798, 2014.
[34] Burschka D., Mair E. Direct pose estimation with a monocular camera. Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics), 4931 LNCS, pp. 440-453, 2008.
[35] Williams B., Klein G., Reid I. Real-time slam relocalisation. IEEE 11th International Conference on Computer Vision, pp. 1-8, Rio de Janeiro, 2007.
[36] Owczarek P., Goslinski J., Rybarczyk D., Kubacki A. Control of an electro-hydraulic manipulator by vision system using central point of a marker estimated via kalman filter. ICA 2016. Advances in Intelligent Systems and Computing, Springer International Publishing, pp. 587-596, 2016.
[37] Stroupe A.W., Martin M.C., Balch T. Merging Gaussian distributions for object localization in multi-robot systems. In Experimental Robotics VII, Lecture Notes in Control and Information Sciences, Vol. 271, pp. 343-352, Springer, Berlin, Heidelberg, 2001.
[38] Fox D., Burgard W., Thrun S. Markov localization for mobile robots in dynamic environments. Journal of Artificial Intelligence Research, 11:391-427, 1999.
[39] O'Kane J. Global localization using odometry. Proceedings 2006 IEEE International Conference on Robotics and Automation (ICRA 2006), pp. 37-42, 2006.
[40] Delahoche L., Pégard C., Marhic B. Vasseur P. Navigation system based on an omnidirectional vision sensor. In Proc. of the Int. Conf. on Intelligent Robotics and Systems 1997 (IROS’97), Vol. 2, pp. 718–724, Grenoble, France, 1997.
[41] Eberhard P., Tang Q. Sensor data fusion for the localization and position control of one kind of omnidirectional mobile robots. Multibody System Dynamics, Robotics and Control, Springer, Chapter 4, pp. 45-73, 2013.
[42] Caruso M., Bratland T., Smith C., Schneider R. Anisotropic magnetoresistive sensors: Theory and applications. Sensors (Peterborough, NH), 16(3):18-26, 1999.
[43] Schott C., Racz R., Betschart F., Popovic R. A new two-axis magnetic position sensor. SENSORS, 2002 IEEE, Vol. 2, pp. 911-915, 2002.
[44] Bosch-Jorge M., Soriano A., Valera A., Sánchez Salmerón A.J. Detección y seguimiento de objetos móviles utilizando un escáner láser mediante filtros de kalman. Comité Español de Automática (CEA-IFAC), Valencia (Spain), 2014.
[45] Woo J., Kim Y., Lee J., Lim M. Localization of mobile robot using particle filter. In 2006 SICE-ICASE International Joint Conference, pp. 3031-3034, 2006.
[46] Vadakkepat P., Jing L. Improved particle filter in sensor fusion for tracking randomly moving object. IEEE Transactions on Instrumentation and Measurement, 55(5):1823-1832, 2006.
[47] Yuan C., Lianchuan M., Ma W. Mobile target tracking based on hybrid open-loop monocular vision motion control strategy. Discrete Dynamics in Nature and Society, 10 pp., 2015.
[48] Boukhers Z., Shirahama K., Li F., Grzegorzek M. Object detection and depth estimation for 3d trajectory extraction. In Proceedings of International Workshop on Content-Based Multimedia Indexing, Vol. 2015, July 2015.
[49] Sakharkar S.S., Kamble S.D., Khobragade A.S. Object detection and tracking using particle filters. International Journal of Computer Trends and Technology(IJCTT), Vol. 22(1), April 2015.
[50] Bouguet J.-Y. Camera calibration toolbox for Matlab, 2015.
[51] Bouguet J.-Y. Camera calibration toolbox for Matlab, 2008.
[52] Y.-I. Curiel-Razo, O. Icasio-Hernández, G. Sepúlveda-Cervantes, J.-B. Hurtado-Ramos, J.-J. González-Barbosa, Leap motion controller three dimensional verification and polynomial correction, Measurement 93 (2016) 258 - 264. doi:http://dx.doi.org/10.1016/j.measurement.2016.07. 017.
[53] D. Elsaesser, Sensor data fusion using a probability density grid, FUSION 2007 - 2007 10th International Conference on Information Fusiondoi:10.1109/ICIF.2007.4408047.
[54] S. Schlegel, N. Korn, G. Scheuermann, On the interpolation of data with normally distributed uncertainty for visualization, IEEE Transactions on Visualization and Computer Graphics 18 (12) (2012) 2305-2314. doi: 10.1109/TVCG.2012.249
[55] S. Thrun, W. Burgard, D. Fox, Probabilistic robotics (intelligent robotics and autonomous agents).
[56] C. Winter, S. Rupp, M. Elter, C. Münzenmayer, H. Gerhäuser, T. Wit tenberg, Automatic adaptive enhancement for images obtained with fiberscopic endoscopes., IEEE transactions on bio-medical engineering 53 (10) (2006) 2035-46. doi:10.1109/TBME.2006.877110.
[57] L. Xu, Q. Yan, Y. Xia, J. Jia, Structure extraction from texture via relative total variation, ACM Transactions on Graphics (TOG).
[58]ISO, 1995. Guide to Expression of Uncertainty in Measurement, 1st ed. International Standard Organization.
[59] Orghidan, R., Salvi, J., Mouaddib, E. M., 2005. Accuracy estimation of a new omnidirectional 3d vision sensor. En: Image Processing, 2005. ICIP 2005. IEEE International Conference on. Vol. 3. IEEE, pp. III–365.
[60] Totani, H., Goto, H., Ikeda, M., Yada, T., 1995. Miniature 3D optical scanning sensor for pipe-inspection robot. En: Microrobotics and Micromechanical Systems. Vol. 2593. International Society for Optics and Photonics, pp. 21–30.
[61] Wakayama, T., Yoshizawa, T., 2013. Uni-axial inner profile measurement. Dimensional Optical Metrology and Inspection for Practical Applications II. Vol. 8839. International Society for Optics and Photonics, p. 88390L.
[62] Matsui, K., Yamashita, A., Kaneko, T., 2010. 3-d shape measurement of pipe by range finder constructed with omni-directional laser and omni-directional camera. Robotics and Automation (ICRA), 2010 IEEE International Conference on. IEEE, pp. 2537–2542.
[63]Sadjadi, H., Hashtrudi-Zaad, K., Fichtinger, G., 2016. Simultaneous localization and calibration for electromagnetic tracking systems. The International Journal of Medical Robotics and Computer Assisted Surgery 12 (2), 189–198.
[64]Liu, Z., Krys, D., 2012. The use of laser range finder on a robotic platform for pipe inspection. Mechanical systems and signal processing 31, 246–257.
[65] Tezerjani, A. D., Mehrandezh, M., Paranjape, R., 2014. 4-dof pose estimation of a pipe crawling robot using a collimated laser, a conic mirror, and a fish-eye camera. En: Image Analysis and Interpretation (SSIAI), 2014 IEEE Southwest Symposium on. IEEE, pp. 45–48.
Published on 08/02/19
Accepted on 03/05/18
Submitted on 12/09/17
Volume 35, Issue 1, 2019
DOI: 10.23967/j.rimni.2018.06.001
Licence: CC BY-NC-SA license