A Three-Dimensional Deformation Monitoring Method: Combining Optical Deformation Monitoring Based on Regression Models and GB-SAR Interferometry

This paper introduces a method for quantifying the three-dimensional deformation of ground targets and outlines the associated process. Initially, ground-based synthetic aperture radar was employed to monitor the radial deformation of targets, and optical equipment monitored pixel-level deformation in the vertical plane of the line of sight. Subsequently, a regression model was established to transform pixel-level deformation into two-dimensional deformation based on a fundamental length unit, and the radar deformation monitoring data were merged with the optical deformation monitoring data. Finally, the fused data underwent deformation, resulting in a comprehensive three-dimensional deformation profile of the target. Through physical data acquisition experiments, the comprehensive three-dimensional deformation of targets was obtained and compared with the actual deformations. The experimental results show that the method has a relative error of less than 10%, and monitoring accuracy is achieved at the millimeter level.


Introduction
In recent years, frequent geological disasters have arisen from the disruption of geological structures triggered by crustal movements such as volcanoes and earthquakes, shifts in climate patterns due to global warming, and human activities such as mineral extraction and deforestation [1,2].Beyond natural calamities, deformation issues affect structures such as buildings.Construction facilities such as bridges, dams, and high-rise buildings commonly undergo specific deformations and settlements throughout their life cycles, sometimes leading to severe incidents, e.g., collapse [3].
Geological disasters typically unfold gradually, with deformation accumulating over the long term rather than occurring suddenly [4].Preceding major disaster events, potential hazards often undergo gradual and incremental deformations, accumulating slowly until instability arises [5,6].Identifying these subtle deformation processes in a timely manner and issuing early warnings can enable targeted disaster mitigation, effectively preventing casualties and economic losses [7].Therefore, it is necessary to carry out the long-term, continuous, and real-time monitoring of slopes and mines, dilapidated buildings, and large buildings with landslide and collapse risks, evaluate their stability, and make early warnings if necessary [8].
Differential interferometric synthetic aperture radar (DInSAR) has been a key player in surface monitoring, boasting high precision, extensive coverage, and all-weather capabilities since its inception.However, challenges persist for space-based SAR in terms of measurement accuracy and resolution, attributed to issues such as time decorrelation, long spatial baselines, and constraints on electromagnetic wave irradiation angles [9,10].Ground-based synthetic aperture radar (GB-SAR), following in the footsteps of space-based and airborne SAR, has gained widespread adoption for high-precision deformation monitoring in areas affected by natural disasters, dams, alpine glaciers, and pier structures.This popularity is attributed to its advantages, including a short space-time baseline, high resolution, a brief measurement period, and ease of operation [11][12][13].
As an emerging technology, GB-SAR provides detailed and accurate measurements of ground deformation [14,15].Despite the significant advancements in GB-SAR research in recent years, challenges persist [16].Notably, GB-SAR's limitation in capturing only the projection component of the target's deformation variable along the radar radial direction can result in underestimated deformations.This limitation becomes particularly evident when a substantial angle exists between the target's deformation direction and the radar radial direction, causing the projection of even significant deformations to appear very small in GB-SAR [17,18].This discrepancy poses the risk of missed alarms and potential hazards [19].In order to address this issue with GB-SAR, some scholars have attempted multi-station monitoring to capture the multi-dimensional deformation of a target.However, this approach introduces increased equipment costs and deployment challenges.Moreover, the current accuracy of this method still remains at the centimeter level [20].
In addition to GB-SAR deformation monitoring, optical imaging deformation technology finds widespread applications in various fields, including materials science and medicine [21].Optical deformation monitoring captures relevant information in the vertical plane of line of sight (LOS).We aim to leverage optical images to complement the third-dimensional information not attainable through radar monitoring.This approach seeks to achieve comprehensive three-dimensional deformation analysis for ground targets, minimizing the likelihood of missed alarms during monitoring.Furthermore, it provides essential technical support for natural disaster warnings and building deformation monitoring [22,23].
In this paper, a technique for obtaining the comprehensive three-dimensional deformation of ground targets is proposed.Firstly, based on GB-SAR and optical images, radar deformation monitoring, and optical deformation monitoring were carried out, respectively, and a nonlinear fitting regression model was established.The pixel-level deformation obtained by optical deformation monitoring was converted into deformation under the basic length unit, and the deformation monitoring data of the two were fused to solve the comprehensive three-dimensional deformation of targets finally.
Key innovations in this paper include the following: 1.
Proposing a method to convert pixel deformation quantity from optical deformation monitoring into deformation quantity under a basic length unit (mm).

2.
Conducting three-dimensional deformation inversion by integrating optical and radar deformation information.

3.
Establishing a three-dimensional solution model and validating the feasibility and superiority of three-dimensional deformation inversion through measured data.In single-station monitoring scenarios, the accuracy of deformation monitoring within tens of meters reaches the millimeter level, with a controlled relative error below 10%.
The paper is organized as follows: Section 2 delves into the principles of deformation monitoring using GB-SAR interferometry measurements.In Section 3, we present a comprehensive overview of the methods proposed in this article.This includes a novel optical deformation monitoring technology based on fitting a regression model, along with a three-dimensional deformation inversion technology that integrates optical and radar images.Section 4 validates the feasibility and effectiveness of the proposed method through experimental testing.In conclusion, Section 5 provides a comprehensive summary and the concluding remarks of the paper.

GB-SAR Deformation Monitoring
The image obtained by GB-SAR monitoring is a single-look-complex (SLC), and its value is distributed in [−π, π].For N radar SLC images in a time series, differential interference processing is performed on each target point at different times: ∆φ n = angle(s i e jφ i (s j e jφ j ) * ) = angle(s i s j e j(φ i −φ j ) ) ( where ∆φ n is the interference phase of the n point, P i = s i e jφ i is the i radar echo image, P j = s j e jφ j is the j radar echo image, and the symbol * is conjugate.The interference phase of the n target point in the differential interferogram can be modeled as where ∆φ dis is the LOS variable, ∆φ atm is the atmospheric delay phase introduced by the change in atmospheric refractive index at different times, ∆φ noise is the phase of incoherent noise, which can be filtered by using a low-pass phase filter.2kπ is the fuzzy phase, and k is an integer, which represents the ambiguity [24].However, it should be noted that in this experiment, we assume that the error phase, such as the atmospheric phase, has little effect, and our deformation monitoring target is a strong scatterer, which is less affected by the error.Therefore, the method proposed in this paper does not consider the influence of atmospheric phase error and random noise error and directly solves the deformation through the interference phase obtained by GB-SAR [25]. GB-SAR deformation monitoring technology uses differential interferometry.This technique fuses and analyzes two SAR complex images obtained by radar in the same target area.By calculating the target phase information obtained by the radar at different times, the millimeter-level precision deformation information of the target is obtained.Each pixel in the image obtained by GB-SAR is a complex number, and its amplitude is usually used to interpret the imaging scene and study the scattering characteristics.When GB-SAR is used for deformation measurements, the radar position is fixed, and the spatial baseline between different images is zero.Differential interferometry can be achieved by using complex conjugate multiplication on the corresponding pixels of two images.GB-SAR is mainly used to monitor the same area on parallel orbits and obtain two (or more) SLC images to form interference [26].
Figure 1 depicts the schematic diagram of GB-SAR differential interference deformation monitoring.Target 1 denotes the initial position of the target, whereas Target 2 corresponds to its position after displacement.The radar performs measurements on these targets both before and after deformation, capturing their respective phase information.Subsequently, the phase difference between the two is determined: By leveraging the correlation between phase difference and distance, we can extract the deformation information of the target projected in the LOS direction of the GB-SAR: where ∆φ is the phase difference, λ is the wavelength of the radar, and ∆d is the deformation information of the target projected in the LOS direction.
Principle of differential interferometry.

Method Flowchart
Since optical deformation monitoring can provide information within the plane perpendicular to the line of sight, we aim to utilize optical images to supplement the three-dimensional deformation information that cannot be obtained through radar monitoring.The workflow of the proposed method is shown in Figure 2. The main steps include radial deformation monitoring based on GB-SAR, optical deformation monitoring, pixel deformation conversion to deformation under the basic length unit, data information fusion, and three-dimensional deformation calculation.

Optical Deformation Monitoring
The deformation monitoring technology based on optical images is a method that can be used to monitor the deformation of objects by using optical measurement principles and imaging technology [21].This method obtains two or more optical images before and after the deformation of the target and then analyzes the information in the image to obtain the deformation information.Deformation monitoring technology based on optical images is widely used in materials science, physics, medicine, and other fields.Optical images are used to measure the expansion, bending, and torsion of materials.The digital image correlation method is one of the most commonly used optical deformation monitoring methods.This method uses two digital images before and after the deformation of the specimen to obtain the deformation information of the region of interest through correlation calculation [27].The basic principle is that the region of interest in the image before deformation is meshed, and the divided sub-regions are regarded as rigid motion.Then, for the selected sub-region, the correlation calculation is carried out according to the pre-defined correlation function by using a certain method, and the region with the maximum cross-correlation coefficient with the sub-region is found in the deformed image, that is, the position of the sub-region after deformation, and then the displacement of the sub-region is obtained.This method has the advantages of a strong anti-interference ability and high measurement accuracy [28].
The goal of correlation calculation is to assess the similarity between the template and the target image, underscoring the significance of choosing an appropriate correlation function in the digital image correlation method.Common correlation functions include direct correlation, normalized correlation, and covariance correlation functions.The expression for the correlation function is defined as follows: where H(X) represents the information entropy of picture X; the expression used for its calculation is where s i represents the total number of pixels in image X with a grayscale value of i, N represents the total number of grayscale levels in image X, and p i represents the probability of grayscale i appearing.H(X, Y) represents the joint entropy of the two images, X and Y: and I(X, Y) represents the mutual information of two pictures: X and Y.Given the suboptimal anti-noise performance of the direct correlation function, the commonly employed alternatives are the normalized correlation function and the covariance correlation function: Utilizing the correlation function calculation to identify the position with the maximum correlation coefficient allowed us to determine the initial location of the target in the optical image prior to any target displacement.This enables the calculation of the pixel-level deformation of the target in the two-dimensional optical plane.

The Proposed Idea
The optical deformation monitoring methods represented by the digital image correlation method have many advantages.However, these optical deformation monitoring methods can only obtain pixel-level deformation, which is widely used in materials science, physics, and other fields.In these fields, it is usually only necessary to qualitatively analyze the deformation of the target or to use the size of the target as a priori known information so as to roughly convert the pixel deformation into the actual deformation.Whether it is a qualitative analysis using pixel deformation quantity or rough conversion using prior information, these methods cannot be applied to physical target deformation monitoring in actual scenes.Therefore, how to apply optical deformation monitoring to the actual deformation monitoring field and convert the pixel deformation measured at any position into actual deformation becomes a difficult problem.
In this paper, a method based on regression fitting is proposed to convert the pixel deformation quantity obtained by optical deformation monitoring into the deformation quantity under the basic length unit.By fitting the optical sensor's optical deformation and monitoring the data of a series of standard samples with a known detection distance, the ratio of pixel deformation at different positions in the fitting range to the deformation under the basic length unit is predicted and obtained.The basic process of the method is shown in Figure 3: Among them, data acquisition is mainly used to collect the diagonal pixel values of standard samples of the same size at different positions.Data preprocessing mainly includes steps such as image interpolation and edge detection.The prediction result is to obtain the predicted diagonal pixel values of standard samples at different positions based on the fitting model so as to calculate the pixel conversion ratio and finally convert the pixel-level optical deformation into the deformation value under the basic length unit.

Construction Model
The nonlinear transformation observed in camera imaging primarily stems from the multiple refractions occurring within the lens system.Among the components of the imaging system, the lens plays a pivotal role, and the optical transformation within a single lens involves phase transformation: It can be seen that the influence of a single lens on optics is an exponential function relationship.Therefore, when nonlinearly fitting camera imaging, an exponential function model should be established: where R is the object distance, and P is the diagonal pixel of the standard sample when the detection distance is R. Since a single lens usually contains multiple sets of lenses in actual imaging, in order to obtain better fitting prediction results, the constructed nonlinear fitting model should contain at least two sets of exponential functions: After obtaining the diagonal pixels of the standard sample at the detection distance R i , the pixel conversion ratio r can be calculated here.The expression of the pixel conversion ratio is defined as where both D and P denote the diagonal dimensions of the plane containing a set of identical standard samples at a detection distance, R i , where D is measured in millimeters (mm), and P is measured in pixels.For the purposes of this study, the standard samples are a set of optical calibration boards with dimensions of 80 mm by 80 mm.Consequently, the diagonal size (D) of this standard sample can be calculated using the formula D = 80 2 + 80 2 mm.For the acquisition of detection distance, the ranging function of GB-SAR can be used.In the actual deformation monitoring scene, GB-SAR can directly read the detection distance of the corresponding position after monitoring the deformation point.When the optical sensor is placed in the radar center of GB-SAR, the detection distance obtained by GB-SAR is R i .
Through the nonlinear fitting model, the pixel P of the standard sample at any R = R i position can be obtained by using Formula (13), and then the pixel conversion ratio can be obtained by using Formula ( 14) so as to realize the function of converting the pixel deformation obtained in Section 3.3.1 based on an optical image under deformation and the basic length.

Image Fusion
Compared with spaceborne SAR imaging, GB-SAR imaging has the problem of a large detection range span, so there is serious spatial variation in radar signal processing, and there is obvious geometric distortion in its imaging [29].Consequently, the conventional approach of fusing spaceborne SAR and optical images is not applicable to GB-SAR and optical image fusion.This paper presents a method for combining GB-SAR and optical images through the utilization of co-ordinate system transformation, achieved by employing transformation matrices and control points [30].Firstly, there is a target point, P 0 , in the actual scene, and its co-ordinates are (x 0 , y 0 , z 0 ).In the optical image, the point corresponding to P 0 is P 1 , and its co-ordinates are (x 1 , y 1 ), as shown in Figure 4.According to the lens imaging principle of optical cameras, it can be seen that where f is the focal length of the camera; the transformation relationship between the actual scene and the optical image can be obtained: where matrix A can be proved to be an invertible matrix.The corresponding point of the target point P 0 on the GB-SAR image is P 2 , and the co-ordinate under the rectangular co-ordinate system is (x 2 , y 2 ).When the azimuth resolution of the GB-SAR is ρ a , and the range resolution is ρ r , the transformation relationship between the actual scene and the optical image can be obtained: It can also be proved that matrix B is invertible; thus Matrix C is the transformation matrix, which can be calculated by control points.In practical data collection scenarios, typically four-six corner reflectors are employed as predefined control points.Following the acquisition of both optical and GB-SAR images, the co-ordinate information of these control points can be manually extracted from the images.

Construct a Three-Dimensional Solution Model
After obtaining the radial deformation and two-dimensional plane deformation of the target, the actual three-dimensional deformation of the target can be solved by establishing a three-dimensional solution model.A basic three-dimensional deformation solution model is shown in Figure 5. where δ is the pitch angle of the target, and ε is the pitch angle of GB-SAR LOS, projecting the target deformation and radar LOS to the plane; in the horizontal plane, α is the angle between the radar radial direction and the target deformation direction, β is the target azimuth, and γ is the azimuth angle of GB-SAR.The three-dimensional deformation solution can be realized by using the following formula: where dx, dz, and dl are the monitoring results for the X-axis, Z-axis, and radar radial direction, respectively, and dt is the real deformation of the target.Finally, the radial deformation measured by GB-SAR and the plane deformation measured by the optical camera are substituted into the model to obtain the complete deformation, dt.

Experiments and Analysis
In order to demonstrate the accuracy and efficacy of the method detailed in this paper, we carried out experiments outdoors to obtain measured data for verification.Specifically, Section 4.1 examines the pixel conversion method's effectiveness, Section 4.2 examines the impact of image fusion, and Section 4.3 shows the obtainment of the comprehensive three-dimensional deformation quantity of the target.

Optical Deformation Monitoring Based on Fitting Regression Models
Firstly, the fitting parameters of Formula (13) were solved by practical experiments.In order to obtain the detection range of the targets, the experimental system includes an X-band GB-SAR system and an optical camera, as shown in Figure 6.The parameters of the system are shown in Table 1.The lens of the optical camera is placed on the center of the radar line of sight of the GB-SAR so that the GB-SAR can obtain information such as the detection range and azimuth angle of the targets.The experimental scene is shown in Figure 7. Identical optical calibration plates are placed in the experimental site as experimental targets.The pattern size of the optical calibration plate is 80 × 80 mm, and its diagonal size under the basic length unit is D = 113.12mm.The detection distance of these optical calibration plates is collected for fitting parameters.The detection distance and diagonal pixels of the six calibration boards are obtained, as shown in Table 2.
Repeat the experiment by changing the position of the calibration plate to obtain the detection distance and diagonal pixels of the six calibration plates and compare them with the fitting prediction data of Equation (24), as shown in Table 3 : The root mean square error (RMSE) and mean relative error (MRE) of the model are calculated, as shown in Table 4.It can be seen that the nonlinear fitting model can be reused after a parameter solution.The root mean square error is 0.9301 mm, and the average relative error can be controlled within 3%.It shows that the prediction effect of the fitting model is good.In order to verify the effect of the pixel conversion method in the actual optical deformation monitoring, a corner reflector is used as the target to move horizontally by 3.0 mm in the experimental site to obtain two images before and after deformation.
After three interpolation runs of the image, the pixel deformation of the target obtained by optical deformation monitoring is six pixels, as shown in Figure 8.The position information of the target obtained by using GB-SAR is = 9.82 m, as shown in Figure 9.
Based on the nonlinear fitting model of Section 3.3, the diagonal pixels of the standard samples at this position are obtained as follows: f (P) = 756.3e (−0.5089R) + 109.3e (−0.04932R) = 72.45(25) Thus, the original pixel conversion ratio here is 113.12/72.45= 1.56.Due to the interpolation of the optical monitoring image by a factor of three, the pixel conversion ratio in the image is r = 1.56/3 = 0.52.Consequently, the pixel deformation quantity obtained from optical deformation monitoring is converted into the deformation quantity under the basic length unit as 6 × 0.52 = 3.12 mm.The relative error between the result and the actual movement of the target is 4%.
In order to better test the performance of the method, we changed the deformation distance of the target, repeated the experiment, and obtained the results shown in Table 5: Azimuth / m It can be seen that the accuracy of optical deformation monitoring using this method can reach a sub-millimeter level within 10 m, and the relative error can be controlled within 6%.In order to better verify the applicability and effectiveness of the proposed method, the detection distance was increased, and the results obtained by repeating the above experiments are shown in Table 6: It can be seen that the accuracy of optical deformation monitoring using this method can reach the millimeter level at a detection distance of tens of meters, and the relative error can be controlled within 7%.
The experiments show that the method proposed in this paper realizes the real-time pixel conversion function at any position in the range of several meters to dozens of meters.The method does not need to solve the parameters again after solving the model parameters once and can be reused in subsequent experiments.It solves the problem of the pixel deformation quantity is difficult to convert into deformation quantity under a basic length unit in optical deformation monitoring, and it greatly improves the applicability and practicability of the whole deformation monitoring system.

Image Fusion
In order to verify the effectiveness of the proposed algorithm in Section 3.4, we carried out multi-objective image fusion verification experiments outdoors.Firstly, some corner reflector targets are randomly arranged in the experimental site, as shown in Figure 10.The GB-SAR image and an optical image of these corner reflectors are depicted in Figure 11, where C n represents the corresponding position of the nth corner reflector in the optical and radar image.Based on the algorithm in Section 3.4, the images are merged, as shown in Figure 12.The enhanced fusion result is evident, displaying improved alignment between the targets in both the GB-SAR and optical images.No significant dislocation is observed, resulting in a more natural and clear visualization effect.The fusion of GB-SAR and optical images allows for the identification of potential hazards through radar, followed by the selection of any corresponding areas in the optical image for optical deformation monitoring.Subsequently, the fusion data from GB-SAR radial deformation monitoring and optical image two-dimensional deformation monitoring are utilized to calculate the comprehensive three-dimensional deformation data of the target.

Obtaining Three-Dimensional Deformation Quantity
In order to the feasibility of the method proposed in this paper, we carried out data acquisition experiments in actual outdoor scenes to obtain the three-dimensional deformation monitoring amount of the target; the experimental scene is shown in Figure 13.By moving the target along the X, Y, and Z axes by 2.0 mm, 2.0 mm, and 2.5 mm, respectively, the actual three-dimensional deformation of the target is 3.77 mm.First, the GB-SAR was used to monitor the radial deformation of the target, and the radial deformation, dl, of the target is 2.13 mm, as shown in Figure 14.
After that, the optical camera is used to obtain the optical two-dimensional deformation monitoring of the target, and the two-dimensional deformation, dx = 3 pixel and dz = 4 pixel, of the target in pixel units is obtained, as shown in Figure 15.
Using the ranging function of GB-SAR to obtain the detection range of the target imaging results in R = 13.73m, as shown in Figure 16.
Based on the pixel conversion algorithm in Section 3.3, the diagonal pixels of the calibration board at this detection distance are P = 56.23 pixels.Therefore, the pixel conversion ratio is r = 2.01.Since the actual image of optical monitoring is interpolated three times, the pixel conversion ratio in the image is r = 2.01 / 3 = 0.67, and then the optical deformation monitoring result-dx = 2.01 mm, dz = 2.68 mm-under the basic length unit is obtained.By substituting the optical deformation monitoring data and GB-SAR monitoring data into the three-dimensional solution model shown in Figure 5, the final comprehensive three-dimensional deformation monitoring result is 3.91 mm, and the relative error is 3.71%.
In order to better show the effect of the experiment, the target was moved the same distance along the X-axis, Y-axis, and Z-axis by 3 mm, 4 mm, and 5 mm, respectively.The comparative experimental results are illustrated in Figure 17.The term "GB-SAR Monitoring" denotes the outcomes obtained through the exclusive use of GB-SAR for deformation monitoring, and "3D Monitoring" signifies the results obtained by applying the method proposed in this paper to deformation monitoring.When observing the results, it becomes evident that when there is a significant angle between the actual deformation direction of the target and the radial direction of the radar, the radial deformation monitored by the radar is notably underestimated compared to the real deformation.However, through the supplementation of optical image information, the three-dimensional monitoring of deformation aligns more closely with the actual deformation.This approach proves effective in preventing missed alarms resulting from substantial angles between the deformation direction and the radial direction of the radar.

2.13
In order to better verify the applicability and effectiveness of the proposed method, the detection distance of the target was changed.The radar imaging is shown in Figure 18, and the detection distance is R = 50.64m.The mean relative error (MRE) and root mean square error (RMSE) of the experimental results of the proposed method under different detection distances were calculated, as shown in Table 7. Through the conducted experiments, it can be proved that the three-dimensional deformation monitoring method proposed in this paper can effectively carry out comprehensive three-dimensional deformation monitoring for targets within a monitoring range of tens of meters.Among them, sub-millimeter-level monitoring accuracy can be achieved at a detection distance of 10 meters, and millimeter-level monitoring accuracy can be achieved at a detection distance of tens of meters.The average relative error of monitoring can be controlled within 10%, and the root mean square error can be controlled within 0.7 mm, indicating that the monitoring effect of this method is good and the comprehensive three-dimensional deformation of the monitoring target can be effectively applied.

Comparison of Conventional Methods and the Proposed Method
Predominantly, the research in three-dimensional deformation monitoring is centered on spaceborne and airborne SAR, with limited methodologies tailored for ground targets.Current research into the three-dimensional deformation monitoring of ground targets predominantly focuses on two main approaches.The primary research directions include multi-station three-dimensional deformation monitoring and the use of prior terrain information to predict deformations in various directions.
In the context of multi-station three-dimensional deformation monitoring, the predominant method entails employing multiple GB-SARs for multi-view interferometry measurements.This technique typically necessitates a minimum of two GB-SARs, and for enhanced accuracy, three GB-SARs are often concurrently deployed.However, the current accuracy level is confined to the centimeter range, accompanied by substantial costs and deployment challenges, significantly impeding the widespread adoption of GB-SARs.
In contrast, the three-dimensional deformation method proposed in this article requires only single-station monitoring.This significantly diminishes the cost and deployment complexities associated with three-dimensional deformation monitoring when compared to the aforementioned method.
Additionally, in contrast to predicting deformation information using prior terrain information, which relies on prior knowledge of the terrain and time series analysis combined with geographic reference models, the proposed method in this article operates without the need for geographic reference models or similar data.This increased flexibility enhances its adaptability across diverse scenarios, effectively expanding the application range of GB-SAR.

Conflict between Accuracy and Monitoring Distance
The accuracy of the three-dimensional deformation monitoring method proposed in this article decreases with the increase in monitoring distance, as demonstrated by a series of experiments in Section 4.3.The main reason for this problem is limited by the current constraints of optical imaging equipment.In theory, the current method can achieve millimeter-level accuracy in three-dimensional deformation monitoring at tens of meters and can control the relative error within 10%, as illustrated in Figures 17 and 19 and Table 7.However, as the detection distance increases to the 100-meter range, the accuracy of the three-dimensional deformation monitoring method will decrease to the centimeter level.With the rapid development of optical devices, long-distance three-dimensional deformation monitoring will gradually be carried out.At that time, the atmospheric phase correction work mentioned in Formula (2) will also come into the picture, further enhancing the applicability of the method in this article.

Conclusions
This paper introduces a new method to obtain the three-dimensional deformation of a target based on ground-based synthetic aperture radar deformation monitoring and optical deformation monitoring.This method only needs single-station monitoring, which avoids the high cost and difficult deployment of dual-station monitoring.It can obtain a comprehensive three-dimensional deformation of the ground-based target.Compared with the traditional use of GB-SAR alone for deformation monitoring, it can effectively reduce the monitoring error and avoid potential safety hazards such as missed alarms.Through outdoor physical acquisition experiments, the proposed method demonstrates sub-millimeter monitoring accuracy within a 10-meter range and millimeter-level accuracy at a range of tens of meters.The average relative error can be controlled within 10%, and the root mean square error can be controlled within 0.7 mm.It can be effectively applied to the three-dimensional deformation monitoring of foundation targets field and provides a new idea for natural disaster warnings and building deformation monitoring.

Figure 2 .
Figure 2. The process of three-dimensional deformation monitoring.

Figure 3 .
Figure 3.The basic process of pixel conversion based on fitting regression.

By using the nonlinear model in 3 . 3 .
2 to fit the parameters a = 756.3,b = −0.5089,c = 109.3,d = −0.04932, the nonlinear regression model in the fitting range is obtained:

Figure 10 .
Figure 10.Optical image of some corner reflectors.

Figure 11 .
Figure 11.GB-SAR image and an optical image of the corner reflectors.

Figure 14 .
Figure 14.GB-SAR acquires the radial deformation of the target.

Figure 17 .
Figure 17.Comparison of monitoring effects at R = 13.73 m.

Figure 18 .
Figure 18.GB-SAR image of the target.The outcomes of the repeated experiments showcasing three-dimensional deformation monitoring are illustrated in Figure 19.

Figure 19 .
Figure 19.Comparison of monitoring effects at R = 50.64m.

Table 1 .
Index parameters of the system.

Table 3 .
Comparison of fitting prediction effect.

Table 4 .
Comparison of fitting prediction effect.

Table 6 .
Optical deformation monitoring results with a detection distance of 55.62 m.

Table 7 .
The monitoring effect under different detection distances.