¿Georreferenciar archivos LAS en coordenadas geográficas (lat / lon) o cartesianas (proyectadas)?

En primer lugar, soy un novato en SIG.

Tengo una nube de puntos lidar con GPS y datos de telemetría asociados. Luego, transformo las coordenadas de cada punto en un sistema de coordenadas ECEF y finalmente transformo las coordenadas angulares en un sistema de coordenadas cartesianas, utilizando la aproximación del plano de tangente local (LTP). (Esta aproximación utiliza las coordenadas gps del primer punto como origen del sistema de referencia y luego, hace referencia al resto de puntos con respecto al primero).

Estoy escribiendo el archivo LAS usando liblas, y realmente no sé cómo configurar la información de georreferencia en el archivo LAS. Estoy escribiendo las coordenadas de cada punto en la aproximación LTP, así que tengo 2 preguntas:

1. ¿Es correcto escribir las coordenadas de los puntos con respecto a la aproximación LTP o debo escribir las coordenadas en unidades lat / lon (gps)? Si esto es correcto, ¿cómo puedo configurar la georreferenciación para respetar el LTP?

2. Si '1' no es correcto y debo escribir las coordenadas en unidades lat / lon, ¿cómo puedo configurar la georreferenciación? ¿Es suficiente escribir un VLR con valor WGS84?

Depende de qué versión de la especificación LAS esté utilizando. Si es 1.3 o menos, entonces las especificaciones definen la información de georreferenciación utilizando registros de longitud variable (VLR) predefinidos (ver especificaciones) usando el mismo formato que el GeoTIFF:

La georreferenciación para el formato LAS utilizará el mismo mecanismo robusto que se desarrolló para el estándar GeoTIFF.

Este formato, aunque a veces es difícil de asimilar, es notablemente flexible. Se basa en tres etiquetas definidas llamadas etiqueta GeoKeyDirectoryTag, que es como una tabla de contenido para datos georef, la etiqueta GeoDoubleParamsTag, que es como un almacén de todos los valores de doble precisión mencionados en GeoKeyDirectoryTag, y la etiqueta GeoAsciiParamsTag, que de manera similar se usa para almacenar todos los valores ASCII (texto). Este sitio proporciona una buena explicación y un ejemplo.

Sin embargo, a partir de LAS v. 1.4 , este método de almacenamiento de información de georreferenciación se cambió para favorecer el formato de texto conocido (WKT), también almacenado en VLR definidos, aunque el formato GeoTIFF todavía se utiliza para el legado:

La información del Sistema de referencia de coordenadas (CRS) para los datos de puntos es necesaria para todos los datos. La información de CRS se colocará en registros de longitud variable o registros de longitud variable extendida (tenga en cuenta que si el escritor desea mantener la compatibilidad heredada, entonces se debe usar GeoTIFF en VLR). El CRS está representado por GeoTIFF o Well Know Text como lo indica el bit de codificación global WKT. Los formatos de registro de puntos 0-5 pueden usar GeoTIFF o WKT (pero no ambos simultáneamente). Los formatos de registro de puntos 6-10 deben usar WKT.

No veo ninguna razón dados estos formatos flexibles por los que no puede almacenar información de puntos en coordenadas geográficas (lat / long), pero esto sería bastante inusual para los datos de LAS, ya que nunca los he visto hacer antes. Me imagino que la razón es que los conjuntos de datos LiDAR tienden a ser de gran escala (pequeña extensión espacial) y, por lo tanto, se prefieren los sistemas de coordenadas proyectados. Facilita el cálculo de las distancias entre puntos, lo cual es importante para algunos algoritmos (por ejemplo, clasificación o filtrado de puntos).