Naar de hoofdinhoud

Voeg context toe aan uw projecten met vrij toegankelijke IGN-puntenwolken (HD LiDAR)

Geschreven door Sjaak Velthoven

Dit artikel is volledig automatisch vertaald.

 

Als onderdeel van het nationale LiDAR HD-programma produceert en distribueert IGN 3D-kaarten van de gehele grond en oppervlakte van Frankrijk in LiDAR-gegevens.

De verspreide gegevens zijn met name gekalibreerde puntenwolken, onverwerkt of geclassificeerd, en 3D-digitale modellen (DEM, DSM, MNH, enz.).

Het is mogelijk met gratis en open-source tools om IGN-puntenwolken (gepubliceerd in .laz- en Lambert 93-formaten) om te zetten in .las- of .e57-formaat (leesbaar in onze 3D-viewer) en in het referentiecoördinatensysteem (RCS) dat overeenkomt met uw project.

 

 

Volg deze stappen om een contextuele puntenwolk aan uw project toe te voegen!

 

De volgende thema's worden in dit artikel beschreven:

 

 

1. Upload uw puntenwolk

Dit is de gemakkelijkste stap!

Download de puntenwolktegel die overeenkomt met de context van uw project via deze interface.

Voor Zwitserland biedt Swisstopo ook puntenwolken aan die het hele grondgebied in .las-formaat bestrijken.

 

 

2. Download de open-source software die nodig is om uw puntenwolk om te zetten

Dit wordt lastig!

Om uw puntenwolk om te zetten, moet u deze twee open-source softwareprogramma's installeren:

  1. Miniconda is een lichtgewicht versie van de gratis en open-source Anaconda distributie van de Python- en R-programmeertalen, toegepast op de ontwikkeling van data science-toepassingen.

  2. PDAL is een open-source bibliotheek voor verwerking van puntenwolkgegevens. Het is een beetje zoals een VLC-speler voor puntenwolken ;)

 

Eigenlijk hebben we de Anaconda Prompt nodig om PDAL te gebruiken, wat zorg draagt voor de transformatie.

Laten we gaan 👇

 

De volgende thema's worden in dit gedeelte beschreven:

 

 

2.1 Miniconda installeren

PDAL is afhankelijk van andere bibliotheken, dus de gemakkelijkste manier is om Miniconda (een lichtgewicht versie van Anaconda) te gebruiken.

  1. Ga naar de officiële Miniconda-downloadpagina:
    👉 https://docs.conda.io/en/latest/miniconda.html

  2. Download de Miniconda-installatieprogramma voor Windows (64-bits, Python 3.x).

  3. Start het installatieprogramma:
    - Accepteer de licentieovereenkomst
    - Kies "Alleen voor mij" (aanbevolen)
    - Houd de standaard installatielocatie aan
    - Controleer "Miniconda3 aan PATH toevoegen" als de optie wordt aangeboden

  4. Nadat u het hebt geïnstalleerd, opent u de Anaconda Prompt (dit is het venster dat u in plaats van CMD/PowerShell voor PDAL zult gebruiken).

 

2.2 Maak een Conda-omgeving voor PDAL

Het wordt aanbevolen om PDAL in zijn eigen omgeving te isoleren.

In Anaconda Prompt kopieert en voert u de volgende opdracht uit:

 

conda create -n pdal-env -c conda-forge pdal python=3.10

 

Dit doet drie dingen:

  • Maak een omgeving met de naam pdal-env

  • Installeer PDAL vanuit de conda-forge-opslagplaats (meest recente versie)

  • Python installeren (nuttig als u PDAL ook in uw Python-scripts wilt gebruiken)

 

2.3 Activeer de PDAL-omgeving

Nog steeds in de Anaconda Prompt, activeer de omgeving alvorens PDAL te gebruiken door het volgende uit te voeren:

 

conda activate pdal-env

 

Nu zou elke pdal-opdracht deze omgeving moeten gebruiken.

(Als u de terminal sluit, voert u conda activate pdal-env opnieuw uit wanneer u deze opnieuw opent.)

 

2.4 Controleer de installatie

Voer uit:

 

pdal --version pdal --drivers

 

Als een versie wordt weergegeven en een lijst met "readers.* / filters.* / writers.*" verschijnt, is PDAL geïnstalleerd. ✅ Controleer of writers.las in de lijst staat als u naar .las wilt exporteren.

Hetzelfde voor het .e57-formaat dat ook door Catenda Hub wordt ondersteund ;)

 

 

3. Transform uw puntenwolk

We zijn er bijna!

Nog maar een paar stappen...

 

De volgende thema's worden in dit gedeelte beschreven:

 

 

3.1 Identificeer de EPSG-code van uw coördinaatreferentiesysteem (CRS)

Om uw puntenwolk om te zetten in PDAL, moet u de EPSG-code kennen die overeenkomt met het coördinatensysteem van uw IFC-model.

Hier is een niet-uitputtende lijst van SCR's die het metropolitaanse Frankrijk en hun respectieve EPSG-codes bestrijken.

 

Naam

EPSG

Opmerking

RGF93 / Lambert-93

2154

De standaard voor metropolitaan Frankrijk, gebruikt voor IGN-puntenwolken. Dekt het hele grondgebied.

RGF93 / Lambert CC42

3942

Voor zone 1 (Corsica).

RGF93 / Lambert CC43

3943

Voor zone 2 (zuiden bij 44° N).

RGF93 / Lambert CC44

3944

Voor zone 3 (43° N tot 45° N).

RGF93 / Lambert CC45

3945

Voor zone 4 (44° N tot 46° N).

RGF93 / Lambert CC46

3946

Voor zone 5 (45° N tot 47° N).

RGF93 / Lambert CC47

3947

Voor zone 6 (46° N tot 48° N).

RGF93 / Lambert CC48

3948

Voor zone 7 (47° N tot 49° N).

RGF93 / Lambert CC49

3949

Voor zone 8 (48° N tot 50° N).

RGF93 / Lambert CC50

3950

Voor zone 9 (noord bij 49° N).

 

 

3.2 Transform uw puntenwolk met Anaconda en PDAL

Nu we de beroemde EPSG-code kennen die overeenkomt met de SCR van ons project, kunnen we EINDELIJK onze puntenwolk transformeren

 

Hier is een voorbeeld van een opdracht waarmee u uw puntenwolk kunt transformeren:

 

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"

 

🤓 Laten we deze opdracht analyseren:

  1. pdal translate is de hoofdopdracht voor het uitvoeren van de transformatie.

  2. "C:\Users\USERNAME\Downloads\LHD_FXX_0766_6282_PTS_LAMB93_IGN69.copc.laz" is het pad naar het .laz-bestand dat uit de IGN-database is gedownload.

  3. "C:\Users\USERNAME\Downloads\LHD_FXX_0766_6282_PTS_LAMB93_IGN69.copc.laz" is het pad waar het toekomstige .las-bestand wordt gemaakt.

  4. de reprojection command --filters.reprojection.in_srs="EPSG:2154" --filters.reprojection.out_srs="EPSG:3943" stelt u in staat om de puntenwolk van het Lambert-93 SCR (EPSG:2154) naar het CC43 SCR (3943) om te projecteren

 

Kopieer de opdracht en vervang eenvoudig de EPSG-paden en codes in invoer en uitvoer.

Zodra de opdracht is uitgevoerd, wordt uw nieuwe .las-bestand gegenereerd op de locatie die is opgegeven in het doelpad.

 

IGN-puntenwolken zijn standaard niet gekleurd, maar er zijn opdrachten in PDAL om kleuren toe te voegen op basis van puntverhoging of door gebruik te maken van een geogerefereerde .tif-orthofoto met dezelfde CRS.

De open-source software QGIS maakt het gemakkelijk om geogerefereerde orthofoto's te genereren.

 

Hier is dezelfde opdracht met een colorization filter die verwijst naar een orthofoto:

 

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" ^ colorization ^ --filters.colorization.raster="C:\Users\USERNAME\Downloads\Orthophoto.tif" ^ --filters.colorization.dimensions="Red:1,Green:2,Blue:3"

 

3.3 Upload uw puntenwolk naar uw Catenda Hub-project

Zodra de transformatie is voltooid (naar de CRS van uw project en in .las- of .e57-formaat), hoeft u alleen maar uw puntenwolk naar uw Catenda Hub-project te uploaden!

Als u alle stappen correct hebt gevolgd, moet deze perfect aansluiten bij uw project :)

 

 

Was dit een antwoord op uw vraag?