CN-115951353-B - Multi-track InSAR extraction high-precision three-dimensional deformation method

CN115951353BCN 115951353 BCN115951353 BCN 115951353BCN-115951353-B

Abstract

The invention discloses a multi-track InSAR extraction high-precision three-dimensional deformation method which comprises the steps of obtaining multi-track LOS deformation, establishing an LOS deformation and three-dimensional deformation relation equation set, establishing an LOS deformation and two-dimensional deformation relation equation set, obtaining three-dimensional deformation, vertical deformation and east-west deformation, establishing a three-dimensional and two-dimensional deformation error transfer model, establishing an additional system error transfer model, establishing an objective function, adding regularization constraint terms to the objective function, obtaining an objective function matrix based on the objective function with constraint terms in combination with a set error transfer model, obtaining LOS deformation errors meeting minimum values of the objective function based on the objective function matrix, substituting the LOS deformation errors into the three-dimensional and two-dimensional deformation error transfer model to obtain three-dimensional deformation errors, and obtaining a three-dimensional deformation result according to the three-dimensional deformation errors. The purpose of improving the three-dimensional deformation precision is achieved, and particularly the deformation precision in the north-south direction is remarkably improved.

Inventors

ZHANG LELE
KONG XIAOLAN

Assignees

内蒙古工业大学

Dates

Publication Date: 20260505
Application Date: 20230113

Claims (10)

1. The multi-track InSAR extraction high-precision three-dimensional deformation method is characterized by comprising the following steps of: the LOS deformation of each track is obtained from the SAR data of multiple tracks, and the LOS deformation displacement of each track is obtained by processing the SAR data of three tracks by adopting an InSAR technology respectively; According to the side view imaging geometrical relationship of the SAR sensor, establishing a relationship equation set of LOS deformation and three-dimensional deformation, and establishing a relationship equation set of LOS deformation and two-dimensional deformation neglecting north-south deformation; Obtaining a three-dimensional deformation result and vertical and east-west deformation by using a least square method according to the relation equation set of LOS deformation and three-dimensional deformation and the relation equation set of LOS deformation and two-dimensional deformation; According to the relation between the LOS deformation errors and the three-dimensional deformation results, a three-dimensional deformation error transfer model is established, and according to the relation between the LOS deformation errors and the vertical deformation errors and the east-west deformation errors, an LOS deformation error and a two-dimensional deformation error transfer model are established; establishing an additional system error transfer model according to the relation between the vertical direction of the two-dimensional deformation error transfer model and the northeast and west systematic errors and the north-south deformation components; Establishing an objective function based on an additional system error transfer model, three-dimensionally resolving north-south deformation, ignoring the relationship between north-south resolved vertical deformation and east-west deformation and LOS deformation, and adding a regularization constraint term to the objective function; Obtaining an objective function matrix based on an objective function with constraint terms and a set error transfer model; obtaining LOS directional deformation errors meeting minimum values of the objective function based on the objective function matrix, substituting the LOS directional deformation errors into a three-dimensional deformation error transfer model, and obtaining three-dimensional deformation errors; And obtaining a three-dimensional deformation result according to the three-dimensional deformation error.
2. The multi-orbit InSAR extraction high-precision three-dimensional deformation method according to claim 1, wherein the establishing LOS-direction deformation and three-dimensional-direction deformation relation equation set is as follows: , In the formula, Indicating i the displacement of the track LOS towards deformation, ; , Indicating the angle of incidence of the radar satellite and the angle of north to the clockwise direction of the radar satellite's flight, , , , The three-dimensional deformation is respectively vertical, east-west and north-south.
3. The multi-orbit InSAR extraction high-precision three-dimensional deformation method according to claim 2, wherein a two-dimensional deformation relation equation set of LOS deformation and neglecting north-south deformation is established, and the equation set is as follows: , In the formula, , To ignore north-south direction vertical deformation in the east-west direction.
4. The multi-orbit InSAR extraction high-precision three-dimensional deformation method according to claim 3, wherein the formula is as follows in the three-dimensional deformation result and the vertical and east-west deformation obtained by the least square method according to the LOS deformation and three-dimensional deformation relation equation set and the LOS deformation and two-dimensional deformation relation equation set: , In the middle of The method comprises the following steps: , The deformation of the radar vision line is respectively: , W is an LOS deformation weight matrix, an equivalent weight matrix is adopted, To solve the resulting three-dimensional deformation or to ignore the two-dimensional deformation from north to south, the resulting two-dimensional deformation is: 。
5. The multi-orbit InSAR extraction high-precision three-dimensional deformation method according to claim 4, wherein the establishing a three-dimensional deformation error transfer model according to the relationship between the LOS deformation error and the three-dimensional deformation result error is: , In the formula, , In the middle of For i track LOS-direction deformation errors, , , The deformation coefficients of the track in the vertical direction, the east-west direction and the north-south direction are respectively i, , , Respectively, three-dimensional direction deformation errors.
6. The multi-orbit InSAR extraction high-precision three-dimensional deformation method according to claim 5, wherein the establishment of the LOS-oriented deformation error and two-dimensional deformation error transfer model according to the relationship between the LOS-oriented deformation error and the vertical and east-west deformation errors is as follows: , In the formula, , In the formula, , And the deformation errors in the vertical direction and the northeast and west directions are ignored respectively.
7. The multi-orbit InSAR extraction high-precision three-dimensional deformation method according to claim 6, wherein the establishing additional system error transfer model according to the relationship between the vertical direction and the northeast and west system error of the two-dimensional deformation error transfer model and the north and south deformation component is as follows: , In the formula, Is a systematic error in the vertical deformation, For the systematic error of the east-west deformation, Ignoring the positive correlation between the vertical systematic error, the east-west systematic error and the north-south deformation component caused by the north-south deformation calculation, wherein the proportionality constant is as follows: , , In the middle of , 。
8. The multi-track InSAR extraction high precision three-dimensional deformation method of claim 7, wherein creating an objective function based on an additional systematic error transfer model, three-dimensionally resolving north-south deformations and ignoring the relationship of north-south resolved vertical and east-west deformations with LOS-directed deformations, and adding regularization constraint terms to the objective function, comprising: adding a stable generalization function to an objective function The objective function form of adding constraint terms is: In the formula, Is the European 2-norm, , Is a constraint term coefficient.
9. The multi-orbit InSAR extraction high-precision three-dimensional deformation method according to claim 8, wherein the obtaining the LOS-direction deformation error satisfying the minimum value of the objective function based on the objective function matrix, substituting the LOS-direction deformation error into the three-dimensional deformation error transfer model, and obtaining the three-dimensional deformation error, comprises: establishing a derivative equation set for the objective function: , the system of derivative equations is represented as a matrix: , solving the equation: , the above is satisfied LOS directional deformation error of (C) Substituting the three-dimensional deformation error obtained by the error transfer model is as follows: , the three-dimensional deformation error is utilized to obtain the optimized north-south deformation And the deformation in the east-west direction And vertical deformation And combining to obtain the optimized three-dimensional deformation.
10. An electronic device, the electronic device comprising: at least one processor, and A memory communicatively coupled to the at least one processor, wherein, The memory stores instructions executable by the at least one processor to enable the at least one processor to perform the multi-track InSAR extraction high-precision three-dimensional morphing method of any one of claims 1-9.

Description

Multi-track InSAR extraction high-precision three-dimensional deformation method Technical Field The invention relates to the field of microwave remote sensing measurement, in particular to a multi-track InSAR extraction high-precision three-dimensional deformation method. Background Interferometric synthetic aperture radar (InSAR) technology has been an effective method for monitoring landslide, earthquakes, volcanic eruptions, mining subsidence and other geological disasters in the past two decades. However, the InSAR technique only provides one-dimensional deformation of the radar sight Line (LOS for short). In practical applications, LOS-direction distortion is not suitable for evaluating disaster ranges in many cases. In recent years, as the number of SAR satellites increases, it becomes possible to acquire multiple sets of LOS-direction deformation data of different heading and incidence angles (orbit up and orbit down data, or left-view and right-view data). And 3 groups of LOS deformation data are utilized to establish an observation equation set, so that the vertical, east-west and north-south deformation data can be solved. But limited by the near-north-south flight orbit of SAR satellites, the accuracy of the north-south solution is far lower than that of the vertical, east-west solution. Disclosure of Invention In order to solve the technical problems in the prior art, the invention provides a multi-track InSAR extraction high-precision three-dimensional deformation method. The method comprises the following steps: In a first aspect of the present invention, a method for extracting high-precision three-dimensional deformation by multi-track InSAR is provided, comprising: the LOS deformation of each track is obtained from multi-track SAR data, which comprises the steps of respectively processing the SAR data of three tracks by adopting an InSAR technology to obtain LOS deformation displacement of each track; According to the side view imaging geometrical relationship of the SAR sensor, establishing a relationship equation set of LOS deformation and three-dimensional deformation, and establishing a relationship equation set of LOS deformation and two-dimensional deformation neglecting north-south deformation; Obtaining a three-dimensional deformation result and vertical and east-west deformation by using a least square method according to the relation equation set of LOS deformation and three-dimensional deformation and the relation equation set of LOS deformation and two-dimensional deformation; According to the relation between the LOS deformation errors and the three-dimensional deformation results, a three-dimensional deformation error transfer model is established, and according to the relation between the LOS deformation errors and the vertical deformation errors and the east-west deformation errors, an LOS deformation error and a two-dimensional deformation error transfer model are established; establishing an additional system error transfer model according to the relation between the vertical direction of the two-dimensional deformation error transfer model and the northeast and west systematic errors and the north-south deformation components; Establishing an objective function based on an additional system error transfer model, three-dimensionally resolving north-south deformation, ignoring the relationship between north-south resolved vertical deformation and east-west deformation and LOS deformation, and adding a regularization constraint term to the objective function; Obtaining an objective function matrix based on an objective function with constraint terms and a set error transfer model; obtaining LOS directional deformation errors meeting minimum values of the objective function based on the objective function matrix, substituting the LOS directional deformation errors into a three-dimensional deformation error transfer model, and obtaining three-dimensional deformation errors; And obtaining a three-dimensional deformation result according to the three-dimensional deformation error. Optionally, in the establishing a relationship equation set between LOS deformation and three-dimensional deformation, the equation set is: , In the formula, Indicating i the displacement of the track LOS towards deformation,; ,Indicating the angle of incidence of the radar satellite and the angle of north to the clockwise direction of the radar satellite's flight,。 ,,The three-dimensional deformation is respectively vertical, east-west and north-south. Optionally, in a two-dimensional direction deformation relation equation set of LOS deformation and neglecting north-south deformation, the equation set is: , In the formula, ,To ignore north-south direction vertical deformation in the east-west direction. Optionally, the formula is that, according to the relationship equation set of the deformation in the LOS direction and the three-dimensional direction and the relationship equation set of the deform