laxpy. Let it handle your spatial queries for you!
laxpy is a way to read and write
.lax files generated from
lasindex in Python. I have ambitions to merge this into
laspy codebase, but for now this standalone package can handle spatial queries in conjuction with
.lax files that match the name of the
.las file you are interested in querying.
Indexing LiDAR data is only relevant if the entire file is not needed. This most frequently happens when LiDAR files need to be clipped using some geometry smaller than the original file extent. Indexing allows loading smaller chunks of the file into memory, rather than the whole thing, at the cost of some initial overhead.
Whether or not indexing ultimately increases computational efficiency depends on the clipping geometry, and the structure and size of the LiDAR file.
Current Release: 0.1.9
laxpy is in a stable condition.
# via pip pip install laxpy # via GitHub git clone https://github.com/brycefrank/laxpy.git cd laxpy pip install .
Clipping a Polygon
laxpy clips a
shapely polygon by adjusting the memory map (an internal mechanism in
laspy) of the original reference object.
Long story short, this is done in place in the following manner:
from laxpy import IndexedLAS my_las = IndexedLAS('my_las.las') my_las.map_polygon(my_shapely_polygon) # Print the points within the polygon print(my_las.points)
Note that an
IndexedLAS object inherits directly from
laspy.file.File, and all dimensions are available for querying,
.y, etc. etc.