Géoréférencer des fichiers LAS en coordonnées géographiques (lat / lon) ou cartésiennes (projetées)?

Tout d'abord, je suis un débutant dans le SIG.

J'ai un nuage de points lidar avec des données GPS et de télémétrie associées. Ensuite, je transforme les coordonnées de chaque point en un système de coordonnées ECEF et enfin je transforme les coordonnées angulaires en un système de coordonnées cartésiennes, en utilisant l'approximation du plan tangent local (LTP). (Cette approximation utilise les coordonnées GPS du premier point comme origine du système de référence, puis référence le reste des points par rapport au premier).

J'écris le fichier LAS à l'aide de liblas et je ne sais pas vraiment comment définir les informations de géoréférence dans le fichier LAS. J'écris les coordonnées de chaque point dans l'approximation LTP donc j'ai 2 questions:

1. Est-il correct d'écrire les coordonnées des points par rapport à l'approximation LTP ou dois-je écrire les coordonnées en unités lat / lon (gps)?. Si cela est correct, comment puis-je configurer le géoréférencement pour respecter le LTP?.

2. Si «1» n'est pas correct et que je dois écrire les coordonnées en unités lat / lon, comment puis-je définir le géoréférencement? L'écriture d'un VLR avec la valeur WGS84 est-elle suffisante?

Cela dépend de la version de la spécification LAS que vous utilisez. S'il est de 1,3 ou moins, les spécifications définissent les informations de géoréférencement à l'aide d'enregistrements de longueur variable (VLR) prédéfinis (voir spécifications) en utilisant le même format que GeoTIFF:

Le géoréférencement pour le format LAS utilisera le même mécanisme robuste qui a été développé pour la norme GeoTIFF.

Ce format, bien que quelque peu difficile à gérer parfois, est remarquablement flexible. Il s'appuie sur trois balises définies appelées balise GeoKeyDirectoryTag, qui est comme une table des matières pour les données géoréférencées, la balise GeoDoubleParamsTag, qui est comme un magasin de toutes les valeurs de double précision référencées dans GeoKeyDirectoryTag, et la balise GeoAsciiParamsTag, qui de même est utilisé pour stocker toutes les valeurs ASCII (texte). Ce site fournit une bonne explication et un exemple.

Cependant, à partir de LAS v. 1.4 , cette méthode de stockage des informations de géoréférencement a été modifiée pour favoriser le format de texte bien connu (WKT), également stocké dans des VLR définis, bien que le format GeoTIFF soit toujours utilisé pour l'héritage:

Les informations du système de référence de coordonnées (CRS) pour les données ponctuelles sont requises pour toutes les données. Les informations CRS seront placées dans des enregistrements de longueur variable ou des enregistrements de longueur variable étendue (notez que si le rédacteur souhaite conserver la compatibilité héritée, alors GeoTIFF dans les VLR doit être utilisé). Le CRS est représenté par GeoTIFF ou Well Know Text comme indiqué par le bit WKT Global Encoding. Les formats d'enregistrement de points 0-5 peuvent utiliser GeoTIFF ou WKT (mais pas les deux simultanément). Les formats d'enregistrement de points 6-10 doivent utiliser WKT.

Je ne vois aucune raison, étant donné ces formats flexibles, pourquoi vous ne pouviez pas stocker les informations de point dans les coordonnées géographiques (lat / long), mais ce serait assez inhabituel pour les données LAS dans la mesure où je ne l'ai jamais vu auparavant. J'imagine que la raison en est que les ensembles de données LiDAR ont tendance à être à assez grande échelle (petite étendue spatiale) et les systèmes de coordonnées projetés sont donc préférés. Il facilite le calcul des distances entre les points, ce qui est important pour certains algorithmes (par exemple, la classification ou le filtrage des points).