Forest canopy height is an important indicator of the forest ecosystem, and an accurate assessment of forest canopy height on a large scale is of great significance for forest resource quantification and carbon sequestration. The retrieval of canopy height based on remote sensing provides a possibility for studying forest ecosystems. This study proposes a new method for estimating forest canopy height based on remote sensing. In this method, the GEDI satellite and ICESat-2 satellite, which are different types of space-borne lidar products, are used to cooperate with the Landsat 9 image and SRTM terrain data, respectively. Two forest canopy height-retrieval models based on multi-source remote sensing integration are obtained using a random forest regression (RFR) algorithm. The study, conducted at a forest site in the northeastern United States, synthesized various remote sensing data sets to produce a robust canopy height model. First, we extracted relative canopy height products, multispectral features, and topographic data from GEDI, ICESat-2, Landsat 9, and SRTM images, respectively. The importance of each variable was assessed, and the random forest algorithm was used to analyze each variable statistically. Then, the random forest regression algorithm was used to combine these variables and construct the forest canopy height model. Validation with airborne laser scanning (ALS) data shows that the GEDI and ICESat-2 models using a single data source achieve better accuracy than the Landsat 9 model. Notably, the combination of GEDI, Landsat 9, and SRTM data (R = 0.92, MAE = 1.91 m, RMSE = 2.78 m, and rRMSE = 12.64%) and a combination of ICESat-2, Landsat 9, and SRTM data (R = 0.89, MAE = 1.84 m, RMSE = 2.54 m, and rRMSE = 10.75%). Compared with the least accurate Landsat 9 model, R increased by 29.58%, 93.48%, MAE by 44.64%, 46.20%, RMSE by 42.80%, 49.40%, and the rRMSE was increased by 42.86% and 49.32%, respectively. These results fully evaluate and discuss the practical performance and benefits of multi-source data retrieval of forest canopy height by combining space-borne lidar data with Landsat 9 data, which is of great significance for understanding forest structure and dynamics. The study provides a reliable methodology for estimating forest canopy height and valuable insights into forest resource management and its contribution to global climate change.

Canopy height is a crucial indicator for assessing the structure and function of the forest ecosystems. It plays a significant role in carbon sequestration, sink enhancement, and promoting green development. This study aimed to evaluate the accuracy of GEDI L2A version 2 data in estimating ground elevation and canopy height by comparing it with airborne laser scanning (ALS) data. Among the six algorithms provided by the GEDI L2A data, algorithm a2 demonstrated higher accuracy than the others in detecting ground elevation and canopy height. Additionally, a relatively strong correlation (R-squared = 0.35) was observed between rh95 for GEDI L2A and RH90 for ALS. To enhance the accuracy of canopy height estimation, this study proposed three backpropagation (BP) neural network inversion models based on GEDI, Landsat 8 OLI, and Landsat 9 OLI-2 data. Multiple sets of relative heights and vegetation indices were extracted from the GEDI and Landsat datasets. The random forest (RF) algorithm was employed to select feature variables with a cumulative importance score of 90% for training the BP neural network inversion models. Validation against RH90 of ALS revealed that the GEDI model outperformed the OLI or OLI-2 data models in terms of accuracy. Moreover, the quality improvement of OLI-2 data relative to OLI data contributed to enhanced inversion accuracy. Overall, the models based on a single dataset exhibited relatively low accuracy. Hence, this study proposed the GEDI and OLI and GEDI and OLI-2 models, which combine the two types of data. The results demonstrated that the combined model integrating GEDI and OLI-2 data exhibited the highest performance. Compared to the weakest OLI data model, the inversion accuracy R-squared improved from 0.38 to 0.74, and the MAE, RMSE, and rRMSE decreased by 1.21 m, 1.81 m, and 8.09%, respectively. These findings offer valuable insights for the remote sensing monitoring of forest sustainability.