A Multi-Level Auto-Adaptive Noise-Filtering Algorithm for Land ICESat-2 Photon-Counting Data

Due to atmospheric scattering, solar radiation, and other factors, the Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) product data suffer from a substantial amount of background noise. This poses a significant challenge when attempting to directly utilize the raw data. Consequently, data denoising becomes an indispensable preprocessing step for its subsequent applications, such as the extraction of forest structure parameters and ground elevation data. While the Density-Based Spatial Clustering of Applications with Noise (DBSCAN) algorithm is currently the most widely used method, it remains susceptible to complexities arising from terrain, low signal-to-noise ratio (SNR), and input parameter variations. This paper proposes an efficient Multi-Level Auto-Adaptive Noise Filter (MLANF) algorithm based on photon spatial density. Its purpose is to extract signal photons from ICESat-2 terrestrial data of different ground cover types. The algorithm follows a two-step process. Firstly, random noise photons are removed from the upper and lower regions of the signal photons through a coarse denoising process. Secondly, in the fine denoising step, the K-Nearest Neighbor (KNN) algorithm selects the K photons to calculate the slope along the track. The calculated slope is then used to rotate the direction of the searching neighborhood in the DBSCAN algorithm. The proposed algorithm was tested in eight datasets of four surface types: forest, grassland, desert, and urban, and the extraction results were compared with those from the ATL08 datasets and the DBSCAN algorithm. Based on the ground-truth signal photons obtained by visual inspection, the classification precision, recall, and F-score of our algorithm, as well as two other algorithms, were calculated. The MLANF could achieve a good balance between classification precision (97.48% averaged) and recall (97.96% averaged). Its F-score (97.69% averaged) was higher than that of the other two methods. This demonstrates that the MLANF algorithm successfully obtained a continuous surface profile from ICESat-2 datasets with different surface cover types, significant topographic relief, and low SNR.