A surface three-dimensional deformation inversion method, device, equipment and medium

Abstract

The application discloses a kind of ground three-dimensional deformation inversion method, device, equipment and medium, involve ground monitoring technical field, only using the monitoring data of surface multi-track InSAR, without high-density monitoring point is arranged on ground to obtain external data, and by obtaining the terrain slope angle of surface from terrain data, by the north-south deformation of slope angle constraint surface, to construct the relationship equation of east-west deformation and north-south deformation in surface three-dimensional deformation, simultaneously, the stress model of surface is coupled into surface three-dimensional deformation inversion model to carry out surface slope stability constraint, the terrain slope constraint of surface can make up the deficiency of InSAR north-south monitoring capability, and the influence of single-point noise is reduced by excess observation in surface slope stability constraint, so as to need external measurement data, so that the result precision of three-dimensional deformation inversion is improved.

Description

Technical Field

[0001] The present invention relates to the field of ground monitoring technology, and in particular to a method, device, equipment and medium for inverting three-dimensional deformation of the ground surface. Background Art

[0002] Interferometric Synthetic Aperture Radar (InSAR) technology has become an important tool for geological disaster early warning, urban subsidence monitoring, and crustal movement research due to its large-scale, high-precision, and all-weather surface deformation monitoring capabilities. Traditional InSAR technology obtains one-dimensional surface deformation information through phase changes in the radar's line-of-sight (LOS) direction. However, its main drawback is that it cannot directly separate the three-dimensional deformation components (eastward, northward, and vertical).

[0003] In order to solve and break through the limitations of the above problems, multi-orbit InSAR joint solution methods have gradually become a research hotspot. The main multi-orbit InSAR joint solution methods are: one is to use multi-source InSAR data fusion, mainly using differential interferometry technology to obtain LOS deformation results, and using multi-aperture InSAR (MAI) or offset tracking technology to obtain line-of-sight deformation results. By fusing ascending and descending InSAR data and different perspectives, the inversion of three-dimensional deformation fields can be achieved; another is to supplement InSAR monitoring through the assumption of prior information, such as surface flow models; and another is to refer to external data, such as GPS, leveling, etc.; however, InSAR satellites usually fly in polar orbits, and the LOS results are less sensitive to north-south deformations. The azimuth results are limited by terrain factors and the influence of InSAR coherence, resulting in poor measurement accuracy, thereby reducing the accuracy of multi-orbit InSAR three-dimensional deformation inversion results.

[0004] Therefore, in order to compensate for the problem of poor measurement accuracy of the surface north-south orientation results due to the limitations of terrain factors and the influence of InSAR coherence, the measurement accuracy is currently improved by utilizing external data and setting prior information. However, when using external data, a high density of monitoring points needs to be deployed on the ground. At the same time, the setting of prior information is hypothetical and unclear, resulting in low accuracy of the three-dimensional deformation inversion results. Summary of the Invention

[0005] The embodiments of the present invention provide a method, device, equipment and medium for inverting three-dimensional surface deformation, which can solve the problem in the existing technology that, at this stage, measurement accuracy is improved by utilizing external data and setting prior information, but when using external data, a high density of monitoring points needs to be deployed on the ground. At the same time, the setting of prior information is hypothetical and unclear, resulting in low accuracy of the three-dimensional deformation inversion results.

[0006] The embodiment of the present application provides a surface three-dimensional deformation inversion method, comprising the following steps:

[0007] Obtaining multi-track InSAR deformation data of the surface;

[0008] According to the digital elevation information of the surface in the multi-track InSAR deformation data of the surface, obtaining the terrain slope angle of the surface;

[0009] According to the LOS deformation and azimuth deformation in the multi-track InSAR deformation data, obtaining the vertical deformation and east-west deformation of the surface in the three-dimensional direction; using the terrain slope angle of the surface, constraining the deformation field of the LOS deformation and azimuth deformation in the north-south component in the multi-track InSAR deformation data, and obtaining the north-south deformation of the surface in the three-dimensional direction after the constraint according to the constrained LOS deformation and azimuth deformation;

[0010] According to the east-west deformation and north-south deformation of the surface in the three-dimensional direction, constructing the east-west deformation and north-south deformation relationship model of the surface three-dimensional deformation; establishing the crustal stress model of the surface according to the vertical deformation of the surface in the three-dimensional direction, and coupling the east-west deformation and north-south deformation relationship model and the crustal stress model into a slope crustal stress inversion model;

[0011] Using the slope crustal stress inversion model to perform surface deformation inversion, and obtaining the three-dimensional deformation of the surface.

[0012] Preferably, the relationship among the terrain slope angle of the surface, the LOS deformation and the azimuth deformation and the surface three-dimensional deformation is expressed as:

[0013] ;

[0014] ;

[0015] ;

[0016] Wherein: represents the LOS deformation of each satellite; represents the LOS deformation of each satellite; represents the vertical deformation; represents the east-west deformation; represents the north-south deformation; represents the satellite incidence angle; represents the satellite heading angle; represents the terrain slope of the surface.

[0017] Preferably, the slope crustal stress inversion model is expressed as:

[0018] ;

[0019] wherein: represents the observation value column vector; N represents the number of adjacent points in the ground strain model; represents the ascending track line of view InSAR monitoring result of the Nth adjacent point; represents the descending track line of view InSAR monitoring result of the Nth adjacent point; represents the residual error; represents the to-be-estimated parameter; respectively represent the vertical, east-west and north-south deformation monitoring results in the ground surface deformation monitoring; respectively represent the strain parameters and rotation parameters in the ground stress model;

[0020] observation matrix is the SAR imaging geometry matrix of the multi-directional deformation relationship model and the matrix of the ground stress model is the product of the two; , and respectively represent the east-west, north-south and vertical distances between the to-be-solved point and the adjacent points.

[0021] Preferably, the SAR imaging geometry matrix of the multi-directional deformation relationship model and the matrix of the ground stress model are respectively represented as:

[0022] ;

[0023] ;

[0024] ;

[0025] wherein: , and respectively represent the vertical, east-west and north-south distances between the to-be-solved point and the adjacent points.

[0026] Preferably, in the ground surface deformation inversion using the slope ground stress inversion model, the error in the deformation inversion is corrected by using the joint phase standard deviation and the reference point distance error, comprising:

[0027] total variance in the ground surface deformation inversion comprising the InSAR observation phase standard deviation and the distance error between the reference points , and the total variance is represented as:

[0028] ;

[0029] ;

[0030] The distance error between the reference points is the inverse of the distance between the two points. The standard deviation of the InSAR observation phase is obtained by opening a regular window with the estimated observation point as the center and counting the residuals of the InSAR observation values ​​within the regular window, which is expressed as:

[0031] ;

[0032] in: Indicates the number of pixels in the window; Represents the observation value of the corresponding observation point in the window; Represents the average value of deformation observations of all observation points in the window;

[0033] Get the weight of the total variance, expressed as:

[0034] ;

[0035] The inverse of the total variance is used as the weight, and the observation error is reduced by controlling the weights of different observations.

[0036] Preferably, obtaining the three-dimensional deformation of the ground surface includes:

[0037] According to the slope stress inversion model, weight matrix and the observation values ​​of each surface observation point, the least squares adjustment algorithm is used to solve the optimal estimate of the three-dimensional deformation rate point by point, which is expressed as:

[0038] ;

[0039] in: represents the three-dimensional deformation rate vector; B The observation matrix representing the three-dimensional surface deformation inversion model; B T represents the transpose of the observation matrix; ω represents the weight matrix; L represents the observation vector.

[0040] An embodiment of the present invention further provides a device for inverting three-dimensional surface deformation, comprising:

[0041] Data module, used to obtain multi-orbit InSAR deformation data of the surface;

[0042] The deformation monitoring module is used to obtain the surface terrain slope angle based on the surface digital elevation information in the multi-track InSAR deformation data of the surface;

[0043] Based on the LOS and azimuth deformations in the multi-orbit InSAR deformation data, the vertical deformation and east-west deformation of the surface in the three-dimensional direction are obtained. The terrain slope angle of the surface is used to constrain the deformation field of the LOS and azimuth deformations in the multi-orbit InSAR deformation data on the north-south component, and based on the constrained LOS and azimuth deformations, the constrained north-south deformation of the surface in the three-dimensional direction is obtained.

[0044] The inversion module is used to construct a relationship model between the east-west and north-south deformation of the three-dimensional surface deformation based on the east-west deformation of the surface and the constrained north-south deformation. The inversion module is used to establish a ground stress model based on the vertical deformation of the surface in the three-dimensional direction, and couple the east-west and north-south deformation relationship models with the ground stress model to form a slope ground stress inversion model.

[0045] The surface deformation inversion model is used to obtain the three-dimensional deformation of the surface.

[0046] An embodiment of the present invention further provides an electronic device, including a memory and a processor;

[0047] The memory is used to store computer programs;

[0048] The processor is used to implement the steps of the above-mentioned method for inverting three-dimensional surface deformation when executing the computer program stored in the memory.

[0049] An embodiment of the present invention further provides a computer-readable storage medium for storing a computer program, which, when executed by a processor, implements the steps of the above-mentioned method for inverting three-dimensional surface deformation.

[0050] The embodiments of the present invention provide a method, apparatus, device, and medium for inverting three-dimensional surface deformation. Compared with the prior art, the methods and devices have the following advantages:

[0051] The present invention only uses multi-orbit InSAR monitoring data on the surface, eliminating the need for high-density monitoring points on the ground to obtain external data. Instead, the surface terrain aspect angle is obtained from the terrain data, and the north-south deformation of the surface is constrained by the aspect angle to construct a relationship model between the east-west deformation and the north-south deformation in the three-dimensional surface deformation, thereby constraining the north-south deformation. At the same time, the surface ground stress model is coupled into the slope aspect ground stress inversion model to constrain the surface slope stability. The surface terrain aspect constraint can compensate for the shortcomings of InSAR's north-south monitoring capabilities, and the surface slope stability constraint reduces the influence of single-point noise through redundant observations, eliminating the need for external measurement data and improving the accuracy of the three-dimensional deformation inversion results. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] Figure 1A schematic diagram of the main process of a method for inverting three-dimensional surface deformation provided by an embodiment of the present invention;

[0053] Figure 2 A schematic diagram of surface elevation simulated by a method for inverting three-dimensional surface deformation provided by an embodiment of the present invention;

[0054] Figure 3 Schematic diagram of the original simulated surface deformation of a method for inverting three-dimensional surface deformation provided by an embodiment of the present invention; (a) is the simulated vertical deformation; (b) is the simulated east-west deformation; and (c) is the simulated north-south deformation.

[0055] Figure 4 Schematic diagram of the results of simulation of noise-added LOS and azimuth deformation data for ascending and descending orbits, according to a method for inverting three-dimensional surface deformation provided by an embodiment of the present invention; (a) shows the LOS deformation for ascending orbits; (b) shows the LOS deformation for descending orbits; (c) shows the line-of-sight deformation for ascending orbits; and (d) shows the line-of-sight deformation for descending orbits.

[0056] Figure 5 A schematic diagram comparing the three-dimensional surface deformation inversion method provided by an embodiment of the present invention using the original least squares method, the ground strain model method and the method of the present invention to invert the three-dimensional surface deformation; (a) is the vertical deformation solved by the least squares method; (b) is the east-west deformation solved by the least squares method; (c) is the north-south deformation solved by the least squares method; (d) is the vertical deformation solved by the ground strain model method; (e) is the east-west deformation solved by the ground strain model method; (f) is the north-south deformation solved by the ground strain model method; (g) is the vertical deformation solved by the method of the present invention; (h) is the east-west deformation solved by the method of the present invention; and (i) is the north-south deformation solved by the method of the present invention. DETAILED DESCRIPTION

[0057] To make the above-mentioned objects, features, and advantages of the present invention more readily apparent, specific embodiments of the present invention are described in detail below with reference to the accompanying drawings. The following description sets forth numerous specific details to facilitate a full understanding of the present invention. However, the present invention can be implemented in many other ways than those described herein, and those skilled in the art may make similar modifications without departing from the scope of the present invention. Therefore, the present invention is not limited to the specific embodiments disclosed below.

[0058] See also Figure 1 The embodiment of the present invention provides a method for inverting three-dimensional surface deformation, comprising the following steps:

[0059] Step S1: Acquisition of multi-track InSAR data.

[0060] Step S2: Obtain InSAR monitoring results through different means such as D-InSAR, POT and MAI.

[0061] The multi-orbit InSAR data is extracted and the SLC file required for the data processing process is generated; based on the spatial geometric relationship between the DEM's geographic projection coordinate system and the SAR's range-Doppler coordinate system, a conversion lookup table between the DEM coordinate system and the SAR coordinate system is calculated and generated, and the DEM is converted to the main image coordinate system using a polynomial fitting algorithm to assist in image registration between SAR images; a certain time and space interval is used as a threshold to generate a time and space baseline distribution map, and images with a time interval less than the threshold are subjected to pairwise differential interferometry processing in chronological order; the baseline combination should ensure the continuity between images and there should also be a certain amount of redundant observations to weaken the influence of noise, and finally the multi-orbit InSAR deformation monitoring results are obtained through filtering, unwrapping, and atmospheric delay correction.

[0062] Step S3: Select adjacent points around the reference point to be determined.

[0063] The purpose of this step is to select the deformation of the adjacent points of the estimated point, and use the deformation of the adjacent points to carry out the estimation of the three-dimensional deformation of the surface in the subsequent steps. The accuracy of the three-dimensional deformation estimation of the surface of the target point is improved through multi-dimensional observation. First, the search window is determined to be 5*5 pixels. The search is performed according to the point selection criteria. If no suitable point is found, the search window length is increased and the search is continued. The selection criteria of the adjacent points include:

[0064] Coherence is an important indicator for measuring InSAR monitoring accuracy. In areas with low coherence, InSAR monitoring accuracy is low. Therefore, the coherence of adjacent points cannot be too low. For the field of geological disasters, the coherence of adjacent points should theoretically be greater than 0.2. InSAR overlapping shadow areas are caused by the influence of the SAR satellite's flight direction and incident angle, resulting in distortion of SAR images. The selection of adjacent points should also avoid these areas.

[0065] Step S4: Establish the observation equation of the slope aspect-ground stress model.

[0066] First, the terrain aspect angle is calculated based on the digital elevation model (DEM) and a regional aspect distribution map is constructed.

[0067] The following equation can be established based on the relationship between InSAR-LOS deformation, InSAR azimuth deformation, terrain slope and surface three-dimensional deformation:

[0068] .

[0069] .

[0070] .

[0071] in: Indicates the LOS deformation of each satellite; Indicates the deformation of each satellite's line of sight; represents vertical deformation; Indicates east-west deformation; Indicates north-south deformation; represents the satellite incidence angle; Indicates the satellite heading angle; Indicates the terrain slope.

[0072] Assume that there exists a The relevant points are used as the relevant points of the research model to establish the slope constraint and ground stress model. The observation equation of multi-track InSAR observations and three-dimensional deformation is expressed as:

[0073] .

[0074] in: represents the column vector of observation values; N represents the number of adjacent points in the ground strain model; It represents the InSAR monitoring result of the Nth adjacent point in the ascending orbit line of sight; Indicates the InSAR monitoring result of the descending azimuth of the Nth adjacent point; represents the residual; represents the parameter to be estimated; Respectively represent the vertical, east-west, and north-south deformation monitoring results in surface deformation monitoring; Respectively represent the strain parameters and rotation parameters in the ground stress model; the observation matrix is the SAR imaging geometric matrix of the multi-directional deformation relationship model and the matrix of the ground stress model The product of 、 and Respectively represent the east-west distance, north-south distance, and vertical distance between the point to be found and the related adjacent points.

[0075] .

[0076] .

[0077]

[0078] Step S5: Determine the weights jointly using the joint phase standard deviation and the reference point distance.

[0079] The total variance is composed of the InSAR observation error and the distance error between the reference points (independent of each other) ; first, the normalization coefficient is calculated, which is used to make the average variance contribution of the two equal through a scale factor, expressed as:

[0080] .

[0081] The total variance is obtained by combining the two variance terms, expressed as:

[0082] .

[0083] Where the reference point spacing error is the inverse of the distance between the two points, and the standard deviation of the InSAR observation phase is estimated by opening a regular window centered on the estimated observation point, and using the following formula, expressed as:

[0084] .

[0085] Where: The number of pixels in the window is represented by N ; The observation of the corresponding observation point in the window is represented by Y i ; The average deformation observation of all observation points in the window is represented by X i.

[0086] The total variance is the sum of the two terms, and the weight is the inverse of it, and the weight of the equation is obtained by normalization, expressed as:

[0087] .

[0088] Step S6: point-by-point solution of the three-dimensional deformation of the ground surface.

[0089] According to the observation matrix, weight matrix and observation value in step S4 and step S5, the least square difference is used to obtain the optimal estimation value of the three-dimensional ground surface deformation rate, and the equation is solved, expressed as:

[0090] .

[0091] For the multi-track InSAR deformation monitoring result in step S2, repeat steps S3-S6 for different image elements until the three-dimensional deformation of all InSAR pixels in the study area is solved.

[0092] In the simulation experiment, a certain area is simulated to simulate the terrain, as shown in Figure 2 The three-dimensional deformation field in the east-west direction, north-south direction and vertical direction in the area, where the east-west and north-south deformation is about-20cm~20cm, and the vertical deformation is about 25cm, as shown in Figure 3 .

[0093] Assuming that the azimuth angle of the ascending satellite is 349.55° and the mean angle of incidence is 37.45°; the azimuth angle of the descending satellite is 190.34° and the mean angle of incidence is 43.98°, the LOS and azimuth projection deformations of the ascending and descending orbits are calculated through the geometric projection relationship between the three-dimensional surface deformation and the imaging between satellites; Gaussian noise is added to the obtained LOS and azimuth projection deformations of the ascending and descending orbits, and its standard deviation is 5mm (ascending orbit LOS direction), 100mm (ascending orbit azimuth direction), 7mm (descending orbit LOS direction), and 140mm (descending orbit azimuth direction), respectively; the original data used for the simulation experiment can be obtained, as shown in Figure 2. Figure 4 As shown, they are respectively the ascending orbit LOS direction, the ascending orbit azimuth direction, the descending orbit LOS direction and the descending orbit azimuth direction.

[0094] At the same time, the three-dimensional surface deformation field was solved using the traditional least squares method, the method based only on the ground strain model, and the method of the present invention. The root mean square error (RMS) between the three-dimensional deformation field and the original simulation was calculated as the accuracy standard, as shown in Table 1. Figure 5 The three-dimensional surface deformation inverted by the three methods is also shown and compared.

[0095] Table 1 Root mean square residuals of three-dimensional surface deformation inversion using different methods

[0096]

[0097] The present invention does not require external data and only uses InSAR monitoring results to estimate the three-dimensional surface deformation by fusing slope constraints with the ground stress model. At the same time, in order to improve the accuracy of the results, the study adopts a combined phase and distance method to determine the weights and improve the accuracy of three-dimensional deformation monitoring.

[0098] The above-described embodiments merely illustrate several implementations of the present invention, and while their descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the patent. It should be noted that a person skilled in the art would be able to make numerous variations and improvements without departing from the spirit of the present invention, all of which fall within the scope of protection of the present invention. Therefore, the scope of protection of the patent for this invention shall be determined by the appended claims.

Claims

1. A method for inverting three-dimensional surface deformation, characterized in that: The following steps are involved: Acquire multi-orbit InSAR deformation data of the earth's surface; Obtain the surface terrain slope angle based on the surface digital elevation information in the multi-orbit InSAR deformation data; Based on the LOS and azimuth deformations in the multi-orbit InSAR deformation data, the vertical deformation and east-west deformation of the surface in the three-dimensional direction are obtained. The terrain slope angle of the surface is used to constrain the deformation field of the LOS and azimuth deformations in the multi-orbit InSAR deformation data on the north-south component, and based on the constrained LOS and azimuth deformations, the constrained north-south deformation of the surface in the three-dimensional direction is obtained. Based on the east-west deformation of the surface in the three-dimensional direction and the constrained north-south deformation, a relationship model of the east-west deformation and the north-south deformation of the surface in the three-dimensional direction is constructed; Based on the vertical deformation of the ground surface in the three-dimensional direction, a ground stress model of the ground surface is established, and the east-west deformation relationship model and the north-south deformation relationship model are coupled with the ground stress model to form a slope ground stress inversion model; The surface deformation inversion model is used to obtain the three-dimensional deformation of the surface.

2. A method for inverting three-dimensional surface deformation according to claim 1, characterized in that: The relationship between the terrain aspect angle, LOS deformation and azimuth deformation of the surface and the three-dimensional deformation of the surface is expressed as: d2-d3tan(π-α aspect )=0; Where: D LOS Indicates the LOS deformation of each satellite; D AZI represents the deformation of each satellite's line of sight; d1 represents the vertical deformation; d2 represents the east-west deformation; d3 represents the north-south deformation; θ inc represents the satellite incident angle; α heading represents the satellite heading angle; α aspect Indicates the terrain slope of the surface.

3. The method for inverting three-dimensional surface deformation according to claim 2, characterized in that: The slope stress inversion model is expressed as: L = BX + V; in: represents the column vector of observation values; N represents the number of adjacent points in the ground strain model; It represents the InSAR monitoring result of the Nth adjacent point in the ascending orbit line of sight; represents the InSAR monitoring result of the Nth adjacent point in the descending orbit azimuth direction; V represents the residual; Represents the parameters to be estimated; d1, d2, and d3 represent the vertical, east-west, and north-south deformation monitoring results in surface deformation monitoring, respectively; ε and ω represent the strain parameter and rotation parameter in the ground stress model, respectively; The observation matrix B is the SAR imaging geometric matrix B of the multi-directional deformation relationship model gm and the matrix B of the ground stress model sm The product of .

4. The method for inverting three-dimensional surface deformation according to claim 3, characterized in that: The SAR imaging geometric matrix B of the multi-directional deformation relationship model gm and the matrix B of the ground stress model sm Respectively expressed as: B=B gm ·B sm ; Among them: Δx1, Δx2 and Δx3 represent the vertical distance, east-west distance and north-south distance between the point to be determined and the relevant adjacent points respectively.

5. The method for inverting three-dimensional surface deformation according to claim 4, characterized in that: In the surface deformation inversion using the slope stress inversion model, the error in the deformation inversion is corrected using the joint phase standard deviation and the reference point spacing error, including: Total variance σ in surface deformation inversion 2 total Including InSAR observation phase standard deviation σ 2 phase The distance error σ between the reference point 2 distance , the total variance is expressed as: s 2 total =k 2 ·s 2 phase +s 2 distance ; The distance error between the reference points is the inverse of the distance between the two points. The standard deviation of the InSAR observation phase is obtained by opening a regular window with the estimated observation point as the center and counting the residuals of the InSAR observation values ​​within the regular window, which is expressed as: Where: num represents the number of pixels in the window; d(i,j) represents the observation value of the corresponding observation point in the window; Represents the average value of deformation observations of all observation points in the window; Get the weight of the total variance, expressed as: The inverse of the total variance is used as the weight, and the observation error is reduced by controlling the weights of different observations.

6. The method for inverting three-dimensional surface deformation according to claim 5, characterized in that: The obtaining of three-dimensional deformation of the earth's surface includes: According to the slope stress inversion model, weight matrix and the observation values ​​of each surface observation point, the least squares adjustment algorithm is used to solve the optimal estimate of the three-dimensional deformation rate point by point, which is expressed as: in: represents the three-dimensional deformation rate vector; B represents the observation matrix of the three-dimensional surface deformation inversion model; B T represents the transpose of the observation matrix; ω represents the weight matrix; L represents the observation value vector.

7. A device for inverting three-dimensional surface deformation, characterized in that: include: Data module, used to obtain multi-orbit InSAR deformation data of the surface; The deformation monitoring module is used to obtain the surface terrain slope angle based on the surface digital elevation information in the multi-track InSAR deformation data of the surface; Based on the LOS and azimuth deformations in the multi-orbit InSAR deformation data, the vertical deformation and east-west deformation of the surface in the three-dimensional direction are obtained. The terrain slope angle of the surface is used to constrain the deformation field of the LOS and azimuth deformations in the multi-orbit InSAR deformation data on the north-south component, and based on the constrained LOS and azimuth deformations, the constrained north-south deformation of the surface in the three-dimensional direction is obtained. The inversion module is used to construct a relationship model between the east-west and north-south deformation of the three-dimensional surface deformation based on the east-west deformation of the surface in the three-dimensional direction and the constrained north-south deformation; Based on the vertical deformation of the ground surface in the three-dimensional direction, a ground stress model of the ground surface is established, and the east-west deformation relationship model and the north-south deformation relationship model are coupled with the ground stress model to form a slope ground stress inversion model; The surface deformation inversion model is used to obtain the three-dimensional deformation of the surface.

8. An electronic device, characterized in that: include: memory and processor; The memory is used to store computer programs; The processor is configured to implement the steps of a method for inverting three-dimensional surface deformation according to any one of claims 1 to 6 when executing the computer program stored in the memory.

9. A computer-readable storage medium, characterized in that Used to store a computer program, which, when executed by a processor, implements the steps of a method for inverting three-dimensional surface deformation as claimed in any one of claims 1 to 6.

Images