Library for common LAS files manipulation with PDAL
Project description
ign-pdal-tools
This repo contains various python tools based on PDAL that are used to work on LiDAR data in the LidarHD project at IGN (Institut National de l'Information Géographique et Forestière / French National Institute of Geographic and Forest Information) to work on LiDAR data.
We've decided to make them available because think that they may be useful to others, but this repo is NOT meant to be substantially modified from community input, and may be amended/completed depending on the fonctionalities that our team needs.
Content
This library contains pdal-based tools to:
- colorize a point cloud using images from Geoplateforme / cartes.gouv.fr (a portal from French government providing access to aerial imagery)
- stitch together LAS files using their location
- standardize LAS files
- unlock LAS files generated by TerraSolid
Installation / Usage
This library can be used in different ways:
- directly from sources:
make install
creates a mamba environment with the required dependencies - installed with
pip
from pypi:pip install ign-pdal-tools
- used in a docker container: see documentation Dockerfile
More details on the contained tools
Colorization
- color.py: Colorize a point cloud from Geoplateforme data
Las infos
Misc tools to get information on a las file, eg. retrieve metadata, find epsg value, find bounds, get parameters to pass to a writer. They are intended to be used from the pdaltools module, for example:
from pdaltools import las_infos
filename = ...
las_infos.las_info_metadata(filename)
Point cloud infos
Misc tools to get information on a point cloud (numpy array). Eg. get expected origin of a point cloud based on a square tiling:
from pdaltools import pcd_infos
points = ...
pcd_infos.get_pointcloud_origin_from_tile_width(points, tile_width=1000)
Stitching
- las_clip.py: crop a LAS file using 2d bounding box
- las_merge.py: merge a LAS file with its neighbors according to their filenames
- las_add_buffer.py: add points to a LAS file from a buffer (border) from its neighbors (using filenames to locate neighbors)
WARNING: In las_merge.py
and las_add_buffer.py
, filenames are used to get the LAS files extents
and to find neighbors.
The naming convention is {prefix1}_{prefix2}_{xcoord}_{ycoord}_{postfix})}
(eg. Semis_2021_0770_6278_LA93_IGN69.laz
).
By default, xcoord
and ycoord
are given in kilometers and the shape of the tile is 1 km * 1 km
Standardization
- standardize_format.py: re-write a LAS file in a standard format (see code for details)
- count_occurences: count occurences for each value of a given attribute in a set of LAS files (initially used for classification)
- count_occurences_for_attribute.py: count occurences in one or several files and save the result in a json file.
- merge_occurences_counts.py: merge counts from several results of count_occurences_for_attribute (json files) into a single json file (used for parallelization)
- replace_attribute_in_las.py: using a json file containing a correspondance map, replace the occurences of each value in a LAS file by its corresponding value from the map.
Unlock
unlock_file.py: overwrite a LAS file in case PDAL raises this error:
readers.las: Global encoding WKT flag not set for point format 6 - 10.
which is due to TerraSolid
malformed LAS output for LAS1.4 files with point format 6 to 10.
Add points in Las
add_points_in_las.py: add points from some vector files (ex: shp, geojson, ...) inside Las. New points will have X,Y and Z coordinates. Other attributes values given by the initial las file are null (ex: classification at 0). These others attributes could be forced by using the '--dimensions/-d' option in the command line (ex : 'add_points_in_las.py -i myLas.las -g myPoints.json -d classification=64' - points will have their classification set to 64). The dimension should be present in the initial las ; this is not allowed to add new dimension.
Dev / Build
Contribute
Every time the code is changed, think of updating the version file: pdaltools/_version.py
Please log your changes in CHANGELOG.md
Before committing your changes, run the precommit hooks. They can be installed to run automatically with make install-precommit
Tests
Create the conda environment: make install
Run unit tests: make testing
Pip package
To generate a pip package and deploy it on pypi, use the Makefile at the root of the repo:
make build
: build the librarymake install
: install the library in an editable way (pip -e
)make deploy
: deploy it on pypi
Docker image
To build a docker image with the library installed: make docker-build
To test the docker image: make docker-test
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
File details
Details for the file ign_pdal_tools-1.7.5.tar.gz
.
File metadata
- Download URL: ign_pdal_tools-1.7.5.tar.gz
- Upload date:
- Size: 27.1 kB
- Tags: Source
- Uploaded using Trusted Publishing? Yes
- Uploaded via: twine/4.0.1 CPython/3.11.10
File hashes
Algorithm | Hash digest | |
---|---|---|
SHA256 | 6b4a2e06f53f2d4affba9915f52ed9e719264db7f63e571e49c7e7d566ae9319 |
|
MD5 | 260c81cd47e64cb61f4d0ca32931007c |
|
BLAKE2b-256 | 21569ce6eb9422667b89c488041730d48fcef0c15e1ec7ad6596427790fdb610 |
File details
Details for the file ign_pdal_tools-1.7.5-py3-none-any.whl
.
File metadata
- Download URL: ign_pdal_tools-1.7.5-py3-none-any.whl
- Upload date:
- Size: 23.8 kB
- Tags: Python 3
- Uploaded using Trusted Publishing? Yes
- Uploaded via: twine/4.0.1 CPython/3.11.10
File hashes
Algorithm | Hash digest | |
---|---|---|
SHA256 | 5aece5e9d53fe0fb9606d6830815efaa74a228862e08c3ef3aaa611d579e8a67 |
|
MD5 | 58476304ab038f56d068d04cbdc18525 |
|
BLAKE2b-256 | bb3a59c2d4573ddc2216545bd0d7e0623a93645e15674f6d218f38d8607dd273 |