Airborne laser scanning (ALS) and terrestrial laser scanning (TLS) are two ways to obtain forest three-dimensional (3D) spatial information. Due to canopy occlusion and the features of different scanning methods, some of the forest point clouds acquired by a single scanning platform may be missing, resulting in an inaccurate estimation of forest structure parameters. Hence, the registration of ALS and TLS point clouds is an alternative for improving the estimation accuracy of forest structure parameters. Currently, forest point cloud registration is mainly conducted based on individual tree attributes (e.g., location, diameter at breast height, and tree height), but the registration is affected by individual tree segmentation and is inefficient. In this study, we proposed a method to automatically fuse ALS and TLS point clouds by using feature points of canopy gap shapes. First, the ALS and TLS canopy gap boundary vectors were extracted by the canopy point cloud density model, and the turning or feature points were obtained from the canopy gap vectors using the weighted effective area (WEA) algorithm. The feature points were then aligned, the transformation parameters were solved using the coherent point drift (CPD) algorithm, and the TLS point clouds were further aligned using the recovery transformation matrix and refined by utilizing the iterative closest point (ICP) algorithm. Finally, individual tree segmentations were performed to estimate tree parameters using the TLS and fusion point clouds, respectively. The results show that the proposed method achieved more accurate registration of ALS and TLS point clouds in four plots, with the average distance residuals of coarse and fine registration of 194.83 cm and 2.14 cm being much smaller compared with those from the widely used crown feature point-based method. Using the fused point cloud data led to more accurate estimates of tree height than using the TLS point cloud data alone. Thus, the proposed method has the potential to improve the registration of ALS and TLS point cloud data and the accuracy of tree height estimation.
Forest canopy height (FCH) is a vital indicator for assessing forest health and ecosystem service capacity. Over the past two decades, full-waveform (FW) LiDAR has been widely employed for estimating forest biophysical variables due to its high precision in measuring vertical forest structures. However, the impact of terrain undulations on forest parameter estimation remains challenging. To address this issue, this study proposes a bisection approximation decomposition (BAD) method for processing GEDI L1B data and FCH estimation. The BAD method analyzes the energy composition of simplified echo signals and determines the fitting parameters by integrating overall signal energy, the differences in unresolved signals, and the similarity of inter-forest signal characteristics. FCH is subsequently estimated based on waveform peak positions. By dynamically adjusting segmentation points and Gaussian fitting parameters, the BAD method achieved precise separation of mixed canopy and ground signals, substantially enhancing the physical realism and applicability of decomposition results. The effectiveness and robustness of the BAD method for FCH estimation were evaluated using 2049 footprints across varying slope conditions in the Harvard Forest region of Petersham, Massachusetts. The results demonstrated that digital terrain models (DTMs) extracted using the GEDI data and the BAD method exhibited high consistency with the DTMs derived using airborne laser scanning (ALS) data (coefficient of determination R2 > 0.99). Compared with traditional Gaussian decomposition (GD), wavelet decomposition (WD), and deconvolution decomposition (DD) methods, the BAD method showed significant advantages in FCH estimation, achieved the smallest relative root mean square error (rRMSE) of 17.19% and greatest mean estimation accuracy of 84.57%, and reduced the rRMSE by 10.74%, 21.49%, and 28.93% compared to GD, WD, and DD methods, respectively. Moreover, the BAD method exhibited a significantly stronger correlation with ALS-derived canopy height mode data than the relative height metrics from GEDI L2A products (r = 0.84, p < 0.01). The robustness and adaptability of the BAD method to complex terrain conditions provide great potential for forest parameters using GEDI data.
Tree height is a critical indicator for estimating forest stock and can be effectively acquired by UAV-LiDAR. Ground filtering works to classify ground points and non-ground points and can impact the tree height extraction results, while the points classification quality obtained by ordinary filtering methods is limited in complex forest conditions. A partitioned cloth simulation filtering (PCSF) method based on different vegetation cover was proposed in this study to improve the classification accuracy, and tree heights were extracted to demonstrate the effectiveness of the proposed method. UAV-LiDAR data and field measurements collected from the Lutou experimental forest farm in the southern subtropical forest region of China were used for validation, and the slope-based filtering, progressive triangulated irregular network densification filtering (PTD), moving surface fitting filtering (MSFF), and CSF were adopted for comparisons. The results showed that the proposed method yielded the best ground filtering effect, reducing the filtering total error by 2.12%–4.22% compared with other methods, and the relative root mean squared error (rRMSE) of extracted tree heights was reduced by 1.24%–3.84%, respectively. The proposed method can achieve a satisfactory filtering effect and tree height extraction result, which provides a methodological basis to precisely extract tree heights in large-scale forests.