Direct Estimation of Tree Volume and Aboveground Biomass Using Deep Regression with Synthetic Lidar Data
Organizations: The University of Melbourne · Private Forestry Service Victoria, Australia · Nicon Rural Services, Australia
Abstract
Accurate estimation of forest biomass is crucial for monitoring carbon sequestration and informing climate change mitigation strategies. Existing methods often rely on allometric models, which estimate individual tree biomass by relating it to measurable biophysical parameters, e.g., trunk diameter and height. This indirect approach is limited in accuracy due to measurement uncertainties and the inherently approximate nature of allometric equations, which may not fully account for the variability in tree characteristics and forest conditions. This study proposes a direct approach that leverages synthetic point cloud data to train a deep regression network, which is then applied to real point clouds for plot-level wood volume and aboveground biomass (AGB) estimation. We created synthetic 3D forest plots with ground truth volume, which were then converted into point cloud data using a lidar simulator. These point clouds were subsequently used to train deep regression networks based on PointNet, PointNet++, DGCNN, and PointConv. When applied to synthetic data, the deep regression networks achieved mean absolute percentage error (MAPE) values ranging from 1.69% to 8.11%. The trained networks were then applied to real lidar data to estimate volume and AGB. When compared against field measurements, our direct approach showed discrepancies of 2% to 20%. In contrast, indirect approaches based on individual tree segmentation followed by allometric conversion, as well as FullCAM, exhibited substantially large underestimation, with discrepancies ranging from 27% to 85%. Our results highlight the potential of integrating synthetic data with deep learning for efficient and scalable forest carbon estimation at plot level.
Figures & tables
| Model (domain in | Diameter | Correction | ||
|---|---|---|---|---|
| parentheses) | measurement | factor (c) | ||
| Eucalypt ( cm) | -2.016 | 2.375 | 1.0668 | |
| Multi ( cm) | -2.757 | 2.474 | 1.0775 | |
| Shrub ( cm) | -3.007 | 2.428 | 1.1281 | |
| Other ( cm) | -1.693 | 2.220 | 1.0436 |
| Setting | Riegl VUX-1UAV22 | HOVERMAP ST |
|---|---|---|
| Altitude (m) | 35 | 30–40 |
| Flight speed (m/s) | 2 | Not available |
| Pulse repetition frequency (kHz) | 200 | 300 |
| Scan frequency (Hz) | 80 | 20 |
| Max. returns per pulse | 5 | 2 |
| Scan angle (°) | 180 | 360 |
| Hyperparameter | Value |
|---|---|
| Initial learning rate | 0.001 |
| Batch size | 16 |
| Number of epochs | 200 |
| Optimizer | AdamW |
| Weight decay | |
| Learning rate scheduler | Cosine annealing warm restarts ( Loshchilov and Hutter, 2016 ) |
| Metric | Synthetic | MF1995.1 | MF1999.2 | MF2007.2 | HP2010.2 | Knewleaves20 | ||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| RS | FPS | RS | FPS | RS | FPS | RS | FPS | RS | FPS | RS | FPS | |
| Plot area (m 2 ) | 300.32 50.38 | 315.15 52.51 | 301.74 53.89 | 313.33 55.40 | 279.73 44.02 | 288.20 44.99 | 339.92 69.24 | 353.64 75.94 | 277.85 38.52 | 284.60 39.26 | 335.20 49.72 | 341.85 42.10 |
| Point density (m -2 ) | 7.03 1.28 | 6.69 1.21 | 7.05 1.53 | 6.79 1.48 | 7.53 1.37 | 7.30 1.32 | 6.45 2.35 | 6.23 2.30 | 7.50 0.97 | 7.32 0.94 | 6.33 1.58 | 6.12 1.07 |
| Point density (m -3 ) | 0.59 0.14 | 0.45 0.08 | 0.45 0.12 | 0.42 0.12 | 0.50 0.15 | 0.47 0.14 | 0.46 0.30 | 0.43 0.28 | 0.50 0.13 | 0.48 0.12 | 0.41 0.17 | 0.38 0.13 |
| Average point spacing (cm) | 0.33 0.03 | 0.71 0.06 | 0.26 0.04 | 0.59 0.10 | 0.28 0.04 | 0.60 0.10 | 0.32 0.05 | 0.72 0.12 | 0.32 0.04 | 0.71 0.10 | 0.41 0.08 | 1 0.19 |
| Model | Random Sampling | Farthest Point Sampling | ||
|---|---|---|---|---|
| Validation Error (m 3 ) | MAPE (%) | Validation Error (m 3 ) | MAPE (%) | |
| PointNet | 1.034 0.019 | 8.999 0.189 | 0.959 0.029 | 8.111 0.204 |
| PointNet++ | 0.151 0.016 | 1.438 0.178 | 0.190 0.010 | 1.797 0.149 |
| DGCNN | 0.195 0.010 | 1.627 0.060 | 0.202 0.015 | 1.691 0.146 |
| PointConv | 0.202 0.019 | 1.705 0.138 | 0.286 0.040 | 2.307 0.254 |
| Model | MF1995.1 | MF1999.2 | MF2007.2 | HP2010.2 | Knewleave20 | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| AGB (t/ha) | C (t/ha) | AGB (t/ha) | C (t/ha) | AGB (t/ha) | C (t/ha) | AGB (t/ha) | C (t/ha) | AGB (t/ha) | C (t/ha) | |
| PointNet | 215 20 | 107 10 | 279 30 | 139 15 | 195 13 | 97 7 | 277 34 | 139 17 | 206 27 | 103 14 |
| PointNet++ | 200 2 | 100 1 | 237 2 | 119 1 | 181 1 | 91 1 | 225 4 | 113 2 | 205 7 | 103 3 |
| DGCNN | 154 12 | 77 6 | 200 12 | 100 6 | 149 8 | 75 4 | 205 15 | 103 7 | 170 11 | 85 5 |
| PointConv | 186 26 | 93 13 | 201 17 | 100 9 | 150 18 | 75 9 | 188 28 | 94 14 | 222 89 | 111 44 |
| Model | MF1995.1 | MF1999.2 | MF2007.2 | HP2010.2 | Knewleave20 | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| AGB (t/ha) | C (t/ha) | AGB (t/ha) | C (t/ha) | AGB (t/ha) | C (t/ha) | AGB (t/ha) | C (t/ha) | AGB (t/ha) | C (t/ha) | |
| PointNet | 130 15 | 65 7 | 193 33 | 97 16 | 125 19 | 63 10 | 179 37 | 89 19 | 114 5 | 57 3 |
| PointNet++ | 135 2 | 67 1 | 188 3 | 94 2 | 119 1 | 60 1 | 170 4 | 85 2 | 113 2 | 57 1 |
| DGCNN | 125 5 | 63 3 | 170 9 | 85 4 | 112 5 | 56 3 | 158 8 | 79 4 | 104 4 | 52 2 |
| PointConv | 130 7 | 65 3 | 181 8 | 91 4 | 117 5 | 58 2 | 164 11 | 82 6 | 115 4 | 58 2 |
| Method | MF1995.1 | MF1999.2 | MF2007.2 | HP2010.2 | Knewleave20 | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| AGB | C | AGB | C | AGB | C | AGB | C | AGB | C | |
| (t/ha) | (t/ha) | (t/ha) | (t/ha) | (t/ha) | (t/ha) | (t/ha) | (t/ha) | (t/ha) | (t/ha) | |
| CHM-segmentation | 79 | 40 | 70 | 35 | 50 | 25 | 55 | 28 | 67 | 34 |
| ( Khoshelham et al., 2023 ) | ||||||||||
| TreeLearn | 143 | 72 | 132 | 66 | 129 | 65 | 121 | 61 | 162 | 81 |
| ( Henrich et al., 2024 ) | ||||||||||