Process to convert laz file to DTM

100%

Process to convert laz file to DTM

✓ Published0🌍 Public

Last edited Feb 13, 2022

Created on Feb 13, 2022

This example demonstrates a command-line workflow for converting a LIDAR point cloud in LAZ format into a digital terrain model (DTM) as a GeoTIFF raster. It shows two approaches using the PDAL library, both relying on its `pipeline` command and `translate` subcommand, paired with the `writers.gdal` writer to output a 0.5-meter resolution TIFF. The process ingests a sample LAZ file from the French IGN LIDAR HD dataset, then uses PDAL’s GDAL writer to generate the DTM, with a failed attempt at using the `writers.ogr` output to shapefile also noted.

AI-generated description

Similar vizzes

Loading thumbnail…

Gist ab738ee0d721842fbea501cc96296160

This example demonstrates how to load and manage LiDAR point cloud layers in QGIS using Python, showing both automatic and manual index generation for LAZ files from the French IGN LIDAR HD dataset. It uses the QgsPointCloudLayer class with the PDAL provider, QgsProject for layer management, and QgsPointCloudDataProvider to check indexing state, generate indexes, and retrieve polygon bounds, point counts, and renderer information. The script also prints coordinate reference system details and renderer type, illustrating direct provider-level access to point cloud metadata.

TThomasG77

63% match

Loading thumbnail…

Gist 0bdc45cce93d15fce0d1991d2b4b9e9d

This example demonstrates how to convert a compressed LiDAR point-cloud dataset from IGN’s French national archive into a ZIP archive for easier distribution. The shell script downloads a 7z file via `wget`, decompresses it with `unp`, recompresses the contents into a ZIP using `zip`, and then reports the final disk usage with `du -cksh`. It relies on standard command-line utilities rather than JavaScript libraries, showing a practical data-preparation workflow for large geospatial files.

TThomasG77

62% match

Loading thumbnail…

Load GTFS as GPKG with GDAL 3.7

This example demonstrates converting an IDFM GTFS transit feed into a GeoPackage database using GDAL’s new GTFS driver. The Bash script first checks that GDAL 3.7 is installed, then uses `ogrinfo` with `-json` to list all available layers. For each layer, it runs `ogr2ogr` with the GeoPackage format, appending layers incrementally with the `-update` flag. The source data is loaded directly from a URL via `/vsicurl/`, and the script allows customizing input and output names through environment variables.

TThomasG77

58% match

Loading thumbnail…

Gist c32a8597373bd00c10029d4f21f2400e

This example demonstrates how to convert a GDAL VRT (Virtual Dataset) into a single raster file using the `ogr2vrt` utility. It shows a command-line workflow where a shell script first installs the `gdal-utils` package via `pip`, then locates the installed samples directory to make the utility scripts executable, and finally runs `ogr2vrt.py` to perform the conversion. The code relies on Python’s `osgeo_utils` module and standard shell commands to manage the GDAL sample tools.

TThomasG77

57% match

Loading thumbnail…

Generate vector tiles for Planet with tilemaker

This example demonstrates a command-line workflow for generating vector tiles from the entire OpenStreetMap planet dataset using tilemaker. It shows how to automate the process of downloading planet-latest.osm.pbf, compiling tilemaker from source, and producing a planet.mbtiles output using the OpenMapTiles configuration and Lua processing scripts. The script also periodically checks the dataset’s last-modified header via curl to trigger regeneration only when new updates exist. External resources like Natural Earth landcover and coastline polygons are fetched and merged into the tile pipeline.

TThomasG77

57% match

Loading thumbnail…

mapshaper-info-recreated-via-gdal

This example recreates the command-line tool mapshaper-info using the gdal-async library to inspect a GeoJSON TopoJSON file containing German administrative regions. It downloads the dataset from GitHub, then prints layer metadata including geometry type, record count, bounding box, and coordinate reference system, followed by a formatted table of attribute names and their first values. The script uses the got library for streaming downloads and the pipeline function from stream/promises to handle the file writing asynchronously.

TThomasG77

55% match

Loading thumbnail…

Export DVF geolocalisées en Geoparquet

This example converts French land-registry transaction data (DVF) from compressed CSV files into a single spatially-indexed GeoParquet file. It downloads yearly datasets from data.gouv.fr, loads them into a DuckDB database, and creates a point geometry for each transaction before exporting with a Hilbert-curve spatial ordering. The script leverages DuckDB’s spatial extension, the `ST_Point` and `ST_Hilbert` functions, and parquet compression via zstd.

TThomasG77

54% match

Loading thumbnail…

GDAL VSICURL tests to remote BDTOPO Geopackage France

This example tests GDAL’s VSICURL virtual file system by streaming a remote GeoPackage of French administrative data, specifically focusing on mairie (town hall) geometries. It demonstrates two ogr2ogr commands that query the BDTOPO dataset using SQL joins to extract mairie footprints as multipolygons and their centroids as points, converting them to GeoJSON with EPSG:4326. The visualization shows how both polygon boundaries and point locations are derived from the same remote source, highlighting the ability to perform spatial SQL operations directly on cloud-hosted vector data without local downloads.

TThomasG77

53% match