• Mapping forest height by combining photon-counting LiDAR data and Landsat OLI data. • Douglas-Peucker algorithm can solve the large turns in photon-counting LiDAR data. • The filtering of forest heights can improve the accuracy of forest height map. Large-scale, accurate and detailed forest height map is worthwhile and necessary for understanding and assessing global carbon cycle and biodiversity. The Ice, Cloud and Land Elevation Satellite-2 (ICESat-2) mission, employing a photon-counting LiDAR (PCL) system, offers an opportunity to map global forest height with high resolution. This study aimed to map forest height at a spatial resolution of 30 m by combining Multiple Altimeter Beam Experimental Lidar (MABEL) data with Landsat 8 Operational Land Imager (OLI) data. There are four key steps to accomplish this goal. First, a segmentation method based on the Douglas-Peucker algorithm was proposed to solve the problem of large turns or calibration maneuvers in MABEL data. Second, we estimated forest height and selected the forest height samples with high accuracy and reliability by developing three filters including signal-to-noise ratio (SNR) filter, slope filter, and canopy photons density (CPD) filter. Third, forest height models based on both random forest (RF) and stepwise multiple regression algorithms were developed to establish relationships between the selected forest height samples and predicator variables of Landsat-derived spectral indices, topographic variables and geographic coordinates. Finally, a wall-to-wall forest height map was generated by applying the developed forest height models to predicator variables, and the accuracy of forest height map was validated using airborne LiDAR-derived forest heights. An area of 160, 000 km 2 in southeast Virginia and east North Carolina was chosen for testing the methods proposed in this study. The results demonstrated that the Douglas-Peucker algorithm can effectively solve the MABEL data overlapping issues caused by large turns and calibration maneuvers in flight lines. The suitable filters for selecting forest heights are SNR > 6, terrain slope < 25°, and 40 < CPD < 170. These developed filters can substantially increase the accuracies of forest height models. The results also indicated that RF-derived forest height models achieved higher modeling accuracy than stepwise regression-derived forest height models. RF-derived forest height models yielded coefficient of determination (R 2 ) values of 0.59, 0.68 and 0.62 and RMSE values of 4.55 m, 3.41 m and 4.08 m for deciduous forest, evergreen forest and mixed forest, respectively. Compared to LiDAR-derived forest height, the forest height map produced in this study has a R 2 value of 0.54 and a RMSE value of 6.85 m, which demonstrates that combination of MABEL data and Landsat 8 OLI data can be used to generate forest height maps with a spatial resolution of 30 m.