Географическая привязка файлов LAS в географических координатах (широта / долгота) или в декартовых (прогнозируемых)?

Прежде всего, я новичок в ГИС.

У меня есть облако точек лидара с данными GPS и телеметрии. Затем я преобразую координаты каждой точки в систему координат ECEF и, наконец, я преобразую угловые координаты в декартову систему координат, используя приближение локальной касательной плоскости (LTP). (Это приближение использует координаты GPS первой точки в качестве источника системы отсчета, а затем ссылается на остальные точки относительно первой).

Я пишу файл LAS с использованием liblas, и я не знаю, как установить информацию о географической привязке для файла LAS. Я пишу координаты каждой точки в приближении LTP, поэтому у меня есть 2 вопроса:

1. Правильно ли записывать координаты точек относительно приближения LTP или я должен записывать координаты в единицах широта / долгота (gps) ?. Если это правильно, как я могу установить привязку к LTP?

2. Если '1' неверно, и я должен записать координаты в широтах / долях, как я могу установить географическую привязку? Достаточно ли писать VLR со значением WGS84?

Это зависит от того, какую версию спецификации LAS вы используете. Если он равен 1,3 или меньше, то спецификации определяют информацию о географической привязке, используя предварительно определенные (см. Спецификации) записи переменной длины (VLR), используя тот же формат, что и GeoTIFF:

При географической привязке для формата LAS будет использоваться тот же надежный механизм, который был разработан для стандарта GeoTIFF.

Этот формат, хотя порой и трудный для восприятия, довольно удивителен. Он опирается на три определенных тега, называемых тегом GeoKeyDirectoryTag, который походит на оглавление для данных georef, тег GeoDoubleParamsTag, который походит на хранилище всех значений двойной точности, указанных в теге GeoKeyDirectoryTag, и тег GeoAsciiParamsTag, которые аналогично используется для хранения всех ASCII (текстовых) значений. Этот сайт дает хорошее объяснение и пример.

Начиная с LAS v. 1.4, однако, этот метод хранения информации о географической привязке был изменен для поддержки формата общеизвестного текста (WKT), также хранящегося в определенных VLR, хотя формат GeoTIFF все еще используется для прежних версий:

Информация о системе координат (CRS) для точечных данных требуется для всех данных. Информация CRS будет помещена в записи переменной длины или записи расширенной длины переменной (обратите внимание, что если автор желает поддерживать унаследованную совместимость, необходимо использовать GeoTIFF в VLR). CRS представлен либо GeoTIFF, либо хорошо известным текстом, как указано битом глобального кодирования WKT. Форматы записи точек 0-5 могут использовать GeoTIFF или WKT (но не оба одновременно). Форматы записи точек 6-10 должны использовать WKT.

Я не вижу причины, учитывая эти гибкие форматы, почему вы не можете хранить информацию о точке в географических координатах (широта / долгота), но это было бы довольно необычно для данных LAS, поскольку я никогда раньше не видел, чтобы это было сделано. Я предполагаю, что причина в том, что наборы данных LiDAR, как правило, имеют довольно большой масштаб (небольшой пространственный экстент), и поэтому предпочтительными являются проекционные системы координат. Это значительно упрощает расчет расстояний между точками, что важно для некоторых алгоритмов (например, классификация точек или фильтрация).