Pure Julia reader of lasindex .lax files
Switch branches/tags
Nothing to show
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Permalink
Failed to load latest commit information.
example
sample
src
test
.codecov.yml
.gitignore
.travis.yml
LICENSE
README.md
REQUIRE
appveyor.yml

README.md

Build Status Build status codecov

LASindex.jl

Pure Julia reader of lasindex .lax files. Translates the .lax file to a quadtree from the RegionTrees package and provides a bounding box intersect method to retrieve a vector of ranges.

julia> using LASindex
julia> using FileIO

julia> laxfile = joinpath("sample", "Palm Beach Post Hurricane.lax");
julia> qt = load(laxfile)
INFO: Processed Palm Beach Post Hurricane.lax with 1924631 points.
Cell: RegionTrees.HyperRectangle{2,Float64}([955000.0, 885000.0], [32000.0, 32000.0])

julia> using StaticArrays
julia> bbox = SVector(955000.0, 985002.0, 885000.0, 988800.0);
julia> r = LASindex.intersect(qt, bbox)
1-element Array{UnitRange{Integer},1}:
 1:1924631

The resulting ranges can be used in combination with LasIO or LazIO to stream a small subset of a (larger than memory) dataset.

using LasIO
header, points = load("Palm Beach Post Hurricane.las", mmap=true)

intersected_points = points[LASindex.intersect(qt, bbox)]

Note that all points inside the bounding box are given, but not all points given are inside the bounding box. In other words, because of how lasindex groups ranges together, some ranges will include points outside the bounding box.

Background

LAX files are quadtree indexes used by the LASTools suite [1] if present. You can generate them with lasindex -i *.laz[2]. There's a good introduction to lasindex here [3].

Future

I'd rather use Cxx to call the lasindex c++ code directly, but at the moment the Cxx package is harder to install on Windows than compiling the lastools shared library itself.

An example how to interface with lastools directly is provided in example/cxx.jl