Improving Tropical Forest Canopy Height Mapping by Fusion of Sentinel-1/2 and Bias-Corrected ICESat-2–GEDI Data

Accurately estimating the forest canopy height is essential for quantifying forest biomass and carbon storage. Recently, the ICESat-2 and GEDI spaceborne LiDAR missions have significantly advanced global canopy height mapping. However, due to inherent sensor limitations, their footprint-level estimates often show systematic bias. Tall forests tend to be underestimated, while short forests are often overestimated. To address this issue, we used coincident G-LiHT airborne LiDAR measurements to correct footprint-level canopy heights from both ICESat-2 and GEDI, aiming to improve the canopy height retrieval accuracy across Puerto Rico’s tropical forests. The bias-corrected LiDAR dataset was then combined with multi-source predictors derived from Sentinel-1/2 and the 3DEP DEM. Using these inputs, we trained a canopy height inversion model based on the AutoGluon stacking ensemble method. Accuracy assessments show that, compared to models trained on uncorrected single-source LiDAR data, the new model built on the bias-corrected ICESat-2–GEDI fusion outperformed in both overall accuracy and consistency across canopy height gradients. The final model achieved a correlation coefficient (R) of 0.80, with a root mean square error (RMSE) of 3.72 m and a relative RMSE of 0.22. The proposed approach offers a robust and transferable approach for high-resolution canopy structure mapping and provides valuable support for carbon accounting and tropical forest management.