Pliki LAS z georeferencją we współrzędnych geograficznych (lat / lon) czy kartezjańskich (rzutowanych)?

Przede wszystkim jestem nowicjuszem w GIS.

Mam chmurę punktów lidar z danymi GPS i telemetrii. Następnie przekształcam współrzędne każdego punktu na układ współrzędnych ECEF, a na koniec przekształcam współrzędne kątowe na kartezjański układ współrzędnych, stosując przybliżenie Lokalnej płaszczyzny stycznej (LTP). (To przybliżenie wykorzystuje współrzędne GPS pierwszego punktu jako początek układu odniesienia, a następnie odnosi się do reszty punktów w odniesieniu do pierwszego).

Piszę plik LAS za pomocą liblas i tak naprawdę nie wiem, jak ustawić informacje georeferencyjne w pliku LAS. Piszę współrzędne każdego punktu w przybliżeniu LTP, więc mam 2 pytania:

1. Czy poprawne jest zapisywanie współrzędnych punktów w odniesieniu do aproksymacji LTP, czy powinienem zapisywać współrzędne w jednostkach lat / lon (gps) ?. Jeśli jest to poprawne, jak mogę ustawić szacunek georeferencji dla LTP ?.

2. Jeśli „1” jest niepoprawne i muszę wpisać współrzędne w jednostkach lat / lon, jak mogę ustawić georeferencję? Czy wystarczy napisanie VLR o wartości WGS84?

Zależy to od używanej wersji specyfikacji LAS. Jeśli jest to wartość 1,3 lub mniejsza, to specyfikacje definiują informacje georeferencyjne przy użyciu predefiniowanych (patrz specyfikacje) rekordów o zmiennej długości (VLR) przy użyciu tego samego formatu co GeoTIFF:

Georeferencje dla formatu LAS będą korzystać z tego samego solidnego mechanizmu, który został opracowany dla standardu GeoTIFF.

Ten format, choć czasem trudny do opanowania, jest niezwykle elastyczny. Opiera się na trzech zdefiniowanych znacznikach o nazwie GeoKeyDirectoryTag, które są jak spis treści dla danych georef, GeoDoubleParamsTag, który jest jak magazyn wszystkich wartości podwójnej precyzji, o których mowa w GeoKeyDirectoryTag, i GeoAsciiParamsTag, które podobnie służy do przechowywania wszystkich wartości ASCII (tekstowych). Ta strona zawiera dobre wyjaśnienie i przykład.

Jednak od wersji LAS v. 1.4 ta metoda przechowywania informacji georeferencyjnych została zmieniona, aby faworyzować dobrze znany format tekstowy (WKT), również przechowywany w zdefiniowanych VLR, chociaż format GeoTIFF jest nadal używany dla starszych wersji:

Informacje o układzie współrzędnych (CRS) dla danych punktów są wymagane dla wszystkich danych. Informacje CRS zostaną umieszczone w rekordach o zmiennej długości lub w rekordach o zmiennej długości (pamiętaj, że jeśli autor chce zachować starszą kompatybilność, należy użyć GeoTIFF w VLR). CRS jest reprezentowany przez GeoTIFF lub Well Know Text, jak wskazuje bit WKT Global Encoding. Formaty rekordów punktowych 0-5 mogą używać GeoTIFF lub WKT (ale nie oba jednocześnie). Formaty rekordów punktowych 6-10 muszą używać WKT.

Nie widzę powodu, biorąc pod uwagę te elastyczne formaty, dlaczego nie można przechowywać informacji o punktach we współrzędnych geograficznych (szerokość / długość), ale byłoby to dość niezwykłe w przypadku danych LAS, ponieważ nigdy wcześniej tego nie widziałem. Wyobrażam sobie, że powodem jest to, że zestawy danych LiDAR mają raczej dużą skalę (mały zasięg przestrzenny), dlatego też preferowane są rzutowane układy współrzędnych. Ułatwia to obliczanie odległości między punktami, co jest ważne dla niektórych algorytmów (np. Klasyfikacja lub filtrowanie punktów).