A Deep Learning Application for Deformation Prediction from Ground-Based InSAR

: Ground-based synthetic aperture radar interferometry (GB-InSAR) has the characteristics of high precision, high temporal resolution, and high spatial resolution, and is widely used in highwall deformation monitoring. The traditional GB-InSAR real-time processing method is to process the whole data set or group in time sequence. This type of method takes up a lot of computer memory, has low efﬁciency, cannot meet the timeliness of slope monitoring, and cannot perform deformation prediction and disaster warning forecasting. In response to this problem, this paper proposes a GB-InSAR time series processing method based on the LSTM (long short-term memory) model. First, according to the early monitoring data of GBSAR equipment, the time series InSAR method (PS-InSAR, SBAS, etc.) is used to obtain the initial deformation information. According to the deformation calculated in the previous stage and the atmospheric environmental parameters monitored, the LSTM model is used to predict the deformation and atmospheric delay at the next time. The phase is removed from the interference phase, and ﬁnally the residual phase is unwrapped using the spatial domain unwrapping algorithm to solve the residual deformation. The predicted deformation and the residual deformation are added to obtain the deformation amount at the current moment. This method only needs to process the difference map at the current moment, which greatly saves time series processing time and can realize the prediction of deformation variables. The reliability of the proposed method is veriﬁed by ground-based SAR monitoring data of the Guangyuan landslide in Sichuan Province.


Introduction
After 20 years of continuous development of ground-based SAR, many scholars and experts have conducted research on its application in different fields, so that ground-based SAR can continue to develop and be suitable for monitoring in different scenarios [1][2][3][4][5]. In 2018, S. Kuraoka successfully monitored the eruption process of Hakone volcano in Japan using the IBIS-L system, proving that GB-InSAR technology is a feasible method for monitoring volcanic activity [6]. In 2018, in order to study the relationship between the evolution of volcanic landforms and slope instability, Di Traglia et al. used the joint monitoring method of COSMO satellite and ground-based SAR to monitor volcanoes in a certain area of Italy. This method emphasized ground-based SAR equipment and spaceborne SAR. The feasibility of joint monitoring volcanic activity provides a feasible method for the subsequent study of volcanic evolution [7]. In 2019, T.carla et al. adopted the method of joint GNSS, spaceborne SAR and ground-based SAR for a landslide in the alpine region of Italy; the experiment proved the feasibility of the joint measurement of GNSS technology and InSAR technology, and has more advantages for GNSS technology in snow-covered areas, which provides ideas and methods for subsequent joint data and InSAR monitoring in alpine regions [8]. In 2021, T. Carla et al. used GB-InSAR technology to process and analyze the monitoring data of a large landslide in the Italian Alps for ten years, and is located on the top of the mountain in the north, with an altitude of 756 m, and the lowest point is located at the bottom of the mountain, with an altitude of 540m, and the height difference between the two is 216 m. The upper part of the landslide is steep, with a gradient of about 15 • to 30 • , and the lower part is relatively gentle, with a gradient of about 5 • to 20 • , and the overall main sliding direction is about 223 • . The landslide starts from the upper part of the hillside, the trailing edge of the landslide is bounded by bedrock, and the west and east sides are bounded by a little exposed bedrock. The possibility of the landslide continuing to expand cannot be ruled out. The landslide is about 851 m long, 120 m to 160 m wide, about 2 m thick on average, and has a distribution area of about 77,297 square meters. It is a large soil landslide. The landslide extends in the southwest direction as a whole, and the lateral extension is short. At present, the deformation of the landslide is mainly concentrated in the top trailing edge, local slow displacement occurs at the soil layer accumulation from the middle to the bottom, and the overall structure tends to be stable. If the landslide deforms significantly and this causes the landslide to become unstable, then the landslide will pose a threat to the residents in the danger zone and the three high-voltage power line towers. Contact monitoring is dangerous, and the satellite time baseline is long, which is not suitable for emergency monitoring. Therefore, GB-InSAR with short time baseline, high resolution and flexible observation angle is used for continuous monitoring. The location of the study area and the field equipment is shown in Figure 1.
Remote Sensing 2022, 14, x FOR PEER REVIEW 3 of 20 On October 7, 2021, a landslide occurred near Yanwo Village, Guangyuan City, Sichuan Province. The landslide was located near Dashiyan, Yanwo Village, Guangyuan City, specifically at 106°0′41.53″ east longitude, 32°22′28.449″ north latitude. The highest point is located on the top of the mountain in the north, with an altitude of 756m, and the lowest point is located at the bottom of the mountain, with an altitude of 540m, and the height difference between the two is 216m. The upper part of the landslide is steep, with a gradient of about 15° to 30°, and the lower part is relatively gentle, with a gradient of about 5° to 20°, and the overall main sliding direction is about 223°. The landslide starts from the upper part of the hillside, the trailing edge of the landslide is bounded by bedrock, and the west and east sides are bounded by a little exposed bedrock. The possibility of the landslide continuing to expand cannot be ruled out. The landslide is about 851m long, 120m to 160m wide, about 2m thick on average, and has a distribution area of about 77,297 square meters. It is a large soil landslide. The landslide extends in the southwest direction as a whole, and the lateral extension is short. At present, the deformation of the landslide is mainly concentrated in the top trailing edge, local slow displacement occurs at the soil layer accumulation from the middle to the bottom, and the overall structure tends to be stable. If the landslide deforms significantly and this causes the landslide to become unstable, then the landslide will pose a threat to the residents in the danger zone and the three high-voltage power line towers. Contact monitoring is dangerous, and the satellite time baseline is long, which is not suitable for emergency monitoring. Therefore, GB-InSAR with short time baseline, high resolution and flexible observation angle is used for continuous monitoring. The location of the study area and the field equipment is shown in Figure 1.  Fig. 1(a) is an orthophoto image processed by the DJI Royal 2 drone on October 30, 2021, showing the geographic location of the study area. P1, P2 and P3 are three study points selected from different locations of the landslide. P1 is the rear wall of the landslide, P2 is the landslide platform, and P3 is the secondary step. The Fig. 1(b) is a photo of the GPRI-II equipment and the landslide mass.
The GAMMA Portable Radar Interferometer (GPRI-II) was set up 500 meters south of the landslide. At the same time, the FSR-8 automatic weather station with the same time baseline was used to collect the necessary meteorological parameters for a total of 16 days, of which there was no interruption in the first 14 days. Day 15 was suspended for 8 hours due to power issues. The main parameters during measurement are show in table 1: The GAMMA Portable Radar Interferometer (GPRI-II) was set up 500 m south of the landslide. At the same time, the FSR-8 automatic weather station with the same time baseline was used to collect the necessary meteorological parameters for a total of 16 days, of which there was no interruption in the first 14 days. Day 15 was suspended for 8 h due to power issues. The main parameters during measurement are show in Table 1:

GB-InSAR Phase Model
After GBSAR was installed, the monitoring target was observed at a certain frequency (minute level). Due to the stable observation attitude, a zero baseline could be achieved. Compared with the spaceborne interference phase model, there is no reference plane phase and elevation phase, and the corresponding interference phase model can be represented by the following formula: Among them, ϕ is the interferometric phase of the wrapping, ϕ de f is the deformation corresponding to the radar line of sight, ϕ atm is the atmospheric delay phase in the radar line of sight, ϕ noise is the interference phase noise, and k is the number of integer cycles [19]. Among them, the atmospheric delay phase is the most important error source affecting the monitoring accuracy of GB-InSAR. The research shows that at a distance of 1 km from the GBSAR device, when the temperature is 20 • C, a 1% humidity change can lead to a measurement error of 2 mm [20,21].
According to the propagation mode of electromagnetic waves in the air, a mathematical model and a physical model were constructed for the atmospheric delay phase. The mathematical model is based on the propagation principle of electromagnetic waves. The radar emits electromagnetic waves with a wavelength of 1 at all times, and the electromagnetic waves reach the target pixel after the propagation distance return. The echo phase model for this pixel is [22]: n is the refractive index related to distance and time.
The physical model assumes that in the ground-based SAR monitoring, the spatial distribution of the atmospheric medium is assumed to be uniform within the monitoring area of 300-500 m, and the refractive index is independent of distance and only related to time, so the atmospheric phase of the pixel can be expressed as: In the formula, the relationship between the refractive index n and the refractive index N can be expressed as: In atmospheric refraction studies, the refractive index is usually expressed by temperature, air pressure, and relative humidity [23]: In the formula, h is the atmospheric pressure and e is the water vapor pressure, which can be expressed as [24] according to the Magnus saturated water vapor pressure formula: e sat = 6.1016 · 10 ( 7.5×T Through the above analysis, the range-to-homogeneous atmospheric phase during GBSAR monitoring can be modeled and eliminated from the interference phase, while reducing the difficulty of phase unwrapping in the space domain.

Long Short-Term Memory for Sequence Modeling
For the phenomenon of gradient disappearance and gradient explosion in RNN, Hochreiter proposed a network structure to improve the gradient disappearance problem, which works in a similar way to RNN. The difference is that a memory storage mechanism is established in some units within RNN to achieve more effective information storage and updates. Figure 2 shows the internal structure of the LSTM. By analyzing Figure 2, it can be seen that LSTM is obviously different from RNN in that it adds a gating unit, and the network transmits the gating unit information and the cell state at the last moment to the next unit in the network. 6.1016 10 Through the above analysis, the range-to-homogeneous atmospheric phase during GBSAR monitoring can be modeled and eliminated from the interference phase, while reducing the difficulty of phase unwrapping in the space domain.

Long Short-Term Memory for Sequence Modeling
For the phenomenon of gradient disappearance and gradient explosion in RNN, Hochreiter proposed a network structure to improve the gradient disappearance problem, which works in a similar way to RNN. The difference is that a memory storage mechanism is established in some units within RNN to achieve more effective information storage and updates. Figure 2 shows the internal structure of the LSTM. By analyzing Figure 2, it can be seen that LSTM is obviously different from RNN in that it adds a gating unit, and the network transmits the gating unit information and the cell state at the last moment to the next unit in the network. After the input data are processed by the LSTM model, the predicted value at the next moment is obtained. The LSTM model can be expressed as:

Proposed Methodology
Firstly, this paper uses the time-series InSAR technology to obtain the deformation variable data set, and organizes and models the meteorological parameter data to obtain the atmospheric phase data set. After that, the LSTM prediction model is trained with two After the input data are processed by the LSTM model, the predicted value at the next moment is obtained. The LSTM model can be expressed as:

Proposed Methodology
Firstly, this paper uses the time-series InSAR technology to obtain the deformation variable data set, and organizes and models the meteorological parameter data to obtain the atmospheric phase data set. After that, the LSTM prediction model is trained with two kinds of data sets, and the prediction accuracy of the model is evaluated by the prediction effect of the two kinds of data. Finally, the prediction model is used to realize the real-time processing of GB-InSAR monitoring and obtain the real-time processing results. The accuracy of the method is evaluated by the precision of inner coincidence and a computational efficiency evaluation. The technical roadmap for this paper is shown in the Figure 3 below: kinds of data sets, and the prediction accuracy of the model is evaluated by the prediction effect of the two kinds of data. Finally, the prediction model is used to realize the realtime processing of GB-InSAR monitoring and obtain the real-time processing results. The accuracy of the method is evaluated by the precision of inner coincidence and a computational efficiency evaluation. The technical roadmap for this paper is shown in the figure 3 below:

Data Normalization
In order to make the neural network model converge faster and obtain higher prediction accuracy, the variable values entering the neural network model should belong to the same range; that is, the value range of the data set should be consistent, and it is necessary to manually adjust the range of variables according to a certain ratio, the variable is scaled, and the following are two common processing methods: Normalization refers to the scaling of variable values in the range zero to one or negative one to one [24]. First, make the difference between the maximum value Xmax and the minimum value of the variable Xmin to obtain the variation range of the variable X, and then divide each observation Xi value in the variable X with the minimum value Xmin of the variable X and divide it by the variation range of the variable X to scale the observed value Xi to zero to one; the process is called maximum and minimum normalization, as shown in the following formula:

Data Normalization
In order to make the neural network model converge faster and obtain higher prediction accuracy, the variable values entering the neural network model should belong to the same range; that is, the value range of the data set should be consistent, and it is necessary to manually adjust the range of variables according to a certain ratio, the variable is scaled, and the following are two common processing methods: Normalization refers to the scaling of variable values in the range zero to one or negative one to one [24]. First, make the difference between the maximum value X max and the minimum value of the variable X min to obtain the variation range of the variable X, and then divide each observation X i value in the variable X with the minimum value X min of the variable X and divide it by the variation range of the variable X to scale the observed value X i to zero to one; the process is called maximum and minimum normalization, as shown in the following formula: The normalized atmospheric phase value and deformation phase value are selected and shown in Tables 2 and 3.

LSTM Network Design
The data set is divided into a training set, validation set and test set according to the ratio of 6:2:2. After the network design is completed, the training set can be input to train the model. In the training process of the model, due to the complex network structure of the neural network, as the number of iterations increases, the model is prone to overfitting. Therefore, in order to have a high generalization ability of the model and avoid overfitting, the model's dropout layer needs to be added in the LSTM layer. The main function of this layer is to discard the neurons in the LSTM layer according to a certain proportion in each training [25].
After experimental testing, the increase in the number of LSTM layers does not change the prediction accuracy of the model, and the selection of timesteps is equally important. Selecting timesteps is three, that is, predicting the deformation variables at the next moment according to the deformation variables of the first three moments, which corresponds to the experiment in this paper. That is, it is recommended to use the one-hour deformation data to predict the deformation data after 20 min as the output. The Units parameter in the LSTM layer represents the number of neurons in this layer, which is the tensor output dimension of the LSTM layer. Activation represents the activation function method. The ReLU activation function is selected in the LSTM layer, and the sigmoid activation function is selected in the Dense layer. This combination method enables the model to obtain higher accuracy. The dropout layer is set to zero point two, that is, twenty percent of neurons are randomly dropped during each training [26][27][28].

Hyperparameter Selection
For neural network models, there are many hyperparameters in the model, which control the structure, function, and efficiency of the model, and are the key to obtaining high-precision results for the model. The hyperparameter selection of the model in the experiment in this paper is adjusted using the grid search method. By setting the range of hyperparameters, the model is evaluated according to a certain growth rate of hyperparameters, so as to select the combination of hyperparameters with the best model performance.
In the evaluation index of the model, we use the loss rate to measure; the loss rate is mainly the calculation method of the MSE (mean square error), and its calculation method is shown in the following: In the training process of the model, the fitting degree of the model is judged by measuring the size of the loss function. When the loss rates of the training set and the validation set finally reach the same level, it proves that the model has achieved a good Remote Sens. 2022, 14, 5067 8 of 20 fitting effect. The combination of hyperparameters of the model determined using the grid search method is shown in Table 4: For deformable data, this paper conducts first-order difference to remove the autocorrelation of the data itself and then trains the model. For atmospheric phase data, since the data itself have periodic components, this paper uses the STL (seasonal trend loss) method to decompose the periodic term, trend term and residual term. After that, the three components are modeled separately. The predicted values of the three components are then added to obtain the overall predicted value. The training set and validation set loss rate curves of the model are shown in Figures 4 and 5: Remote Sensing 2022, 14, x FOR PEER REVIEW 9 of 20 validation set finally reach the same level, it proves that the model has achieved a good fitting effect. The combination of hyperparameters of the model determined using the grid search method is shown in Table 4: For deformable data, this paper conducts first-order difference to remove the autocorrelation of the data itself and then trains the model. For atmospheric phase data, since the data itself have periodic components, this paper uses the STL (seasonal trend loss) method to decompose the periodic term, trend term and residual term. After that, the three components are modeled separately. The predicted values of the three components are then added to obtain the overall predicted value. The training set and validation set loss rate curves of the model are shown in Figure 4 and Figure 5:   validation set finally reach the same level, it proves that the model has achieved a good fitting effect. The combination of hyperparameters of the model determined using the grid search method is shown in Table 4: For deformable data, this paper conducts first-order difference to remove the autocorrelation of the data itself and then trains the model. For atmospheric phase data, since the data itself have periodic components, this paper uses the STL (seasonal trend loss) method to decompose the periodic term, trend term and residual term. After that, the three components are modeled separately. The predicted values of the three components are then added to obtain the overall predicted value. The training set and validation set loss rate curves of the model are shown in Figure 4 and Figure 5:   and finally reaches a stable value near 0.000002, indicating that the prediction accuracy of the deformation prediction model is very high.

Obtaining the Initial Deformation Value
PS-InSAR technology is a method of selecting PS (Permanent Scatterer) point targets for time series analysis by statistically analyzing the SAR time series data in the study area according to the amplitude and phase. Different from the D-InSAR method [29], it only models and analyzes point targets with high coherence, strong scattering and stability, which effectively reduces the impact of decoherence on data processing. Sequence analysis can clearly obtain the change in the deformation of the target area over time, and more intuitively understand the change trend of the surface, which is helpful for the investigation and prevention of geological disasters. The technical process of PS-InSAR is shown in Figure 6 below: In the figure 5, the horizontal axis represents the epoch, and the vertical axis represents the size of the model loss. When the epoch reaches 500 times, the loss rate converges well, and finally reaches a stable value near 0.000002, indicating that the prediction accuracy of the deformation prediction model is very high.

Obtaining the Initial Deformation Value
PS-InSAR technology is a method of selecting PS (Permanent Scatterer) point targets for time series analysis by statistically analyzing the SAR time series data in the study area according to the amplitude and phase. Different from the D-InSAR method [29], it only models and analyzes point targets with high coherence, strong scattering and stability, which effectively reduces the impact of decoherence on data processing. Sequence analysis can clearly obtain the change in the deformation of the target area over time, and more intuitively understand the change trend of the surface, which is helpful for the investigation and prevention of geological disasters. The technical process of PS-InSAR is shown in Figure 6 below: In Figure 6, in the monitoring of GBSAR, the radar transmitter moves along a fixed orbit in each measurement. After the main image is selected, PS points are selected by setting thresholds for the phase dispersion index and coherence coefficient, and then the auxiliary image and the main image are subjected to interference processing, respectively, to construct a phase model function. Then, build a Delaunay triangulation for all adjacent PS points, use the solution space search method for the phase difference of PS points on the same baseline, iteratively solve the linear deformation rate, analyze the residual phase, and estimate the atmospheric phase component. After the atmospheric phase is separated from the phase, the PS point is reselected and the above operation is repeated to obtain the deformation time series results of the monitoring area. In Figure 6, in the monitoring of GBSAR, the radar transmitter moves along a fixed orbit in each measurement. After the main image is selected, PS points are selected by setting thresholds for the phase dispersion index and coherence coefficient, and then the auxiliary image and the main image are subjected to interference processing, respectively, to construct a phase model function. Then, build a Delaunay triangulation for all adjacent PS points, use the solution space search method for the phase difference of PS points on the same baseline, iteratively solve the linear deformation rate, analyze the residual phase, and estimate the atmospheric phase component. After the atmospheric phase is separated from the phase, the PS point is reselected and the above operation is repeated to obtain the deformation time series results of the monitoring area.

Real-Time Processing
For time t − 1, use the scene data before t − 1 to analyze the time series of the monitoring area, and obtain the deformation variable data within this time range, and use these data as the data set of the LSTM model to predict the deformation variable at time t. The meteorological parameters [30] collected by the weather station at the previous time t − 1 are used to fit and predict the atmospheric phase at time t, which can be: The formula represents the predicted phase at time t, including the predicted value of the atmospheric phase and the predicted value of the deformation phase. At time t, if you want to perform time series analysis, you need to remodel and solve the data at the previous time t, which will waste more time and resources. Interference is performed at time t, and the phase includes: The atmospheric phase change and deformation phase change at time t relative to time t − 1 are subtracted from the difference result between time t − 1 and time t. At this time, there are only residual atmospheric phase and residual deformation phase in the phase. bit and noise phase. At this time, the spatial domain unwrapping of the interference phase is performed after the spatial domain filtering, which is easier than the original differential phase, and a better unwrapping result can be obtained. After the unwrapped phase is obtained, according to the deformation result at time t − 1, a mask is generated to cover up the area with a large deformation amount. For the stable area, the residual phase is the atmospheric phase of the area. The residual atmospheric phase of the entire monitoring area is obtained by interpolation, and the residual deformation phase can be obtained by subtracting the phase from the previous step. After being converted into a deformation variable and added to the predicted deformation variable at time t, the predicted deformation variable at time t can be corrected, and the time for time series analysis can be saved in order to achieve the effect of real-time monitoring.

Results
In this section, the GBSAR monitoring data of Guangyuan City, Sichuan Province were used for experiments, and 20% of the total data were used to predict the atmospheric phase and shape variables. Finally, 40 views of GBSAR data from 14:39 on October 28 to 03:59 on October 29 were used for the real-time processing experiment. Figures 7 and 8, respectively, show the interferogram and the selected PS point diagram at 16:00 on the 17th and 16:00 on the 19th. Analysis of Figure 7 shows that the landslide body has undergone large deformation, and the coherence is poor in the vegetation coverage area.
The formula represents the predicted phase at time t, including the p of the atmospheric phase and the predicted value of the deformation pha you want to perform time series analysis, you need to remodel and solve previous time t, which will waste more time and resources. Interference i time t, and the phase includes: The atmospheric phase change and deformation phase change at tim time t-1 are subtracted from the difference result between time t-1 and time there are only residual atmospheric phase and residual deformation phas bit and noise phase. At this time, the spatial domain unwrapping of the inte is performed after the spatial domain filtering, which is easier than the orig phase, and a better unwrapping result can be obtained. After the unwra obtained, according to the deformation result at time t-1, a mask is genera the area with a large deformation amount. For the stable area, the residu atmospheric phase of the area. The residual atmospheric phase of the en area is obtained by interpolation, and the residual deformation phase can subtracting the phase from the previous step. After being converted into variable and added to the predicted deformation variable at time t, deformation variable at time t can be corrected, and the time for time seri be saved in order to achieve the effect of real-time monitoring.

Results
In this section, the GBSAR monitoring data of Guangyuan City, Sic were used for experiments, and 20% of the total data were used to predict t phase and shape variables. Finally, 40 views of GBSAR data from 14:39 on 03:59 on October 29 were used for the real-time processing experiment.

Atmosphere Prediction
In this paper, the monitoring area was monitored for 16 days of meteo parameters. By ensuring that the sampling frequency of GBSAR was consistent of GBSAR, the changes in atmospheric parameters during the observation per obtained. According to the atmospheric phase model of GBSAR, the atmosphe based on the distance direction was obtained. The data set was constructed using and the LSTM prediction result of the monitoring area was obtained by app LSTM network model, which is plotted in Figure 9.

Deformation Prediction
In this paper, a total of 1,398,120 PS points and the deformation amount of in the entire monitoring area were predicted, and the deformation prediction r partially shown in Figure 10 below according to a certain time interval. To difference, we drew the prediction residuals as shown in Figure 11 and Figur drew the prediction residuals in Figures 12 and 14. By analyzing the above resu be seen that the LSTM network model has high prediction accuracy in t

Atmosphere Prediction
In this paper, the monitoring area was monitored for 16 days of meteorological parameters. By ensuring that the sampling frequency of GBSAR was consistent with that of GBSAR, the changes in atmospheric parameters during the observation period were obtained. According to the atmospheric phase model of GBSAR, the atmospheric phase based on the distance direction was obtained. The data set was constructed using this data, and the LSTM prediction result of the monitoring area was obtained by applying the LSTM network model, which is plotted in Figure 9.

Atmosphere Prediction
In this paper, the monitoring area was monitored for 16 days of meteorological parameters. By ensuring that the sampling frequency of GBSAR was consistent with that of GBSAR, the changes in atmospheric parameters during the observation period were obtained. According to the atmospheric phase model of GBSAR, the atmospheric phase based on the distance direction was obtained. The data set was constructed using this data, and the LSTM prediction result of the monitoring area was obtained by applying the LSTM network model, which is plotted in Figure 9.

Deformation Prediction
In this paper, a total of 1,398,120 PS points and the deformation amount of 59 hours in the entire monitoring area were predicted, and the deformation prediction results are partially shown in Figure 10 below according to a certain time interval. To make a difference, we drew the prediction residuals as shown in Figure 11 and Figure 13, and drew the prediction residuals in Figures 12 and 14. By analyzing the above results, it can be seen that the LSTM network model has high prediction accuracy in the entire monitoring area, and the prediction results of the key deformation area and the stable area are in line with the facts.

Deformation Prediction
In this paper, a total of 1,398,120 PS points and the deformation amount of 59 h in the entire monitoring area were predicted, and the deformation prediction results are partially shown in Figure 10 below according to a certain time interval. To make a difference, we drew the prediction residuals as shown in Figures 11 and 12, and drew the prediction residuals in Figures 13 and 14. By analyzing the above results, it can be seen that the LSTM network model has high prediction accuracy in the entire monitoring area, and the prediction results of the key deformation area and the stable area are in line with the facts.      was generated between the date that is marked in the low right corner and 1 acquisition date of GBSAR image). The deformation of the landslide is obvio from -104mm to 135mm.    Remote Sensing 2022, 14, x FOR PEER REVIEW Figure 10. Historical results of deformation prediction along LOS direction was generated between the date that is marked in the low right corner and 1 acquisition date of GBSAR image). The deformation of the landslide is obvio from -104mm to 135mm.    Figure 10. Historical results of deformation prediction along LOS direction. The displacement map was generated between the date that is marked in the low right corner and 17 October 2021 (the first acquisition date of GBSAR image). The deformation of the landslide is obvious. The color bar ranges from -104mm to 135mm.
plot of real value and predicted lue.   igure 10. Historical results of deformation prediction along LOS direction. The displacement map as generated between the date that is marked in the low right corner and 17 October 2021 (the first cquisition date of GBSAR image). The deformation of the landslide is obvious. The color bar ranges om -104mm to 135mm.
lot of real value and predicted e.

Real-Time Processing
This experiment uses the ground-based SAR data of 40 scenes from 14:39 on 28 October to 03:59 on 29 October. First, the difference results between time t − 1 and time t are obtained as shown in Figure 15. From the difference map, it can be clearly seen that there is an atmospheric phase with a trend term, and there is relatively large atmospheric disturbance in the area farthest from the radar center. The atmospheric phase and deformation phase predicted by the LSTM model are subtracted from the interference phase, and the results are shown in Figure 16. Analysis of the results shows that when fitting the atmospheric phase in the horizontal direction in this area, there is still atmospheric turbulence in the elevation direction and in a small area. For this part of the atmospheric phase, the method of mask interpolation after unwrapping is used to remove it. For the phase residuals of adjacent moments, spatial domain filtering is first performed, and the filtering results are shown in Figure 17. The filtering will remove the noise interference in the residuals. After that, the minimum cost flow method is used to unwrap, and the unwrapping results are shown in Figure 18. After unwrapping, the atmospheric phase of the trend term can be clearly seen, and then the stable region is interpolated after masking, and the fitted atmospheric phase is shown in Figure 19. After subtracting the residual atmospheric phase from the residual, the residual deformation is obtained as shown in Figure 20. Analysis of Figure 18 shows that the residual deformation of most regions is predicted to be zero, which proves that the LSTM model can be used to complete large regions with high accuracy. In the deformation prediction, the real-time processing deformation variable at time t is obtained by adding the deformation residual to the predicted value at the end.
shown in Figure 20. Analysis of Figure 18 shows that the regions is predicted to be zero, which proves that the LSTM large regions with high accuracy. In the deformation predi deformation variable at time t is obtained by adding th predicted value at the end.  mospheric phase with a trend term, and there is relatively large atmospheric nce in the area farthest from the radar center. The atmospheric phase and tion phase predicted by the LSTM model are subtracted from the interference nd the results are shown in Figure 16. Analysis of the results shows that when the atmospheric phase in the horizontal direction in this area, there is still heric turbulence in the elevation direction and in a small area. For this part of the heric phase, the method of mask interpolation after unwrapping is used to remove e phase residuals of adjacent moments, spatial domain filtering is first performed, filtering results are shown in Figure 17. The filtering will remove the noise ence in the residuals. After that, the minimum cost flow method is used to unwrap, unwrapping results are shown in Figure 18. After unwrapping, the atmospheric f the trend term can be clearly seen, and then the stable region is interpolated after , and the fitted atmospheric phase is shown in Figure 19. After subtracting the l atmospheric phase from the residual, the residual deformation is obtained as in Figure 20. Analysis of Figure 18 shows that the residual deformation of most is predicted to be zero, which proves that the LSTM model can be used to complete gions with high accuracy. In the deformation prediction, the real-time processing tion variable at time t is obtained by adding the deformation residual to the d value at the end.
ges. The interference , atmospheric phase, le.     diagram. In the study ing the shape variable he stable region is ospheric delay phase. processing 40 scene data to simulate real-time processing in the monitoring ment, the final deformation rate graph is obtained as shown in Figure 21:   Interpolated atmospheric phase diagram. In the study area, the stable region is selected by taking the shape variable as the threshold, and the phase of the stable region is interpolated to obtain the overall atmospheric delay phase. By processing 40 scene data to simulate real-time environment, the final deformation rate graph is obtained Figure 19. Interpolated atmospheric phase diagram. In the study area, the stable region is selected by taking the shape variable as the threshold, and the phase of the stable region is interpolated to obtain the overall atmospheric delay phase.

of 20
ing. The residual phase ing algorithm, and the ined. e diagram. In the study king the shape variable the stable region is ospheric delay phase. processing 40 scene data to simulate real-time processing in the monitoring By processing 40 scene data to simulate real-time processing in the monitoring environment, the final deformation rate graph is obtained as shown in Figure 21: the threshold, and the phase of the stable region is olated to obtain the overall atmospheric delay phase. deformation phase can be interpolated atmospher By processing 40 scene data to simulate real-time environment, the final deformation rate graph is obtained  Figure 22. Real deformation curves of three research points. P1 is located at the crown crack located at the head, and P3 is located at the minor scarp. According to the analysis cur deformation rate of the three points is different, but all of them keep rising. Figure 22. Real deformation curves of three research points. P1 is located at the crown cracks, P2 is located at the head, and P3 is located at the minor scarp. According to the analysis curve, the deformation rate of the three points is different, but all of them keep rising. Figure 22. Real deformation curves of three research points. P1 is located at the crown cracks, P2 located at the head, and P3 is located at the minor scarp. According to the analysis curve, t deformation rate of the three points is different, but all of them keep rising. By analyzing the change curves of the three points in Figure 23, the time series resu obtained by this real-time processing method are consistent with the GB-InSAR tim series analysis results, which are in line with the real deformation variable changes, whi proves that this method can realize real-time deformation processing. The process is follows: select the deformation variable in the stable region for statistical analysis, a draw the frequency statistical histogram of the deformation variable in this region shown in Figure 24: By analyzing the change curves of the three points in Figure 23, the time series results obtained by this real-time processing method are consistent with the GB-InSAR time series analysis results, which are in line with the real deformation variable changes, which proves that this method can realize real-time deformation processing. The process is as follows: select the deformation variable in the stable region for statistical analysis, and draw the frequency statistical histogram of the deformation variable in this region as shown in    Generally speaking, the GB-INSAR deformation results need to be obtained through timing processing. Timing processing requires modeling and analysis of the entire data, using PS-InSAR or SBAS-InSAR methods. Due to the large amount of data in this method, the processing needs a lot of time, and the data set needs to be recalculated every time an additional scene is added. In this paper, the time of sequence analysis of 40 scenes of GBSAR data is counted, and the computer configuration used in the experiment is shown in Table 5. The time consumed by the real-time processing method proposed in this paper is plotted in Table 6. According to Table 5, GB-InSAR real-time monitoring based on LSTM network model has high processing efficiency, which can realize real-time monitoring and obtain highprecision results in slope monitoring. The average time of obtain deformation results by the real-time processing method is 4.73 min, which is much faster than the ordinary time series processing method. It is proved that the method can obtain results quickly and play a greater role in site landslide monitoring.

Discussion
The results show that the application of the LSTM model can effectively predict the deformation variables of the landslide in Guangyuan City, Sichuan Province. By evaluating the accuracy of the prediction results, the results show that the LSTM model can provide accurate prediction results. Compared with the linear prediction model, using the LSTM model for prediction can mine the correlation between historical data and make more accurate predictions for future data; compared with the machine learning model, it has the advantages of faster model convergence and higher prediction accuracy. At the same time, the distance from the study area to the atmosphere was predicted, and good prediction results were also obtained. GB-InSAR technology has the advantages of high precision, high spatial and temporal resolution, and flexible observation. The combination of the LSTM model and GB-InSAR provides more means for ground observation.
Based on the LSTM model, the GB-InSAR real-time processing method proposed in this paper can complete the extraction of deformation information in a short time. In the case of fewer computer resources and shorter time consumption, GB-INSAR time series processing is completed, and the time series results of deformation variables are obtained. The method proposed in this paper can provide case reference and verification for subsequent GB-InSAR monitoring technology and data-processing methods.

Conclusions
As a non-contact measurement technology, GB-InSAR technology can measure landslides with high precision and high spatial and temporal resolution. Compared with SAR satellites, ground-based SAR has a short time baseline, no spatial baseline, high resolution, and high-resolution observation. The advantages of flexible angle are widely used in landslide monitoring. Based on this, this paper proposes a GB-InSAR timing processing method based on the LSTM network model, aiming at the importance of deformation prediction in GBSAR monitoring, and the shortcomings of GB-InSAR timing analysis, such as slow processing efficiency, high memory consumption and inability to carry out realtime monitoring. The experimental results show that the method can quickly obtain the deformation prediction results in the monitoring area, and the inner coincidence accuracy evaluation proves that the method has high deformation prediction accuracy. At the same time, the real-time processing module can obtain the shape variables in real time, and the hardware occupancy rate is low, which proves that the method has high-precision real-time processing ability. The following conclusions are drawn through the experiment: (1) The LSTM model has high deformation prediction accuracy for the deformation results of GB-InSAR, and the standard deviation of the deformation prediction residual is 0.46 mm, which has the advantages of high prediction accuracy and fast model convergence speed.
(2) The proposed method can quickly obtain deformation results without occupying a large amount of memory, and the real-time processing time of a single scene is better than 5 min.
(3) Compared with the traditional GB-InSAR time series processing method, the standard deviation of the proposed method in the slow deformation region is between 0.1-0.3 mm, and the average deviation from the traditional method is within 0.13 mm.