1. Introduction
Interferometric synthetic aperture radar (InSAR) is a remote sensing technique that combines electromagnetic wave interferometry with synthetic aperture radar (SAR) and extracts information on ground elevation or surface deformation, facilitating ground surveillance and measurement, by processing SAR image data through interferometry [
1,
2]. In recent years, with the rapid advancement of remote sensing technology and the widespread use of high-resolution satellites, spaceborne InSAR has increasingly become crucial in applications such as ground elevation measurement, surface deformation studies, and atmospheric research [
3,
4,
5,
6]. InSAR has many advantages in terrain mapping and surface deformation monitoring, such as being robust to time and weather, but the InSAR processing process can be affected by many factors, such as interference decorrelation, atmospheric phase, and phase unwrapping. These factors could cause reduced measurement accuracy or measurement failure, especially for single-baseline InSAR systems.
Among these factors, phase unwrapping is one of the key steps to obtaining high-precision terrain and surface deformation using InSAR [
7,
8,
9]. The phase unwrapping method used in InSAR processing is directly related to the accuracy of the solved terrain and surface deformation information [
10,
11]. There still does not exist a suitable phase unwrapping method for all cases, as especially in the case of large phase noise and steep terrain, many phase unwrapping methods will produce large errors. So far, many phase unwrapping methods have been proposed, which can mainly be divided into the following three categories, path-following methods [
12,
13,
14], optimal estimation methods [
15], and minimum norm methods [
16,
17,
18,
19], as well as other methods [
20,
21,
22,
23,
24,
25], such as the Kalman filtering methods [
20] and machine learning methods [
21,
22].
The first category, path-following methods, is represented by the Goldstein branch-cutting method [
12,
14] and the quality map guidance method [
13]. The Goldstein branch-cutting method identifies positive and negative residual points, connects the residual points according to the nearest neighbor principle, and generates “branch cuts”. It prevents the propagation of errors by following the principle that the integral path does not cross the “branch cuts”. The Goldstein branch-cutting method has the obvious advantages of fast speed and high accuracy under the condition of high signal-to-noise ratio (SNR) and low phase gradient. However, when the noise is large and the phase gradient is steep, it is difficult for this method to accurately connect branch tangents, which often leads to errors in integration path selection and introduces errors into the unwrapping stage [
14]. The quality map guidance method does not identify residual points and does not set branch tangents. It is a method that defines the wrapping phase quality through the phase quality map, controls the integration path to advance along the high-quality pixels to the low-quality pixels, and sequentially unwraps the wrapping phase. However, its accuracy depends largely on the reliability of the quality map; when the noise is large and the phase gradient is steep, the residuals will be distributed in the high-quality area, and the quality-map guidance method will not work properly [
13].
The second category, optimal estimation methods, such as the network flow method, is also commonly used for phase unwrapping. Unlike the methods mentioned earlier, the network flow method transforms the phase problem into a network flow problem with the concept of minimum cost flow, and it limits the phase error transmission in the low-mass region and obtains a globally optimal solution by minimizing the difference between the phase gradient of the unwrapped phase and the discrete partial derivatives of the unwrapped phase. However, for low SNR and steep phase gradient, the error can still be significant.
Compared to the methods mentioned above, the third category, the minimum norm-based methods, has the fastest unwrapping speed and smoothest phase unwrapping results [
15]. The least squares (LS) phase unwrapping method is one of the most commonly used methods of the minimum norm-based methods. The LS phase unwrapping method can be understood as finding the solution to the Poisson equation under Neumann boundary conditions, which can be achieved through iterative methods, discrete cosine transforms (DCT), fast Fourier transform (FFT), and an unweighted multi-level grid method, but it also produces significant unwrapping errors with large phase noise and high phase gradient [
17]. In recent years, some scholars have improved the LS phase unwrapping method based on its rapidity and smoothness, making the LS phase unwrapping method suitable for high phase noise and high phase gradient. For example, Xia et al. [
26] proposed a phase calibration unwrapping algorithm for phase data corrupted by strong decorrelation speckle noise (PCUA). This algorithm uses the average phase gradient to replace the phase gradient that is greater than the standard deviation when solving the Poisson equation. The advantages of the LS phase unwrapping method can be exploited by reducing the phase gradient, and then the error phase is corrected by the iterative method. However, the phase unwrapping effect is unstable when the phase noise and phase gradient change. After this, to make the phase unwrapping results more stable, Zong et al. [
27] proposed a robust phase unwrapping algorithm for noisy and segmented phase measurements (RPUA). When solving the Poisson equation, it uses the
e function to correct the phase gradient, which further improves the phase unwrapping accuracy under high noise and high phase gradient. However, in true interference data phase unwrapping, the results obtained are also unstable.
Based on the above analysis, there are certain difficulties in phase unwrapping. Using traditional category I and category II methods, phase unwrapping under high phase noise and steep phase gradient will produce large errors [
24,
25]. Xia et al. and Zong et al. show that using iterative LS phase unwrapping can reduce the phase unwrapping error under high phase noise and steep phase gradient, but when calculating the equivalent phase gradient, it is too dependent on the overall steepness of the phase gradient to change the steep phase gradient into a sufficiently large number of flat phase gradients. However, the calculation of an equivalent flat phase gradient is too dependent on the overall steepness of the phase gradient. It is difficult to change a steep phase gradient into a sufficiently flat phase gradient.
Therefore, to make the LS phase unwrapping method more stable under high phase noise and high phase gradient, this paper proposes an improved LS phase unwrapping method based on the Chebyshev filter [
28] and designs its adaptive filtering threshold so that the high phase gradient is changed into flatter phase gradient as much as possible. Finally, an iterative method of unwrapping phase errors is used to eliminate errors under high phase noise and steep phase gradient.
This paper is organized as follows:
Section 2 describes the principle of LS phase unwrapping and the improvement method of this paper.
Section 3 is the experimental results and analysis part.
Section 4 is a discussion of phase unwrapping methods.
Section 5 is the conclusion of this paper.
3. Results
The simulated data and the real InSAR phase data are used in this section to verify the effectiveness of the proposed method of unwrapping. The simulated wrapped phase data were generated using the PEAKS function in MATLAB (Version 9.8.0.1323502, R2020a), and the real InSAR phase data were obtained by InSAR processing of the Sentinel-1 satellite [
29].
To verify the superiority of the improved LS phase unwrapping method in this paper, we compared the proposed method with traditional LS, as well as the recently improved LS phase unwrapping methods, PCUA and RPUA. The proposed method is an improvement on the LS method, so the LS method was used in the comparison. PCUA and RPUA are also improved methods of the LS method in recent years and have good phase unwrapping effects. Therefore, comparison with LS, PCUA, and RPUA can better demonstrate the effectiveness of our proposed method.
The iterative method is also used in PCUA and RPUA, and
and
are set in the experiment. Root mean square error (RMSE) between the unwrapped phase and the true phase is used to evaluate the quality of phase unwrapping [
30]:
3.1. The Simulated Data Results
The simulated wrapped phase was generated using the PEAKS function in MATLAB with a peak to valley of 10, with an adding Gaussian noise of 0 mean and 1 standard deviation. The simulated data have a horizontal and vertical size of 512 × 512, as shown in
Figure 2. The color bar on the right side of the figures represents the phase value, the unit is rad.
In
Figure 2a, the highest phase is 83.84 rad and the lowest phase is −68.23 rad. The wrapped phase map is shown in
Figure 2b. Since the wrapped phase is distributed between −π and π, noise with a standard deviation of 1 is added, which satisfies the high-noise situation. When the Y direction is at the 250th pixel, the phase gradient map of
Figure 2b is shown in
Figure 3, which shows that the simulated phase has a steep phase gradient. Therefore, the simulated phase used has a large noise and steep phase gradient, which has a positive effect on verifying the method proposed in this paper.
Then, LS, PCUA, RPUA, and the proposed method are used for phase unwrapping. The unwrapping error maps are shown in
Figure 4.
Table 1 presents the time and RMSE values for each of the four methods. As shown in
Figure 4a, the LS phase unwrapping method has a lot of unwrapped areas. As shown in
Table 1, the advantage of the LS phase unwrapping method is that it can perform phase unwrapping fastest. But it has the largest RMSE and the worst unwrapping effect. As shown in
Figure 4b,c, PCUA and RPUA have greatly improved the LS phase unwrapping method, and have good phase unwrapping results, but there are still a few unwrapped areas. As shown in
Table 1, PCUA and RPUA have smaller RMSE compared to the LS phase unwrapping method but require more unwrapping time, which has room for improvement. As shown in
Figure 4d, the proposed method has the least unwrapped areas, and as shown in
Table 1, the proposed method has the smallest RMSE, which is reduced by 94.53% compared with LS, 31.50% compared with PCUA and 15.11% compared with RPUA. The proposed method has a faster unwrapping speed compared with PCUA and RPUA, which is improved by 55.29% and 10.91%. Therefore, this simulation experiment shows that the proposed method has a better effect than LS, PCUA, and RPUA.
Furthermore, to test the robustness of the proposed method, phase unwrapping tests were performed by increasing the phase gradient while keeping the noise constant. Based on the above simulated experimental data, the same MATLAB PEAKS function was used to generate experimental data. While the noise standard deviation remained unchanged at 1, the peak to the valley gradually increased from 10 to 15. The RMSE of the unwrapping results with different phase gradients is shown in
Figure 5. Among them, the RMSE of LS phase unwrapping is the largest and gradually increases as the phase gradient increases. The unwrapping effect is the worst and is not suitable for steep phase unwrapping. PCUA, RPUA, and the proposed method can all achieve lower RMSE than LS and can perform stable phase unwrapping under a high phase gradient compared to LS. However, when the peaks to the valley is 13.5, the RMSE of the unwrapping results of PCUA and RPUA will mutate and be unstable. For proposed the method, compared with the other three methods, the RMSE of the unwrapping results has always been at the lowest, and there is no RMSE mutation when the phase gradient changes. Therefore, the proposed method is very little affected by the phase gradient can almost be ignored, and can perform stable phase unwrapping at steep phases.
Then, with the peak to the valley constant at 10, the standard deviation of the noise is gradually increased from 1 to 1.5, and the unwrapped resulting RMSE is shown in
Figure 6. Like the above phase unwrapping results, the RMSE of LS phase unwrapping is the largest and gradually increases with the increase in phase noise. The phase unwrapping effect of LS is the worst and is not suitable for phase unwrapping under high noise. PCUA, RPUA, and the proposed method all improve the effect of LS phase unwrapping under the influence of noise. From the RMSE of the phase unwrapping results in
Figure 6, it can be seen that RPUA is affected by phase noise the most compared to PCUA and the proposed method. The phase unwrapping results of the proposed method have the lowest RMSE and are least affected by phase noise. Therefore, the proposed method has the best phase unwrapping effect under high noise compared with LS, PCUA, and RPUA.
According to the above simulation data analysis, the method proposed in this paper has the best robustness compared with LS, PCUA, and RPUA under different phase gradients and different phase noises.
3.2. The Real InSAR Phase Unwrapping Results of Sentinel-1
We selected two sentinel-1 IW data from 12 January 2023 and 24 January 2023 covering the position (30.43 N, 103.02 E), and the experimental data are distributed over strip 6 of IW2. The Sentinel-1 satellite is a C-band SAR launched by ESA in 2014; its data are currently available for free on the Copernicus Data Space Ecosystem. It has been widely used in marine environment monitoring, land change detection, earthquake detection, and landslide detection. Using real Sentinel-1 InSAR phase data for phase unwrapping experiments can better reflect the advantages and disadvantages of the phase unwrapping algorithm.
The terrain data located in Yaan City, Sichuan Province, China, with a maximum and minimum elevation difference of about 500 m was used. Before phase unwrapping, 1 × 5 multi-looking was applied to the data in azimuth and range. Phase interference was performed after image registration. Then the reference ellipsoid phase was removed, and then the Goldstein filter [
31] with a coefficient of 0.5 was selected for phase filtering, and finally, a 500 × 500 wrapped phase map in the horizontal and vertical was obtained. As shown in
Figure 7a, it can be seen that the phase interference fringes are relatively dense and there is a large amount of noise residue, which provides favorable conditions for verifying the effectiveness of this method.
Figure 7b shows the real-terrain simulated phase, which uses the Digital Elevation Model of Shuttle Radar Topography Mission [
32] information and was generated under the same conditions as the InSAR phase data.
Using the same unwrapping methods as the simulated test, the unwrapping error results are shown in
Figure 8 and
Table 2.
It can be seen from
Figure 8a that the incompletely unwrapped phase range in the error map obtained by the LS phase unwrapping method is the highest. From
Table 2, it can be seen that the RMSE obtained by the LS phase unwrapping method is the largest and the error range is the widest; its biggest advantage is that it can be quickly unwrapping. It can be seen from
Figure 8b,c that PCUA and RPUA reduce the incompletely unwrapped area and phase range compared to the LS phase unwrapping method, but there will still be more incompletely unwrapped areas and larger phase error. As can be seen from
Table 2, PCUA takes the most phase unwrapping time but has a lower RMSE than LS and RPUA. Although RPUA reduces the processing time, it has a larger RMSE and a wider error range than PCUA. As can be seen from
Figure 8d and
Table 2, the proposed method has the smallest incomplete unwrapping area and the smallest incomplete unwrapping phase range, and the unwrapping time is about 62.86% and 11.64% shorter than PCUA and RPUA, respectively. Compared with LS, PCUA, and RPUA, the proposed method has the smallest RMSE, which is reduced by 63.91%, 35.38%, and 54.39%, respectively.
4. Discussion
InSAR has many advantages in terrain height measurement and surface deformation measurement, but InSAR processing will be affected by many factors, especially in terms of phase unwrapping, which will directly affect the accuracy of the measurement results. Currently, existing phase unwrapping methods have certain flaws, and there is currently no perfect phase unwrapping method suitable for all situations. Especially when the phase noise is large and the phase gradient is steep, the mainstream phase unwrapping methods will produce large unwrapping errors.
In recent years, many researchers have made improvements to the LS phase unwrapping method to make it suitable for high phase noise and steep phase gradient, such as PCUA and RPUA, which are experimentally compared in this paper. However, there will still be a large phase unwrapping error, and when the phase noise and phase gradient change, PCUA and RPUA are unstable, and the generated RMSE is sudden.
Therefore, this paper makes further improvements based on LS phase unwrapping. In experiments, it was found that the factor that affects the accuracy of phase unwrapping is the phase gradient, and noise is also one of the factors that affect the change of the phase gradient. When the phase gradient is relatively flat, the LS phase unwrapping method can perform phase unwrapping quickly and accurately. However, when the phase gradient is steep, the LS phase unwrapping method can only perform phase unwrapping quickly and cannot accurately obtain phase unwrapping results. Therefore, this paper uses a low-pass filter to change the steep phase gradient into a flat phase gradient for phase unwrapping. Based on the properties of many low-pass filters, the Chebyshev low-pass filter is selected and improved, which is added in phase unwrapping. To maintain the phase characteristics, an iterative method is used to add all the phase unwrapped results after low-pass filtering to obtain the final unwrapped phase.
5. Conclusions
This paper proposes an improved least squares phase unwrapping method based on the Chebyshev filter. The accuracy and robustness of the method were verified by comparing the simulated phase and the real InSAR phase. Simulation data indicate that the proposed method’s RMSE is 94.53%, 31.50%, and 15.11% lower than that of LS, PCUA, and RPUA, respectively. The proposed method has a faster unwrapping speed compared with PCUA and RPUA, which is improved by 55.29% and 10.91%. The robustness of the proposed method was validated through its performance in phase unwrapping under varying phase gradients and phase noise. Under the real InSAR phase unwrapping, the phase unwrapping RMSE of the proposed method is reduced by 63.91%, 35.38%, and 54.39% compared to LS, PCUA, and RPUA. The phase unwrapping time is reduced by 62.86% and 11.64% compared to PCUA and RPUA. Based on the analysis of the above results, the proposed method has the most accurate unwrapping effect and best robustness compared to LS, PCUA, and RPUA.
However, in the simulated and real phase unwrapping results obtained, there is still some incompletely unwrapped phase, and the proposed method speed is slower than the LS phase unwrapping speed. In the future, the speed and efficiency of the phase unwrapping can be improved by adding more appropriate filters to divide the phase of the steep gradient into more flat gradient phases or adding Chebyshev filters to other phase unwrapping methods to improve the phase unwrapping effect of other methods.