Finding the boundary¶
This exercise uses PDAL to find a tight-fitting boundary of an aerial scan.
Printing the coordinates of the boundary for the file is quite simple using a
pdal info call, but visualizing the boundary is more complicated. To
complete this exercise, we are going to use QGIS to view the
boundary, which means we must first install it on our system.
We are going to run using the Uncompahgre data in the
pdal info ./exercises/analysis/density/uncompahgre.laz --boundary
… a giant blizzard of coordinate output scrolls across our terminal. Not very useful.
Instead, let’s generate some kind of vector output we can visualize with
pdal tindex is the “tile index” command, and it outputs a
vector geometry file for each point cloud file it reads. It generates this
boundary using the same mechanism we invoked above – filters.hexbin.
We can leverage this capability to output a contiguous boundary of the
1 2 3
pdal tindex create --tindex ./exercises/analysis/boundary/boundary.sqlite \ --filespec ./exercises/analysis/density/uncompahgre.laz \ -f SQLite
1 2 3
pdal tindex create --tindex ./exercises/analysis/boundary/boundary.sqlite ^ --filespec ./exercises/analysis/density/uncompahgre.laz ^ -f SQLite
Once we’ve run the tindex, we can now visualize our output:
Open QGIS and select Add Vector Layer:
Navigate to the
exercises/analysis/boundary directory and
then open the
The PDAL boundary computation is an approximation based on a hexagon tessellation. It uses the software at http://github.com/hobu/hexer to do this task.
The tindex can be used to generate boundaries for large collections of data. A boundary-based indexing scheme is commonly used in LiDAR processing, and PDAL supports it through the
tindexapplication. You can also use this command to merge data together (query across boundaries, for example).