
There are a number of methods available, check this
Jupyter notebook for an introduction.
Functions include reading .las files to numpy array, writing to .las as well as a number of methods to dice, slice and tile
LiDAR data.
The main set of functions found in forestlas.canopyComplexity.
These allow you to derive metrics of vertical canopy structure such as Pgap and also estimate number of canopy layers.
More information can be found in this paper Wilkes, P. et al. (2016). Using discrete-return airborne laser scanning to
quantify number of canopy strata across diverse forest types. Methods in Ecology and Evolution, 7(6), 700–712.
This research was funded by the Australian Postgraduate Award, Cooperative Research Centre for Spatial Information under Project 2.07, TERN/AusCover and Commonwealth Scientific and IndustrialResearch Organisation (CSIRO) Postgraduate Scholarship.