Point cloud data processing
Project description
PDAL Python support allows you to process data with PDAL into Numpy arrays. It supports embedding Python in PDAL pipelines with the readers.numpy and filters.python stages, and it provides a PDAL extension module to control Python interaction with PDAL.
Additionally, you can use it to fetch schema and metadata from PDAL operations.
Installation
PyPI
PDAL Python support is installable via PyPI:
pip install PDALGitHub
The repository for PDAL’s Python extension is available at https://github.com/PDAL/python
Python support released independently from PDAL itself as of PDAL 1.7.
Usage
Simple
Given the following pipeline, which simply reads an ASPRS LAS file and sorts it by the X dimension:
json = """
{
"pipeline": [
"1.2-with-color.las",
{
"type": "filters.sort",
"dimension": "X"
}
]
}"""
import pdal
pipeline = pdal.Pipeline(json)
count = pipeline.execute()
arrays = pipeline.arrays
metadata = pipeline.metadata
log = pipeline.logReading using Numpy Arrays
The following more complex scenario demonstrates the full cycling between PDAL and Python:
Read a small testfile from GitHub into a Numpy array
Filters those arrays with Numpy for Intensity
Pass the filtered array to PDAL to be filtered again
Write the filtered array to an LAS file.
data = "https://github.com/PDAL/PDAL/blob/master/test/data/las/1.2-with-color.las?raw=true"
json = """
{
"pipeline": [
{
"type": "readers.las",
"filename": "%s"
}
]
}"""
import pdal
import numpy as np
pipeline = pdal.Pipeline(json % data)
count = pipeline.execute()
# get the data from the first array
# [array([(637012.24, 849028.31, 431.66, 143, 1,
# 1, 1, 0, 1, -9., 132, 7326, 245380.78254963, 68, 77, 88),
# dtype=[('X', '<f8'), ('Y', '<f8'), ('Z', '<f8'), ('Intensity', '<u2'),
# ('ReturnNumber', 'u1'), ('NumberOfReturns', 'u1'), ('ScanDirectionFlag', 'u1'),
# ('EdgeOfFlightLine', 'u1'), ('Classification', 'u1'), ('ScanAngleRank', '<f4'),
# ('UserData', 'u1'), ('PointSourceId', '<u2'),
# ('GpsTime', '<f8'), ('Red', '<u2'), ('Green', '<u2'), ('Blue', '<u2')])
arr = pipeline.arrays[0]
print (len(arr)) # 1065 points
# Filter out entries that have intensity < 50
intensity = arr[arr['Intensity'] > 30]
print (len(intensity)) # 704 points
# Now use pdal to clamp points that have intensity
# 100 <= v < 300, and there are 387
clamp =u"""{
"pipeline":[
{
"type":"filters.range",
"limits":"Intensity[100:300)"
}
]
}"""
p = pdal.Pipeline(clamp, [intensity])
count = p.execute()
clamped = p.arrays[0]
print (count)
# Write our intensity data to an LAS file
output =u"""{
"pipeline":[
{
"type":"writers.las",
"filename":"clamped.las",
"offset_x":"auto",
"offset_y":"auto",
"offset_z":"auto",
"scale_x":0.01,
"scale_y":0.01,
"scale_z":0.01
}
]
}"""
p = pdal.Pipeline(output, [clamped])
count = p.execute()
print (count)
Requirements
PDAL 2.2+
Python >=3.6
Cython (eg
pip install cython)Numpy (eg
pip install numpy)Packaging (eg
pip install packaging)scikit-build (eg
pip install scikit-build)
Changes
2.3.5
Fix memory leak https://github.com/PDAL/python/pull/74
Handle metadata with invalid unicode by erroring https://github.com/PDAL/python/pull/74
2.3.0
PDAL Python support 2.3.0 requires PDAL 2.1+. Older PDAL base libraries likely will not work.
Python support built using scikit-build
readers.numpy and filters.python are installed along with the extension.
Pipeline can take in a list of arrays that are passed to readers.numpy
readers.numpy now supports functions that return arrays. See https://pdal.io/stages/readers.numpy.html for more detail.
2.0.0
PDAL Python extension is now in its own repository on its own release schedule at https://github.com/PDAL/python
Extension now builds and works under PDAL OSGeo4W64 on Windows.
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
File details
Details for the file PDAL-2.3.6.tar.gz.
File metadata
- Download URL: PDAL-2.3.6.tar.gz
- Upload date:
- Size: 942.8 kB
- Tags: Source
- Uploaded using Trusted Publishing? No
- Uploaded via: twine/3.2.0 pkginfo/1.6.1 requests/2.25.0 setuptools/50.3.2 requests-toolbelt/0.9.1 tqdm/4.54.1 CPython/3.8.6
File hashes
| Algorithm | Hash digest | |
|---|---|---|
| SHA256 | b5e22c79043081e6b60d948af1f1036e4a5b1a503f843c0cb49cb669eff59c5e |
|
| MD5 | 997c229bdde3659a715604b7c12031bb |
|
| BLAKE2b-256 | 1eb31e736d42e0579048b1a4f9b686cd252a18d230603de2985220d9bc9cda4a |
