Forest canopy height is a fundamental parameter of forest structure, and plays a pivotal role in understanding forest biomass allocation, carbon stock, forest productivity, and biodiversity. Spaceborne LiDAR (Light Detection and Ranging) systems, such as GEDI (Global Ecosystem Dynamics Investigation), provide large-scale estimation of ground elevation, canopy height, and other forest parameters. However, these measurements may have uncertainties influenced by topographic factors. This study focuses on the calibration of GEDI L2A and L1B data using an airborne LiDAR point cloud, and the combination of Sentinel-2 multispectral imagery, 1D convolutional neural network (CNN), artificial neural network (ANN), and random forest (RF) for upscaling estimated forest height in the Guangxi Gaofeng Forest Farm. First, various environmental (i.e., slope, solar elevation, etc.) and acquisition parameters (i.e., beam type, Solar elevation, etc.) were used to select and optimize the L2A footprint. Second, pseudo-waveforms were simulated from the airborne LiDAR point cloud and were combined with a 1D CNN model to calibrate the L1B waveform data. Third, the forest height extracted from the calibrated L1B waveforms and selected L2A footprints were compared and assessed, utilizing the CHM derived from the airborne LiDAR point cloud. Finally, the forest height data with higher accuracy were combined with Sentinel-2 multispectral imagery for an upscaling estimation of forest height. The results indicate that through optimization using environmental and acquisition parameters, the ground elevation and forest canopy height extracted from the L2A footprint are generally consistent with airborne LiDAR data (ground elevation: R2 = 0.99, RMSE = 4.99 m; canopy height: R2 = 0.42, RMSE = 5.16 m). Through optimizing, ground elevation extraction error was reduced by 45.5% (RMSE), and the canopy height extraction error was reduced by 30.3% (RMSE). After training a 1D CNN model to calibrate the forest height, the forest height information extracted using L1B has a high accuracy (R2 = 0.84, RMSE = 3.13 m). Compared to the optimized L2A data, the RMSE was reduced by 2.03 m. Combining the more accurate L1B forest height data with Sentinel-2 multispectral imagery and using RF and ANN for the upscaled estimation of the forest height, the RF model has the highest accuracy (R2 = 0.64, RMSE = 4.59 m). The results show that the extrapolation and inversion of GEDI, combined with multispectral remote sensing data, serve as effective tools for obtaining forest height distribution on a large scale. [ABSTRACT FROM AUTHOR]