Arquivos LAS de georreferenciamento em coordenadas geográficas (lat / lon) ou cartesianas (projetadas)?

Primeiro de tudo, eu sou um novato em GIS.

Eu tenho uma nuvem de pontos lidar com dados de GPS e telemetria associados. Em seguida, transformo as coordenadas de todos os pontos em um sistema de coordenadas ECEF e, finalmente, transformo as coordenadas angulares em um sistema de coordenadas cartesianas, usando a aproximação do plano tangente local (LTP). (Essa aproximação usa as coordenadas gps do primeiro ponto como uma origem do sistema de referência e, em seguida, referencia o restante dos pontos em relação ao primeiro).

Estou escrevendo o arquivo LAS usando liblas e realmente não sei como definir as informações de georreferência no arquivo LAS. Estou escrevendo as coordenadas de cada ponto na aproximação LTP, então eu tenho 2 perguntas:

1. É correto escrever as coordenadas dos pontos em relação à aproximação LTP ou devo escrever coordenadas em unidades lat / lon (gps) ?. Se isso estiver correto, como posso definir o georreferenciamento para respeitar o LTP?

2. Se '1' não estiver correto e precisar escrever as coordenadas em unidades lat / lon, como posso definir o georreferenciamento? Escrever um VLR com o valor WGS84 é suficiente?

Depende da versão da especificação LAS que você está usando. Se for 1,3 ou menos, as especificações definirão informações de georreferenciamento usando registros de comprimento variável (VLRs) predefinidos (consulte as especificações) usando o mesmo formato que o GeoTIFF:

A georreferenciação para o formato LAS usará o mesmo mecanismo robusto desenvolvido para o padrão GeoTIFF.

Esse formato, apesar de um pouco desafiador para se destacar às vezes, é notavelmente flexível. Ele se baseia em três tags definidas, denominadas GeoKeyDirectoryTag, que são como um sumário de dados georef, a tag GeoDoubleParamsTag, que é como um armazenamento de todos os valores de precisão dupla referidos na GeoKeyDirectoryTag e a tag GeoAsciiParamsTag, da mesma forma é usado para armazenar todos os valores ASCII (texto). Este site fornece uma boa explicação e um exemplo.

No entanto, a partir do LAS v. 1.4 , esse método de armazenamento de informações de georreferenciamento foi alterado para favorecer o formato de texto conhecido (WKT), também armazenado em VLRs definidas, embora o formato GeoTIFF ainda seja usado para legado:

As informações do Sistema de referência de coordenadas (CRS) para os dados do ponto são necessárias para todos os dados. As informações do CRS serão colocadas em Registros de comprimento variável ou Registros de comprimento variável estendidos (observe que, se o gravador desejar manter a compatibilidade herdada, o GeoTIFF nos VLRs deverá ser usado). O CRS é representado pelo GeoTIFF ou pelo Well Know Text, conforme indicado pelo bit WKT Global Encoding. Os formatos de registro de ponto 0-5 podem usar GeoTIFF ou WKT (mas não os dois simultaneamente). Os formatos de registro de ponto 6-10 devem usar o WKT.

Não vejo nenhuma razão, devido a esses formatos flexíveis, por que você não pode armazenar informações de pontos em coordenadas geográficas (lat / long), mas isso seria bastante incomum para os dados do LAS, pois nunca vi isso anteriormente. Imagino que a razão é que os conjuntos de dados LiDAR tendem a ser de grande escala (pequena extensão espacial) e, portanto, os sistemas de coordenadas projetados são os preferidos. Isso facilita muito o cálculo das distâncias entre pontos, o que é importante para alguns algoritmos (por exemplo, classificação ou filtragem de pontos).