LiDAR-HD Downloader

Description

Utility to download France IGN LiDAR HD to LAZ file based on AOI (Area of Interest) defined in a GeoPandas GeoDataframe.

Installation

Create a new Conda environment and install PDAL (see PDAL documentation) :

conda update -n base -c conda-forge conda
conda create --yes --name pdal_env --channel conda-forge pdal
conda activate pdal_env
pdal --version # Check PDAL version

Clone this repository and install dependencies :

git clone https://github.com/thiboyyy/lidarhd.git
cd lidarhd && pip install -r requirements.txt

Usage

CLI example

python lidarhd.py --lat 47.34 --lon 2.85 --radius 100 --name "Example Area Name" --output ./lidarhd_downloads/

Python example

from lhd import LiDARHD
from shapely.geometry import LineString
import geopandas as gpd

# Download LiDARHD tiling database
lidarhd = LiDARHD(folder_path='./lidarhd/', overwrite=False)

# Create a GeoDataframe as input (here a simple segment in Paris)
linestring = LineString([(2.3522, 48.8566), (2.3639, 48.8667)])
gdf = gpd.GeoDataFrame({
    'name': ['Paris Hotel de Ville - Paris Place de la République'],
    'geometry': [linestring]}, crs=4326).to_crs(2154)
gdf['geometry'] = gdf.geometry.buffer(30)

# Download LiDARHD data covering the GeoDataframe area using PDAL and COPC polygon reading capabilities
lidarhd.download(gdf, "./lidarhd/points.laz")