The Ice, Cloud and Land Elevation Satellite-2 (ICESat-2) mission employs a micro-pulse photon-counting LiDAR system for mapping and monitoring the biomass and carbon of terrestrial ecosystems over large areas. In preparation for ICESat-2 data processing and applications, this paper aimed to develop and validate an effective algorithm for better estimating ground elevation and vegetation height from photon-counting LiDAR data. Our new proposed algorithm consists of three key steps. Firstly, the noise photons were filtered out using a noise removal algorithm based on localized statistical analysis. Secondly, we classified the signal photons into canopy photons and ground photons by conducting a series of operations, including elevation frequency histogram building, empirical mode decomposition (EMD), and progressive densification. At the same time, we also identified the top of canopy (TOC) photons from canopy photons by percentile statistics method. Thereafter, the ground and TOC surfaces were generated from ground photons and TOC photons by cubic spline interpolation, respectively. Finally, the ground elevation and vegetation height were estimated by retrieved ground and TOC surfaces. The results indicate that the noise removal algorithm is effective in identifying background noise and preserving signal photons. The retrieved ground elevation is more accurate than the retrieved vegetation height, and the results of nighttime data are better than those of the corresponding daytime data. Specifically, the root-mean-square error (RMSE) values of ground elevation estimates range from 2.25 to 6.45 m for daytime data and 2.03 to 6.03 m for nighttime data. The RMSE values of vegetation height estimates range from 4.63 to 8.92 m for daytime data and 4.55 to 8.65 m for nighttime data. Our algorithm performs better than the previous algorithms in estimating ground elevation and vegetation height due to lower RMSE values. Additionally, the results also illuminate that the photon classification algorithm effectively reduces the negative effects of slope and vegetation coverage. Overall, our paper provides an effective solution for estimating ground elevation and vegetation height from micro-pulse photon-counting LiDAR data.

The terrain slope is one of the most important surface characteristics for quantifying the Earth surface processes. Space-borne LiDAR sensors have produced high-accuracy and large-area terrain measurement within the footprint. However, rigorous procedures are required to accurately estimate the terrain slope especially within the large footprint since the estimated slope is likely affected by footprint size, shape, orientation, and terrain aspect. Therefore, based on multiple available datasets, we explored the performance of a proposed terrain slope estimation model over several study sites and various footprint shapes. The terrain slopes were derived from the ICESAT/GLAS waveform data by the proposed method and five other methods in this study. Compared with five other methods, the proposed method considered the influence of footprint shape, orientation, and terrain aspect on the terrain slope estimation. Validation against the airborne LiDAR measurements showed that the proposed method performed better than five other methods (R² = 0.829, increased by ~0.07, RMSE = 3.596°, reduced by ~0.6°, n = 858). In addition, more statistics indicated that the proposed method significantly improved the terrain slope estimation accuracy in high-relief region (RMSE = 5.180°, reduced by ~1.8°, n = 218) or in the footprint with a great eccentricity (RMSE = 3.421°, reduced by ~1.1°, n = 313). Therefore, from these experiments, we concluded that this terrain slope estimation approach was beneficial for different terrains and various footprint shapes in practice and the improvement of estimated accuracy was distinctly related with the terrain slope and footprint eccentricity.