Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) shows great potential for forest height retrieval. However, there are abundant noise photons in the ICESat-2 data, which make the accurate extraction of global forest heights challenging. In this letter, a novel algorithm based on the clustering method of ordering points to identify the clustering structure (OPTICS) was proposed to remove noise photons. First, we modified the circular shape of the search area in the OPTICS algorithm to an elliptical shape. Second, a distance ordering of all photons was generated using the modified OPTICS algorithm. Finally, signal photons were effectively detected using distance thresholds set by the Otsu method. To evaluate the algorithm performance, both the simulated and real ICESat-2 data were applied to our proposed algorithm. In addition, we compared our algorithm with another noise removal algorithm based on the modified density-based spatial clustering of applications with noise (DBSCAN). The results show that our algorithm works well in distinguishing the signal and noise photons as indicated by high $F$ values. Compared with the modified DBSCAN, our algorithm performs better in filtering out noise photons regardless of the simulated or real ICESat-2 data sets. In addition, the results also indicate that our algorithm is robust because it is insensitive to the clustering parameters. Overall, the new proposed algorithm is effective for removing noise photons in the ICESat-2 data.

Accurate estimation of ground elevation on a large scale is essential and worthwhile in topography, geomorphology, and ecology. The Ice, Cloud and Land Elevation Satellite-2 (ICESat-2) mission, launched in September 2018, offers an opportunity to obtain global elevation data over the earth’s surface. This paper aimed to evaluate the performance of ICESat-2 data for ground elevation retrieval. To fulfill this objective, our study first tested the availability of existing noise removal and ground photon identification algorithms on ICESat-2 data. Second, the accuracy of ground elevation data retrieved from ICESat-2 data was validated using airborne LiDAR data. Finally, we explored the influence of various factors (e.g., the signal-to-noise ratio (SNR), slope, vegetation height and vegetation cover) on the estimation accuracy of ground elevation over forest, tundra and bare land areas in interior Alaska. The results indicate that the existing noise removal and ground photon identification algorithms for simulated ICESat-2 data also work well for ICESat-2 data. The overall mean difference and RMSE values between the ground elevations retrieved from the ICESat-2 data and the airborne LiDAR-derived ground elevations are −0.61 m and 1.96 m, respectively. In forest, tundra and bare land scenarios, the mean differences are −0.64 m, −0.61 m and −0.59 m, with RMSE values of 1.89 m, 2.05 m, and 1.76 m, respectively. By analyzing the influence of four error factors on the elevation accuracy, we found that the slope is the most important factor affecting the accuracy of ICESat-2 elevation data. The elevation errors increase rapidly with increasing slope, especially when the slope is greater than 20°. The elevation errors decrease with increasing SNR, but this decrease varies little once the SNR is greater than 10. In forest and tundra areas, the errors in the ground elevation also increase with increasing vegetation height and the amount of vegetation cover.

The upcoming space-borne LiDAR satellite Ice, Cloud and land Elevation Satellite-2 (ICESat-2) is scheduled to launch in 2018. Different from the waveform LiDAR system onboard the ICESat, ICESat-2 will use a micro-pulse photon-counting LiDAR system. Thus new data processing algorithms are required to retrieve vegetation canopy height from photon-counting LiDAR data. The objective of this paper is to develop and validate an automated approach for better estimating vegetation canopy height. The new proposed method consists of three key steps: 1) filtering out the noise photons by an effective noise removal algorithm based on localized statistical analysis; 2) separating ground returns from canopy returns using an iterative photon classification algorithm, and then determining ground surface; 3) generating canopy-top surface and calculating vegetation canopy height based on canopy-top and ground surfaces. This automatic vegetation height estimation approach was tested to the simulated ICESat-2 data produced from Sigma Space LiDAR data and Multiple Altimeter Beam Experimental LiDAR (MABEL) data, and the retrieved vegetation canopy heights were validated by canopy height models (CHMs) derived from airborne discrete-return LiDAR data. Results indicated that the estimated vegetation canopy heights have a relatively strong correlation with the reference vegetation heights derived from airborne discrete-return LiDAR data (R² and RMSE values ranging from 0.639 to 0.810 and 4.08 m to 4.56 m respectively). This means our new proposed approach is appropriate for retrieving vegetation canopy height from micro-pulse photon-counting LiDAR data.

The accurate estimation of crop biomass during the growing season is very important for crop growth monitoring and yield estimation. The objective of this paper was to explore the potential of hyperspectral and light detection and ranging (LiDAR) data for better estimation of the biomass of maize. First, we investigated the relationship between field-observed biomass with each metric, including vegetation indices (VIs) derived from hyperspectral data and LiDAR-derived metrics. Second, the partial least squares (PLS) regression was used to estimate the biomass of maize using VIs (only) and LiDAR-derived metrics (only), respectively. Third, the fusion of hyperspectral and LiDAR data was evaluated in estimating the biomass of maize. Finally, the biomass estimates were validated by a leave-one-out cross-validation (LOOCV) method. Results indicated that all VIs showed weak correlation with field-observed biomass and the highest correlation occurred when using the red edge-modified simple ratio index (ReMSR). Among all LiDAR-derived metrics, the strongest relationship was observed between coefficient of variation (${\mathrm{H}}_{CV}$of digital terrain model (DTM) normalized point elevations with field-observed biomass. The combination of VIs through PLS regression could not improve the biomass estimation accuracy of maize due to the high correlation between VIs. In contrast, the ${\mathrm{H}}_{CV}$ combined with ${\mathrm{H}}_{mean}$ performed better than one LiDAR-derived metric alone in biomass estimation (R² = 0.835, RMSE = 374.655 g/m², RMSE_CV = 393.573 g/m²). Additionally, our findings indicated that the fusion of hyperspectral and LiDAR data can provide better biomass estimates of maize (R² = 0.883, RMSE = 321.092 g/m², RMSE_CV = 337.653 g/m²) compared with LiDAR or hyperspectral data alone.

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.