Estimating Leaf Area Index (LAI), leaf area density, and leaf inclination angles using Terrestrial Laser Scanning (Boise Center Aerospace Lab and SPRUCE)
Welcome!
The main goal of this project is to estimate annual canopy structural metrics for SPRUCE plots from August 2015 - August 2022. To do so, we use a voxel-based contact frequency model. In this model, the occupied voxels are quantified for each canopy layer, providing the leaf area density of each layer. Summing these layers results in the LAI. The vegetation structural data were measured from twelve plots (12 m in diameter) in the Marcell Experimental Forest in northern Minnesota, USA 2015 - 2022. The lidar data comes from a Riegl VZ-1000 TLS with a 1550 nm laser. The original data is publically available:
The scripts provided in this project perform the following steps:
- Separate leaf and wood points.
- Voxelize the point cloud data.
- Compute normal vectors and estimate leaf angles.
- Slice the tree plots and estimate contact frequencies for each slice to get leaf area density.
- Sum leaf area densities to estimate the LAI of each plot.
Below you can find the project organization.
├── LICENSE
├── README.md <- The top-level README file overview of the project.
├── LAI_SPRUCE_dest_0325.csv <- Destructive tree LAI estimates used in validation regression.
├── Canopy_Structure.ipynb <- Runs full workflow to estimate all structural metrics using functions from utils_final.py.
├── Leaf_Volume.py <- Contains functions to estimate the voxel volume of plots.
├── SPRUCE_RF.ipynb <- Random Forest Classifier to separate leaves and wood.
├── allLAI_plots.ipynb <- All LAI results (10 model iterations) scaled by convex hull.
├── destructive_trees.zip <- Destructively harvested tree files used as model validation data.
├── leaf_volume_output.ipynb <- Notebook containing Leaf_Volume.py functions applied.
├── utils_final.py <- Contains functions to estimate all structural metrics.
├── RF_training_txtfiles.zip <- Contains training files for SPRUCE_RF.
├── environment.yml <- The conda environment generated with 'conda env export > environment.yml', name for activation: canopy_structure_env