Modelo para evaluación de la exactitud de escáner láser terrestre - TLS
Model for evaluation of the accuracy of terrestrial laser scanner - TLS
Modelo para avaliação da precisão do scanner a laser terrestre - TLS
William Barragán Zaque* | Lorena Lizeth Herrera de la Rosa | Lilia Edith Aparicio Pico
Universidad Distrital Francisco José de Caldas, Bogotá, Colombia
Resumen
El presente trabajo determina la exactitud de la nube de puntos capturada por un láser escáner terrestre (TLS, por sus siglas en inglés), con el objetivo de desarrollar un modelo estadístico que permita la evaluación de la exactitud del escáner terrestre. Este modelo se basa en la experimentación, con variación de parámetros de referencia en mediciones directas. Los datos se obtienen por medio de procedimientos que incluyen capturar información en diferentes tipos de distancias, ángulos y
superficies (target). Para dar contundencia a la validez del modelo obtenido, los datos conseguidos son comparados con respecto a otros adquiridos con una estación topográfica total. El desarrollo matemático y estadístico del modelo utiliza la teoría del diseño del experimento, donde las mediciones realizadas en cada uno de los escaneos son independientes entre sí, al igual que cada tipo target. Proponiendo la implementación del análisis de modelo lineal general (MLG), para ajustar una expresión matemática al error medio (distancias euclidianas entre las coordenadas observadas y teóricas). El modelo desarrollado permitió determinar la exactitud de los datos adquiridos por el láser terrestre, estableciendo que por cada metro adicional en la distancia del escáner con relación a una superficie de estudio, el error medio esperado aumentará entre 0,033% y 1,5%.
Palabras clave: escáner láser terrestre, evaluación de precisión, análisis de precisión, exactitud TLS
Summary
The present work determines the accuracy of the point cloud captured by a terrestrial laser scanner (TLS), with the objective of developing a statistical model that allows the evaluation of the accuracy of the terrestrial scanner. This model is based on experimentation, with a variation of reference parameters in direct measurements. The data is obtained through procedures that include capturing information on different types of distances, angles and surfaces (target). To give force to the validity of the model obtained, the data obtained are compared with respect to that which was acquired with a topographic station total. The mathematical and statistical development of the model uses the theory of the design of the experiment, where the measurements made in each of the scans are independent of each other, as is each target type. Proposing the implementation of the general linear model analysis (GLM) to adjust a mathematical expression to the average error (Euclidean distances between the observed and theoretical coordinates). The developed model allowed for the determination of the accuracy of the data acquired by the terrestrial laser, establishing that for each additional meter in the distance of the scanner in relation to a study surface, the expected average error will increase by between 0.033% and 1.5%.
Keywords: terrestrial laser scanner, precision evaluation, precision analysis, TLS accuracy
Sumário
O presente trabalho determina a precisão da nuvem de pontos capturada por um scanner a laser terrestre (TLS), com o objetivo de desenvolver um modelo estatístico que permita avaliar a precisão do scanner terrestre. Este modelo é baseado em experimentação, com variação de parâmetros de referência em medições diretas. Os dados são obtidos através de procedimentos que incluem captar informações sobre diferentes tipos de distâncias, ângulos e superfícies (alvo). Para dar força à validade do modelo obtido, os dados captados são comparados com outros adquiridos com uma estação topográfica total. O desenvolvimento matemático e estatístico do modelo usa a teoria do desenho do experimento, onde as medidas feitas em cada uma das varreduras são independentes uma da outra, assim como cada tipo de alvo. Propondo a implementação da análise do modelo linear geral (MLG), para ajustar uma expressão matemática ao erro médio (distâncias euclidianas entre as coordenadas observadas e teóricas). O modelo desenvolvido permitiu determinar a precisão dos dados adquiridos pelo laser terrestre, estabelecendo que para cada metro adicional à distância do scanner em relação à superfície de estudo, o erro médio esperado aumentará entre 0,033% e 1,5%.
Palavras-chave: : Scanner a laser terrestre, avaliação de precisão, análise de precisão, precisão TLS
El Grupo de Investigación GEOTOPO1 realizó el diseño metodológico de esta investigación, que sometió a prueba de campo y el propósito fue estudiar la exactitud de los datos capturados por un escáner láser terrestre (TLS) marca y referencia Faro Focus 3D (FARO, 2018), perteneciente al laboratorio de topografía de la Universidad Distrital Francisco José de Caldas, en Bogotá, Colombia. El estudio permitió la construcción de un modelo soportado estadísticamente, que comprobó la exactitud de los datos capturados, todo ello en razón a que las especificaciones técnicas del fabricante del escáner mencionado se encuentran delimitadas a unos rangos de distancias entre los 0,6 m y 120 m para los cuales se establece un error sistemático de ± 2 mm por cada intervalo entre 10 m y 25 m. Razón por la cual se presenta la necesidad de hallar un modelo estadístico que permita establecer los rangos de error de los datos adquiridos teniendo en cuenta los ángulos de incidencia (Ramos, Marchamalo & Martínez, 2017; Voegtle & Ruoming, 2009), distancias de captura (Voegtle & Schwab, 2008), y targets u objetivos de captu-ra para el escaneo, para evidenciar de manera estadística la influencia que estas variables tienen en los errores de medición de los datos capturados por un TLS Faro Focus 3D (Kaasalainen, Jaakkola, Kaasalainen, Krooks & Kukko, 2011; Soudarissanane, Lindenberg & Menenti, 2009).
De esta manera, cuando se capturan datos sobre un
escenario y se genera un modelo en 3D de las superficies (Al-Durgham, Habib & Kwak, 2013), estas mediciones son vulnerables a los efectos de ruido que se presentan en el montaje de un proyecto (Gressin, Mallet, Demantké & Da, 2013), debido a muchos factores que pueden distorsionar las dimensiones de los objetos. Por ello, se propone el desarrollo de un modelo estadístico que incluya las variables de objetos o targets (Zhou, Cui & Yang 2009), distancia y ángulos de incidencia que pueden afectar la exactitud de los datos capturados dentro de la generación de modelos 3D (Lichti & Licht, 2006; Ramos et al., 2015), a partir de datos obtenidos de un TLS.
Esta investigación tiene como objetivo principal determinar un modelo estadístico que represente la precisión de los datos adquiridos por un TLS (Lohani & Sasidhaan, 2017) mediante la realización de pruebas de campo en diferentes escenarios donde deben interactuar las variables de distancia, ángulo de incidencia y objetivos o targets a escanear (Chiba & Masuda, 2016; López & Ruiz, 2011). Dicho objetivo se plantea con la intención de establecer los errores en medición que se tienen al implementar estas nuevas tecnologías que son de gran avance para el campo de la topografía (Barragán & Escobar, 2017). En la primera parte de este artículo se establece la metodología de la investigación en donde se define el diseño del experimento y los modelos estadísticos a implementar para obtener el
mejor ajuste posible de datos que permita establecer la exactitud de las mediciones realizadas con el TLS Faro Focus 3D frente a las obtenidas por topografía convencional (Aoki, Aoki & Miyamoto, 2017; Masuda & Tanaka, 2010), a partir de los cuales se presentan los resultados y conclusiones obtenidas de la investigación, con la obtención del modelo estadístico estimado, validado e interpretado que se ajusta al objetivo principal de este estudio (Date, Yokoyama & Sugawara, 2018; Takai, Date & Kanai, 2013).
Construcción de la metodología
Diseño del experimento
Se plantea el desarrollo de un escenario, donde se pueda escanear una superficie en un rango de distancia visible de 0 m a 120 m, teniendo en cuenta las especificaciones técnicas del TLS Faro Focus 3D, para lo cual se definió una locación que cumplía con dicha condición dentro del campus de la Universidad Nacional de Colombia, en un costado del Estadio Alfonso López (foto 1).
Luego de establecer la zona de estudio (foto 2), se procede a diseñar la localización de cada uno de los puntos, sobre los cuales se tendrán en cuenta la variación de ángulos y distancias con relación a la superficie de captura. Para cumplir estas condiciones, se hace el diseño del esquema en oficina, de cada uno de los tres ejes con ángulos de 30°, 60° y 90° y distancias conocidas (figura 1). El diseño establecido se materializa en campo por medio de técnicas de levantamiento topográfico convencional, con una estación topográfica total, que hoy día es considerada
una medición extremadamente confiable (Barragán & Martínez, 2015).
Cada uno de los puntos materializados en campo, es utilizado para ubicar sobre cada uno de ellos el TLS Faro Focus 3D, y realizar cada uno de los escaneos propuestos.
|
|
|
Figura 1. Diagrama de diseño de experimento en planta. |
Diseño de objetos o targets de reconocimiento
Para la identificación de objetos dentro de los escaneos fue necesario diseñar una serie de targets u objetivos de reconocimiento para los escaneos (figura 2), los cuales fueron pensados para que se distinguieran en diversos colores, materiales, tamaños (véase tabla 1) y formas, que además fueran de fácil elaboración, accesibilidad a sus materias primas y representativos en los entornos naturales del espacio físico que tradicionalmente se documenta por medio del escáner terrestre.
Tabla 1. Dimensiones de los tipos de target a escanear
|
TIPO DE TARGET |
DIMENSIONES(m) |
|
Tipo 1 |
0,08*0,08 |
|
Tipo 2 |
0,08*0,08 |
|
Tipo 3 |
0,05*0,085 |
|
Tipo 4 |
0,1*0,1 |
|
Tipo 5 |
0,1*0,1 |
|
Tipo 6 |
0,08*0,08 |
|
Tipo 7 |
Diámetro: 0,07 |
|
Tipo 8 |
Contorno: 1*0,7. Polígonos internos: 2. Polígonos de 0,40*0,15. Un polígono de 0,40*0,10 y dos polígonos 0,15*0,12. |
Implementación del diseño
En este proceso, el diseño del experimento definido en oficina se materializa en la práctica por medio del
levantamiento topográfico convencional. Seguidamente, se realiza el escaneado de la superficie de estudio desde diferentes puntos localizados en campo (Theiler, Wegner & Schindler, 2013).
Levantamiento topográfico convencional
Este se realiza según lo presentado en la figura 1, donde a partir del diagrama de diseño se muestra la localización de cada uno de los puntos, donde se hará el levantamiento con el escáner terrestre. Puntos colocados sobre los tres alineamientos definidos (en ángulos de 30°, 60° y 90° con respecto a la superficie de escaneo). La localización precisa de dichos puntos en el terreno se realiza por medio de métodos de topografía convencional, con estación topográfica total (foto 3).
Tras la demarcación en campo con estacas de cada uno de los puntos sobre los cuales se hacen los escaneos, se elige el punto ubicado a 15.595 m de distancia en un ángulo perpendicular a la superficie de escaneo. Precisamente, como base para realizar el levantamiento topográfico de los elementos dentro de esa superficie, incluidos la cantidad de objetivos de escaneo por cada uno de los tipos targets, se distribuyen aleatoriamente a lo largo de la superficie, teniendo en cuenta que se cubra de manera equitativa la totalidad del área.
Este levantamiento topográfico sirve de insumo en el proceso de identificación y extracción de coordenadas de los targets, para realizar un análisis diferencial de la exactitud de los datos capturados por medio del TLS con relación a las coordenadas reales levantadas por topografía convencional.
Levantamiento TLS
El TLS Faro Focus 3D es un equipo que captura gran cantidad de puntos a alta velocidad, que facilitan la medición y la documentación aplicables a diferentes campos de ejecución. Por otro lado, la información resultante de un escaneo está formada por nubes de puntos, que representan una reproducción digital exacta de la realidad capturada en un corto lapso de tiempo (foto 4).
Para la captura de datos por medio del TLS Faro Focus 3D, enmarcado en el diseño del experimento planteado, se procede con la ubicación del equipo, en cada uno de los puntos y desde allí hacer el levantamiento con los siguientes parámetros de configuración por escaneo.
Si se incrementa el factor de calidad, se reduce el ruido en los datos y, por tanto, aumenta la calidad del escaneo, lo que se traduce en un aumento del tiempo de adquisición (véase tabla 2). Respecto al nivel de calidad aplicado a un escaneo, se hace con tasas de medición o aplicando compresión de ruido adicional (Wujanz, Burger &
Mettenleiter, 2017).
Tabla 2. Parámetros de escaneo
|
ESCANEO |
TIEMPO |
RESOLUCIÓN |
CALIDAD |
|
1 |
00:15:06 |
1/2 |
4X |
|
2 |
00:10:21 |
1/2 |
4X |
|
3 |
00:15:06 |
1/2 |
4X |
|
4 |
00:08:48 |
1/2 |
4X |
|
5 |
00:07:56 |
1/2 |
4X |
|
6 |
00:07:56 |
1/2 |
4X |
|
7 |
00:07:56 |
1/2 |
4X |
|
8 |
00:04:47 |
1/2 |
4X |
Luego de la adquisición de datos capturados por el TLS Faro Focus 3D, estos deben administrarse y procesarse mediante el reconocimiento de objetos o targets establecidos en el experimento.
Resultados
Definición del modelo estadístico
Dentro de la especificación del modelo teórico se incluyen las variables protagonistas del desarrollo del mismo, las cuales se establecen como la distancia de escaneo (Ozendi, Akca & Topan, 2017), los targets u objetivos de reconocimiento del escaneo y los ángulo de incidencia propuestos (Soudarissanane, Lindenbergh & Menenti, 2011), que para el caso puntual de la modelación se van a tomar en sentido horario con respecto a la superficie de captura, por lo que equivaldría a hablar de ángulos de 60°, 90° y 150°, respectivamente, teniendo en cuenta el diagrama del diseño del experimento (Soudarissanane, Lindenbergh & Menenti, 2009).
En la búsqueda de un modelo estadístico, que explique el comportamiento de las observaciones obtenidas en campo de los targets levantados por medio de la topografía convencional y el TLS (Poreba & Goulette, 2015), se acudió a diferentes áreas de la estadística, entre las cuales se encuentran la geoestadística y los modelos lineales generalizados (MLG), siendo esta última la solución óptima para el fin del estudio. Por un lado, la geoestadística prometía ser una herramienta útil en la predicción del error obtenido al realizar las mediciones con el escáner láser, puesto que considera la correlación espacial que pueda existir entre las observaciones en un mismo escaneo; sin embargo, dichas predicciones tenían un horizonte de tan solo 20 metros y estaban limitadas al espacio donde se realizaron las medidas, el cual es bastante ideal; contrario al sinnúmero de configuraciones que pueden encontrarse en la práctica (Sumi, Date & Kanai, 2018).
A partir de los datos obtenidos en los diferentes escaneos realizados con el TLS y el del levantamiento por medio de topografía convencional con estación total, se hace un análisis estadístico diferencial que toma como reales las coordenadas obtenidas por este último levantamiento; es decir, como variable de referencia para hallar la exactitud de
los errores en los datos obtenidos con el TLS (Yang, Dong & Liang, 2016). De lo anterior, se hace el cálculo de la distancia euclidiana (véase ecuación 1) donde se obtienen los valores de error para cada uno de los targets.
(1)
Donde la distancia euclidiana entre x, y, z está dada por la raíz de la suma del cuadrado de las diferencias de las variables de la función. El resultado de esta ecuación representa el error de los datos en la medición con el TLS para cada uno de los targets, y a partir de este valor se
calculó la media por cada target dentro de cada uno
de los escaneos, que de ahora en adelante será entendida como la variable de respuesta para el desarrollo del modelo estadístico por medio de la implementación de MLG.
Haciendo uso del postulado donde las mediciones obtenidas entre targets son independientes y, asimismo,
entre los diferentes tipos de escaneo, se utilizó el análisis de MLG (López & Ruiz, 2011) para ajustar un modelo al error medio (media de las distancias euclidianas entre las coordenadas observadas y las teóricas) obtenido para cada escaneo en cada target (variable respuesta), teniendo como covariables el ángulo, la distancia y, claro está, el target (Yoshimura, Date & Kanai, 2016).
El proceso de estimación de parámetros se realiza calculando el valor de los coeficientes del modelo examinado a partir del conjunto de datos de las observaciones del error medio, para determinar si el modelo teórico
propuesto es aceptable como una representación aproximada que se ajusta de la mejor forma a los datos (Dong, Cai & Du, 2016; Xu, Boerner & Yao, 2017).
Durante la evaluación del modelo, se examinan y se comprueban que las observaciones individuales cumplan los supuestos de normalidad, linealidad (Cohen, 1962), homocedasticidad (Little, 2013) e independencia, características del MLG (Enderlein, McCullagh & Nelder, 1987). A su vez, la validación del modelo se realiza a partir de los gráficos de bandas y el criterio de la distancia de Cook, donde se verifica el ajuste del modelo con los valores del error medio de las distancias euclidianas por target en cada uno de los escaneos.
Para la implementación del análisis por medio de MLG se usa el software libre R donde se realiza el desarrollo del modelo estadístico propuesto en la investigación. Mediante la teoría del diseño del experimento, las mediciones realizadas en cada uno de los escaneos son independientes entre sí (ejemplo: escaneo 2 es independiente de los escaneos 1, 3, 4, 5, 6, 7, 8 y 9), ya que fueron ejecutadas en diferentes ángulos y distancias con respecto a la superficie de captura. Adicional a esto, los datos obtenidos de los tipos de targets también cumplen con esta condición de independencia, pues las mediciones obtenidas para cada tipo target son independientes según el tipo de target (ejemplo: las medidas del target tipo 1, son independientes a las medidas de los targets tipos 2, 3, 4, 5, 6, 7 y 8). Por lo que se propuso la implementación del análisis de MLG para ajustar un modelo al error medio (media de las distancias euclidianas entre las coordenadas observadas y las teóricas) obtenido para cada escaneo en cada target (variable respuesta), teniendo como covariables el ángulo, la distancia y, claro está, el target. Para la definición de la variable respuesta dentro del modelo, fue necesario descartar los escaneos 7 y 8, puesto que la resolución de estos era demasiada pobre y no se logró extraer información de coordenadas de los targets (Lu, Wang & Guo, 2018).
A partir de lo anterior se comenzó con la implementación del modelo en el software estadístico R, para lo cual se inició con la premisa de considerar que el error es una variable cuyo dominio son los reales positivos. Por lo que se propone un modelo de la familia de distribuciones Gamma para explicar el comportamiento de esta variable (Forstner
& Khoshelham, 2017). Para verificar si la familia de distribuciones Gamma es adecuada para el ajuste del modelo,
se emplean las estadísticas descriptivas (véase tabla 3).
De las estadísticas descriptivas del error medio, puede decirse que la variable respuesta es sesgada a la derecha (la media es mayor que la mediana). Para verificar este hecho e identificar la intensidad de dicho sesgo, se ilustra la densidad de la dicha variable respuesta (figura 3).
Tabla 3. Estadística descriptiva
|
Mínimo |
0,01203 |
|
Cuartil 1 |
0,02463 |
|
Mediana |
0,04068 |
|
Media |
0,05990 |
|
Cuartil 3 |
0,07245 |
|
Máximo |
0,27843 |
En la figura 3 se observa, a la derecha, un sesgo bastante pronunciado, razón por la cual se establece la familia de distribuciones Gamma para la estimación del modelo.
De los gráficos que se ilustran en las figuras 4 y 5, puede decirse que la dispersión no es constante entre las observaciones. Nótese que el error medio toma valores entre 0 y 0,15 aproximadamente, para todas las distancias medidas excepto para los 80 m donde se observa un posible dato atípico con un valor superior a 0,25 y una concentración de los puntos de cerca de 0 y 0,05. De manera análoga, en el gráfico de dispersión de los ángulos se nota la presencia de un posible valor atípico para los ángulos de 90° y 150°, y de nuevo (Ge, 2017) la dispersión de los puntos parece no ser homogénea entre ángulos. Por último, del box-plot en la figura 6, se observa que las cajas poseen anchos diferentes, lo cual indica que no existe la misma dispersión de las observaciones entre targets. Con base en lo anterior, se concluye que no es adecuado asumir una dispersión constante; por tal motivo, se decide ajustar un modelo lineal generalizado doble.
Considerando que la dispersión de la respuesta es más notoria entre targets, se propone el siguiente modelo estadístico (véase ecuación 2).
(2)
Donde:
Para la identificación de las variables cuyo efecto es estadísticamente significativo sobre el error medio esperado
se acudió al estadístico de Wald (Lemonte & Vanegas, 2005), el cual está dado por la siguiente expresión (véase ecuación 3).
(3)
Donde 
y 
corresponden a las estimaciones de 
y su error estándar asociado 
, esto para j = 2,…,10. El estadístico de Wald está asociado al sistema de hipótesis (véase ecuación 4).
(4)
Donde 
es equivalente a decir que el efecto de la variable 
no es estadísticamente significativo sobre el error medio esperado.
La regla de decisión del sistema de hipótesis es:
Donde 
n siendo el tamaño de muestra y p el número de parámetros del modelo.
De manera análoga, se tienen los sistemas de hipótesis asociados a los efectos de los 
sobre la dispersión, dados por la siguiente expresión (véase ecuación 5).
(5)
Esto para j = 2,…,8.
Se obtiene que para una significancia del 1% la distancia es la única variable que tiene un efecto estadísticamente significativo sobre el error medio esperado con un p-valor = 0,000974. Por otro lado, se concluyó que el parámetro 
es constante para todo k; es decir, las covariables asociadas a los targets no tienen un efecto estadísticamente significativo sobre la dispersión. Por lo cual, el modelo
final teórico queda definido como sigue (véase ecuación 6).
(6)
En la tabla 4 se ilustran los parámetros estimados, su error estándar y su p-valor asociado.
Tabla 4. Estimación de los parámetros
|
Parámetro |
Estimación |
Error |
Estadística de Wald |
p-valor |
|
|
-3,2142 |
0,2240 |
-14,349 |
5,5752 × 10-19 |
|
|
0,007689 |
0,002257 |
3,407 |
0,000974 |
|
|
-0,7155 |
0,1856 |
-3,8851 |
0,0001197 |
Nótese que el parámetro 
tiene un p-valor mayor al nivel de significancia del 1%; sin embargo, este se conserva por cuestiones de ajuste del modelo. Cabe resaltar que este no tiene interpretación práctica.
Reemplazando el valor de las estimaciones sobre el
modelo propuesto, se obtiene el siguiente modelo estadístico final estimado (véase ecuación 7).
(7)
Validación del modelo
Para la validación del modelo se emplean los gráficos de
bandas de 99% confianza y el criterio de la distancia
de Cook para la identificación de observaciones potencialmente influyentes. Al observar las figuras 7 y 8 que ilustran las bandas de confianza, se encuentra que existen dos puntos fuera de estas (observaciones 19 y 18) para el modelo asociado a la media; sin embargo, esta distancia no es considerable y puede explicarse por la variación aleatoria de los datos. Por otro lado, para el modelo asociado a la dispersión se observa que no existe punto alguno fuera de las bandas, a partir de esto se concluye que el modelo se ajusta de forma adecuada. En los gráficos asociados a la distancia de Cook que se ilustran en las figuras 9 y 10, se tiene de nuevo que la observación 18 presenta valores posiblemente atípicos y en este caso especial potencialmente influyentes; a esto se le suma que las observaciones 15 y 23 también se postulan como potencialmente influyentes; sin embargo, al extraer estas observaciones y realizar el ajuste del modelo no se generan cambios estadísticamente significativos en el efecto de la distancia sobre el error medio esperado.
Residuo componente del desvío para dispersión
Finalmente, como 
puede concluirse que el error medio esperado aumentará un 
por cada metro adicional en la distancia del escáner. Teniendo en cuenta que un intervalo del 99% de confianza para 
está dado por [
es decir [0,00033, 0,015], al sacar el exponencial a los límites se tiene que el intervalo queda dado por [1,00033, 1,01515], a partir del cual se concluye que por cada metro adicional en la distancia del escáner el error medio esperado aumentará entre 0,033% y 1,5%.
Conclusiones
Es posible obtener un modelo con base en fundamentación estadística para determinar la exactitud de los datos adquiridos por medio de un TLS Faro Focus 3D, a partir del cual se estableció que por cada metro adicional en la distancia del escáner con relación a una superficie de estudio, el error medio esperado aumentará entre 0,033% y 1,5%.
Dentro de las especificaciones técnicas del equipo Faro Focus 3D, se describe de manera generalizada el error en las medidas del equipo como ± 2 mm por cada intervalo entre 10 m y 25 m limitando dicho error en un rango de distancia máxima de hasta 120 m. A partir del desarrollo del modelo con base en las pruebas de campo se pudo establecer que a una distancia de 120 m el error de medición que se puede tener es calculado entre 52,3 mm y 122,9 mm.
A partir del análisis estadístico en el desarrollo del modelo, se pudo establecer que de las variables o parámetros involucrados inicialmente, la distancia es la que mayor incidencia tiene sobre los errores que pueden presentarse en los datos adquiridos con el TLS. Pero, por otro lado, la variación de ángulos y targets no presenta gran significancia dentro de la exactitud de los datos al compararlos con los obtenidos por topografía convencional. Sin embargo, en las pruebas de campo realizadas para adquirir datos por medio del TLS Faro Focus 3D, se pudo observar que la resolución del escaneo o la cantidad de puntos capturados disminuye para la misma superficie de estudio, al aumentar la distancia y al hacer escaneos desde ángulos entre dos rangos aproximados, el primero de 0° a 60° y el segundo de 120° a 180° si se toman los ángulos en sentido horario con relación a la superficie de captura.
Teniendo en cuenta la variable respuesta introducida en el modelo como la media de las distancias euclidianas entre las coordenadas observadas y las teóricas, se pudo establecer que el tipo de target que arroja mejor exactitud es el tipo 2 que corresponde con madera tipo cedro y, por el contrario, el que más error presenta es el tipo 6, hecho de papel brillante.
En el diseño del experimento se planeó comparar dos escaneos realizados desde un mismo punto mediante el TLS Faro Focus 3D, pero configurado con diferentes parámetros de calidad, como en el caso de los escaneos 8 y 9, tomados a 120 m a un ángulo de 90° con una resolución de 1/2 y variando el factor de calidad por 4x y 6x, respectivamente. Se pudo validar que el parámetro de calidad del TLS está relacionado con la reducción del ruido en los datos de escaneo, lo cual hace que el tiempo de captura se incremente y que exista gran pérdida de información, caso que no fue conveniente en esta investigación. Como consecuencia de ello, el escaneo 9 capturado con factor de calidad 6x, se descartó en la modelación estadística de los datos.
Recomendaciones
Teniendo en cuenta que, en el resultado de la implementación del modelo, se obtuvo que las variables de ángulo de incidencia y target no presentan una afectación significativa dentro de la exactitud de los datos adquiridos con el TLS Faro Focus 3D, sería un gran aporte realizar una investigación donde se pueda establecer la posible relación de estas dos variables con relación a la resolución del escaneo, es decir, si de estas depende la cantidad de puntos capturados en un escaneo.
El modelo estadístico generado aquí, se realizó teniendo en cuenta los errores medios de las distancias euclidianas entre el levantamiento topográfico con estación y el levantamiento por medio del TLS Faro Focus 3D, estableciendo para cada uno de esos errores una configuración de los parámetros de distancia, ángulo de incidencia y target. Si se quieren obtener los rangos de error para otro tipo de TLS, sería recomendable validar, por medio de los procedimientos plasmados en este trabajo, el ajuste de dichos datos al planteamiento de este modelo.
Referencias
Al-Durgham, K., Habib, A., & Kwak, E. (2013). RANSAC approach for automated registration of terrestrial laser scans using linear features. ISPRS Ann. Photogramm. Remote Sens. Spatial Inf. Sci., II-5/W2, 13-18.
Aoki, R., Aoki, S., & Miyamoto, R. (2017). Key point localization for 3d model generation from facial illustrations using SURF and color features (pp. 55-56). In IEEE 7th International Conference on Consumer Electronics-Berlin (ICCE-Berlin).
Barragán, W., & Escobar, K. (2017). Obtención de parámetros óptimos en la clasificación de nubes de puntos LiDAR, a partir de sensores aerotransportados. Revista Avances Investigación en Ingeniería, 14(1). Bogotá: Universidad Libre.
Barragán, W., & Martínez A. (2015). Generación automática de cartografía para edificaciones, utilizando fotografías digitales aéreas verticales de alta resolución espacial y nubes de puntos. Revista de Ingeniería, 42. Bogotá: Universidad de los Andes.
Chiba, A., & Masuda, H. (2016). Reconstruction of polygonal prisms from point-clouds of engineering facilities. Journal of Computational Design and Engineering, 3, 322-329.
Cohen, J. (1962). The statistical power of abnormal-social psychological research: A review. AThe Journal of Abnormal and Social Psychology, 65(3), 145-153.
Date, H., Yokoyama, T., & Sugawara, T. (2018). Efficient registration of laser-scanned point clouds of bridges using linear features. International Journal of Automation Technology, 328-338.
Dong, J., Cai, Z., & Du, S. (2016). Improvement of affine iterative closest point algorithm for partial registration. IET Computer Vision, 11(2), 135-144.
Enderlein, G., McCullagh, P., & Nelder, J. A. (1987). Generalized linear models. Chapman and Hall. London, New York, 1983. Biometrical Journal, 29(2), 206-206 [sic].
Faro Focus. (2018) | FARO SPAIN, S.L.U. [en línea]. Accedido: 24-may-2018. Disponible en https://www.faro.com/es es/productos/construccion-bim-cim/faro-focus.
Forstner, W., & Khoshelham, K. (2017). Efficient and accurate registration of point clouds with plane to plane correspondences (pp. 2165-2173). In Proceedings of the IEEE International Conference on Computer Vision.
Ge, X. (2017). Automatic markerless registration of point clouds with semantic-keypoint-based 4 points congruent sets. ISPRS Journal of Photogrammetry and Remote Sensing, 130, 344-357.
Gressin, A., Mallet, C., Demantké, J., & Da, N. (2013). Towards 3D lidar point cloud registration improvement using optimal neighborhood knowledge. ISPRS Journal of Photogrammetry and Remote Sensing, 79, 240-251.
Kaasalainen, S., Jaakkola, A., Kaasalainen, M., Krooks, A., & Kukko, A. (2011). Analysis of incidence angle and distance effects on terrestrial laser scanner intensity: Search for correction methods. Journal Remote Sensing ISPRS, 59-68.
Lemonte, A. J., & Vanegas, L. H. (2005). Una comparación entre la inferencia basada en las estadísticas de Wald y razón de verosimilitud en los modelos probabilísticos. Revista Colombiana de Estadística, 28, 77-97.
Li, W., & Song, P. (2015). A modified ICP algorithm based on dynamic adjustment factor for registration of point cloud and CAD model. Pattern Recognit., 65, 88-94.
Lichti, D. D., & Licht, M. G. (2006). Experiences with terrestrial laser scanner modelling and accuracy assessment. IAPRS, XXXVI. Germany: Dresden.
Little, T. D. (2013). The Oxford Handbook of Quantitative Methods in Psychology, vol. 2: Statistical analysis. OUP USA.
Lohani, B., & Sasidhaan, S. (2017). An evaluation of intensity augmented ICP for terrestrial LiDAR data registration. Journal of Geomatics, 11(2), 139-148.
López, G. E., & Ruiz, S. M. (2011). Análisis de datos con el modelo lineal generalizado. Una aplicación con R. Revista Española de Pedagogía, 69(248), 59-80.
Lu, J., Wang, W., & Guo, C. (2018). Point cloud registration based on CPD algorithm (pp. 8235-8240). 37th Chinese Control Conference (CCC) IEEE.
Masuda, H., & Tanaka, I. (2010). As-built 3D modeling of large facilities based on interactive feature editing. Computer-Aided Design and Applications, 7(3), 349-360.
Mellado, N., Aiger, D., & Mitra, N. J. (2014). Super 4pcs fast global point cloud registration via smart indexing. Computer Graphics Forum, 33(5), 205-215.
Ozendi, M., Akca, D., & Topan, H. (2017). A generic point error model for TLS derived point clouds. In: Proceedings of SPIE, 10332, 103320J.
Poreba, M., & Goulette, F. (2015). A robust linear featurebased procedure for automated registration of point clouds. Sensors, 15, 1435-1457.
Ramos, L., Marchamalo, M., & Martínez, R. (2017). Estimating and plotting TLS midrange precisions in field conditions: Application to dam monitoring. International Journal of Civil Engineering, 15, 299-307.
Ramos, L., Marchamalo, M., Rejas, J. G., & Martínez, R. (2015). Aplicación del láser escáner terrestre (TLS) a la modelización de estructuras: precisión, exactitud y diseño de la adquisición de datos en casos reales. Informes de la Construcción, 67(538): e074. http://dx.doi.org/10.3989/ic.13.103
Saeedi, S., Trentini, M., Seto, M., & Li, H. (2016). Multiple robot simultaneous localization and mapping: A review. Journal of Field Robotics, 33(1), 3-46.
Soudarissanane, S., Lindenbergh, R., & Menenti, M. (2011). Scanning geometry: Influencing factor on the quality of terrestrial laser scanning points. ISPRS Journal of Photogrammetry and Remote Sensing, 66, 389-399.
Soudarissanane, S., Lindenbergh, R., Menenti, M. (2009). Incidence angle influence on the quality of terrestrial laser scanning points. Proceedings of ISPRS Workshop Laser Scanning, 183-188.
Sumi, T., Date, H., & Kanai, S. (2018). Multiple TLS point cloud registration based on point projection images. Int. Arch. Photogramm. Remote Sens. Spatial Inf. Sci., XLII-2, 1083-1090.
Takai, S., Date, H., & Kanai, S. (2013). Accurate registration of MMS point clouds of urban areas using trajectory. ISPRS Ann. Photogramm. Remote Sens. Spatial Inf. Sci., II-5/W2, 277-282.
Theiler, P. W., Wegner, J. D., & Schindler, K. (2013). Markerless point cloud registration with keypoint-based 4-points congruent sets. ISPRS Ann. Photogramm. Remote Sens. Spatial Inf. Sci., II-5/W2, 283-288.
Voegtle, A., & Ruoming, R. (2009). Effects on the measurements of the terrestrial laser scanner Hds 6000 (Leica) caused by different object materials. Institute of Photogramm. Remote Sensing (IPF). Germany: Universitaet Karlsruhe (TH).
Voegtle, A., & Schwab A. (2008). Influences of different materials on the measurements of a terrestrial laser scanner (TLS). Proc. of the XXI Congress, the International Society for Photogrammetry and Remote Sensing, ISPR.
Wujanz, D., Burger, M., & Mettenleiter, M. (2017). An intensity-based stochastic model for terrestrial laser scanners. ISPRS Journal of Photogrammetry and Remote Sensing, 125, 146-155.
Xu, Y., Boerner, R., & Yao, W. (2017). Automated coarse registration of point clouds in 3D urban scenes using voxel based plane constraint. ISPRS Annals of Photogrammetry, Remote Sensing & Spatial Information Sciences, IV-2/W4, 185-191.
Yang, B., Dong, Z., & Liang, F. (2016). Automatic registration of large-scale urban scene point clouds based on semantic feature points. ISPRS Journal of Photogrammetry and Remote Sensing, 113, 43-58.
Yoshimura, R., Date, H., & Kanai, S. (2016). Automatic registration of MLS point clouds and SfM meshes of urban area. Geo-spatial Information Science, 19(3), 171-181
Zhou, Y., Cui, M., & Yang, L. (2009). Application of 3D laser scanner in topographic change monitor and analysis (pp. 4-382-4-385). In 9th International Conference on Electronic Measurement Instruments. ICEMI.
1 Este grupo de investigación es dirigido por William Barragán Zaque, desde 2007 desarrolla sus actividades en la Universidad Distrital Francisco José de Caldas en Bogotá. El grupo cuenta con el apoyo del proyecto curricular de Ingeniería topográfica y la maestría en Ciencias de la información y las comunicaciones con énfasis en Geomática. Sus áreas de trabajo se han concentrado en los temas desarrollados en los laboratorios de cartografía, fotogrametría y topografía de la misma universidad.
