High-precision 3D deformation extraction method using multi-track InSAR

Abstract

This invention discloses a high-precision three-dimensional deformation extraction method using multi-track InSAR, comprising: acquiring multi-track LOS-direction deformation; establishing a set of equations relating LOS-direction deformation to three-dimensional deformation, and establishing a set of equations relating LOS-direction deformation to two-dimensional deformation; obtaining three-dimensional deformation, as well as vertical and east-west deformation; establishing three-dimensional and two-dimensional deformation error propagation models; establishing an additional systematic error propagation model; establishing an objective function and adding regularization constraints to the objective function; obtaining an objective function matrix based on the objective function with constraints and the established error propagation model; obtaining the LOS-direction deformation error that satisfies the minimum value of the objective function based on the objective function matrix; substituting the LOS-direction deformation error into the three-dimensional and two-dimensional deformation error propagation models to obtain the three-dimensional deformation error; and obtaining the three-dimensional deformation result based on the three-dimensional deformation error. This method achieves the goal of improving the accuracy of three-dimensional deformation, especially significantly improving the accuracy of north-south deformation.

Description

Technical Field

[0001] This invention relates to the field of microwave remote sensing measurement, and more specifically to a method for extracting high-precision three-dimensional deformation using multi-track InSAR. Background Technology

[0002] Interferometric Synthetic Aperture Radar (InSAR) technology has been an effective method for monitoring geological hazards such as landslides, earthquakes, volcanic eruptions, and mining subsidence for the past two decades. However, InSAR technology only provides a one-dimensional deformation along the radar line of sight (LOS). In practical applications, in most cases, the LOS deformation is not suitable for assessing the extent of a hazard.

[0003] In recent years, with the increase in the number of SAR satellites, it has become possible to acquire multiple sets of LOS deformation data with different headings and incident angles (ascending and descending orbit data, or left-looking and right-looking data). Using three sets of LOS deformation data to establish a set of observation equations, vertical, east-west, and north-south deformation data can be obtained. However, due to the limitation of the near-north-south flight orbit of SAR satellites, the accuracy of the north-south solution is much lower than that of the vertical and east-west solutions. Summary of the Invention

[0004] To address the technical problems in existing technologies, this invention provides a multi-track InSAR method for extracting high-precision three-dimensional deformations. It includes the following:

[0005] A first aspect of the present invention provides a method for extracting high-precision three-dimensional deformation using multi-track InSAR, comprising:

[0006] The LOS deformation of each orbit is obtained from multi-orbit SAR data, including processing the SAR data of the three orbits separately using InSAR technology to obtain the LOS deformation displacement of each orbit.

[0007] Based on the side-view imaging geometry of the SAR sensor, a set of equations relating LOS deformation to three-dimensional deformation is established, and a set of equations relating LOS deformation to two-dimensional deformation ignoring north-south deformation is also established.

[0008] Based on the equations relating LOS deformation to 3D deformation and LOS deformation to 2D deformation, the least squares method is used to obtain the 3D deformation results as well as the vertical and east-west deformations.

[0009] Based on the relationship between LOS deformation error and deformation result error in three dimensions, a three-dimensional deformation error propagation model is established. Based on the relationship between LOS deformation error and deformation errors in the vertical and east-west directions, a LOS deformation error and two-dimensional deformation error propagation model is established.

[0010] Based on the relationship between the vertical and east-west systematic errors and the north-south deformation components in the two-dimensional deformation error propagation model, an additional systematic error propagation model is established.

[0011] An objective function is established based on the additional system error propagation model, the relationship between the three-dimensional solution of north-south deformation and the vertical and east-west deformations ignoring the north-south solution and the LOS deformation, and a regularization constraint term is added to the objective function.

[0012] Based on the objective function with constraints and the established error propagation model, the objective function matrix is ​​obtained;

[0013] The LOS-directed deformation error that satisfies the minimum value of the objective function matrix is ​​obtained. The LOS-directed deformation error is then substituted into the three-dimensional deformation error propagation model to obtain the three-dimensional deformation error.

[0014] The three-dimensional deformation result is obtained based on the three-dimensional deformation error.

[0015] Optionally, in establishing the equation set relating LOS deformation to three-dimensional deformation, the equation set is as follows:

[0016] ,

[0017] In the formula, This indicates the LOS-direction deformation displacement of orbit i. ;

[0018] , This indicates the radar satellite's incident angle and the clockwise angle between the northward direction and the radar satellite's flight direction. .

[0019] , , The deformations are three-dimensional, namely vertical, east-west, and north-south.

[0020] Optionally, in establishing a set of equations relating LOS-direction deformation to neglecting north-south deformation in two dimensions, the equation set is as follows:

[0021] ,

[0022] In the formula, , To ignore the vertical and east-west deformations when the north-south direction is ignored.

[0023] Optionally, based on the equations relating LOS deformation to three-dimensional deformation and the equations relating LOS deformation to two-dimensional deformation, the least squares method is used to obtain the deformation results in the three-dimensional direction, as well as the deformations in the vertical and east-west directions. The formulas are as follows:

[0024] ,

[0025] In the formula for:

[0026] ，

[0027] The deformation along the radar line of sight is as follows:

[0028] ,

[0029] W is the LOS-oriented deformation weight matrix, using an equivalent weight matrix.

[0030] The calculated three-dimensional deformation, or the two-dimensional deformation obtained by ignoring the north-south direction, is as follows:

[0031] .

[0032] Optionally, the three-dimensional deformation error propagation model is established based on the relationship between the LOS deformation error and the deformation result error in the three-dimensional direction as follows:

[0033] ,

[0034] In the formula, ,

[0035] In the formula For i-orbital LOS-direction deformation error, , , These are the deformation coefficients of orbit i in the vertical, east-west, and north-south directions, respectively. , , These represent the deformation errors in the three-dimensional direction.

[0036] Optionally, based on the relationship between LOS-direction deformation error and vertical and east-west deformation errors, a LOS-direction deformation error and two-dimensional deformation error propagation model is established as follows:

[0037] ,

[0038] In the formula, ,

[0039] In the formula, , These represent the vertical and east-west deformation errors when the north-south direction is ignored.

[0040] Optionally, based on the relationship between the vertical and east-west systematic errors and the north-south deformation components in the two-dimensional deformation error propagation model, the additional systematic error propagation model is established as follows:

[0041] ,

[0042] In the formula, This is the error in the vertical deformation system. For east-west deformation systematic error,

[0043] Ignoring the vertical systematic error caused by the north-south deformation calculation, the east-west systematic error is positively correlated with the north-south deformation component, and its proportionality constant is:

[0044] ,

[0045] ,

[0046] In the formula ,

[0047] .

[0048] Optionally, an objective function is established based on the additional systematic error propagation model, the relationship between the three-dimensional solution of north-south deformation and the vertical and east-west deformations ignoring the north-south solution and the LOS deformation, and regularization constraints are added to the objective function, including:

[0049] Add a stable generalization function to the objective function The objective function with added constraints is in the form of:

[0050]

[0051] In the formula, It is a Euclidean 2-norm. , These are the coefficients of the constraint terms.

[0052] Optionally, the step of obtaining the LOS-directed deformation error that satisfies the minimum value of the objective function based on the objective function matrix, and substituting the LOS-directed deformation error into the three-dimensional deformation error propagation model to obtain the three-dimensional deformation error includes:

[0053] Establish a system of derivative equations for the objective function:

[0054] ,

[0055] The system of derivative equations can be represented as a matrix:

[0056] ,

[0057] The solution to the equation is:

[0058] ,

[0059] The above is satisfied LOS deformation error Substituting into the error propagation model, the three-dimensional deformation error is obtained as follows:

[0060] ,

[0061] Optimized north-south deformation is obtained by utilizing the three-dimensional directional deformation error. With east-west deformation and vertical deformation The optimized three-dimensional deformation is obtained by combining these factors.

[0062] A second aspect of the present invention provides an electronic device, the electronic device comprising:

[0063] At least one processor; and,

[0064] A memory communicatively connected to the at least one processor; wherein,

[0065] The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the multi-track InSAR high-precision three-dimensional deformation extraction method described in any of the first aspects.

[0066] A third aspect of the present invention provides a computer-readable storage medium storing computer instructions for causing a computer to perform the multi-track InSAR extraction high-precision three-dimensional deformation method described in any of the first aspects.

[0067] The multi-orbit InSAR high-precision three-dimensional deformation extraction method of the present invention is based on multi-orbit LOS deformation data to extract three-dimensional deformation, reducing the strict requirements of SAR satellite parameters for extracting high-precision three-dimensional deformation results, and eliminating the need for given parameter conditions and complex signal processing methods. It obtains high-precision three-dimensional deformation results through mathematical calculation, thereby improving the accuracy of north-south deformation. Attached Figure Description

[0068] Figure 1 This is a flowchart of the multi-track InSAR high-precision three-dimensional deformation extraction method disclosed in this embodiment. Detailed Implementation

[0069] Various exemplary embodiments of the present invention will now be described in detail. This detailed description should not be considered as a limitation of the present invention, but rather as a more detailed description of certain aspects, features, and embodiments of the present invention.

[0070] It should be understood that the terminology used in this invention is merely for describing particular embodiments and is not intended to limit the invention. Furthermore, with respect to numerical ranges in this invention, it should be understood that the upper and lower limits of the range and each intermediate value between them are specifically disclosed. Any stated value or intermediate value within a stated range, as well as each smaller range between any other stated value or intermediate value within said range, are also included in this invention. The upper and lower limits of these smaller ranges may be independently included or excluded from the range.

[0071] Unless otherwise stated, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art. While only preferred methods and materials have been described herein, any methods and materials similar or equivalent to those described herein may be used in the implementation or testing of this invention. All references to this specification are incorporated by way of citation to disclose and describe methods and / or materials associated with those references. In the event of any conflict with any incorporated reference, the content of this specification shall prevail.

[0072] like Figure 1 As shown, a multi-track InSAR method for extracting high-precision three-dimensional deformation includes:

[0073] Step S1: Obtain the LOS deformation of each orbit from the SAR data of the multi-orbit system. Specifically, process the SAR data of the three orbits using InSAR technology to obtain the LOS deformation displacement of each orbit.

[0074] Specifically, geocoded LOS-oriented deformation maps can be obtained based on InSAR technology. Geocoding refers to transforming the radar image coordinate system to the Universal Transverse Mercator (UTM) coordinate system.

[0075] Step S2: Based on the side-view imaging geometry of the SAR sensor, establish a set of equations relating LOS deformation to 3D deformation, and establish a set of equations relating LOS deformation to 2D deformation (ignoring north-south deformation). Specifically, based on the side-view imaging geometry of the SAR sensor, the matrix form of the equations relating LOS deformation to 3D deformation is as follows:

[0076] ,

[0077] In the formula, This represents the deformation displacement along the LOS direction of orbit i, obtained using InSAR technology.

[0078] , This represents the radar satellite incident angle and the clockwise angle between the north direction and the radar satellite's flight direction; its values ​​are obtained from the SAR data file.

[0079] , , It exhibits three-dimensional deformation in the vertical, east-west, and north-south directions;

[0080] Since SAR satellites fly near the north-south direction, north-south deformation contributes little to the LOS deformation displacement. To reduce the impact of north-south deformation on the solution results, it is ignored, and the vertical and east-west deformations are calculated with higher accuracy. The matrix form of the equations relating LOS deformation to two-dimensional deformation is as follows:

[0081] ,

[0082] In the formula, This represents the LOS-direction deformation of the i-orbital, obtained using InSAR technology;

[0083] , This indicates the radar satellite's incident angle and the clockwise angle between the north direction and the radar satellite's flight direction;

[0084] , To ignore the vertical and east-west deformations when the north-south direction is ignored.

[0085] Step S3: Based on the equations relating LOS deformation to 3D deformation and LOS deformation to 2D deformation, the least squares method is used to obtain the 3D deformation result, as well as the vertical and east-west deformations. Specifically, based on the equations relating LOS deformation to 3D deformation and LOS deformation to 2D deformation, the least squares method is used to calculate the 3D deformation, and the more accurate vertical and east-west deformations when the north-south direction is ignored. The calculation formula is:

[0086] ,

[0087] In the formula for:

[0088] ，

[0089] In the formula, The deformation along the radar line of sight is as follows:

[0090] ,

[0091] In the formula, W is the LOS-directed deformation weight matrix. In this method, an equivalent weight matrix is ​​used, and W is:

[0092] ,

[0093] In the formula, The calculated three-dimensional deformation and the two-dimensional deformation obtained by ignoring the north-south direction are respectively:

[0094] .

[0095] Step S4: Based on the relationship between the LOS deformation error and the deformation result error in the three-dimensional direction, establish a three-dimensional deformation error propagation model. Based on the relationship between the LOS deformation error and the vertical and east-west deformation errors, establish a LOS deformation error and two-dimensional deformation error propagation model. Specifically, based on the relationship between the LOS deformation error and the three-dimensional deformation error, establish a LOS deformation error and three-dimensional deformation error propagation model. The LOS deformation error and three-dimensional deformation error propagation model is as follows:

[0096] ,

[0097] In the formula, ,

[0098] In the formula For i-orbital LOS-direction deformation error, , , Let be the deformation coefficients of orbit i in the vertical, east-west, and north-south directions. , , This represents the three-dimensional deformation error. , , ;

[0099] Based on the relationship between LOS deformation error and 2D deformation error, a propagation model for LOS deformation error and 2D deformation error is established. The propagation model for LOS deformation error and 2D deformation error is as follows:

[0100] ,

[0101] In the formula, ,

[0102] In the formula, , To ignore the vertical and east-west deformation errors when the north-south direction is ignored.

[0103] Step S5: Based on the relationship between the vertical and east-west systematic errors and the north-south deformation components in the two-dimensional deformation error propagation model, establish an additional systematic error propagation model. Specifically, based on the relationship between the vertical and east-west systematic errors and the north-south deformation components when the north-south solution is ignored, establish an additional systematic error propagation model:

[0104] ,

[0105] In the formula, This is the error in the vertical deformation system. For east-west deformation systematic error,

[0106] Ignoring the vertical systematic error caused by the north-south deformation calculation, the east-west systematic error is positively correlated with the north-south deformation component, and its proportionality constant is:

[0107] , ,

[0108] In the formula ,

[0109] .

[0110] Step S6: Based on the additional system error propagation model, the objective function is established by calculating the north-south deformation in 3D and ignoring the vertical and east-west deformations in the north-south solution and the LOS deformation. A regularization constraint term is added to the objective function. Specifically, the objective function is established based on the north-south component of the 3D deformation and the relationship between the vertical and east-west components in the north-south deformation solution and the LOS deformation. Based on the ill-conditioned nature of the objective function's coefficient matrix, a stable generalization function is added to the objective function using regularization. The objective function with added constraints is in the form of:

[0111]

[0112] In the formula, It is a Euclidean 2-norm. , These are the coefficients of the constraint terms.

[0113] In the objective function, after eliminating the LOS deformation error, vertical deformation component, and east-west deformation component from the LOS deformation, the result is the north-south deformation component. When the minimum value is taken, it is considered that... It approximates the north-south deformation component.

[0114] Step S7: Based on the objective function with constraints and the established error propagation model, the objective function matrix is ​​obtained. Specifically, according to the relationship between LOS error and deformation error in the error propagation model, the objective function is combined with the error propagation model. The multivariate function consists of six variables. It becomes three variables. And express it in matrix form. The objective function is in the form of:

[0115] ,

[0116] In the formula, For constraint term coefficients,

[0117] ,

[0118] .

[0119] Step S8: Based on the objective function matrix, obtain the LOS-directed deformation error that satisfies the minimum value of the objective function. Substitute the LOS-directed deformation error into the three-dimensional deformation error propagation model to obtain the three-dimensional deformation error. Specifically, establish a system of derivative equations for the objective function, use singular value decomposition to calculate the matrix form of the system of derivative equations, obtain the LOS-directed deformation error that satisfies the minimum value condition of the objective function, combine it with the error propagation model to obtain the three-dimensional deformation error, and obtain a high-precision three-dimensional deformation result based on the three-dimensional deformation error. The specific process is as follows:

[0120] Establish a system of derivative equations for the objective function:

[0121] ,

[0122] The matrix form of its system of derivative equations can be expressed as:

[0123] ,

[0124] Using the SVD method, the equation is solved, and the solution takes the following form:

[0125]

[0126] I represents the identity matrix.

[0127] The above is satisfied LOS deformation error Substituting into the error propagation model, the three-dimensional deformation error is obtained as follows:

[0128] ,

[0129] Step S9: Obtain the three-dimensional deformation result based on the three-dimensional deformation error.

[0130] Optimized north-south deformation is obtained by utilizing the three-dimensional directional deformation error. With east-west deformation and vertical deformation Combined into optimized high-precision three-dimensional deformation.

[0131] This embodiment of the high-precision 3D deformation extraction method using multi-track InSAR utilizes InSAR technology to obtain radar line-of-sight deformation; based on the radar side-looking imaging geometry, an observation equation set is established; using the least squares method, the 3D deformation (east-west, north-south, and vertical) is calculated, while neglecting the 2D deformation in the north-south direction; an error propagation model is established based on the relationship between the 3D deformation error and the line-of-sight deformation error; an additional systematic error propagation model is established based on neglecting the relationship between the vertical and east-west systematic errors and the north-south deformation components; an objective function is established based on the relationship between the 3D deformation components and the line-of-sight deformation, and a regularization constraint term is added; the objective function is combined with the error propagation model and expressed in matrix form; the objective function is differentiated, and the derivative equation is solved using the SVD method to obtain the LOS deformation error, which is then substituted into the error model to obtain the 3D deformation error, and high-precision 3D deformation is obtained using the 3D deformation error.

[0132] This embodiment uses three pairs of InSAR images to extract high-precision three-dimensional deformation of the land surface. It overcomes the dependence of the azimuth deformation component on radar imaging mode and high image coherence in traditional methods, and does not require deformation parameters of the study area and complex signal processing methods. Through mathematical methods, it greatly improves the north-south deformation accuracy and obtains high-precision three-dimensional deformation, providing a new technical solution for extracting high-precision three-dimensional deformation.

[0133] The electronic device disclosed in this embodiment includes a memory and a processor. The memory is used to store non-transitory computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0134] The processor may be a central processing unit (CPU) or other processing unit with data processing and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of this disclosure, the processor is used to execute computer-readable instructions stored in the memory, causing the electronic device to perform all or part of the steps of the multi-track InSAR high-precision three-dimensional deformation extraction method described in the foregoing embodiments of this disclosure.

[0135] Those skilled in the art will understand that, in order to solve the technical problem of how to achieve a good user experience, this embodiment may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included within the protection scope of this disclosure.

[0136] Electronic devices may include processing units (such as central processing units, graphics processing units, etc.) that can perform various appropriate actions and processes based on programs stored in read-only memory (ROM) or programs loaded from storage devices into random access memory (RAM). RAM also stores various programs and data required for the operation of the electronic device. The processing unit, ROM, and RAM are interconnected via a bus. Input / output (I / O) interfaces are also connected to the bus.

[0137] Typically, the following devices can be connected to the I / O interface: input devices, such as sensors or visual information acquisition devices; output devices, such as displays; storage devices, such as magnetic tapes or hard drives; and communication devices. Communication devices allow electronic devices to communicate wirelessly or wiredly with other devices (such as edge computing devices) to exchange data. However, it should be understood that it is not required to implement or possess all of the devices shown. More or fewer devices may be implemented or possessed alternatively.

[0138] In particular, according to embodiments of this disclosure, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, embodiments of this disclosure include a computer program product comprising a computer program carried on a non-transitory computer-readable medium, the computer program containing program code for performing the methods shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via a communication device, or installed from a storage device, or installed from a ROM. When the computer program is executed by a processing device, all or part of the steps of the multi-track InSAR method for extracting high-precision three-dimensional deformation according to embodiments of this disclosure are performed.

[0139] For a detailed description of this embodiment, please refer to the corresponding descriptions in the foregoing embodiments, which will not be repeated here.

[0140] A computer-readable storage medium according to embodiments of the present disclosure stores non-transitory computer-readable instructions. When these non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the multi-track InSAR high-precision three-dimensional deformation extraction method described in the foregoing embodiments of the present disclosure are performed.

[0141] The aforementioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or portable hard drive), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).

[0142] For a detailed description of this embodiment, please refer to the corresponding descriptions in the foregoing embodiments, which will not be repeated here.

[0143] The basic principles of this disclosure have been described above with reference to specific embodiments. However, it should be noted that the advantages, benefits, and effects mentioned in this disclosure are merely examples and not limitations, and should not be considered as essential features of each embodiment of this disclosure. Furthermore, the specific details disclosed above are for illustrative and facilitative purposes only, and are not limitations. These details do not limit the scope of this disclosure to the necessity of employing the aforementioned specific details for implementation.

[0144] In this disclosure, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. The block diagrams of devices, apparatuses, devices, and systems involved in this disclosure are merely illustrative examples and are not intended to require or imply that they must be connected, arranged, or configured in the manner shown in the block diagrams. As those skilled in the art will recognize, these devices, apparatuses, devices, and systems can be connected, arranged, and configured in any manner. Words such as "comprising," "including," "having," etc., are open-ended terms meaning "including but not limited to," and are used interchangeably with them. The terms "or" and "and" as used herein refer to the terms "and / or," and are used interchangeably with them unless the context clearly indicates otherwise. The term "such as" as used herein refers to the phrase "such as but not limited to," and is used interchangeably with it.

[0145] Additionally, as used herein, the "or" used in a list of items beginning with "at least one" indicates a separate list, such that a list of, for example, "at least one of A, B, or C" means A or B or C, or AB or AC or BC, or ABC (i.e., A and B and C). Furthermore, the word "exemplary" does not imply that the described example is preferred or better than other examples.

[0146] It should also be noted that in the systems and methods of this disclosure, the components or steps can be decomposed and / or recombined. These decompositions and / or recombinations should be considered as equivalent solutions to this disclosure.

[0147] Various changes, substitutions, and modifications can be made to the technology described herein without departing from the teachings defined by the appended claims. Furthermore, the scope of the claims of this disclosure is not limited to the specific aspects of the processes, machines, manufactures, events, means, methods, and actions described above. Currently existing or later-developed processes, machines, manufactures, events, means, methods, or actions that perform substantially the same function or achieve substantially the same result as the corresponding aspects described herein can be utilized. Therefore, the appended claims include such processes, machines, manufactures, events, means, methods, or actions within their scope.

[0148] The above description of the disclosed aspects is provided to enable any person skilled in the art to make or use this disclosure. Various modifications to these aspects will be readily apparent to those skilled in the art, and the general principles defined herein may be applied to other aspects without departing from the scope of this disclosure. Therefore, this disclosure is not intended to be limited to the aspects shown herein, but rather to be carried out within the widest scope consistent with the principles and novel features disclosed herein.

[0149] The above description has been given for purposes of illustration and description. Furthermore, this description is not intended to limit the embodiments of this disclosure to the forms disclosed herein. Although numerous exemplary aspects and embodiments have been discussed above, those skilled in the art will recognize certain variations, modifications, alterations, additions, and sub-combinations therein.

Claims

1. A method for extracting high-precision three-dimensional deformation using multi-track InSAR, characterized in that, include: The LOS deformation of each orbit is obtained from the SAR data of the three orbits, including processing the SAR data of the three orbits separately using InSAR technology to obtain the LOS deformation displacement of each orbit. Based on the side-view imaging geometry of the SAR sensor, a set of equations relating LOS deformation to three-dimensional deformation is established, and a set of equations relating LOS deformation to two-dimensional deformation ignoring north-south deformation is also established. Based on the equations relating LOS deformation to 3D deformation and LOS deformation to 2D deformation, the least squares method is used to obtain the 3D deformation results as well as the vertical and east-west deformations. Based on the relationship between LOS deformation error and deformation result error in three dimensions, a three-dimensional deformation error propagation model is established. Based on the relationship between LOS deformation error and deformation errors in the vertical and east-west directions, a LOS deformation error and two-dimensional deformation error propagation model is established. Based on the relationship between the vertical and east-west systematic errors and the north-south deformation components in the two-dimensional deformation error propagation model, an additional systematic error propagation model is established. An objective function is established based on the additional system error propagation model, the relationship between the three-dimensional solution of north-south deformation and the vertical and east-west deformations ignoring the north-south solution and the LOS deformation, and a regularization constraint term is added to the objective function. Based on the objective function with constraints and the established error propagation model, the objective function matrix is ​​obtained; The LOS-directed deformation error that satisfies the minimum value of the objective function matrix is ​​obtained. The LOS-directed deformation error is then substituted into the three-dimensional deformation error propagation model to obtain the three-dimensional deformation error. The three-dimensional deformation result is obtained based on the three-dimensional deformation error.

2. The method for extracting high-precision three-dimensional deformation using multi-track InSAR according to claim 1, characterized in that, The equations establishing the relationship between LOS deformation and three-dimensional deformation are as follows: ， In the formula, This indicates the LOS-direction deformation displacement of orbit i. ; , This indicates the radar satellite's incident angle and the clockwise angle between the north direction and the radar satellite's flight direction. ， ， , , The deformations are three-dimensional, namely vertical, east-west, and north-south.

3. The multi-track InSAR method for extracting high-precision three-dimensional deformation according to claim 2, characterized in that, In establishing the system of equations relating LOS-direction deformation to neglecting north-south deformation in two dimensions, the system of equations is as follows: ， In the formula, , To ignore the vertical and east-west deformation when the north-south direction is ignored.

4. The method for extracting high-precision three-dimensional deformation using multi-track InSAR according to claim 3, characterized in that, Based on the equations relating LOS deformation to 3D deformation and LOS deformation to 2D deformation, the least squares method is used to obtain the deformation results in the 3D direction, as well as the deformation in the vertical and east-west directions. The formulas are as follows: ， In the formula for: ， The deformation along the radar line of sight is as follows: ， W is the LOS-oriented deformation weight matrix, using an equivalent weight matrix. The calculated three-dimensional deformation, or the two-dimensional deformation obtained by ignoring the north-south direction, is as follows: 。 5. The multi-track InSAR method for extracting high-precision three-dimensional deformation according to claim 4, characterized in that, The three-dimensional deformation error propagation model is established based on the relationship between the LOS deformation error and the deformation result error in the three-dimensional direction as follows: ， In the formula, , In the formula For i-orbital LOS-direction deformation error, , , These are the deformation coefficients of orbit i in the vertical, east-west, and north-south directions, respectively. , , These represent the deformation errors in the three-dimensional direction.

6. The multi-track InSAR method for extracting high-precision three-dimensional deformation according to claim 5, characterized in that, Based on the relationship between LOS deformation error and vertical and east-west deformation errors, the LOS deformation error and two-dimensional deformation error propagation model are established as follows: ， In the formula, , In the formula, , These represent the vertical and east-west deformation errors when the north-south direction is ignored.

7. The multi-track InSAR method for extracting high-precision three-dimensional deformation according to claim 6, characterized in that, The additional systematic error propagation model is established based on the relationship between the vertical and east-west systematic errors and the north-south deformation components in the two-dimensional deformation error propagation model as follows: ， In the formula, This is the error in the vertical deformation system. For east-west deformation systematic error, Ignoring the vertical systematic error caused by the north-south deformation calculation, the east-west systematic error is positively correlated with the north-south deformation component, and its proportionality constant is: ， ， In the formula , 。 8. The multi-track InSAR method for extracting high-precision three-dimensional deformation according to claim 7, characterized in that, An objective function is established based on the additional systematic error propagation model, the relationship between the three-dimensional solution of north-south deformation and the vertical and east-west deformations neglecting the north-south solution, and the LOS deformation. Regularization constraints are added to the objective function, including: Add a stable generalization function to the objective function The objective function with added constraints is in the form of: In the formula, It is a Euclidean 2-norm. , These are the coefficients of the constraint terms.

9. The method for extracting high-precision three-dimensional deformation using multi-track InSAR according to claim 8, characterized in that, The LOS-directed deformation error, which satisfies the minimum value of the objective function matrix, is obtained by substituting the LOS-directed deformation error into the three-dimensional deformation error propagation model, resulting in the three-dimensional deformation error, including: Establish a system of derivative equations for the objective function: ， The system of derivative equations can be represented as a matrix: ， The solution to the equation is: ， The above is satisfied LOS deformation error Substituting into the error propagation model, the three-dimensional deformation error is obtained as follows: ， Optimized north-south deformation is obtained by utilizing the three-dimensional deformation error. With east-west deformation and vertical deformation The optimized three-dimensional deformation is obtained by combining these factors.

10. An electronic device, characterized in that, The electronic device includes: At least one processor; and, A memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the multi-track InSAR high-precision three-dimensional deformation extraction method according to any one of claims 1-9.

Images