Removing InSAR Topography-Dependent Atmospheric Effect Based on Deep Learning

: Atmospheric effects are among the primary error sources affecting the accuracy of interferometric synthetic aperture radar (InSAR). The topography-dependent atmospheric effect is particularly noteworthy in reservoir areas for landslide monitoring utilizing InSAR, which must be effectively corrected to complete the InSAR high-accuracy measurement. This paper proposed a topography-dependent atmospheric correction method based on the Multi-Layer Perceptron (MLP) neural network model combined with topography and spatial data information. We used this proposed approach for the atmospheric correction of the interferometric pairs of Sentinel-1 images in the Baihetan dam. We contrasted the outcomes with those obtained using the generic atmospheric correction online service for InSAR (GACOS) correction and the traditional linear model correction. The results indicated that the MLP neural network model correction reduced the phase standard deviation of the Sentinel-1 interferogram by an average of 64% and nearly eliminated the phase-elevation correlation. Both comparisons outperformed the GACOS correction and the linear model correction. Through two real-world examples, we demonstrated how slopes with displacements, which were previously obscured by a signiﬁcant topography-dependent atmospheric delay, could be successfully and clearly identiﬁed in the interferograms following the correction by the MLP neural network. The topography-dependent atmosphere can be better corrected using the MLP neural network model suggested in this paper. Unlike the previous model, this proposed approach could be adjusted to ﬁt each interferogram, regardless of how much of the topography-dependent atmosphere was present. In order to improve the effectiveness of DInSAR and time-series InSAR solutions, it can be applied immediately to the interferogram to retrieve the effective displacement information that cannot be identiﬁed before the correction.


Introduction
InSAR technology is a geodetic measurement technology [1] based on spaceborne Earth observation, which has been developed rapidly in the past three decades. Satellite radar remote sensing has prominent strengths, such as wide-range coverage, high accuracy, and all-day and all-weather applicability [2]. At present, although InSAR technology has been widely used in landslide identification and monitoring [3][4][5][6][7][8], seismic

Data Source
A total of 17 Sentinel-1 images, including 10 images of the descending orbit and 7 images of the ascending orbit, were used in this study. The Sentinel-1 satellite incorporates four modes of operation: stripmap mode (SM), interferometric wide swath mode (IW), extra-wide swath mode (EW), and wave mode (WV). In this paper, we selected singlelook complex (SLC) images of VV polarization in IW mode, spanning from 21 April 2021 to 4 July 2021, all after the impoundment of the Baihetan Reservoir. The Sentinel-1 single satellite uses topographic observation in progressive scans (TOPS) mode to acquire images on a global scale every 12 days [43]; that is, the satellite has a revisit period of 12 days. Other parameters, such as the incidence angle and resolution of the Sentinel-1 satellite, are shown in Table 1.

Data Source
A total of 17 Sentinel-1 images, including 10 images of the descending orbit and 7 images of the ascending orbit, were used in this study. The Sentinel-1 satellite incorporates four modes of operation: stripmap mode (SM), interferometric wide swath mode (IW), extra-wide swath mode (EW), and wave mode (WV). In this paper, we selected single-look complex (SLC) images of VV polarization in IW mode, spanning from 21 April 2021 to 4 July 2021, all after the impoundment of the Baihetan Reservoir. The Sentinel-1 single satellite uses topographic observation in progressive scans (TOPS) mode to acquire images on a global scale every 12 days [43]; that is, the satellite has a revisit period of 12 days. Other parameters, such as the incidence angle and resolution of the Sentinel-1 satellite, are shown in Table 1.

DInSAR Processing
Differential Interferometry Synthetic Aperture Radar (DInSAR) monitor surface displacement was performed using the phase information of 2 SAR images acquired from the SAR sensor at different times in the same area [44]. As shown in Figure 2, SAR images covering the same area were acquired and co-registered with the precise orbit data to form interferograms. Shuttle radar topography mission digital elevation model (SRTM DEM) data, with a 30 m resolution, were used to remove the topographic and flat effect from the interferogram. The Goldstein adaptive filtering method [45] was used for filtering to remove the noise. Finally, a displacement map was acquired after phase unwrapping using the minimum cost flow (MCF) method [46,47].

Date of launch
April 2014 The angle of incidence 29.

DInSAR Processing
Differential Interferometry Synthetic Aperture Radar (DInSAR) monitor surface displacement was performed using the phase information of 2 SAR images acquired from the SAR sensor at different times in the same area [44]. As shown in Figure 2, SAR images covering the same area were acquired and co-registered with the precise orbit data to form interferograms. Shuttle radar topography mission digital elevation model (SRTM DEM) data, with a 30 m resolution, were used to remove the topographic and flat effect from the interferogram. The Goldstein adaptive filtering method [45] was used for filtering to remove the noise. Finally, a displacement map was acquired after phase unwrapping using the minimum cost flow (MCF) method [46,47].

MLP-Based Topography-Dependent Atmospheric Correction Method
In the deep learning community, we have witnessed the rapid development of neural networks. The architecture of neural networks developed from the MLP (Multi-Layer Perceptron), to the CNN (Convolutional Neural Network), to the Transformer, and back to the MLP [48]. In this study, we proposed an MLP with ten layers, referred to as FatMLP, for removing the atmospheric phase in InSAR. FatMLP was designed to learn high-order non-linear interactions of the digital elevation model (DEM), longitude and latitude.
The architecture of FatMLP is shown in Figure 2, which consists of ten layers. In FatMLP, an MLP layer was used to learn the linear correlations of the input variables. We used a wider MLP to capture more linear correlations. We designed deeper layers for higher-order nonlinear interactions of the input variables. The input layer included a fully connected layer with 3 neurons, which represent the DEM, longitude and latitude, respectively. Each pixel of the DEM was treated as a sample. First, the input data were lifted into a high-dimensional space, in which the representation of the original data was limited [49]. To this end, the first two hidden layers were fully connected layers with 4096 neurons. Such a high-dimensional space is useful for learning complex interactions of the input variables. Second, the intermediate sparse data were gradually compressed. The third and fourth hidden layers were two fully connected layers with 2048 neurons. The fifth and sixth layers were two fully connected layers with 1024 neurons. The seventh and eighth layers were two fully connected layers with 512 neurons. Finally, the output layer was fully connected with 1 neuron, representing the atmospheric phase. For each hidden layer, the activation function was the Rectified Linear Unit (ReLU), which was used to learn non-linear interactions of the input variables. The depth of FatMLP was greater than conventional MLPs. Deeper neural networks are useful for learning highorder interactions of the input variables. Equations (1)- (5) show the implementation details of the proposed neural network, where FC represents a fully connected layer. The following hyperparameters and configurations are displayed: 9192 is the batch size; the Adaptive Moment Estimation (Adam) is the optimizer, which is one of the most widely used optimization algorithms in deep learning, with the advantages of fewer hyperparameters and good performance; the rate of learning is 0.001; and the proposed neural network was trained for 50 epochs.

Linear Model Correction
The distribution of atmospheric water vapour is the primary cause of tropospheric delays, with impacts of atmospheric vertical stratification introducing topography-dependent delayed phases and turbulence introducing local delayed phases. In contrast to delays caused by turbulent mixing, vertical stratification effects are regarded as static stratification in a particular location and time scale. In the case of changes in surface elevation, the vertical stratification tropospheric delay increases as the water vapour concentration and pressure to temperature ratio rise [50]. The concentration of water vapour declines exponentially with altitude, and the theoretical delay curve is an exponential function of altitude [29]. Suppose the interferogram is majorly affected by the static tropospheric delay. In this case, the delay signal is the difference between the corresponding delay curves of the primary and secondary SAR images, and the resulting exponential function is a Taylor-expanded and second-order function. Higher-order terms are neglected, and a linear relationship can be obtained as follows [51]: where b is a deviation term, h is the topography, ∆∅ is the phase, and K is the transfer function between the topography and the phase. The correction of tropospheric delay using a linear model is applicable to most weather conditions, except for special cases, such as reverse or non-monotonic troposphere. Since the elevation-phase relationship may have a local trade-off with the elevation displacement relationship, the transfer function K is best determined through global linear regression [52].

GACOS Correction Method
GACOS, an InSAR atmospheric correction model developed by Zhenhong Li' (2006) [22] at Newcastle University, United Kingdom, is a numerical atmospheric product based on  [53]. InSAR users can obtain ZTD products for the corresponding time period by submitting the location of the study area and the time of SAR image acquisition through (http://www.gacos.net) (accessed on 5 January 2022).

Analysis of Correction Effect
As shown by the analysis in Figure 3, the interferograms a-f suffered from more serious atmospheric delay effects associated with elevation, and the atmospheric delay phase was related to the elevation change. The atmospheric simulated interferogram (i.e., MLP-modelled atmospheric correction that needs to be removed from the original interferograms) is shown in Figure 3g,h. The topography-dependent atmospheric delay effects in Figure 3m-r were reduced to different degrees after correction by the MLP neural network method, and the phases of the interferograms were smoother. Overall, the effect of atmospheric delay phase removal was more obvious after correction using the MLP neural network method. Meanwhile, the method is not required to build models, does not introduce any external data, has certain applicability for topography-dependent atmospheric corrections, and the atmospheric correction effect is relatively significant.  Table 2 shows the statistics of the change in the phase standard deviation of the original deconvolution of each interferometric pair. The analysis revealed that the phase standard deviation of all interferometric pairs after correction by the MLP neural network model was significantly reduced compared to the phase standard deviation before correction, with the 20210604-20210610 interferometric pair having the best correction effect and  Table 2 shows the statistics of the change in the phase standard deviation of the original deconvolution of each interferometric pair. The analysis revealed that the phase standard deviation of all interferometric pairs after correction by the MLP neural network model was significantly reduced compared to the phase standard deviation before correction, with the 20210604-20210610 interferometric pair having the best correction effect and the highest change rate of up to 84%.  Figure 4 shows the corrective effect of the Sentinel-1 (21 April 2021 to 4 July 2021) interferogram with the three correction methods. As seen in Figure 4, there was a strong association between the interferogram phase spatial distribution and the topographic relief of the study area, indicating that the delayed vertical stratification of the atmosphere predominates in terms of the overall trend. After correction using the conventional linear model, GACOS, and MLP neural network model, the combination of the second and third rows of Figure 4 demonstrated that the topography-dependent atmospheric delay in the linear model and GACOS correction findings were attenuated to varying degrees. The analysis of the fourth row of Figure 4 revealed, however, that after the correction of the MLP neural network model, the simulated phase of the atmosphere was visibly lowered in comparison to the phase value before correction, and the phase value in the majority of places was close to 0. Thus, the MLP neural network correction had a greater impact than the linear model and GACOS.

Linear Model and GACOS Correction Method
To optimize the use of different interferometric SAR atmospheric correction methods and to avoid potential uncertainty, it is necessary to establish appropriate statistical metrics and evaluate the correction performance [54]. The phase standard deviation was applied to measure the effect of atmospheric delay estimation and correction. In particular, the reduction in the interferometric phase was used to demonstrate the reduction in the corrected atmospheric delay phase [36,37]. Using the three aforementioned techniques, atmospheric effects were removed from 27 interferograms, and the standard deviation of the interferometric phase before and after correction was calculated to evaluate the overall performance of the correction results ( Figure 5). Figure 5 shows the magnitude of the reduction in phase standard deviation of the 27 interferometric pairs in the Baihetan reservoir region after atmospheric correction using linear model correction, external data GACOS and MLP neural network, respectively. The analysis revealed that after linear correction and GACOS correction, the phase standard deviation of some interference pairs decreased, which was not significant, while the phase standard deviation of some interference pairs increased instead. Nevertheless, with MLP neural network correction, the phase standard deviation of all 27 interference pairs was reduced, and the reduction was more significant. rows of Figure 4 demonstrated that the topography-dependent atmospheric delay in the linear model and GACOS correction findings were attenuated to varying degrees. The analysis of the fourth row of Figure 4 revealed, however, that after the correction of the MLP neural network model, the simulated phase of the atmosphere was visibly lowered in comparison to the phase value before correction, and the phase value in the majority of places was close to 0. Thus, the MLP neural network correction had a greater impact than the linear model and GACOS. To optimize the use of different interferometric SAR atmospheric correction methods and to avoid potential uncertainty, it is necessary to establish appropriate statistical metrics and evaluate the correction performance [54]. The phase standard deviation was applied to measure the effect of atmospheric delay estimation and correction. In particular, the reduction in the interferometric phase was used to demonstrate the reduction in the corrected atmospheric delay phase [36,37]. Using the three aforementioned techniques, atmospheric effects were removed from 27 interferograms, and the standard deviation of  Table 3 shows that after correction using the linear model, the interferometric pairs with reduced phase standard deviation accounted for 59% of the 27 interferometric pairs, which is the lowest among the three methods. The GACOS atmospheric correction had 88% of the interferometric pairs with reduced phase standard deviation of the unwinding. The percentage of interferograms with reduced standard deviation was larger than that of the linear model. However, the percentage of interferometric pairs with reduced phase standard deviation after correction using the MLP neural network was 100%, which was significantly higher than that of the linear model and GACOS. Hence, for the topography-dependent atmospheric delay in the Baihetan reservoir area, the atmospheric correction results of the MLP neural network were both superior to those of the linear model and the external data GACOS, and the GACOS was superior to the linear model. Compared with the atmospheric delay correction accuracy, the phase standard deviation of all interferometric pairs was reduced by 11.1% on average after using the linear model for atmospheric delay correction, which was the lowest among the three methods. However, after atmospheric delay correction using GACOS, the average reduction in the standard deviation of the unwrapping phase was 17.4%, which was higher than the linear model and located in the second rank. Finally, after the correction of atmospheric delay using the MLP neural network, the average reduction in phase standard deviation was approximately 64%, which was much higher than the first two correction methods, and the correction effect was more obvious for the topography-dependent atmospheric delay in the Baihetan reservoir area.
Remote Sens. 2022, 14, x FOR PEER REVIEW 9 of 16 the interferometric phase before and after correction was calculated to evaluate the overall performance of the correction results ( Figure 5).  Figure 5 shows the magnitude of the reduction in phase standard deviation of the 27 interferometric pairs in the Baihetan reservoir region after atmospheric correction using linear model correction, external data GACOS and MLP neural network, respectively. The analysis revealed that after linear correction and GACOS correction, the phase standard deviation of some interference pairs decreased, which was not significant, while the phase standard deviation of some interference pairs increased instead. Nevertheless, with MLP neural network correction, the phase standard deviation of all 27 interference pairs was reduced, and the reduction was more significant. Table 3 shows that after correction using the linear model, the interferometric pairs with reduced phase standard deviation accounted for 59% of the 27 interferometric pairs, which is the lowest among the three methods. The GACOS atmospheric correction had 88% of the interferometric pairs with reduced phase standard deviation of the unwinding. The percentage of interferograms with reduced standard deviation was larger than that of the linear model. However, the percentage of interferometric pairs with reduced phase standard deviation after correction using the MLP neural network was 100%, which was significantly higher than that of the linear model and GACOS. Hence, for the topographydependent atmospheric delay in the Baihetan reservoir area, the atmospheric correction results of the MLP neural network were both superior to those of the linear model and the external data GACOS, and the GACOS was superior to the linear model. Compared with the atmospheric delay correction accuracy, the phase standard deviation of all interferometric pairs was reduced by 11.1% on average after using the linear model for atmospheric delay correction, which was the lowest among the three methods. However, after atmospheric delay correction using GACOS, the average reduction in the standard deviation of the unwrapping phase was 17.4%, which was higher than the linear model and located in the second rank. Finally, after the correction of atmospheric delay using the  Table 3. The number of, percentage of, and average reduction in interference pairs with reduced phase standard deviation and average increment of interference pairs with increased phase standard deviation of the interferogram after the correction using the three correction methods.  [36] estimated the correlation coefficient between the interferometric phase and elevation as an alternative to the ratio between turbulent and vertical stratified delays. Here, we used the correlation coefficient to describe the elevation dependence in the interferogram. The smaller difference from zero indicates the lower proportion of signals related to the terrain in the interferometric phases. If the absolute value of the correlation decreases after atmospheric correction, it represents the mitigation of the vertical stratified components [36].

The Interferometric Pairs with a Reduced Standard Deviation
As demonstrated by the example of the analysis of the correlation atmospheric delay correction for the elevation of the Baihetan reservoir area in Figure 6, the correlation between phase and elevation for the original interferograms 20210604-20210610 (Figure 6a) was −0.092. After the linear model and GACOS correction, the correlation between phase and elevation was reduced to −0.082 and −0.002, respectively, as shown in Figure 6b,c. However, there was still some correlation between phase and elevation. Figure 6d illustrates that after rectification using the MLP neural network model, the correlation was reduced to −0.000036, which is quite close to 0. The correlation between phase and elevation was rather weak at this stage. metric phase and elevation as an alternative to the ratio between turbulent and vertica stratified delays. Here, we used the correlation coefficient to describe the elevation de pendence in the interferogram. The smaller difference from zero indicates the lower pro portion of signals related to the terrain in the interferometric phases. If the absolute valu of the correlation decreases after atmospheric correction, it represents the mitigation o the vertical stratified components [36].
As demonstrated by the example of the analysis of the correlation atmospheric delay correction for the elevation of the Baihetan reservoir area in Figure 6, the correlation be tween phase and elevation for the original interferograms 20210604-20210610 (Figure 6a was −0.092. After the linear model and GACOS correction, the correlation between phas and elevation was reduced to −0.082 and −0.002, respectively, as shown in Figure 6b,c However, there was still some correlation between phase and elevation. Figure 6d illus trates that after rectification using the MLP neural network model, the correlation wa reduced to −0.000036, which is quite close to 0. The correlation between phase and eleva tion was rather weak at this stage.

Displacement Identification after MLP Correction
When InSAR technology is applied to early displacement monitoring, some minor displacement information is prone to being obscured due to the atmospheric delay effect in the study area, which poses an invisible threat to human life and property. Therefore, it is essential to effectively identify the deformers in some areas through atmospheric correction. In this paper, the interferometric pair 20210517-20210604 with severe atmospheric delay in the Baihetan reservoir area was atmospherically corrected using the above three methods. The results of the atmospheric correction were processed for phase turn displacement, and it was found that there were signs of displacement in the local area. The specific displacement results are shown in Figure 7. The corrected displacement results of the original displacement results, linear model, GACOS, and MLP neural network model are represented in Figure 7a-d, respectively.
As seen in Figure 7, the original displacement results map can be seen in the slope area without signs of displacement, surrounded by atmospheric phase error. It was seriously impacted, and was unable to determine the exact displacement location and displacement magnitude accurately. Combined with Figure 7b,c, the linear model and GACOS correction were ineffective, and the displacement region was not easily identified. However, from Figure 7d, it is clear that the MLP neural network model had a remarkably effective correction effect. In addition to the obvious displacement in the slope ranges, the displacement volume in the non-displacement area was approximately 0mm, and it was able to analyse the displacement area more accurately. rection. In this paper, the interferometric pair 20210517-20210604 with severe atmospheric delay in the Baihetan reservoir area was atmospherically corrected using the above three methods. The results of the atmospheric correction were processed for phase turn displacement, and it was found that there were signs of displacement in the local area. The specific displacement results are shown in Figure 7. The corrected displacement results of the original displacement results, linear model, GACOS, and MLP neural network model are represented in Figure 7a-d, respectively. As seen in Figure 7, the original displacement results map can be seen in the slope area without signs of displacement, surrounded by atmospheric phase error. It was seriously impacted, and was unable to determine the exact displacement location and displacement magnitude accurately. Combined with Figure 7b,c, the linear model and GA-COS correction were ineffective, and the displacement region was not easily identified. However, from Figure 7d, it is clear that the MLP neural network model had a remarkably effective correction effect. In addition to the obvious displacement in the slope ranges, the displacement volume in the non-displacement area was approximately 0mm, and it was able to analyse the displacement area more accurately.
The displacement area, located near Sanjia Village, Yogu Country, Huidong County, Liangshan Yi Autonomous Prefecture, Sichuan Province, with a central longitude of 103°1′19″ and a central latitude of 26°36′59″, was further analysed. This slope elevation was approximately 772-1028 m, with a differential height of 256 m and a relatively steep slope. Combined with Figure 8a,b, the strong displacement area of this slope was more obvious (as shown by the blue dotted line), which was located in the median-lower part of the slope, with a maximum line-of-sight displacement of up to approximately 28 mm. Meanwhile, the interferometric pairs 20210604-20210610 with severe atmospheric delay in the Baihetan reservoir area were corrected for atmospheric correction using the three methods mentioned above, and the results of the atmospheric correction were processed for phase turn displacement. From the specific displacement results of Figure 9, we found that the local area showed signs of displacement. The original displacement results, linear model, GACOS, and MLP neural network model corrected displacement results are represented in Figure 9a-d, respectively. Meanwhile, the interferometric pairs 20210604-20210610 with severe atmospheric delay in the Baihetan reservoir area were corrected for atmospheric correction using the three methods mentioned above, and the results of the atmospheric correction were processed for phase turn displacement. From the specific displacement results of Figure 9, we found that the local area showed signs of displacement. The original displacement results, linear model, GACOS, and MLP neural network model corrected displacement results are represented in Figure 9a-d, respectively.
Meanwhile, the interferometric pairs 20210604-20210610 with severe atmospheric delay in the Baihetan reservoir area were corrected for atmospheric correction using the three methods mentioned above, and the results of the atmospheric correction were processed for phase turn displacement. From the specific displacement results of Figure 9, we found that the local area showed signs of displacement. The original displacement results, linear model, GACOS, and MLP neural network model corrected displacement results are represented in Figure 9a-d, respectively.  Figure 9a showed that the atmospheric delay effect seriously influenced the original displacement results, and the displacement location and displacement magnitude could not be accurately identified. The linear model correction was not effective, and the atmospheric phase remained serious after the correction (Figure 9b). From Figure 9c, it can be seen that the GACOS correction had a certain effect but failed to accurately identify the location of the landslide. Figure 9d reveals that the MLP neural network model had the optimal correction effect. Moreover, with the exception of the obvious displacement in the  Figure 9a showed that the atmospheric delay effect seriously influenced the original displacement results, and the displacement location and displacement magnitude could not be accurately identified. The linear model correction was not effective, and the atmospheric phase remained serious after the correction (Figure 9b). From Figure 9c, it can be seen that the GACOS correction had a certain effect but failed to accurately identify the location of the landslide. Figure 9d reveals that the MLP neural network model had the optimal correction effect. Moreover, with the exception of the obvious displacement in the slope ranges, the displacement volume in the non-displacement area was approximately 0 mm, allowing for a more accurate analysis of the displacement area.
The further analysis of the displacement area showed that the area is located near Dapingzi village, Nagu Town, Huize County, Qujing City, Yunnan Province, with a central longitude and latitude of 103 • 4 24 and 26 • 31 3 , respectively. The elevation of the slope is approximately 1160-1492 m, with a differential height of 332 m, and the slope is relatively steep. Based on the displacement monitoring results shown in Figure 10a, it is evident that the area of strong displacement in this region was more apparent (as indicated by the blue dotted line), with the maximum line-of-sight displacement reaching approximately 40 mm/year. Figure 10b shows the remote sensing interpretation of the displacement area with clear displacement signs. This active slope with displacement was also validated by the SBAS-InSAR result and our fieldwork.
Among the advantages of GACOS are its versatility and global availability. GACOS has advantages in dealing with long-wave atmospheric delay errors and may have some limitations in local applications with large undulating mountains. Therefore, atmospheric errors can be further and more effectively removed in local applications with the help of more refined models and methods (such as the deep learning methods proposed in this paper). It should be mentioned that in order to maintain a high coherence of interferometry, the temporal baseline range was within one month in this study. The effect of the temporal atmospheric water content variations from wet/dry seasons should be evaluated if interferograms with long temporal baselines were used.
slope is approximately 1160-1492 m, with a differential height of 332 m, and the slope is relatively steep. Based on the displacement monitoring results shown in Figure 10a, it is evident that the area of strong displacement in this region was more apparent (as indicated by the blue dotted line), with the maximum line-of-sight displacement reaching approximately 40 mm/year. Figure 10b shows the remote sensing interpretation of the displacement area with clear displacement signs. This active slope with displacement was also validated by the SBAS-InSAR result and our fieldwork. Among the advantages of GACOS are its versatility and global availability. GACOS has advantages in dealing with long-wave atmospheric delay errors and may have some limitations in local applications with large undulating mountains. Therefore, atmospheric errors can be further and more effectively removed in local applications with the help of more refined models and methods (such as the deep learning methods proposed in this paper). It should be mentioned that in order to maintain a high coherence of interferometry, the temporal baseline range was within one month in this study. The effect of the temporal atmospheric water content variations from wet/dry seasons should be evaluated if interferograms with long temporal baselines were used.

Conclusions
This study proposed a topography-dependent atmospheric correction method based on the MLP neural network model integrated with topographic and spatial data. The method was used to correct the atmosphere of an interferometric pair composed of Sentinel-1 ascending and descending orbit SAR in 17 views of the Baihetan dam. Moreover, the correction results were compared with the traditional linear model and GACOS to evaluate the effect of correction and the rate of change of the phase standard deviation of the three models.
The results demonstrated that the MLP neural network method could effectively estimate and correct the topography-dependent atmospheric delay, and the corrected Sentinel-1 wide-field InSAR interferogram phase standard deviation mean reduction showed a 64% improvement compared to that of the traditional linear model and GACOS. Furthermore, through fitting the phase-elevation linear relationship, it was clear that the linear coefficient of the phase and elevation correlation after the MLP neural network model

Conclusions
This study proposed a topography-dependent atmospheric correction method based on the MLP neural network model integrated with topographic and spatial data. The method was used to correct the atmosphere of an interferometric pair composed of Sentinel-1 ascending and descending orbit SAR in 17 views of the Baihetan dam. Moreover, the correction results were compared with the traditional linear model and GACOS to evaluate the effect of correction and the rate of change of the phase standard deviation of the three models.
The results demonstrated that the MLP neural network method could effectively estimate and correct the topography-dependent atmospheric delay, and the corrected Sentinel-1 wide-field InSAR interferogram phase standard deviation mean reduction showed a 64% improvement compared to that of the traditional linear model and GACOS. Furthermore, through fitting the phase-elevation linear relationship, it was clear that the linear coefficient of the phase and elevation correlation after the MLP neural network model correction was almost 0, which means that the phase and elevation are relatively uncorrelated. Eventually, two practical applications confirmed that the MLP neural network model could successfully identify two landslides in Sanjia and Dapingzi villages in the Baihetan reservoir area. The MLP neural network model proposed in this paper is more credible for topographydependent atmospheric correction results, which can achieve the high-accuracy application of InSAR technology in reservoir area displacement monitoring.