LiDAR-Based Forest Inventory: Methods and Implementation

Twenty years ago, estimating forest biomass across a large property meant sending field crews to measure sample plots and extrapolating statistically to the entire area. The approach worked, but it was expensive, time-consuming, and left considerable uncertainty about what lay between the plots.

Airborne LiDAR has fundamentally changed this equation. By providing detailed three-dimensional structural data across entire landscapes, LiDAR enables more accurate biomass estimates at lower cost per acre, particularly for large areas. But realizing these benefits requires understanding both the technology’s capabilities and its proper implementation.

How LiDAR Sees Forests

An airborne LiDAR system fires laser pulses toward the ground and measures how long the returns take to come back. Because forest canopies aren’t solid, a single pulse often generates multiple returns: one from the top of a tree, others from branches and understory, and finally one from the ground surface. Processing millions of these returns creates a point cloud that captures the three-dimensional structure of the forest.

Modern airborne systems achieve extraordinary detail. Point densities of 8-20 points per square meter are now standard for manned aircraft — the USGS 3DEP program specifies 8 points per square meter as its QL1 standard, and advances in sensor technology have reduced acquisition costs at this density below what 2 points per square meter cost a decade ago. UAV-mounted systems routinely deliver hundreds of points per square meter, enabling even finer structural resolution over smaller areas. Intensity values from each return indicate surface reflectance characteristics. GPS and inertial measurement units track the aircraft position precisely enough to georeference returns within centimeters.

The processing workflow transforms raw laser returns into usable forest information. Ground classification algorithms separate terrain from vegetation, enabling calculation of canopy heights. Normalized point clouds express each return as height above ground rather than absolute elevation. From there, various metrics can be extracted: height percentiles, canopy density at different strata, crown cover percentage, and statistical measures of structural variation.

Two Approaches to Forest Measurement

LiDAR-based inventory follows one of two fundamental approaches, each with distinct strengths.

The area-based method works statistically. Field crews measure sample plots using conventional techniques: diameter, height, species, and condition for every tree above threshold. LiDAR metrics are extracted for the same plot areas. Regression models then relate the field measurements to the LiDAR data. Once calibrated, these models can predict forest attributes wall-to-wall across the entire LiDAR coverage.

This approach handles complex forests well. Mixed species stands, multi-layered canopies, and variable conditions all yield to statistical modeling given adequate field sampling. Importantly, area-based methods produce robust predictions at relatively low point densities — research has shown that pulse densities as low as 0.5-1 points per square meter yield reliable stand-level estimates of basal area, volume, and biomass, with diminishing returns above that threshold (Gobakken & Næsset, 2008). The method provides uncertainty estimates through prediction intervals and integrates naturally with existing inventory frameworks. Most operational forest carbon projects use area-based methods, a pattern confirmed by the widespread adoption documented in the Canadian Forest Service’s best practices guide (White et al., 2013).

Individual tree detection takes a different path. Algorithms identify discrete tree crowns in the point cloud using local maxima detection, watershed segmentation, or increasingly deep learning approaches such as panoptic segmentation networks. Each detected tree gets measured: location, height, crown diameter, and sometimes crown volume. Species classification can combine LiDAR geometric features with multispectral imagery.

This approach yields tree-level data compatible with growth models and individual management prescriptions. It reveals spatial patterns, gap distributions, and structural diversity relevant to habitat assessment. Higher point densities are beneficial here — reliable individual tree detection typically requires 10-20 points per square meter (Xiang et al., 2024), and international benchmarking has shown that the extraction algorithm matters as much as point density (Kaartinen & Hyyppä, 2012). But accuracy drops in dense forests where crowns overlap extensively. Understory trees frequently go undetected beneath closed canopies. Existing segmentation procedures typically detect over 90% of overstory trees but barely 60% of understory trees. Accurate understory segmentation requires point densities around 170 points per square meter — far above standard operational acquisitions (Hamraz et al., 2017). Computational requirements increase substantially for large areas.

The practical implication: lower point densities can robustly predict area-based statistics, but may not be sufficient when the objective requires mapping every individual tree.

From Structure to Biomass

Converting LiDAR measurements to biomass estimates requires allometric relationships, equations that predict tree or stand biomass from measured dimensions.

For individual tree methods, allometric equations typically use diameter and height. Since LiDAR provides height and crown dimensions but not diameter, additional relationships are needed. Crown-based allometry exists for some species but remains less developed than traditional diameter-based equations — DBH remains the most widely validated predictor in allometric databases, though recent work combining crown area and height into composite indices is narrowing the gap (Wang et al., 2021; Pu et al., 2023). Detection errors, measurement uncertainty, and species misclassification all propagate into final biomass estimates.

Area-based methods develop statistical models directly predicting biomass from LiDAR metrics. Height percentiles, canopy density measures, and intensity statistics serve as predictor variables. Model forms range from multiple linear regression through random forest ensembles to neural networks. Random forest remains the most widely used machine learning algorithm for LiDAR-based biomass estimation due to its robustness with high-dimensional data, though recent multi-modal deep learning architectures combining point cloud, spectral, and topographic inputs are showing competitive results. Validation with independent test data quantifies prediction accuracy through familiar metrics: R-squared, root mean square error, and bias.

Either approach requires field reference data. The quality of final biomass estimates depends critically on the quality and representativeness of ground measurements used for model development.

Getting Implementation Right

Flight planning determines what the LiDAR can see. For area-based inventory, point densities as low as 1-2 points per square meter can produce reliable stand-level predictions, but individual tree detection benefits from 8 points per square meter or higher — the USGS 3DEP QL1 specification. In structurally complex stands where tree-level resolution matters, higher densities improve results. Scan angles beyond 15-20 degrees from nadir introduce bias in canopy density metrics and compromise vertical accuracy (Roussel et al., 2018; Liu et al., 2018). Timing matters: leaf-on conditions for deciduous forests, consistent phenology across the acquisition area. Weather constraints exclude rain, fog, and high winds.

Field sampling design drives model quality. Plot sizes of 0.04-0.1 hectares match typical model grid cells. Circular plots simplify LiDAR metric extraction. Sample intensity depends on forest heterogeneity, typically 0.5-2% of the project area. Stratification across forest types, age classes, and density conditions ensures models generalize across the landscape. Field timing should synchronize with LiDAR acquisition, ideally within the same growing season.

Quality assurance catches problems before they contaminate results. Vertical accuracy should be verified against RTK-GPS ground surveys. Overlap areas between flight lines should show consistent values. Data gaps and low-density patches require investigation. Ground classification performance affects canopy height accuracy and deserves explicit validation.

The Cost Question

LiDAR economics favor larger projects. Acquisition costs for forestry-grade airborne data run $20-100 per acre depending on area size, required point density, and terrain complexity. Field sampling adds $500-2000 per plot for complete measurement. Processing and model development add $10-30 per acre with economies of scale for larger areas.

A key economic advantage is reduced field sampling. With LiDAR data available for stratification and model-assisted estimation, you can cut field plots by 50-75% without meaningful accuracy loss (Gobakken & Næsset, 2008). One recent analysis found a 41% reduction in required plots for a fixed 10% uncertainty target when LiDAR-aided stratification replaced simple random sampling (Castaño-Díaz et al., 2023). And when you account for the economic cost of bad decisions from inaccurate inventory — not just data collection costs — LiDAR pays for itself even at higher acquisition prices (Eid et al., 2004).

For manned aircraft acquisitions, properties under 1,000 acres often find conventional inventory more cost-effective given high mobilization costs. UAV-based LiDAR has shifted this threshold downward: drone acquisitions are now viable for properties as small as a few hundred acres, at costs of $30-120 per acre with point densities far exceeding manned aircraft. Between 1,000 and 10,000 acres, airborne LiDAR becomes competitive, particularly when reduced field sampling intensity is accounted. Above 10,000 acres, LiDAR typically reduces cost per acre substantially while improving accuracy — a comparison on the Malheur National Forest found LiDAR-based inventory comparable in accuracy and cost to traditional stand exams while providing wall-to-wall coverage (Hummel et al., 2011).

Repeat measurements amplify the advantage. Once field plots are established and models calibrated, subsequent LiDAR acquisitions require minimal additional ground work. Monitoring programs with periodic resurvey find LiDAR increasingly attractive over time.

Real-World Applications

At Arbos, we’ve implemented LiDAR-based inventory across diverse forest types and regulatory contexts.

On a 180,000-acre California timberland property, we developed aligned area-based and individual tree models using 200 field plots across mixed conifer and hardwood stands. The resulting inventory achieved under 10% relative RMSE, meeting California’s Compliance Offset Protocol requirements. This became the first LiDAR-based inventory to pass field verification under California’s Compliance Offset Protocol. Total inventory cost ran approximately 40% below what conventional sampling would have required at comparable precision — consistent with the field plot reductions documented in the literature for LiDAR-assisted stratified inventory at this scale.

Pacific Northwest old-growth presents different challenges. Individual tree detection works better in structurally open conditions, though complex crown shapes require careful algorithm tuning. Integration with multispectral imagery enables species classification for biomass model stratification. Monitoring protocols for carbon project permanence verification build on change detection between LiDAR acquisitions.

Tropical forest applications face cloud cover constraints that limit acquisition windows. Stratified sampling across degradation gradients captures the full range of forest conditions. LiDAR transects calibrated with intensive field plots provide biomass baselines. Integration with radar enables cloud-free change detection for ongoing monitoring.

Emerging Technologies

Spaceborne LiDAR now provides global forest structure data. NASA’s GEDI mission, which resumed operations in April 2024 after a period of hibernation, samples forest canopy height and biomass using waveform LiDAR from the International Space Station. ICESat-2 adds photon-counting altimetry for vegetation structure, with biomass estimates now computed at 30-meter resolution from over 19 million observations. Fusion of GEDI and ICESat-2 footprints with Sentinel-2 imagery is producing wall-to-wall canopy height maps at 10-meter resolution — though these modeled products complement rather than replace airborne acquisition for project-level applications.

UAV-based LiDAR systems have matured rapidly. Costs have dropped, point densities routinely exceed hundreds of points per square meter, and operational flexibility enables rapid deployment. Applications include intensive inventory of high-value stands, validation of area-based models, and change detection in localized areas. Under-canopy autonomous flight for terrestrial-LiDAR-like measurements is an emerging capability. Coverage limitations and regulatory constraints keep UAV LiDAR complementary to rather than replacing airborne acquisition for large areas.

New sensor architectures are expanding what airborne platforms can achieve. Single-photon LiDAR systems like the Leica SPL100 enable high-density acquisition from higher altitudes, reducing per-acre costs for large-area mapping. Geiger-mode systems achieve comparable accuracy from even greater heights. Multispectral LiDAR systems collecting returns at multiple wavelengths improve species discrimination — a recent study showed multispectral ALS reduced classification error by 65% compared to single-channel geometric features (Wang et al., 2024). Fire risk assessment and forest health monitoring applications are developing around these spectral capabilities.

Sensor fusion approaches combine LiDAR with optical and radar data. Hyperspectral imagery adds species identification and health assessment capabilities. Synthetic aperture radar provides all-weather monitoring and different penetration characteristics. Deep learning methods — including transformer architectures adapted for 3D point clouds and panoptic segmentation networks — are enabling automated individual tree segmentation at accuracy levels that were recently impractical.

Making It Work

Successful LiDAR-based inventory requires attention throughout the project lifecycle. Acquisition specifications must match inventory objectives. Processing should follow documented, reproducible workflows. Model development demands statistical rigor with proper validation and uncertainty quantification. Reporting should transparently document methods, accuracy assessment, and known limitations.

The technology has matured. LiDAR-based forest inventory now represents standard practice for large-area biomass assessment when implemented properly. The question isn’t whether it works, but whether a specific implementation follows the methodological requirements for its intended application.

References

• Castaño-Díaz, M. et al. (2023). Effect of sample size on the estimation of forest inventory attributes using airborne LiDAR data in large-scale subtropical areas. Annals of Forest Science. DOI: 10.1186/s13595-023-01209-4
• Eid, T., Gobakken, T., & Næsset, E. (2004). Comparing stand inventories based on photo interpretation and laser scanning by means of cost-plus-loss analyses. Scandinavian Journal of Forest Research, 19, 512-523. DOI: 10.1080/02827580410019463
• Gobakken, T. & Næsset, E. (2008). Assessing effects of laser point density, ground sampling intensity, and field sample plot size on biophysical stand properties derived from airborne laser scanner data. Canadian Journal of Forest Research, 38(5), 1095-1109.
• Hamraz, H. et al. (2017). Forest understory trees can be segmented accurately within sufficiently dense airborne laser scanning point clouds. Scientific Reports, 7, 6770. DOI: 10.1038/s41598-017-07200-0
• Hummel, S. et al. (2011). A comparison of accuracy and cost of LiDAR versus stand exam data for landscape management on the Malheur National Forest. Journal of Forestry, 109(5), 267-273.
• Kaartinen, H. & Hyyppä, J. (2012). EuroSDR/ISPRS project: benchmarking of single tree detection methods. Remote Sensing, 4(4), 950-974. DOI: 10.3390/rs4040950
• Liu, J. et al. (2018). Large off-nadir scan angle of airborne LiDAR can severely affect the estimates of forest structure metrics. ISPRS Journal of Photogrammetry and Remote Sensing, 136, 13-25.
• Pu, Y. et al. (2023). Precise aboveground biomass estimation of plantation forest trees using the novel allometric model and UAV-borne LiDAR. Frontiers in Forests and Global Change, 6, 1166349.
• Roussel, J.R. et al. (2018). A mathematical framework to describe the effect of beam incidence angle on metrics derived from airborne LiDAR. Remote Sensing of Environment, 224, 492-508.
• Wang, Q. et al. (2021). Lidar biomass index: A novel solution for tree-level biomass estimation using 3D crown information. Forest Ecology and Management, 499, 119542. DOI: 10.1016/j.foreco.2021.119542
• Wang, Q. et al. (2024). Attribute-aware cross-branch transformer for multispectral ALS tree species classification. Remote Sensing of Environment, 315, 114434.
• White, J.C. et al. (2013). A best practices guide for generating forest inventory attributes from airborne laser scanning data using an area-based approach. Canadian Forest Service, Information Report FI-X-010.
• Xiang, B. et al. (2024). ForAINet: deep learning for panoptic segmentation of forest point clouds. arXiv preprint, 2312.15084.

Need support with LiDAR-based forest inventory? Contact Arbos to discuss project-specific methodology development and implementation.