A collection of easy to use command line tools for processing of point cloud data:
- basic management of data (info, translate, merge, tile, thin, clip, boundary, density, to_vector)
- export of raster grids (to_raster, to_raster_tin)
- handle virtual point clouds (build_vpc)
And more algorithms will be added in the future!
The whole suite of tools has been integrated into the QGIS Processing framework in QGIS starting from QGIS 3.32 (release in June 2023). There is no plugin to install - everything is available in QGIS core - just open the Processing toolbox and search for point cloud algorithms.
Most of the tools are multi-threaded, making good use of all available CPUs for fast processing.
All tools are based on PDAL pipelines, also using the other usual geospatial libraries: GDAL/OGR, GEOS and PROJ.
Some may ask why not let users build their pipelines in PDAL directly. There are several reasons:
- ease of use: PDAL is great for advanced users, but not everyone finds it easy to manually craft JSON files with pipelines, study manuals of the many stages and read details about file formats involved
- parallel execution: PDAL runs pipelines in a single thread - only one CPU gets to do the work normally - and users need to implement their own parallelism if they wish so
You will need: PDAL >= 2.5 and GDAL >= 3.0, both with development files.
Building is done with CMake:
mkdir build
cd build
cmake ..
make
If PDAL is not installed in the usual system paths (e.g. in /usr), then add PDAL_DIR parameter when running CMake - for example if your PDAL installation dir is /home/martin/pdal-inst:
cmake -DPDAL_DIR=/home/martin/pdal-inst/lib/cmake/PDAL ..
PDAL runs point cloud pipelines in a single thread and any parallelization is up to users of the library. PDAL wrench has parallel processing built in and tries to run pipelines in parallel. This generally happens in two scenarios:
- input dataset is in COPC or EPT format - these formats allow efficient spatial queries (to access only data in a particular bounding box), so PDAL wrench can split the input into multiple tiles that can be processed independently in parallel
- input dataset is a virtual point cloud (VPC) - such datasets are composed of a number of files, so the whole work can be split into jobs where each parallel job processes one or more input files
If the input is a single LAS/LAZ file, no parallelization is attempted. This may change in the future with introduction of more complex algorithms (where the cost of reading the input is much lower than the cost of the actual algorithm).
Prints basic metadata from the point cloud file:
pdal_wrench info --input=data.las
Convert to a different file format (e.g. create compressed LAZ):
pdal_wrench translate --input=data.las --output=data.laz
Reproject point cloud to a different coordinate reference system:
pdal_wrench translate --input=data.las --output=reprojected.las --transform-crs=EPSG:3857
Assign coordinate reference system (if not present or wrong):
pdal_wrench translate --input=data-with-invalid-crs.las --output=data.las --assign-crs=EPSG:3857
Exports a polygon file containing boundary. It may contain holes and it may be a multi-part polygon.
pdal_wrench boundary --input=data.las --output=boundary.gpkg
Exports a raster file where each cell contains number of points that are in that cell's area.
pdal_wrench density --input=data.las --resolution=1 --output=density.tif
Outputs only points that are inside of the clipping polygons.
pdal_wrench clip --input=data.las --polygon=clip.gpkg --output=data_clipped.las
Merges multiple point cloud files to a single one.
pdal_wrench merge --output=merged.las data1.las data2.las data3.las
Creates tiles from input data. For example to get tiles sized 100x100:
pdal_wrench tile --length=100 --output=/data/tiles data1.las data2.las data3.las
Creates a thinned version of the point cloud by only keeping every N-th point (every-nth mode) or keep points based on their distance (sample mode).
For example, to only keep every 20th point, so only 5% of points will be in the output:
pdal_wrench thin --output=thinned.las --mode=every-nth --step-every-nth=20 data.las
Alternatively, to sample points using Poisson sampling of the input:
pdal_wrench thin --output=thinned.las --mode=sample --step-sample=20 data.las
Exports point cloud data to a 2D raster grid, having cell size of given resolution, writing values from the specified attribute. Uses inverse distance weighting.
pdal_wrench to_raster --output=raster.tif --resolution=1 --attribute=Z data.las
Exports point cloud data to a 2D raster grid like to_raster does, but using a triangulation of points and then interpolating cell values from triangles. It does not produce any "holes" when some data are missing. Only supports output of Z attribute.
pdal_wrench to_raster_tin --output=raster.tif --resolution=1 data.las
Exports point cloud data to a vector layer with 3D points (a GeoPackage), optionally with extra attributes:
pdal_wrench to_vector --output=data.gpkg data.las
This is similar to GDAL's VRT - a single file referring to other files that contain actual data. Software then may handle all data as a single dataset.
Virtual point clouds based on STAC protocol's ItemCollection, which is in fact a GeoJSON feature collection with extra metadata written in a standard way. See VPC spec for more details on the file format.
To create a virtual point cloud:
pdal_wrench build_vpc --output=hello.vpc data1.las data2.las data3.las
Or, if the inputs are listed in a text file:
data1.las
data2.las
data3.las
You can provide that file as input:
pdal_wrench build_vpc --output=hello.vpc --input-file-list=inputs.txt
Afterwards, other algorithms can be applied to a VPC:
pdal_wrench clip --input=hello.vpc --polygon=clip.gpkg --output=hello_clipped.vpc
This will create a grid for each data file separately in parallel and then merge the results to a final GeoTIFF:
pdal_wrench density --input=hello.vpc --resolution=1 --output=density.tif
When algorithms create derived VPCs, by default they use uncompressed LAS, but --output-format=laz option can switch to compressed LAZ.
| Algorithm | VPC | Notes |
|---|---|---|
| info | yes | |
| boundary | multi-threaded | per file |
| density | multi-threaded | spatial tiling |
| clip | multi-threaded | per file |
| merge | not supported | |
| tile | multi-threaded | per file |
| thin | multi-threaded | per file |
| to_raster | multi-threaded | spatial tiling |
| to_raster_tin | multi-threaded | spatial tiling |
| to_vector | multi-threaded | per file |
| translate | multi-threaded | per file |