Continuous mapping of forest canopy height using ICESat-2 data and a weighted kernel integration of multi-temporal multi-source remote sensing data aided by Google Earth Engine.

Forest Canopy Height (FCH) is a crucial parameter that offers valuable insights into forest structure. Spaceborne LiDAR missions provide accurate FCH measurements, but a significant challenge is their point-based measurements lacking spatial continuity. This study integrated ICESat-2's ATL08-derived FCH values with multi-temporal and multi-source remote sensing (RS) datasets to generate continuous FCH maps for northern forests in Iran. Sentinel-1/2, ALOS-2 PALSAR-2, and FABDEM datasets were prepared in Google Earth Engine (GEE) for FCH mapping, each possessing unique spatial and geometrical characteristics that differ from those of the ATL08 product. Given the importance of accurately representing the geometrical characteristics of the ATL08 segments in modeling FCH, a novel Weighted Kernel (WK) approach was proposed in this paper. The WK approach could better represent the RS datasets within the ATL08 ground segments compared to other commonly used resampling approaches. The correlation between all RS data features improved by approximately 6% compared to previously employed approaches, indicating that the RS data features derived after convolving the WK approach are more predictive of FCH values. Furthermore, the WK approach demonstrated superior performance among machine learning models, with random forests outperforming other models, achieving a coefficient of determination (R) of 0.71, root mean square error (RMSE) of 4.92 m, and mean absolute percentage error (MAPE) of 29.95%. Furthermore, in contrast to previous studies using only summer datasets, this study included spring and autumn data from Sentinel-1/2, resulting in a 6% increase in R and a 0.5-m decrease in RMSE. The proposed methodology filled the research gaps and improved the accuracy of FCH estimations.