Dans le cadre du programme national LiDAR HD, l’IGN produit et diffuse une cartographie 3D de l’intégralité du sol et du sursol de la France en données LiDAR. Les données diffusées sont en particulier des nuages de points recalés, bruts ou classifiés, et des modélisations numériques 3D (MNT, MNS, MNH…).
Il est possible avec des outils gratuits et open source de transformer les nuages de points de l'IGN (publiés en format .laz et en Lambert 93) en format .las ou .e57 (lisibles dans notre visionneuse 3D) et vers le Système de Coordonnées de Référence (SCR) correspondant à votre projet.
Suivez les étapes suivantes pour ajouter un nuage de points contextuel à votre projet !
1. Téléchargez votre nuage de points
C'est l'étape la plus simple ! Téléchargez via cette interface la dalle du nuage de points correspondant au contexte votre projet.
2. Téléchargez les logiciels open source requis pour la transformation de votre nuage de points
C'est ici que ça se corse ! Pour transformer votre nuage de points, vous aurez besoin d'installer ces deux logiciels open source :
PDAL, une librairie open source pour le traitement des données de nuages de points. C'est un peu comme un lecteur VLC pour les nuages de points ;)
En fait, on a besoin de l'Anaconda Prompt pour utiliser PDAL, qui s'occupera de la transformation. Allez, c'est parti 👇
2.1. Installez Miniconda
PDAL dépend d’autres bibliothèques, donc la manière la plus simple est d’utiliser Miniconda (une version légère d’Anaconda).
Allez sur la page officielle de téléchargement de Miniconda :
👉 https://docs.conda.io/en/latest/miniconda.htmlTéléchargez l’installateur Miniconda pour Windows (64 bits, Python 3.x).
Lancez l’installateur :
Acceptez la licence
Choisir « Just for me » (recommandé)
Gardez l’emplacement d’installation par défaut
Cochez « Add Miniconda3 to PATH » si l’option est proposée
Une fois installé, ouvrir l’Anaconda Prompt (c’est cette fenêtre que vous utiliserez au lieu de CMD/PowerShell pour PDAL).
2.2. Créez un environnement Conda pour PDAL
Il est recommandé d’isoler PDAL dans son propre environnement.
Dans l’Anaconda Prompt, copiez et exécutez la commande suivante :
conda create -n pdal-env -c conda-forge pdal python=3.10
Cela fait trois choses :
Création d'un environnement appelé
pdal-envInstallation de PDAL depuis le dépôt conda-forge (version la plus récente)
Installation de Python (utile si vous voulez aussi utiliser PDAL dans vos scripts Python)
2.3. Activez l’environnement PDAL
Toujours dans l'Anaconda Prompt, avant d’utiliser PDAL, activez l’environnement en exécutant :
conda activate pdal-env
Maintenant, n'importe quelle commande pdal devrait utiliser cet environment. (Si vous fermez le terminal, tapez conda activate pdal-env une nouvelle fois lors de la réouverture de ce dernier.)
2.4. Vérifier l’installation
Exécutez :
pdal --version pdal --drivers
Si une version s’affiche et qu’une liste de “readers.* / filters.* / writers.*” apparaît, PDAL est bien installé. ✅
Vérifiez que writers.las est dans la liste si vous voulez exporter en .las. Pareil pour le format .e57 également supporté par Catenda Hub ;)
3. Transformez votre nuage de points
On y est presque ! Encore quelques étapes...
3.1. Identifiez le code EPSG de votre Système de Coordonnées de Référence (SCR)
Pour transformer votre nuage de points dans PDAL, vous devez connaître le code EPSG correspondant à votre système de coordonnées. Voici une liste non exhaustive des SCR couvrant la France métropolitaine ainsi que leur code EPSG respectif.
Nom | EPSG | Remarques |
RGF93 / Lambert-93 | 2154 | Le standard pour la France métropolitaine, utilisé pour les nuages de points de l'IGN. Couvre l'ensemble du territoire. |
RGF93 / Lambert CC42 | 3942 | Pour la zone 1 (Corse). |
RGF93 / Lambert CC43 | 3943 | Pour la zone 2 (sud de 44° N). |
RGF93 / Lambert CC44 | 3944 | Pour la zone 3 (43° N à 45° N). |
RGF93 / Lambert CC45 | 3945 | Pour la zone 4 (44° N à 46° N). |
RGF93 / Lambert CC46 | 3946 | Pour la zone 5 (45° N à 47° N). |
RGF93 / Lambert CC47 | 3947 | Pour la zone 6 (46° N à 48° N). |
RGF93 / Lambert CC48 | 3948 | Pour la zone 7 (47° N à 49° N). |
RGF93 / Lambert CC49 | 3949 | Pour la zone 8 (48° N à 50° N). |
RGF93 / Lambert CC50 | 3950 | Pour la zone 9 (nord de 49° N). |
3.2. Transformez votre nuage de points en utilisant Anaconda et PDAL
Maintenant que l'on connait le fameux code EPSG correspondant au SCR de notre projet, on peut ENFIN transformer notre nuage de points !
Voici un exemple de commande qui permet de transformer votre nuage de points :
pdal translate "C:\Users\USERNAME\Downloads\LHD_FXX_0766_6282_PTS_LAMB93_IGN69.copc.laz" "C:\Users\USERNAME\Downloads\LHD_FXX_0766_6282_PTS_LAMB93_IGN69.copc.las" reprojection --filters.reprojection.in_srs="EPSG:2154" --filters.reprojection.out_srs="EPSG:3943"
🤓 Décomposons cette commande ensemble :
pdal translateest la commande principale pour effectuer la transformation."C:\Users\USERNAME\Downloads\LHD_FXX_0766_6282_PTS_LAMB93_IGN69.copc.laz"est le chemin du fichier.laztéléchargé depuis la base de données de l'IGN."C:\Users\USERNAME\Downloads\LHD_FXX_0766_6282_PTS_LAMB93_IGN69.copc.laz"est le chemin vers lequel le futur fichier.lassera créé.la commande
reprojection --filters.reprojection.in_srs="EPSG:2154" --filters.reprojection.out_srs="EPSG:3943"permet de faire la reprojection du nuage de points depuis le SCR Lambert-93 (EPSG:2154) vers le SCR CC43 (3943)
Copiez la commande et remplacez simplement les chemins et codes EPSG en input et output. Une fois la commande exécutée, votre nouveau fichier .las sera généré à l'emplacement spécifié par le chemin cible.
3.3. Téléversez votre nuage de points dans votre projet Catenda Hub
Une fois la transformation effectuée (au SCR de votre projet, et en format .las ou .e57), il ne vous reste plus qu'à téléverser votre nuage de points dans votre projet Catenda Hub ! Si vous avez bien suivi toutes les étapes, il devrait être parfaitement aligné avec votre projet :)
PS : par défaut, les nuages de point de l'IGN ne sont pas colorisés, mais il existe des commandes dans PDAL permettant d'ajouter des couleurs à partir de l'élévation des points ou en se servant d'une orthophoto .tif. Ce sera peut-être l'objet d'un article complémentaire ;)