CN-121739966-B - Method, device, medium and equipment for acquiring three-dimensional deformation field data

CN121739966BCN 121739966 BCN121739966 BCN 121739966BCN-121739966-B

Abstract

The application relates to the technical field of surface measurement, in particular to a method, a device, a medium and equipment for acquiring three-dimensional deformation field data, wherein the method can comprise the steps of constructing a three-dimensional deformation model to be solved based on multi-source data of observation points in any observation window in a target research area; the multi-source data comprises synthetic aperture radar interferometry data, system observation data of a global satellite navigation positioning system and level measurement values, observation residual errors of each type of data are obtained based on a weight array of each type of data in the multi-source data and the three-dimensional deformation model, the observation residual errors are calculated to obtain variance component estimated values, and three-dimensional deformation component data in the three-dimensional deformation model are output under the condition that the variance component estimated values meet preset conditions. The embodiment of the application can efficiently and accurately acquire the three-dimensional deformation data of the earth surface.

Inventors

GUO NANNAN
ZHAN WEI
CHEN CHANGYUN

Assignees

中国地震局第一监测中心

Dates

Publication Date: 20260508
Application Date: 20260227

Claims (10)

1. A method of acquiring three-dimensional deformation field data, comprising: Constructing a three-dimensional deformation model to be solved based on multi-source data of observation points in any observation window in a target research area, wherein the multi-source data comprises synthetic aperture radar interferometry data, system observation data of a global satellite navigation positioning system and level measurement values; acquiring an observation residual error of each type of data based on the weight array of each type of data in the multi-source data and the three-dimensional deformation model; Calculating the observation residual to obtain a variance component estimated value, wherein the variance component estimated value is calculated by a Welch variance estimation method; Outputting three-dimensional deformation component data in the three-dimensional deformation quantity model under the condition that the variance component estimated value meets a preset condition; the three-dimensional deformation model is as follows: Wherein, the Each parameter of the matrix is respectively expressed as the observation data of GNSS east-west direction, south-north direction and vertical direction, the vertical leveling value and the observation values of the ascending track and the descending track vision line in InSAR data; And The angle of incidence and the azimuth angle under the SAR image reference system are respectively expressed, and a and d respectively express the ascending track and the descending track deformation observation of InSAR; three-dimensional deformation component data (namely a data matrix) to be solved, which represents components in the east-west direction, the north-south direction and the vertical direction; the calculation method of the Welch variance component estimation value comprises the following steps: Wherein, the For the variance component estimate value, In order to observe the transposed matrix of the residual, For the number of observation points of each type of data, tr () is the trace of the matrix; the target research area is divided into a plurality of observation windows, and model construction, observation residual acquisition, variance component estimation and three-dimensional deformation component data output are respectively executed in parallel or in series for each observation window.
2. The method of claim 1, wherein prior to constructing the three-dimensional deformation model to be solved based on multi-source data of observation points in any observation window in the target study area, the method further comprises: And determining the reference coordinates and the spatial resolution of the three-dimensional deformation area according to the observation range of the synthetic aperture radar interferometry and the distribution condition of the sites and the level points of the global satellite navigation positioning system, wherein the reference coordinates and the spatial resolution are used for enabling the multi-source data to be under a unified frame.
3. The method of claim 1 or 2, wherein the multi-source data is obtained by: and inserting the system observation data and the leveling value into the synthetic aperture radar interferometry data to obtain the multi-source data.
4. The method according to claim 1 or 2, wherein the constructing a three-dimensional deformation model to be solved based on multi-source data of observation points in any observation window in the target study area includes: Taking an observation value matrix constructed by the multi-source data as an output side, and taking the product of a data matrix of the three-dimensional deformation component data to be solved and an observation design matrix as an input side; And taking the equality relation between the output side and the input side as the three-dimensional deformation model.
5. The method of claim 4, wherein the obtaining an observation residual for each type of data based on the weight matrix for each type of data in the multi-source data and the three-dimensional deformation model comprises: obtaining an observation value of each type of data and the observation design matrix from the three-dimensional deformation model; And calculating the observation value, the observation design matrix and the weight matrix to obtain the observation residual error.
6. The method of claim 4, wherein said calculating the observation residual to obtain the variance component estimate comprises: and solving the transposed matrix of the observation residual, the weight matrix and the ratio of the product of the observation residual to the number of observation points to obtain the variance component estimated value.
7. The method according to claim 1 or 2, wherein the outputting three-dimensional deformation component data in the three-dimensional deformation model in a case where the variance component estimation value satisfies a preset condition includes: The following operations are circularly executed until the variance component estimated value remains unchanged, and the three-dimensional deformation component data is output: updating the weight matrix based on the variance component estimate; Updating the observation residual by using the weight matrix; updating the variance component estimate by the observation residual; confirming that the variance component estimate changes.
8. An apparatus for acquiring three-dimensional deformation field data, comprising: The model construction module is used for constructing a three-dimensional deformation model to be solved based on multi-source data of observation points in any observation window in a target research area, wherein the multi-source data comprises synthetic aperture radar interferometry data, system observation data of a global satellite navigation positioning system and leveling values; the residual error acquisition module is used for acquiring the observation residual error of each type of data based on the weight array of each type of data in the multi-source data and the three-dimensional deformation model; The calculation module is used for calculating the observation residual to obtain a variance component estimation value, wherein the variance component estimation value is calculated by a Welch variance estimation method; The output module is used for outputting three-dimensional deformation component data in the three-dimensional deformation quantity model under the condition that the variance component estimated value meets a preset condition; the three-dimensional deformation model is as follows: Wherein, the Each parameter of the matrix is respectively expressed as the observation data of GNSS east-west direction, south-north direction and vertical direction, the vertical leveling value and the observation values of the ascending track and the descending track vision line in InSAR data; And The angle of incidence and the azimuth angle under the SAR image reference system are respectively expressed, and a and d respectively express the ascending track and the descending track deformation observation of InSAR; three-dimensional deformation component data (namely a data matrix) to be solved, which represents components in the east-west direction, the north-south direction and the vertical direction; The calculating method of the variance component estimated value comprises the following steps: Wherein, the For the variance component estimate value, In order to observe the transposed matrix of the residual, For the number of observation points of each type of data, tr () is the trace of the matrix; the target research area is divided into a plurality of observation windows, and model construction, observation residual acquisition, variance component estimation and three-dimensional deformation component data output are respectively executed in parallel or in series for each observation window.
9. A computer readable storage medium, characterized in that the computer readable storage medium has stored thereon a computer program, wherein the computer program when run by a processor performs the method according to any of claims 1-7.
10. An electronic device comprising a memory, a processor, and a computer program stored on the memory and running on the processor, wherein the computer program when run by the processor performs the method of any one of claims 1-7.

Description

Method, device, medium and equipment for acquiring three-dimensional deformation field data Technical Field The application relates to the technical field of surface measurement, in particular to a method, a device, a medium and equipment for acquiring three-dimensional deformation field data. Background The surface deformation measurement refers to the observation and calculation process of geometrical and physical parameters such as position change, displacement, speed, strain and the like of the earth surface in time and space by using technical means such as geodetic measurement, remote sensing, geophysical and the like, and is a core technical means for monitoring the surface dynamic change. Synthetic aperture radar interferometry (InSAR, interferometric Synthetic Aperture Radar) is a side-looking remote sensing technology capable of acquiring earth surface deformation information in a high spatial resolution and large-scale coverage mode, but the observed value is only one-dimensional projection of radar vision, and horizontal deformation components and vertical deformation components are difficult to directly separate. The global satellite navigation positioning system can acquire three-dimensional deformation with high time resolution and high precision, but has low spatial resolution and higher station building cost. Leveling as a classical geodetic measure can directly reflect high-precision vertical deformations, but it is difficult to achieve large-scale intensive observations. Although each of the above techniques has its own advantages, there are also significant limitations. Therefore, how to provide a more accurate method for acquiring three-dimensional deformation field data becomes a technical problem to be solved. Disclosure of Invention The application aims to provide a method, a device, a medium and equipment for acquiring three-dimensional deformation field data, and the technical scheme of the embodiment of the application can improve the accuracy and efficiency of acquiring the three-dimensional deformation field data of the earth surface and has higher practicability. According to the method, a three-dimensional deformation field data is obtained, the method comprises the steps of constructing a three-dimensional deformation model to be solved based on multi-source data of observation points in any observation window in a target research area, wherein the multi-source data comprise synthetic aperture radar interferometry data, system observation data of a global satellite navigation positioning system and level measurement values, obtaining observation residual errors of each type of data based on a weight array of each type of data in the multi-source data and the three-dimensional deformation model, calculating the observation residual errors to obtain variance component estimated values, and outputting the three-dimensional deformation component data in the three-dimensional deformation model under the condition that the variance component estimated values meet preset conditions. According to some embodiments of the application, a corresponding three-dimensional deformation quantity model is constructed on multi-source data of observation points in any observation window, then an observation residual error calculated by combining weight arrays of each type of data is used for determining a variance component estimated value, and corresponding three-dimensional deformation component data is output under the condition that the variance component estimated value reaches a preset condition. According to the embodiment of the application, the single calculated amount can be reduced in a window division mode, meanwhile, the three-dimensional deformation quantity model is constructed and solved, the accuracy and efficiency of acquiring the data of the ground surface three-dimensional deformation field can be improved, and the practicability is high. In some embodiments, before constructing the three-dimensional deformation model to be solved based on multi-source data of observation points in any observation window in the target research area, the method further comprises determining reference coordinates and spatial resolution of the three-dimensional deformation area according to the observation range of synthetic aperture radar interferometry and the site and level point distribution condition of the global satellite navigation positioning system, wherein the reference coordinates and the spatial resolution are used for enabling the multi-source data to be under a unified frame. According to the embodiments of the application, the consistency of multi-source data acquisition can be ensured and the standardized processing of the data can be realized by constructing uniform reference coordinates and spatial resolution for the synthetic aperture radar interferometry, the global satellite navigation positioning system and the level points. In some embodiments, the multi-source data is ob