Next Article in Journal
Deception Island 1967–1970 Volcano Eruptions from Historical Aerial Frames and Satellite Imagery (Antarctic Peninsula)
Next Article in Special Issue
KaRIn Noise Reduction Using a Convolutional Neural Network for the SWOT Ocean Products
Previous Article in Journal
Global Mangrove Watch: Monthly Alerts of Mangrove Loss for Africa
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Robust Algorithm for Photon Denoising and Bathymetric Estimation Based on ICESat-2 Data

1
Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China
2
School of Geography and Information Engineering, China University of Geosciences, Wuhan 430074, China
3
College of Earth and Planetary Science, University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(8), 2051; https://doi.org/10.3390/rs15082051
Submission received: 25 February 2023 / Revised: 30 March 2023 / Accepted: 4 April 2023 / Published: 13 April 2023
(This article belongs to the Special Issue Applications of Satellite Altimetry in Ocean Observation)

Abstract

:
The Ice, Cloud, and Land Elevation Satellite 2 (ICESat-2) is equipped with an Advanced Terrain Laser Altimeter System (ATLAS) with the capability of penetrating water bodies, making it a widely utilized tool for the bathymetry of various aquatic environments. However, the laser sensor often encounters a significant number of noise photons due to various factors such as sunlight, water quality, and after-pulse effect. These noise photons significantly compromise the accuracy of bathymetry measurements. In an effort to address this issue, this study proposes a two-step method for photon denoising by utilizing a method combining the DBSCAN algorithm and a two-dimensional window filter, achieving an F1 score of 0.94. A robust M-estimation method was employed to estimate the water depth of the denoised and refraction-corrected bathymetric photons, achieving an RMSE of 0.30 m. The method proposed in this paper preserves as much information as possible about signal photons, increases the number of bathymetric points, enhances the resistance to gross error, and guarantees the accuracy of bathymetry measurements while outlining the underwater topography. While the method is not fully automated and requires setting parameters, the fixed parameter values allow for efficient batch denoising of underwater photon points in different environments.

1. Introduction

NASA’s new-generation Ice, Cloud, and Land Elevation Satellite 2 (ICESat-2) laser altimetry mission, launched in 2018, aims to precisely measure changes in ice sheets and sea ice and furnish a deeper understanding of their dynamics [1,2]. It carries a single-micropulse, photon-counting laser named the Advanced Terrain Laser Altimetry System (ATLAS), which works at a frequency of 10 kHz and is combined with the diffraction optical system to collect six independent profiles and six beams from 88°N to 88°S along the orbit over a repeated period of 91 days. The six beams are divided into three pairs in the orbit; each pair has a strong beam and a weak beam, and the intensity ratio of the strong beam to the weak beam is 4:1. ICESat-2’s ATLAS laser pulses 10,000 times per second, firing 300 trillion bright green photons toward the ground [3,4]. The length of time it takes for a single laser photon to bounce off Earth’s surface and travel back to the satellite is used by scientists to measure altitude [4].
With an advancement in research and the application of ICESat-2, the role of ICESat-2 in coastal areas is becoming increasingly important. For example, underwater topography in coastal areas is crucial for understanding coastal climate change and its impact on the nearshore environment [5]. In addition, ICESat-2 satellite data have great potential in coastal zone classification. Therefore, accurate bathymetric data are crucial for marine and coastal ecosystems and related studies [6,7]. The ATLAS carried on ICESat-2 adopts photon-computing lidar at the 532 nm band, which has a specific penetration ability for water bodies [8]. The ATLAS system employs photomultiplier tubes (PMTs) as detectors in photon-counting modes. This allows for detection within the ICESat-2 data acquisition system upon the reflection of a single photon. The ICESat-2 laser emits pulses with an interval of only 1.5 nanoseconds, which may broaden during transmission, but typically results in received pulses of a few nanoseconds. Each laser emission, therefore, results in returned laser signal photons from both the surface and underwater bottom concentrated within a limited time duration [9]. These underwater photons may be affected by factors such as solar background noise, water scattering, and the after-pulse effects of lidar. After-pulse effects occur when a lidar pulse interacts with a surface and causes the emission of secondary photons. These secondary photons can be detected by the lidar receiver as if they are from the original pulse, and multiple surface echoes caused by after-pulse effects are typically seen in photons from very smooth open water surfaces [10]. The presence of noise photons within the original photon data received by the detector necessitates the crucial step of point cloud preprocessing, which aims to remove noise and extract target point clouds from the two-dimensional photon point cloud profile along the orbit [11]. Point cloud denoising, a crucial aspect of photon-counting laser point cloud processing, is based on the principle that noise photons, which are generated by solar radiation and scattering, are generally distributed loosely and have different distribution characteristics compared to signal photons [12]. Through classification based on these distinct distribution features, it is possible to eliminate noise photons and attain the most accurate bathymetric estimates.
For coastal areas, various algorithms have been developed to extract signal photons from the original photon data, and among which the density estimation method was widely used [13,14]. This method usually extracts signal photons according to the distribution characteristics and density difference between signal and noise photons. For example, since the Density-Based Spatial Clustering of Applications with Noise (DBSCAN) [15] algorithm was proposed, it has been widely used in denoising the data of various point cloud types. Duan et al. [16] proposed an LDBSCAN based on the DBSCAN algorithm to reduce the difficulty in parameter selection. However, the signal photons of ICESat-2 usually present a strip-like point cluster distribution, and the circle filtering method is challenging to adapt to the distribution characteristics of signal photons [17]. Zhang et al. [18] proposed the method of a ellipse-based DBSCAN algorithm, which can automatically adjust the size of the ellipse according to the density around the photon and calculate the ellipse, but the algorithm needs to try different combinations of parameters to enhance the denoising effect [19]. Magruder et al. [20] proposed a denoising algorithm based on raster image processing, which requires optimal parameter selection to ensure acceptable noise and signal photon filtering, and needs to rasterize the point cloud, which will cause information loss and reduce the accuracy of filtering [21]. Based on the above analysis, it can be found that the performance of the above two methods depends mainly on the input parameters. Therefore, Zhu et al. [22] proposed an OPTICS method based on an elliptic region photon search, which greatly reduced the sensitivity of the input parameters and improved the algorithm’s self-adaptation. Although the OPTICS algorithm was insensitive to parameters, the accurate extraction of underwater signal photons also requires the secondary filtering of OPTICS-filtered photons, which cannot be automated and therefore still requires parameter adjustment for different data sets [23]. Xie et al. [24] proposed a self-adaptive fine denoising approach that calculates the density value and employs a self-adaptive elliptical LDBSCAN algorithm to effectively eliminate noise around the signal segment. Although these methods perform well in forest height retrieval, there is still a lack of relevant experimental evidence showing that these methods are directly applicable to underwater sparse signal photons and underwater terrains of different complexities. Whether they can be used in large areas, especially complex underwater terrain areas, needs to be further tested [14,17]. Heidi Ranndal et al. [25] proposed a sliding window method to obtain the median value of photon elevation, which is an improvement in the two-dimensional window filter [26]. According to the length of the window set by the ICESat-2 spot size, the photon with the median elevation in the window was selected as the window’s center point, and the window’s width was set to slide the window in order to screen the bathymetric photons. This method can adapt to the distribution features of the signal photons and is robust at the same time; so, it can effectively filter out noise photons to improve the precision of bathymetry. However, this method needs to ensure that the number of signal photons in the window is larger than the number of noise photons, so that the selected median elevation can be the elevation of the signal photons. Due to the complex underwater environment, the laser passes through the water column and is absorbed and scattered by the water column. The photons are very sparse under water, so the signal photon density reflected by the electromagnetic wave that penetrates the water column to reach the bottom of the water is sparse [17]. Especially in the daytime, it is affected by the background noise of the sun, and the density of the noise photons is large; so, the number of signal photons in a window is not guaranteed to be larger than the number of noise photons. In addition, the denoised and corrected signal photons have a certain thickness after being affected by rescattering, therefore depth estimation is needed to obtain accurate bathymetric values. In summary, most previous studies directly used photons as the bathymetric points after photon denoising and refraction correction without further bathymetric estimation.
When calculating depth, the traditional least squares algorithm often does not have the ability of gross error interference resistance. The original outlier detection method mainly uses the hypothesis test theory to construct test statistics between the null hypothesis and the alternative hypothesis, only to make decisions and identify outliers [5]. Robust estimation can effectively detect and eliminate the outliers; at the same time, the equivalent weight function is used to add weak weights for small- and medium-sized outliers, which can more accurately estimate the model parameters [27,28,29,30,31,32,33].
In this paper, we combine the DBSCAN algorithm with a two-dimensional window filter to achieve a two-step method for photon denoising. Firstly, to ensure the maximal retention of signal photon points during the denoising process, a “loose” setting for the DBSCAN denoising parameters was adopted, which only needed to ensure that the number of signal photon points was larger than the number of noise photon points in a step range for the denoising in the second step. Given that the density of signal photon points typically exceeds that of noise photon points, the parameter setting of the first step was straightforward. Secondly, preprocessed bathymetric photons were used for further denoising using the two-dimensional window filter. The parameters of the two-dimensional window filter in the second step were based on the thickness of the bathymetric signal photons. Following denoising, refraction correction was performed on the bathymetric photons, and we proposed a depth value estimation method for ICESat-2 bathymetric protons using robust M-estimation in order to obtain bathymetric valuations at various points along the track. In addition, in order to evaluate the accuracy of the two-step denoising method and the robust M-estimation algorithm, the original underwater photon data were manually denoised, and the manually denoised bathymetric photon point data were used as reference values. Finally, our proposed method was tested on both manual denoising data and in situ CZMIL data, demonstrating its effectiveness across a range of environmental conditions. The accuracy of these methods was evaluated to assess their denoising performance and bathymetric performance. This paper compares several denoising methods, including the two-step method, which retains as many signal photons as possible while achieving a better filtering effect. The findings of our experimental analysis indicate that the utilization of robust M-estimation in comparison to solely relying on denoised photon points as bathymetric values yields a noteworthy enhancement in bathymetric accuracy, and enables a more accurate inversion of underwater topography. In this study, all of the parameters of the proposed methods were fixed and were not adjusted for environmental factors (e.g., sunlight intensity and complexity of the seafloor topography) to enable the batch processing of the underwater photon data.

2. Study Areas and Data

In this study, we choose the northwest coast of Hawaii (155.91°W, 19.90°N) as the study area, and used the ICESat-2 photon altimetry data from the official NASA Satellite Data website (https://search.earthdata.nasa.gov, accessed on 28 October 2021). The ICESat-2 altimetry data of two periods along the coast of Hawaii Island were processed. The study area is shown in Figure 1, and the lines in the figure are photon tracks of the two orbits.
The data acquisition time and corresponding terrain for the two orbits are shown in Table 1. The profile of the row data in the region is shown in Figure 2. To compare the denoising effect, the day and night flat and undulating areas were selected, respectively, for denoising comparison.

3. Methods

3.1. Two-Step Method

Here, we used the combination of DBSCAN density clustering and a two-dimensional window filter to classify signal photons and noise photons and then denoise them. DBSCAN is a clustering method based on machine learning [34]. To search for signal photons, the DBSCAN algorithm uses a circular search area with a radius specified by the variable “eps” (epsilon) and a minimum number of points, “MinPts”, which should be included within the circle. To identify signal photons, we used the DBSCAN algorithm, which searches for clusters of photon points in the lidar signal data. As shown in Figure 3, MinPts = 3. The algorithm starts by randomly selecting a photon point, “A”, as the center of a circular search area with a radius of “eps”. Point A and other blue points are core points because their neighborhood (blue circles in the figure) contains at least 3 points (including themselves). Yellow photon points “B” and “C” fall within the neighborhood of a core point; so, they are also included as members of that core point’s cluster. Points “B” and “C” are not core points but fall within the neighborhood of a core point; so, they are considered to be boundary points. The red point “N” is not a core point and does not fall within the neighborhood of any core point; so, it is considered to be a noise point. The DBSCAN algorithm iteratively identifies core points and boundary points until all such points have been identified. The points identified as core points or boundary points are considered to be signal photon points. The points identified as noise points are rejected. In the ATL03 observation data of ICESat-2, the density of signal photons returned from the surface was relatively large; so, DBSCAN was suitable for processing the large number of photon points returned from the surface.
DBSCAN denoising can achieve good results for photon information with good data quality returned from simple topography on Earth’s surface, but the denoising effect is not ideal when the noise point density increases. Since the DBSCAN domain search is conducted according to the circular region of a clustering method, the group is suitable for the classification of point clusters. However, ICESat-2 data acquisition is a strip scan, and the collected signal photon points form strip point clusters; therefore, to improve this algorithm to apply it to the photon data of ATL03 products, it is necessary to make an improvement based on the DBSCAN method by using a two-dimensional window filter. The parameters of DBSCAN are fixed to eps = 6 and Minpts = 3 for primary denoising. The median method has strong robustness; therefore, given the photon numerical characteristics after denoising via DBSCAN, the two-dimensional window filter is used for secondary denoising to achieve a better denoising effect.
  • Since the size of the ICESat-2 spot is about 17 m, set a window 17 m long with a sliding step of 17 m to slide along the underwater photon track;
  • In this window, arrange the photon elevation values within the 17 m window to find the median value, H m , and set the half-height h of the window to 0.7 m (this is an empirical value);
  • When the photon elevation H in the window meets H m h H H m + h , the photon is classified as a signal photon;
  • Slide to the next window and repeat Step 3.

3.2. Robust M-Estimation

Robust estimation is divided into L-estimation, R-estimation, and M-estimation, and the most widely used in the field of surveying and mapping is robust M-estimation, i.e., the maximum likelihood estimate. Robust estimation is the core of the steps in the right-to-choose iteration process [35,36]. In short, the optimal valuation of the parameters is calculated by continuously reducing the weight ratio of the outliers through weight selection iterations [37,38]. Most previous studies directly used the photons as the bathymetric points after photon denoising and refraction correction without further bathymetric estimation. Therefore, this paper uses the robust M-estimation to improve the accuracy of bathymetry and the number of bathymetric points. Each photon emitted by ATLAS can be regarded as an independent observation value. Therefore, the weight ratio of noise photons can be reduced by the weight selection iteration method of independent observation with equal weight to estimate the underwater depth parameters.
Assuming that the observation vector is L and the unknown parameter vector is β, the Gauss–Markov model is constructed, and the error equation can be expressed as follows [33].
V = A β L D L = σ 0 2 P 1
where V is the residual vector, A is the coefficient matrix, L is the observation vector, P is the weight of the observation, and σ 0 2 is the unknown variance component of the unit weight.
Combining Equation (1) with M-estimation, the extreme value function can be obtained as follows:
i = 1 n p i ρ V i = i = 1 n p i ρ a i β ^ L i = m i n
where a i is the ith component of A , L i is the ith component of L , β ^ is the estimated value of the true parameter value β, and ρ V i is the weight function. To determine the partial derivative of β ^ and assign it a value of 0, ψ ( V i ) = ρ / V i is inserted into Equation (2) to obtain Equation (3):
i = 1 n p i ψ ( V i ) a i = 0
Let ψ V i V i = ω i and p i = p i ω i ; p i is the equivalent weight factor. Then, the following equation is obtained:
A T P V = 0
P = d i a g ( p 1 , p 2 , p 3 , . . , p n ) is the equivalent weight array and n is the number of observations. It can be obtained from Equation (1):
A T P A β ^ = A T P L
Therefore, robust estimation solutions of model parameters can be obtained as follows:
β ^ = ( A T P A ) 1 A T P L
The difference between the above model and the least squares estimation is that the weight matrix of the least squares estimation is replaced by the equivalent weight matrix ( P ) shown in Equation (6), which is generally solved via iterative calculation. The iterative solution of the k + 1 step can be expressed as follows.
β ^ k + 1 = ( A T P k A ) 1 A T P k L
Therefore, usually, a ρ function defines the M-estimation; the determination of the ρ function with a set of observations of variance factor σ 0 , obtained through the ρ function, calculates the equivalent weights of observations.
Considering the prior measurement, the IGG3 [30,32,33,39] scheme constructs a robust equivalent weight function using the following formula:
P i = P i V ~ i k 0 P i k 0 V ~ i ( k 1 V ~ i k 1 k 0 ) 2 k 0 < V ~ i k 1 0 V ~ i > k 1
where V ~ i is the ith component of the standardized residuals based on the residuals V , and k 0 and k 1 are constants ( k 0 = 1.0 1.5 , k 1 = 2.5 3.0 ). The choice in these two constants is an empirical model determined according to the size of the data. The error in the post-test unit weight was used as the variance factor:
σ 0 = V T P V n t
σ 0 k = V T P V n , t
where n is the number of observations, n , is the number of P i = 0 corresponding to the elimination of the detected outliers, t is the necessary number of observations, P is the weight of the original observation, and k is the number of iterative calculations.
Since the ICESat-2 spot size is 17 m and each sequential ICESat-2 laser pulse is only 0.7 m away from the preceding pulse, set the estimated window size length to 17 m, the sliding step to 0.7 m, and the median elevation of the depth finder photon for each window to the initial value β ^ . As shown in Figure 4, the green point can be regarded as the bathymetric photon point of each window, the red point is the initial bathymetric valuation point, the solid black line is the contour line where the initial point is located, and the length of the blue dashed line is the vertical distance, V , between one of the bathymetric photon points and the initial point. The specific calculation process is as follows:
  • Calculate the initial value β ^ , residual V , and variance factor σ 0 of the parameter estimate based on the least squares estimation of Equations (1) and (9);
  • Calculate the corresponding standardized residuals, V ~ = V / σ 0 , from P i = P i ρ , using Equation (8) IGG3, and calculate the corresponding equivalent weights P ;
  • The equivalent weight is obtained according to the above steps, and the parameters β ^ 1 , variance factor σ 0 1 , and residual V of the first iteration model are calculated in combination with Equations (6) and (9b);
  • The residual V and variance factor σ 0 , obtained in Step 3, are substituted into Equation (8) to calculate the model parameter β ^ 2 of the second iteration;
  • Repeat Steps 3 and 4 until meeting the conditions β ^ k + 1 β ^ k < 10 6 . The final estimated model parameters can be calculated using Equation (7).

3.3. Manual Denoising

There are no direct real results for validation against photon point cloud data; so, this paper conducts an accuracy evaluation of the two-step denoising and robust M-estimation techniques, using manually denoised bathymetric photon points as reference values. To achieve this, interactive mapping software capable of importing and exporting photon point geographic information was required. The southern CASS system, developed on the AutoCAD platform, was utilized for this purpose. The coordinate information of the underwater photon row points, extracted after classification, was imported into CASS. To facilitate visual interpretation, the latitude and elevation of the photon points were designated as the horizontal and vertical coordinates, respectively, while the aspect ratio of the coordinates was adjusted and the horizontal coordinate was enlarged 10,000 times. Subsequently, manual interaction and visual interpretation were used to remove noise photon points, allowing for the extraction and exportation of signal photon point data. Although this method has a good denoising effect, it is inefficient and cannot remove the noise photons of all data in a batch.

4. Results

4.1. Photon Denoising Based on DBSCAN

To accurately determine the water depth along the flight path of ICESat-2, it is important to distinguish between surface and bottom photons, as underwater photons are affected by the refraction effect of the water column. Thus, relative photon classification is required before photon denoising. As shown in Figure 5, the classified photon points are the brown dots above the water, the blue ones are on the water, and the green ones are underwater, and the red dashed line is the dividing line between the surface photons and the underwater photons. Considering ICESat-2’s limited bathymetric capability and its restriction to shallow water measurements, the dividing line can be considered to be approximately horizontal at short distances offshore. The determination of the dividing line is based on the calculation of d = lm−Δ, where lm represents the local mean sea level, and Δ equals 3σ or 2σ, with σ being the root mean square wave height. The lm and σ are calculated by the mean and standard deviation of the detected sea surface photon points. In this paper, Δ is taken as 3σ, to fully distinguish the surface photons from the bottom photons [6].
This paper set different parameters, and DBSCAN was used for denoising. Due to the limitations of DBSCAN and the sparse density of underwater signal photons, when the parameter setting is “loose” (eps = 6, MintPts = 3), noise photons cannot be removed, and many noise points are retained. When the parameter setting is “strict” (eps = 3, MintPts = 5), some signal photons are removed simultaneously. Figure 6 shows the denoising contrast diagram of 1036gt1r of the DBSCAN denoising method with different parameter conditions. It is apparent that denoising leads to a loss in signal photons and causes discontinuity in the underwater terrain. These two parameters were adjusted to achieve the “best” denoising effect; Figure 7 shows the underwater photon denoising map of each beam via the DBSCAN denoising method after parameter adjustment.

4.2. Photon Denoising Based on the Two-Dimensional Window Filter

Since the spot size of ICESat-2 is 17 m and a maximum of 15 geolocation photons can be observed per meter in water depth less than 1 m [25], the length of the window was set to 17 m, the photon elevation values within the 17 m were arranged to take the median value, and the photon position of the median elevation value was taken as the center point of the window. The window height was set to 1.4 m (an empirical value), and the window was slid away from the coastline. The photon points passing through the window were regarded as signal photon points after denoising, as shown in Figure 8.

4.3. Photon Denoising Based on the Two-Step Method

In the case of the DBSCAN algorithm, to maximize the retention of signal photons, a “loose” parameter setting was used to perform the initial denoising step. However, to further enhance the denoising effect, a two-dimensional window filter was applied to the denoised photon data obtained from the previous step. The filter’s parameters were set to fixed values of eps = 6 and MinPts = 3, which eliminated the need for parameter adjustments and ensured consistent parameters across all data sets.
The two-step denoising method resulted in a significant improvement in the quality of the photon data, as demonstrated by the denoising effect shown in Figure 9. The combination of the DBSCAN algorithm with a two-dimensional window filter effectively removed the noise present in the original data set while retaining as many signal photons as possible, resulting in a clear and accurate representation of the underlying signal.

4.4. Draw Bathymetry Profiles

As seen in Figure 10, after photon denoising, the refraction correction of the de-noised underwater photons was carried out using the method proposed by Parrish et al. [40]. The blue points are the surface photon points, the brown points are the coastal photon points, the red points are the underwater bathymetry photon points before refraction correction, and the green points are the underwater bathymetry photon points after refraction correction. Combined with this robust M-estimation algorithm, the underwater bathymetry profile of each orbit was obtained, as shown in Figure 11.

4.5. Accuracy Assessment

Using the manual denoised underwater photon dots as real signal photon points, this paper compares and verifies the accuracy of the signal photon points after manual denoising with those after denoising via each denoising method. The effect of various denoising methods on underwater photon denoising was evaluated using the following four metrics.
A c c u r a c y = T P + T N T P + F P + T N + F N
P r e c i s i o n = T P T P + F P
R e c a l l = T P T P + F N
F 1 s c o r e = 2 * P r e c i s i o n * R e c a l l P r e c i s i o n + R e c a l l
TP (true positive) is the number of photons judged to be signal photons which are in fact signal photons, FP (false positive) is the number of photons judged to be signal photons which are in fact noise photons, TN (true negative) is the number of photons judged to be noise photons which are in fact noise photons, and FN (false negative) is the number of photons judged to be noise photons which are in fact signal photons. The metrics of each denoising method are shown in Table 2:
In order to verify the accuracy of robust M-estimation, the estimated values were verified by the photon bathymetry points after manual denoising, and since the manual denoising was based on the denoising of the row data without refraction correction for the underwater photons, the bathymetry valuation here is the result of the robust M-estimation without refraction correction, as shown in Figure 12. The red line is the y = x line, the blue line is the regression line, and N is the number of bathymetry points via M-estimation. The bathymetry points obtained via M-estimation were matched with the nearest bathymetry photon points, and scatter plots of four regions were obtained, which reflect the correlation between the M-estimated bathymetry points and the manually denoised photon points.
In order to ensure the reliability of the estimated bathymetric values, we used in situ data collected on the Big Island of Hawaii by the Coastal Zone Mapping and Imaging Lidar (CZMIL) system provided by NOAA (https://coast.noaa.gov/dataviewer/, accessed on 8 November 2021). The accuracy of the denoised bathymetric photon points and the accuracy of the M-estimated bathymetric points were verified separately. The CZMIL vertical positional accuracy of the shallow bathymetry is 0.25 2 + 0.0075 2 d meters at a 95% confidence level, and d is the depth. The CZMIL in situ data horizontal positions were referenced from the North American Datum of 1983 National Adjustment of 2011 (NAD83 PA11), and vertical positions were referenced from the (NAD83 PA11) ellipsoid. The photon event data in ATL03 are presented in geodetic spherical coordinates as geodetic latitude, longitude, and photon event height above the WGS-84 ellipsoid (or ellipsoidal height) [41]. The NAD83 PA11 ellipsoid heights of the CZMIL points are converted to WGS-84 ellipsoid heights via the Vdatum tool (https://vdatum.noaa.gov/vdatumweb/, accessed on 29 November 2022). In addition, to compensate for the tidal effect, the difference between the regional mean seawater surface lm from ICESat-2 data and the WGS-84 ellipsoidal heights of each bathymetric datum were used in the validation process. Since there are no artificial structures in the area, the topographic variation between the CZMIL and ICESat-2 bathymetric data of the scanned area are not considered, due to the large time span. As shown in Figure 13, the bathymetry of both flat and undulating areas was well estimated under different degrees of solar interference during the day and night, and the bathymetry estimation points were increased under the premise of ensuring the bathymetric accuracy.

5. Discussion

The results of the three compared denoising algorithms, DBSCAN, the two-dimensional window filter, and the two-step method are shown in Table 2. The performance metrics include the F1 score, accuracy, precision, and recall. Based on these metrics, the two-step method shows the best performance, as it has the highest F1 score, accuracy, and precision, and a relatively high recall. However, it is worth noting that the best algorithm for denoising data depends on the specific requirements and characteristics of the data being processed. The two-step method appears to have the highest F1 score compared to DBSCAN or the two-dimensional window filter alone, in both flat and undulating terrain, both during the night and day. In flat terrain during the night, the two-step method had an F1 score of 0.9435, which is an improvement of 4.8% and 5.9% compared to using DBSCAN and the two-dimensional window filter alone, respectively. In undulating terrain during the night, the two-step method had an F1 score of 0.8343, which is an improvement of 15.1% and 12.3% compared to using DBSCAN or the two-dimensional window filter alone, respectively. In flat terrain during the day, the two-step method had an F1 score of 0.7038, which is an improvement of 4.5% compared to using DBSCAN or the two-dimensional window filter alone. In undulating terrain during the day, the two-step method had an F1 score of 0.7780, which is an improvement of 25.7% and 32.1% compared to using DBSCAN or the two-dimensional window filter alone, respectively. It can be seen that DBSCAN has a better denoising effect only when processing data from flat areas at night, which results in good data quality, but the denoising effect is not good when dealing with other types of data. In conclusion, the combination of DBSCAN and the two-dimensional window filter method appears to have the greatest advantage in terms of F1 score, achieving the highest F1 score compared to DBSCAN or the two-dimensional window filter alone, in both flat and undulating terrain, both during the night and day.
If only the denoising photon point is used as the bathymetric estimation point, the number of bathymetric points obtained is very small, and the noise photon points of the uneven area data are difficult to remove; so, this paper uses M-estimation to further process the denoising photon point.
Figure 12 shows the accuracy evaluation of the bathymetric M-estimation points with manual denoising points; the RMSE values range from 0.4 to 1.2 and the R2 values range from 0.91 to 0.98. Figure 13 provides performance information for depth estimation before and after M-estimation in different conditions. The conditions include the time of day, the terrain type, the algorithm used, and the number of bathymetric points obtained. As can be seen from the table, in all cases, using the “M-estimation” algorithm is superior to using only the “two-step method” in terms of RMSE and R2. For each depth estimation algorithm, RMSE is generally lower for “flat” terrain data compared to “undulating” terrain data, and low for night data compared to day data.
The RMSE values for the two-step method range from 0.39 to 0.87, while the RMSE values for the M-estimation algorithm range from 0.30 to 0.79. This indicates that the use of the M-estimation algorithm provides more accurate depth estimation results. The M-estimation algorithm also performs better than the two-step method in terms of R2. The R2 values for the M-estimation algorithm range from 0.89 to 0.98, while the R2 values for the two-step method range from 0.86 to 0.96. A higher R2 value indicates a better fit of the data to the model, which suggests that the M-estimation algorithm is better at modeling the relationship between the bathymetric points and the depth values. The M-estimation algorithm also substantially increases the number of bathymetric points. Therefore, the analysis shows that a single denoising method often fails to remove the noise photons cleanly in the case of high solar background noise and complex underwater terrain, or requires multiple parameter adjustments for different data sets to achieve better denoising effects. The DBSCAN denoising method requires constant adjustment of the parameters to achieve the best denoising effect, and is not suitable for strip cluster denoising. The two-dimensional window filter needs to ensure that the number of signal photons in a window is larger than the number of noise photons, which is not suitable for environments with high noise point density. Therefore, after the first step of rough denoising, a second step of fine denoising is needed. In the case of sparse underwater noise photon points, the two-step method tends to avoid too many adjustments of parameters and can achieve better denoising effects.

6. Conclusions

Regarding photon denoising, a two-step method is proposed to address the poor denoising performance of DBSCAN in complex terrains and the need for constant parameter adjustments. In addition, the high noise density of ATL03 data, caused by increased solar background noise during the day, presents challenges for DBSCAN and the two-dimensional window filter in eliminating noisy photons effectively. As shown in Table 1, the DBSCAN night method demonstrates good denoising performance in processing flat terrain data. However, when dealing with complex terrain or noise photons with increasing data, the two-dimensional window filter provides greater robustness but requires ensuring that the number of signal photons in a window is greater than the number of noise photons. When the density of noise points in a window exceeds that of signal points, the denoising performance deteriorates, particularly when processing data affected by solar background noise during the day. To address this issue, this study combined the DBSCAN algorithm with a two-dimensional window filter to retain more signal photons while eliminating noise photons and fixing the parameters of DBSCAN. This method overcomes the limitation of DBSCAN’s poor denoising performance when handling large noise photon density and avoids the need for constant parameter adjustments.
In addition, this study employed robust M-estimation in calculating bathymetry, which effectively estimates the depth of underwater topography while demonstrating resistance to high levels of noise photons in the ATLAS of underwater laser bathymetry application. This approach improves depth estimation performance and facilitates the acquisition of a greater number of bathymetry points, resulting in a clearer representation of the underwater topography. This enhancement in data density is unattainable only using the photon-denoising algorithm. Although the proposed approach is not fully automated and requires setting parameters, the parameters used in this study are fixed and not adjusted for environmental factors, and they still achieve good denoising and bathymetric estimation, avoiding too much parameter adjustment and allowing the batch processing of underwater photon data, improving the efficiency. Moreover, this algorithm can be applied to various types of laser point cloud data, particularly in scenarios involving significant outlier interference and requiring high-precision monitoring, offering broad potential applications. Since laser energy is substantially weakened when passing through water, especially in environments with poor water quality, the signal photons underwater are sparser, and both denoising and bathymetric estimation provide good results. Therefore, the approach also has broad application potential in other areas, and better results may be achieved on Earth’s surfaces with better data quality, such as ground, water, and glaciers. However, further research is needed to investigate the performance of the method in different regions and under different environmental conditions.

Author Contributions

Conceptualization, X.L. and J.Z.; methodology, J.Z.; validation, X.S. and X.L.; writing—original draft preparation, J.Z.; writing—review and editing, L.J.; visualization, J.Z.; supervision, X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key R&D Program Projects of Hubei Province (grant no. 2021BID011) and the National Key R&D Program of China (grant no. 2018YFC1406102).

Data Availability Statement

The open data used in this article, the access to the data is listed in the article, ICESat-2 data can be obtained by accessing website https://search.earthdata.nasa.gov (accessed on 25 February 2023), and CZMIL data can be obtained by accessing website https://coast.noaa.gov/dataviewer/ (accessed on 25 February 2023).

Acknowledgments

We sincerely thank the Goddard Space Flight Center for distributing the ICESat-2 data and NOAA for distributing the bathymetric data of the CZMIL.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abdalati, W.; Zwally, H.J.; Bindschadler, R.; Csatho, B.; Farrell, S.L.; Fricker, H.A.; Harding, D.; Kwok, R.; Lefsky, M.; Markus, T.; et al. The ICESat-2 Laser Altimetry Mission. Proc. IEEE 2010, 98, 735–751. [Google Scholar] [CrossRef]
  2. Markus, T.; Neumann, T.; Martino, A.; Abdalati, W.; Brunt, K.; Csatho, B.; Farrell, S.; Fricker, H.; Gardner, A.; Harding, D.; et al. The Ice, Cloud, and land Elevation Satellite-2 (ICESat-2): Science requirements, concept, and implementation. Remote Sens. Environ. 2017, 190, 260–273. [Google Scholar] [CrossRef]
  3. Hagopian, J.; Bolcar, M.; Chambers, J.; Crane, A.; Eegholm, B.; Evans, T.; Hetherington, S.; Mentzell, E.; Thompson, P.L.; Ramos-Izquierdo, L.; et al. Advanced topographic laser altimeter system (ATLAS) receiver telescope assembly (RTA) and transmitter alignment and test. In Proceedings of the Earth Observing Systems XXI, San Diego, CA, USA, 28 August–1 September 2016; pp. 55–67. [Google Scholar] [CrossRef]
  4. Neumann, T.A.; Martino, A.J.; Markus, T.; Bae, S.; Bock, M.R.; Brenner, A.C.; Brunt, K.M.; Cavanaugh, J.; Fernandes, S.T.; Hancock, D.W.; et al. The Ice, Cloud, and Land Elevation Satellite—2 mission: A global geolocated photon product derived from the Advanced Topographic Laser Altimeter System. Remote Sens. Environ. 2019, 233, 111325. [Google Scholar] [CrossRef] [PubMed]
  5. Urban, T.J.; Schutz, B.E.; Neuenschwander, A.L. A Survey of ICESat Coastal Altimetry Applications: Continental Coast, Open Ocean Island, and Inland River. Terr. Atmos. Ocean. Sci. 2008, 19, 1–19. [Google Scholar] [CrossRef] [Green Version]
  6. Ma, Y.; Xu, N.; Liu, Z.; Yang, B.; Yang, F.; Wang, X.H.; Li, S. Satellite-derived bathymetry using the ICESat-2 lidar and Sentinel-2 imagery datasets. Remote Sens. Environ. 2020, 250, 112047. [Google Scholar] [CrossRef]
  7. Hsu, H.-J.; Huang, C.-Y.; Jasinski, M.; Li, Y.; Gao, H.; Yamanokuchi, T.; Wang, C.-G.; Chang, T.-M.; Ren, H.; Kuo, C.-Y.; et al. A semi-empirical scheme for bathymetric mapping in shallow water by ICESat-2 and Sentinel-2: A case study in the South China Sea. ISPRS J. Photogramm. Remote Sens. 2021, 178, 1–19. [Google Scholar] [CrossRef]
  8. Nguyen, V.-A.; Ren, H.; Huang, C.-Y.; Tseng, K.-H. Bathymetry derivation in shallow water of the South China Sea with ICESat-2 and Sentinel-2 data. J. Appl. Remote Sens. 2021, 15, 044513. [Google Scholar] [CrossRef]
  9. Jasinski, M.F.; Stoll, J.D.; Cook, W.B.; Ondrusek, M.; Stengel, E.; Brunt, K. Inland and Near-Shore Water Profiles Derived from the High-Altitude Multiple Altimeter Beam Experimental Lidar (MABEL). J. Coast. Res. 2016, 76, 44–55. [Google Scholar] [CrossRef]
  10. Lu, X.; Hu, Y.; Yang, Y.; Vaughan, M.; Palm, S.; Trepte, C.; Omar, A.; Lucker, P.; Baize, R. Enabling value added scientific applications of ICESat-2 data with effective removal of afterpulses. Earth Space Sci. 2021, 8, e2021EA001729. [Google Scholar] [CrossRef]
  11. Cao, B.; Fang, Y.; Gao, L.; Hu, H.; Jiang, Z.; Sun, B.; Lou, L. An active-passive fusion strategy and accuracy evaluation for shallow water bathymetry based on ICESat-2 ATLAS laser point cloud and satellite remote sensing imagery. Int. J. Remote Sens. 2021, 42, 2783–2806. [Google Scholar] [CrossRef]
  12. Zhang, Z.; Liu, X.; Ma, Y.; Xu, N.; Zhang, W.; Li, S. Signal Photon Extraction Method for Weak Beam Data of ICESat-2 Using Information Provided by Strong Beam Data in Mountainous Areas. Remote Sens. 2021, 13, 863. [Google Scholar] [CrossRef]
  13. Yujia, L.; Xiaoqing, Z.; Guoyuan, L.; Jinquan, G.; Yue, M.; Yifu, C. Progress and prospect of space-borne photon-counting lidar shallow water bathymetry technology. Infrared Laser Eng. 2022, 51, 20220003. [Google Scholar]
  14. Xiaoxiao, Z.; Cheng, W.; Xiaohuan, X.; Sheng, N.; Xuebo, Y.; Dong, L. Research progress of ICESat-2/ATLAS data processing and applications. Infrared Laser Eng. 2020, 49, 20200259. [Google Scholar] [CrossRef]
  15. Ester, M.; Kriegel, H.-P.; Sander, J.; Xu, X. A density-based algorithm for discovering clusters in large spatial databases with noise. In Proceedings of the 2nd International Conference on Knowledge Discovery and Data Mining, IAAI, Portland, OR, USA, 2–4 August 1996; pp. 226–231. [Google Scholar]
  16. Duan, L.; Xu, L.; Guo, F.; Lee, J.; Yan, B. A local-density based spatial clustering algorithm with noise. Inf. Syst. 2007, 32, 978–986. [Google Scholar] [CrossRef]
  17. Leng, Z.; Zhang, J.; Ma, Y.; Zhang, J.; Zhu, H. A novel bathymetry signal photon extraction algorithm for photon-counting LiDAR based on adaptive elliptical neighborhood. Int. J. Appl. Earth Obs. Geoinf. 2022, 115, 103080. [Google Scholar] [CrossRef]
  18. Zhang, J.; Kerekes, J.; Csathó, B.; Schenk, T.; Wheelwright, R. A clustering approach for detection of ground in micropulse photon-counting LiDAR altimeter data. In Proceedings of the IEEE Geoscience and Remote Sensing Symposium, Quebec City, QC, Canada, 13–18 July 2014; pp. 177–180. [Google Scholar] [CrossRef]
  19. Wang, C.; Wang, A.; Rong, W.; Tao, Y.; Fu, R. Adaptive Denoising Algorithm for Photon-Counting LiDAR Point Clouds. Laser Optoelectron. Prog. 2021, 58, 1428001. [Google Scholar]
  20. Magruder, L.A.; Wharton III, M.E.; Stout, K.D.; Neuenschwander, A.L. Noise filtering techniques for photon-counting ladar data. In Proceedings of the Laser Radar Technology and Applications XVII, Baltimore, MD, USA, 23–27 April 2012; pp. 237–245. [Google Scholar]
  21. Wang, Z.; Chen, S.; Kong, W.; Liu, X. Comparison and Analysis of Denoising for Photon-Counting LiDAR Data. Laser Optoelectron. Prog. 2023, 60, 0628001. [Google Scholar]
  22. Zhu, X.; Nie, S.; Wang, C.; Xi, X.; Wang, J.; Li, D.; Zhou, H. A noise removal algorithm based on OPTICS for photon-counting LiDAR data. IEEE Geosci. Remote Sens. Lett. 2020, 18, 1471–1475. [Google Scholar] [CrossRef]
  23. Wang, Z.; Xi, X.; Nie, S.; Wang, C. Bathymetric Method of Nearshore Based on ICESat-2/ATLAS Data—A Case Study of the Islands and Reefs in The South China Sea. In Proceedings of the IGARSS 2022-2022 IEEE International Geoscience and Remote Sensing Symposium, Kuala Lumpur, Malaysia, 17–22 July 2022; pp. 2868–2871. [Google Scholar]
  24. Xie, D.; Li, G.; Wang, J.; Wang, Z.; Ye, F.; Yang, X. A novel denoising algorithm for photon-counting laser data based on LDBSCAN. In Proceedings of the AOPC 2019: Advanced Laser Materials and Laser Technology, Beijing, China, 7–9 July 2019; pp. 267–280. [Google Scholar]
  25. Ranndal, H.; Christiansen, P.S.; Kliving, P.; Andersen, O.B.; Nielsen, K. Evaluation of a Statistical Approach for Extracting Shallow Water Bathymetry Signals from ICESat-2 ATL03 Photon Data. Remote Sens. 2021, 13, 3548. [Google Scholar] [CrossRef]
  26. Huang, T.; Yang, G.; Tang, G. A fast two-dimensional median filtering algorithm. IEEE Trans. Acoust. Speech Signal Process. 1979, 27, 13–18. [Google Scholar] [CrossRef] [Green Version]
  27. Huber, P.J. Robust estimation of a location parameter. In Breakthroughs in Statistics: Methodology and Distribution; Springer: Berlin /Heidelberg, Germany, 1992; pp. 492–518. [Google Scholar]
  28. Huber, P.J. Robust statistics. In International Encyclopedia of Statistical Science; Springer: Berlin /Heidelberg, Germany, 2011; pp. 1248–1251. [Google Scholar]
  29. Wieser, A.; Brunner, F.K. Short Static GPS Sessions: Robust Estimation Results. GPS Solutions 2002, 5, 70–79. [Google Scholar] [CrossRef]
  30. Yuanxi, Y. Robust estimation for dependent observations. Manuscr. Geod. 1994, 19, 10–17. [Google Scholar]
  31. Yang, Y.; He, H.; Xu, G. Adaptively robust filtering for kinematic geodetic positioning. J. Geodesy 2001, 75, 109–116. [Google Scholar] [CrossRef]
  32. Yang, Y.; Song, L.; Xu, T. Robust estimator for correlated observations based on bifactor equivalent weights. J. Geodesy 2002, 76, 353–358. [Google Scholar] [CrossRef]
  33. Qu, W.; Chen, H.; Zhang, Q.; Gao, Y.; Wang, Q.; Hao, M. A robust estimation algorithm for the increasing breakdown point based on quasi-accurate detection and its application to parameter estimation of the GNSS crustal deformation model. J. Geodesy 2021, 95, 1–17. [Google Scholar] [CrossRef]
  34. Deng, D. DBSCAN clustering algorithm based on density. In Proceedings of the 2020 7th International Forum on Electrical Engineering and Automation (IFEEA), Hefei, China, 25–27 September 2020; pp. 949–953. [Google Scholar]
  35. Zanella, G.; Roberts, G. Scalable Importance Tempering and Bayesian Variable Selection. J. R. Stat. Soc. Ser. B Stat. Methodol. 2019, 81, 489–517. [Google Scholar] [CrossRef] [Green Version]
  36. Weining, Q. Method for selecting weight iteration with robust initial value. Geomat. Inf. Sci. Wuhan Univ. 2003, 28, 452–454. [Google Scholar]
  37. Angiulli, F.; Basta, S.; Pizzuti, C. Distance-based detection and prediction of outliers. IEEE Trans. Knowl. Data Eng. 2005, 18, 145–160. [Google Scholar] [CrossRef]
  38. Green, P.J. Iteratively reweighted least squares for maximum likelihood estimation, and some robust and resistant alternatives. J. R. Stat. Soc. Ser. B 1984, 46, 149–170. [Google Scholar] [CrossRef]
  39. Yang, Y. Robust estimation of geodetic datum transformation. J. Geodesy 1999, 73, 268–274. [Google Scholar] [CrossRef]
  40. Parrish, C.E.; Magruder, L.A.; Neuenschwander, A.L.; Forfinski-Sarkozi, N.; Alonzo, M.; Jasinski, M. Validation of ICESat-2 ATLAS Bathymetry and Analysis of ATLAS’s Bathymetric Mapping Performance. Remote Sens. 2019, 11, 1634. [Google Scholar] [CrossRef] [Green Version]
  41. Neumann, T.; Brenner, A.; Hancock, D.; Robbins, J.; Saba, J.; Harbeck, K.; Gibbons, A.; Lee, J.; Luthcke, S.; Rebold, T. ATLAS/ICESat-2 L2A Global Geolocated Photon Data, version 3; NASA National Snow and Ice Data Center Distributed Active Archive Center: Boulder, CO, USA, 2020; Volume 10. [Google Scholar]
Figure 1. (a) Photon tracks of 1036 orbit: the blue, green, and red lines are the photon tracks of the gt3r, gt2r, and gt1r beams; (b) photon tracks of 480 orbit: the blue, green, and red lines are the photon tracks of the gt1r, gt2r, and gt3r beams.
Figure 1. (a) Photon tracks of 1036 orbit: the blue, green, and red lines are the photon tracks of the gt3r, gt2r, and gt1r beams; (b) photon tracks of 480 orbit: the blue, green, and red lines are the photon tracks of the gt1r, gt2r, and gt3r beams.
Remotesensing 15 02051 g001
Figure 2. (a) Photon profile of the ATL03 raw data for the 480 orbit gt1r beam; (b) photon profile of the ATL03 raw data for the 480 orbit gt2r beam; (c) photon profile of the ATL03 raw data for the 1036 orbit gt2r beam; (d) photon profile of the ATL03 raw data for the 1036 orbit gt1r beam.
Figure 2. (a) Photon profile of the ATL03 raw data for the 480 orbit gt1r beam; (b) photon profile of the ATL03 raw data for the 480 orbit gt2r beam; (c) photon profile of the ATL03 raw data for the 1036 orbit gt2r beam; (d) photon profile of the ATL03 raw data for the 1036 orbit gt1r beam.
Remotesensing 15 02051 g002
Figure 3. Data point types of the DBSCAN algorithm.
Figure 3. Data point types of the DBSCAN algorithm.
Remotesensing 15 02051 g003
Figure 4. Schematic of robust M-estimation of water depth.
Figure 4. Schematic of robust M-estimation of water depth.
Remotesensing 15 02051 g004
Figure 5. Photon classification diagram of 480 orbit gt1r beam.
Figure 5. Photon classification diagram of 480 orbit gt1r beam.
Remotesensing 15 02051 g005
Figure 6. (a) DBSCAN denoising graph with “loose” parameter settings (eps = 6, MintPts = 3); (b) DBSCAN denoising graph with “strict” parameter settings (eps = 3, MintPts = 5).
Figure 6. (a) DBSCAN denoising graph with “loose” parameter settings (eps = 6, MintPts = 3); (b) DBSCAN denoising graph with “strict” parameter settings (eps = 3, MintPts = 5).
Remotesensing 15 02051 g006
Figure 7. (a) DBSCAN denoised 480 orbit gt1r beam after parameter adjustment (eps = 4, Mintpts = 3); (b) DBSCAN denoised 480 orbit gt2r beam after parameter adjustment (eps = 3, Mintpts = 2); (c) DBSCAN denoised 1036 orbit gt2r beam after parameter adjustment (eps = 4, Mintpts = 4); (d) DBSCAN denoised 1036 orbit gt1r beam after parameter adjustment (eps = 3, Mintpts = 4).
Figure 7. (a) DBSCAN denoised 480 orbit gt1r beam after parameter adjustment (eps = 4, Mintpts = 3); (b) DBSCAN denoised 480 orbit gt2r beam after parameter adjustment (eps = 3, Mintpts = 2); (c) DBSCAN denoised 1036 orbit gt2r beam after parameter adjustment (eps = 4, Mintpts = 4); (d) DBSCAN denoised 1036 orbit gt1r beam after parameter adjustment (eps = 3, Mintpts = 4).
Remotesensing 15 02051 g007aRemotesensing 15 02051 g007b
Figure 8. (a) Two-dimensional window filter denoised 480 orbit gt1r beam; (b) two-dimensional window filter denoised 480 orbit gt2r beam; (c) two-dimensional window filter denoised 1036 orbit gt2r beam; (d) two-dimensional window filter denoised 1036 orbit gt1r beam.
Figure 8. (a) Two-dimensional window filter denoised 480 orbit gt1r beam; (b) two-dimensional window filter denoised 480 orbit gt2r beam; (c) two-dimensional window filter denoised 1036 orbit gt2r beam; (d) two-dimensional window filter denoised 1036 orbit gt1r beam.
Remotesensing 15 02051 g008aRemotesensing 15 02051 g008b
Figure 9. (a) The two-step method denoised 480 orbit gt1r beam; (b) the two-step method denoised 480 orbit gt2r beam; (c) the two-step method denoised 1036 orbit gt2r beam; (d) the two-step method denoised 1036 orbit gt1r beam.
Figure 9. (a) The two-step method denoised 480 orbit gt1r beam; (b) the two-step method denoised 480 orbit gt2r beam; (c) the two-step method denoised 1036 orbit gt2r beam; (d) the two-step method denoised 1036 orbit gt1r beam.
Remotesensing 15 02051 g009aRemotesensing 15 02051 g009b
Figure 10. Correction of gt1r beam refraction for orbit 480.
Figure 10. Correction of gt1r beam refraction for orbit 480.
Remotesensing 15 02051 g010
Figure 11. (a) Point profile of bathymetric estimation of 480 orbit gt1r beam; (b) point profile of bathymetric estimation of 480 orbit gt2r beam; (c) point profile of bathymetric estimation of 1036 orbit gt2r beam; (d) point profile of bathymetric estimation of 1036 orbit gt1r beam.
Figure 11. (a) Point profile of bathymetric estimation of 480 orbit gt1r beam; (b) point profile of bathymetric estimation of 480 orbit gt2r beam; (c) point profile of bathymetric estimation of 1036 orbit gt2r beam; (d) point profile of bathymetric estimation of 1036 orbit gt1r beam.
Remotesensing 15 02051 g011
Figure 12. Comparisons of the robust M-estimation depth and manual denoising depth with a (a) 480 orbit gt1r beam; (b) 480 orbit gt2r beam; (c) 1036 orbit gt2r beam; and (d) 1036 orbit gt1r beam.
Figure 12. Comparisons of the robust M-estimation depth and manual denoising depth with a (a) 480 orbit gt1r beam; (b) 480 orbit gt2r beam; (c) 1036 orbit gt2r beam; and (d) 1036 orbit gt1r beam.
Remotesensing 15 02051 g012
Figure 13. Error scatter diagrams of the two methods in different environments. (a) Two-step method for 480 orbit gt1r beam, (b) 480 orbit gt2r beam, (c) 1036 orbit gt2r beam, and (d) 1036 orbit gt1r beam; (e) M-estimation for 480 orbit gt1r beam, (f) 480 orbit gt2r beam, (g) 1036 orbit gt2r beam, and (h) 1036 orbit gt1r beam.
Figure 13. Error scatter diagrams of the two methods in different environments. (a) Two-step method for 480 orbit gt1r beam, (b) 480 orbit gt2r beam, (c) 1036 orbit gt2r beam, and (d) 1036 orbit gt1r beam; (e) M-estimation for 480 orbit gt1r beam, (f) 480 orbit gt2r beam, (g) 1036 orbit gt2r beam, and (h) 1036 orbit gt1r beam.
Remotesensing 15 02051 g013
Table 1. The ICESat-2 data acquisition time and terrain details.
Table 1. The ICESat-2 data acquisition time and terrain details.
Local TimeDataTerrain
19:42:52–19:49:54480-20210425-gt1rflat
(night)480-20210425-gt2rundulating
10:28:06–10:35:161036-20210301-gt2rflat
(day)1036-20210301-gt1rundulating
Table 2. Table of denoising accuracy of different methods in different time periods and different terrains.
Table 2. Table of denoising accuracy of different methods in different time periods and different terrains.
Local Time DataTerrainDenoising AlgorithmF1 ScoreAccuracyPrecisionRecall
19:42:52–19:49:54
(night)
480-20210425-gt1rflatDBSCAN0.90460.83100.86800.9443
Two-dimensional window filter0.89960.84140.99560.8205
The two-step method0.94350.90870.99210.8995
480-20210425-gt2rundulatingDBSCAN0.71870.56880.66960.7755
Two-dimensional window filter0.76650.71740.92750.6531
The two-step method0.83430.83190.90960.7704
10:28:06–10:35:16
(day)
1036-20210301-gt2rflatDBSCAN0.67570.72960.58480.8000
Two-dimensional window filter0.46890.72080.71010.3500
The two-step method0.70380.83270.93490.5643
1036-20210301-gt1rundulatingDBSCAN0.67020.77450.58760.7798
Two-dimensional window filter0.49300.78110.77190.3621
The two-step method0.77800.89230.98730.6420
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhong, J.; Liu, X.; Shen, X.; Jiang, L. A Robust Algorithm for Photon Denoising and Bathymetric Estimation Based on ICESat-2 Data. Remote Sens. 2023, 15, 2051. https://doi.org/10.3390/rs15082051

AMA Style

Zhong J, Liu X, Shen X, Jiang L. A Robust Algorithm for Photon Denoising and Bathymetric Estimation Based on ICESat-2 Data. Remote Sensing. 2023; 15(8):2051. https://doi.org/10.3390/rs15082051

Chicago/Turabian Style

Zhong, Junsheng, Xiuguo Liu, Xiang Shen, and Liming Jiang. 2023. "A Robust Algorithm for Photon Denoising and Bathymetric Estimation Based on ICESat-2 Data" Remote Sensing 15, no. 8: 2051. https://doi.org/10.3390/rs15082051

APA Style

Zhong, J., Liu, X., Shen, X., & Jiang, L. (2023). A Robust Algorithm for Photon Denoising and Bathymetric Estimation Based on ICESat-2 Data. Remote Sensing, 15(8), 2051. https://doi.org/10.3390/rs15082051

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop