Forest canopy height is an important parameter for studying biodiversity and the carbon cycle. A variety of techniques for mapping forest height using remote sensing data have been successfully developed in recent years. However, the demands for forest height mapping in practical applications are often not met, due to the lack of corresponding remote sensing data. In such cases, it would be useful to exploit the latest, cheaper datasets and combine them with free datasets for the mapping of forest canopy height. In this study, we proposed a method that combined ZiYuan-3 (ZY-3) stereo images, Shuttle Radar Topography Mission global 1 arc second data (SRTMGL1), and Landsat 8 Operational Land Imager (OLI) surface reflectance data. The method consisted of three procedures: First, we extracted a digital surface model (DSM) from the ZY-3, using photogrammetry methods and subtracted the SRTMGL1 to obtain a crude canopy height model (CHM). Second, we refined the crude CHM and correlated it with the topographically corrected Landsat 8 surface reflectance data, the vegetation indices, and the forest types through a Random Forest model. Third, we extrapolated the model to the entire study area covered by the Landsat data, and obtained a wall-to-wall forest canopy height product with 30 m × 30 m spatial resolution. The performance of the model was evaluated by the Random Forest's out-of-bag estimation, which yielded a coefficient of determination (R 2 ) of 0.53 and a root mean square error (RMSE) of 3.28 m. We validated the predicted forest canopy height using the mean forest height measured in the field survey plots. The validation result showed an R 2 of 0.62 and a RMSE of 2.64 m.Moderate Resolution Imaging Spectroradiometer (MODIS) data [2,3,[14][15][16]. The SAR is well known for its high temporal resolution, because it is less affected by weather and illumination conditions [17]. The SAR-based techniques for forest canopy height mapping include radargrammetry, interferometry SAR (InSAR), and polarimetric interferometric SAR (PolInSAR). The radargrammetry is based on SAR stereo images [5,18,19]. The InSAR is based on the phase differences between two complex SAR images [6,20]. The PolInSAR includes phase difference methods and the model-based methods, which rely on the complex coherent coefficient model [7,21]. In these methods, the TerraSAR-X (X-band), Radarsat-2 (C-band), and ALOS-2 (L-band) are commonly used SAR data sources. It is expected that the techniques for mapping forest canopy height using the LiDAR and SAR will be further developed when data from the ICESat-2 [22], the GEDI [23], the Biomass (P-band) (due for launch in 2021) [24], and the Tandem-L (L-band) (due for launch in 2022) [25] missions become available. The airborne DP systems have a large field of view, high cruising altitude, fast data acquisition speed, and easy flight planning [26]. The spaceborne DP systems can perform continuous and repetitive observations. Several spaceborne systems, such as the WorldView series, provide sub-meter resolution ...