Optimizing the temporal resolution of BeiDou-assisted InSAR surface deformation series

Abstract

This application relates to the field of high-precision deformation monitoring technology, and in particular to a method for optimizing the temporal resolution of Beidou-assisted InSAR surface deformation sequences. The method comprises: extracting one-dimensional surface deformations based on SBAS-InSAR; extracting three-dimensional deformations based on Beidou / GNSS monitoring points; converting Beidou three-dimensional deformation sequences into one-dimensional deformations along the InSAR radar line of sight; constructing a high-quality Beidou and InSAR monitoring dataset; constructing a deformation law fusion model driven by Beidou and InSAR monitoring data; and obtaining high-temporal-resolution InSAR surface deformations. This method solves the problems of existing surface deformation monitoring technologies, such as high monitoring costs, long station deployment cycles, and low spatiotemporal resolution.

Description

Technical Field

[0001] The present application relates to the technical field of high-precision deformation monitoring, and in particular to a method for optimizing the temporal resolution of a BeiDou-assisted InSAR surface deformation sequence. Background Art

[0002] Frequent geological disasters pose a serious threat to people's lives and property. To monitor surface deformation and prevent geological disasters, monitoring technologies such as Interferometric Synthetic Aperture Radar (InSAR) and Global Navigation Satellite System (GNSS) have emerged as emerging technologies.

[0003] InSAR technology offers advantages such as 24 / 7 operation and high spatial resolution. However, it can only observe deformation in one-dimensional line of sight and has low temporal resolution. GNSS technology is an automated monitoring technology that can obtain high-precision three-dimensional surface deformation in real time, with a maximum accuracy of millimeters. However, compared to InSAR technology, GNSS ground monitoring equipment is more expensive and requires a longer deployment cycle, resulting in lower spatial resolution.

[0004] In summary, the existing surface deformation monitoring technology has high monitoring costs, long station deployment cycles, and low temporal and spatial resolution, which urgently need to be addressed. Summary of the Invention

[0005] This application provides a method for optimizing the temporal resolution of BeiDou-assisted InSAR surface deformation series to solve the problems of existing surface deformation monitoring technologies, such as high monitoring cost, long station deployment cycle, and low temporal and spatial resolution.

[0006] The first aspect of the present application provides a method for optimizing the temporal resolution of a Beidou-assisted InSAR surface deformation sequence, comprising the following steps: extracting one-dimensional surface deformation information of a target monitoring area based on a preset SBAS-InSAR strategy; deploying a Beidou / GNSS receiver in the target monitoring area to obtain Beidou three-dimensional point deformation information of the target monitoring area, and converting the Beidou three-dimensional point deformation information into one-dimensional data in the LOS direction, and constructing a Beidou monitoring dataset and an InSAR monitoring dataset of the target detection area based on the one-dimensional data in the LOS direction and the one-dimensional surface deformation information; constructing a deformation law fusion model corresponding to Beidou and InSAR at the same time based on the Beidou monitoring dataset and the InSAR monitoring dataset, so as to obtain the InSAR surface deformation amount of the target monitoring area based on the deformation law fusion model and the Beidou monitoring dataset.

[0007] Optionally, in one embodiment of the present application, the Beidou three-dimensional point deformation information is converted into LOS one-dimensional data, and the Beidou monitoring data set and the InSAR monitoring data set of the target detection area are constructed based on the LOS one-dimensional data and the one-dimensional surface deformation information, including: obtaining the azimuth and altitude of the InSAR remote sensing satellite image data corresponding to the target detection area, and converting the Beidou three-dimensional point deformation information into the LOS one-dimensional data based on the azimuth and the altitude; respectively determining the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference corresponding to the LOS one-dimensional data and the one-dimensional surface deformation information; based on the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference, performing a spatial reference unification operation on the LOS one-dimensional data and the one-dimensional surface deformation information to obtain the LOS one-dimensional data and the one-dimensional surface deformation information under a unified coordinate reference, and establishing the Beidou monitoring data set and the InSAR monitoring data set based on the LOS one-dimensional data and the one-dimensional surface deformation information under a unified coordinate reference.

[0008] Optionally, in one embodiment of the present application, based on the BeiDou monitoring data set and the InSAR monitoring data set, a deformation law fusion model corresponding to the BeiDou and InSAR at the same time is constructed to obtain the InSAR surface deformation of the target monitoring area according to the deformation law fusion model and the BeiDou monitoring data set, including: obtaining the deformation law fusion model of the target monitoring area; m Beidou deformation monitoring points M Monitor deformation over time, and n InSAR monitoring pixels N Monitor the deformation periodically and determine m and n The first quantitative relationship between N and M The second quantitative relationship between n 、 m 、 N 、 M are all positive integers, and n > m , N < M ; Based on the first quantitative relationship, the second quantitative relationship, the BeiDou monitoring data set and the InSAR monitoring data set, a monitoring training data set is constructed; a pre-constructed neural network model is trained by the monitoring training data set to generate the deformation law fusion model; the deformation law fusion model is interpolated according to the deformation law fusion model n InSAR monitoring pixels ( M - N) period of monitoring deformation, and combining the n InSAR monitoring pixels ( M - N ) period of monitoring deformation and the N The deformation is monitored regularly to generate the InSAR surface deformation of the target monitoring area.

[0009] Optionally, in one embodiment of the present application, the mathematical expression of the deformation law fusion model is:

[0010]

[0011] in, f DNN Representing the pre-built neural network model; w Represents the weight parameter of the deformation law fusion model; b represents the bias parameter of the deformation law fusion model; d BDS,m×1 Representing the LOS as one-dimensional data; d InSAR,n×1 Indicates the one-dimensional surface deformation information.

[0012] The second aspect of the present application provides an optimization device for the time resolution of a Beidou-assisted InSAR surface deformation sequence, including: an extraction module for extracting one-dimensional surface deformation information of a target monitoring area based on a preset SBAS-InSAR strategy; a conversion module for deploying Beidou / GNSS receivers in the target monitoring area to obtain Beidou three-dimensional point deformation information of the target monitoring area, and converting the Beidou three-dimensional point deformation information into LOS-direction one-dimensional data, and constructing a Beidou monitoring data set and an InSAR monitoring data set of the target detection area based on the LOS-direction one-dimensional data and the one-dimensional surface deformation information; a fusion module for constructing a deformation law fusion model corresponding to Beidou and InSAR at the same time based on the Beidou monitoring data set and the InSAR monitoring data set, so as to obtain the InSAR surface deformation amount of the target monitoring area based on the deformation law fusion model and the Beidou monitoring data set.

[0013] Optionally, in one embodiment of the present application, the conversion module includes: a first acquisition unit, used to acquire the azimuth and altitude of the InSAR remote sensing satellite image data corresponding to the target detection area, and based on the azimuth and altitude, convert the Beidou three-dimensional point deformation information into the LOS one-dimensional data; a determination unit, used to respectively determine the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference corresponding to the LOS one-dimensional data and the one-dimensional surface deformation information; a spatial reference unification unit, used to perform a spatial reference unification operation on the LOS one-dimensional data and the one-dimensional surface deformation information based on the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference, so as to obtain the LOS one-dimensional data and the one-dimensional surface deformation information under a unified coordinate reference, and establish the Beidou monitoring data set and the InSAR monitoring data set based on the LOS one-dimensional data and the one-dimensional surface deformation information under a unified coordinate reference.

[0014] Optionally, in one embodiment of the present application, the fusion module includes: a second acquisition unit for acquiring the target detection area m Beidou deformation monitoring points M Monitor deformation over time, and n InSAR monitoring pixels N Monitor the deformation periodically and determine m and n The first quantitative relationship between N and M The second quantitative relationship between n 、 m 、 N 、 M are all positive integers, and n > m , N < M A construction unit is used to construct a monitoring training data set based on the first quantitative relationship, the second quantitative relationship, the BeiDou monitoring data set and the InSAR monitoring data set; a training unit is used to train a pre-constructed neural network model through the monitoring training data set to generate the deformation law fusion model; a combination unit is used to interpolate the deformation law fusion model according to the deformation law fusion model. n InSAR monitoring pixels ( M - N ) period of monitoring deformation, and combining the n InSAR monitoring pixels ( M - N ) period of monitoring deformation and the N The deformation is monitored regularly to generate the InSAR surface deformation of the target monitoring area.

[0015] Optionally, in one embodiment of the present application, the mathematical expression of the deformation law fusion model is:

[0016]

[0017] in, f DNN Representing the pre-built neural network model; w Represents the weight parameter of the deformation law fusion model; b represents the bias parameter of the deformation law fusion model; d BDS,m×1 Representing the LOS as one-dimensional data; d InSAR,n×1 Indicates the one-dimensional surface deformation information.

[0018] The third aspect of the present application provides an electronic device, comprising: a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor executes the program to implement the method for optimizing the temporal resolution of the Beidou-assisted InSAR surface deformation series as described in the above embodiment.

[0019] The fourth aspect of the present application provides a computer-readable storage medium, which stores a computer program. When the program is executed by a processor, it implements the above-mentioned method for optimizing the time resolution of the Beidou-assisted InSAR surface deformation series.

[0020] The fifth aspect of the present application provides a computer program product, including a computer program, which is executed to implement the above-mentioned method for optimizing the time resolution of the Beidou-assisted InSAR surface deformation series.

[0021] Therefore, the embodiments of the present application have the following beneficial effects:

[0022] The embodiments of the present application can extract the one-dimensional surface deformation information of the target monitoring area based on a preset SBAS-InSAR strategy; deploy BeiDou / GNSS receivers in the target monitoring area to obtain BeiDou three-dimensional point deformation information of the target monitoring area, and convert the BeiDou three-dimensional point deformation information into LOS one-dimensional data, and construct a BeiDou monitoring data set and an InSAR monitoring data set of the target detection area based on the LOS one-dimensional data and the one-dimensional surface deformation information; based on the BeiDou monitoring data set and the InSAR monitoring data set, construct a deformation law fusion model corresponding to BeiDou and InSAR at the same time, so as to obtain the InSAR surface deformation of the target monitoring area based on the deformation law fusion model and the BeiDou monitoring data set. The present application uses the respective advantages of InSAR and BeiDou / GNSS to integrate and invert the three-dimensional surface deformation of areas such as large landslides, thereby obtaining three-dimensional deformation with high spatiotemporal resolution to better assist in the monitoring and early warning of geological disasters. Thus, the problems of high monitoring cost, long station deployment cycle, and low spatiotemporal resolution of existing surface deformation monitoring technologies are solved.

[0023] Additional aspects and advantages of the present application will be given in part in the description below, and in part will become apparent from the description below, or will be learned through practice of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] The above and / or additional aspects and advantages of the present application will become apparent and easily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0025] Figure 1 A flowchart of a method for optimizing the temporal resolution of a BeiDou-assisted InSAR surface deformation series according to an embodiment of the present application is provided;

[0026] Figure 2 A schematic diagram of an SBAS-InSAR processing flow provided for one embodiment of the present application;

[0027] Figure 3 A schematic diagram of a BeiDou data processing flow provided for one embodiment of the present application;

[0028] Figure 4 A schematic diagram of a deformation law fusion model construction process provided in one embodiment of the present application;

[0029] Figure 5 This is an example diagram of a device for optimizing the time resolution of a BeiDou-assisted InSAR surface deformation series according to an embodiment of the present application;

[0030] Figure 6 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application.

[0031] Among them, 10 is a device for optimizing the time resolution of Beidou-assisted InSAR surface deformation series; 100 is an extraction module, 200 is a conversion module, 300 is a fusion module; 601 is a memory, 602 is a processor, and 603 is a communication interface. DETAILED DESCRIPTION

[0032] The following describes in detail embodiments of the present application, examples of which are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present application, and should not be construed as limiting the present application.

[0033] The following describes a method for optimizing the temporal resolution of a BeiDou-assisted InSAR surface deformation sequence according to an embodiment of the present application with reference to the accompanying drawings. In response to the problems mentioned in the above background technology, the present application provides a method for optimizing the temporal resolution of a BeiDou-assisted InSAR surface deformation sequence, wherein one-dimensional surface deformation information of a target monitoring area is extracted based on a preset SBAS-InSAR strategy; BeiDou / GNSS receivers are deployed in the target monitoring area to obtain BeiDou three-dimensional point deformation information of the target monitoring area, and the BeiDou three-dimensional point deformation information is converted into one-dimensional data in the LOS (Line Of Sight) direction, and a BeiDou monitoring dataset and an InSAR monitoring dataset of the target detection area are constructed based on the one-dimensional LOS direction data and the one-dimensional surface deformation information; based on the BeiDou monitoring dataset and the InSAR monitoring dataset, a deformation law fusion model corresponding to the BeiDou and InSAR at the same time is constructed, so as to obtain the InSAR surface deformation amount of the target monitoring area based on the deformation law fusion model and the BeiDou monitoring dataset. This application leverages the respective strengths of InSAR and BeiDou / GNSS to invert 3D surface deformation in areas such as large landslides, generating high-resolution 3D deformation data to better aid in geological disaster monitoring and early warning. This addresses the challenges of existing surface deformation monitoring technologies, such as high monitoring costs, long deployment cycles, and low spatial and temporal resolution.

[0034] Specifically, Figure 1 This is a flowchart of a method for optimizing the temporal resolution of a BeiDou-assisted InSAR surface deformation series provided in an embodiment of the present application.

[0035] like Figure 1 As shown in FIG, the method for optimizing the temporal resolution of the BeiDou-assisted InSAR surface deformation series includes the following steps:

[0036] In step S101 , one-dimensional surface deformation information of the target monitoring area is extracted based on a preset SBAS-InSAR strategy.

[0037] It should be noted that, based on the surface features of natural scenes, the embodiment of the present application may use SBAS-InSAR (Small Baseline Subset InSAR, differential interferometry short baseline set temporal analysis) technology to obtain temporal deformation results. The radar interferometry phase model of the embodiment of the present application is as follows:

[0038]

[0039] It can be seen from formula (1) that the radar interference phase in the embodiment of the present application is Depend on is the deformation phase, For terrain phase, is the atmospheric phase, Orbital phase, Error phase component.

[0040] In addition, if Figure 2 As shown, the embodiment of the present application sets appropriate time baseline and space baseline thresholds to perform multi-main image interference pair combination, fully utilizes the advantages of interference subset technology, and selects points based on the backscattering coefficient, thereby effectively supplementing the shortcomings of the coherence point selection method in the low coherence area.

[0041] After selecting the target point, the target point is unwrapped based on the minimum cost flow algorithm of the irregular grid; then, the linear deformation, DEM (Digital Elevation Model) error, and residual phase are separated through the regression model. Multiple iterations should be performed according to the actual situation, and the two-dimensional regression model should be continuously updated to minimize the difference between the two-dimensional linear phase model and the differential phase. Therefore, the embodiment of the present application can adopt a strategy of step-by-step correction of the model estimation parameters. Through multiple iterative regression analysis, the elevation error, deformation rate, baseline error, and atmospheric phase are improved, thereby obtaining terrain elevation error and line-of-sight deformation sequence information (i.e., one-dimensional surface deformation information). .

[0042] The embodiment of the present application extracts the one-dimensional surface deformation (i.e., one-dimensional surface deformation information) of the target monitoring area based on SBAS-InSAR technology, thereby providing reliable data support for improving the temporal resolution of subsequent Beidou-assisted InSAR surface deformation sequences.

[0043] In step S102, Beidou / GNSS receivers are deployed in the target monitoring area to obtain Beidou three-dimensional point deformation information of the target monitoring area, and the Beidou three-dimensional point deformation information is converted into LOS one-dimensional data, and the Beidou monitoring dataset and InSAR monitoring dataset of the target detection area are constructed based on the LOS one-dimensional data and the one-dimensional surface deformation information.

[0044] Furthermore, the embodiments of the present application deploy Beidou / GNSS receivers in the target monitoring area, extract Beidou three-dimensional point deformation information of Beidou / GNSS monitoring points, and convert the Beidou three-dimensional deformation sequence (i.e., Beidou three-dimensional point deformation information) into InSAR radar line of sight one-dimensional deformation (i.e., LOS one-dimensional data), thereby constructing Beidou and InSAR high-quality monitoring data sets.

[0045] It should be noted that the embodiments of the present application can use Beidou / GNSS double-difference carrier phase observations to achieve high-precision positioning and obtain time series deformation results. The double-difference observation equation is as follows:

[0046]

[0047] in, 、 are pseudorange and carrier phase observations respectively; 、 are the receiver clock error and satellite clock error respectively; 、 are the tropospheric delay error and the ionospheric delay error respectively; is the wavelength of the carrier phase; is the whole-week ambiguity; 、 are the time when the station receiver receives the satellite signal and the time when the satellite transmits the signal; 、 are the multipath effects corresponding to pseudorange and carrier phase respectively; 、 The measurement accuracy of pseudorange and carrier phase respectively; is the constant speed of light; It is the spatial geometric distance between the satellite and the ground station receiver.

[0048] Specifically, the process of obtaining Beidou three-dimensional point deformation information of the deformation monitoring area (i.e., the target monitoring area) in the embodiment of the present application is as follows:

[0049] 1. Satellite observation acquisition

[0050] Deploy BeiDou / GNSS receivers in areas at risk of large-scale landslide disasters to obtain BeiDou / GNSS satellite signals;

[0051] 2. Preparation of observation values ​​and ephemeris files

[0052] The receiver observation data is uploaded to the receiving platform via the 4G network. The solution platform obtains the original observation data from the receiving platform via the FTP protocol and downloads the ephemeris data from various domestic and international GNSS research institutions.

[0053] 3. Unified data format

[0054] By decoding and converting the BeiDou / GNSS satellite signals, receiver observation data and ephemeris data, the data formats are unified to obtain standard satellite signals, standard observation data and standard ephemeris data.

[0055] 4. Data quality analysis

[0056] Data quality analysis software such as teqc is used to perform statistical analysis on the original observation data in terms of effective time, data integrity rate, multipath error, etc.

[0057] 5. Baseline solution

[0058] Use professional baseline calculation software to calculate the baseline of the observation data in the relevant area and obtain the relative position information of the monitoring point relative to the benchmark point;

[0059] 6. Network Adjustment

[0060] The results of baseline solution are imported into the universal adjustment software package COSAGPS to perform three-dimensional vector network adjustment and obtain the three-dimensional position information of each monitoring station;

[0061] 7. Time Series Analysis

[0062] By analyzing the deformation time series of each monitoring station after adjustment, the three-dimensional point deformation information of the deformation monitoring area is obtained. .

[0063] The following describes the process of extracting the three-dimensional deformation of monitoring points based on Beidou / GNSS according to an embodiment of the present application through a specific embodiment and in combination with the accompanying drawings.

[0064] Figure 3 The figure is a schematic diagram of the extraction process of three-dimensional deformation of monitoring points based on BeiDou / GNSS. Figure 3 As shown, in the embodiment of the present application, the process of extracting the three-dimensional deformation of the monitoring point based on Beidou / GNSS is as follows:

[0065] 1. Beidou observation data collation:

[0066] (1) Obtaining original observation data from BeiDou / GNSS stations in large landslide areas;

[0067] (2) Obtain the synchronous observation time, multipath error, signal-to-noise ratio, and satellite visibility corresponding to the original observation data;

[0068] (3) Effectively screen and eliminate gross errors from the original observation data based on the synchronous observation time, multipath error, signal-to-noise ratio, and satellite visibility;

[0069] (4) Obtaining “clean” BeiDou observation data from all BeiDou / GNSS stations after data screening and gross error removal;

[0070] 2. Establishment of a unified stability benchmark:

[0071] (1) Obtaining “clean” BeiDou observation data from IGS (International GNSS Service) stations and reference stations;

[0072] (2) Obtain the three-dimensional coordinates of the reference station corresponding to the Beidou observation data of the IGS station and the Beidou observation data of the reference station;

[0073] (3) Based on the InSAR geographic coordinate system, a unified stability benchmark for large landslide areas corresponding to the three-dimensional coordinates of the base station was established;

[0074] 3. 3D deformation information extraction:

[0075] (1) Obtaining “clean” BeiDou observation data from the base station and monitoring points;

[0076] (2) Perform high-precision BeiDou data processing on the “clean” BeiDou observation data from the base stations and monitoring points;

[0077] (3) Obtaining the time series of monitoring points corresponding to the “clean” BeiDou observation data of the base station and monitoring points after high-precision BeiDou data processing;

[0078] (4) Perform time series filtering and data extraction operations such as periodic items and trend items on the monitoring point time series;

[0079] (5) Based on the data of periodic items, trend items and the time series of monitoring points after time series filtering, the three-dimensional surface point deformation information of large landslide areas (i.e., the three-dimensional deformation of monitoring points) is extracted.

[0080] Therefore, the embodiments of the present application effectively ensure that InSAR and BeiDou / GNSS complement each other's strengths by extracting the three-dimensional deformation of monitoring points based on BeiDou / GNSS.

[0081] Optionally, in one embodiment of the present application, Beidou three-dimensional point deformation information is converted into LOS one-dimensional data, and a Beidou monitoring data set and an InSAR monitoring data set of the target detection area are constructed based on the LOS one-dimensional data and the one-dimensional surface deformation information, including: obtaining the azimuth and altitude of the InSAR remote sensing satellite image data corresponding to the target detection area, and based on the azimuth and altitude, converting the Beidou three-dimensional point deformation information into LOS one-dimensional data; determining the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference corresponding to the LOS one-dimensional data and the one-dimensional surface deformation information respectively; performing a spatial reference unification operation on the LOS one-dimensional data and the one-dimensional surface deformation information based on the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference to obtain the LOS one-dimensional data and the one-dimensional surface deformation information under a unified coordinate reference, and establishing a Beidou monitoring data set and an InSAR monitoring data set based on the LOS one-dimensional data and the one-dimensional surface deformation information under the unified coordinate reference.

[0082] After that, the embodiment of the present application also needs to obtain the azimuth of the InSAR remote sensing satellite image data and altitude angle ; Secondly, obtain three-dimensional data of Beidou monitoring points ; Finally, the three-dimensional data of the Beidou monitoring points in the large landslide area are converted into one-dimensional data in the LOS direction.

[0083] InSAR radar line of sight (LOS) observation vector and three-dimensional deformation components The functional relationship between them can be expressed as

[0084]

[0085] The BeiDou and InSAR deformations of the monitoring area obtained above are d BDS and d InSAR The coordinate reference for Beidou deformation monitoring is the CGCS2000 coordinate system; the coordinate reference for InSAR monitoring data is determined by the reference point coordinate system, which can be CGCS2000, WGS84, or ITRF (International Terrestrial Reference Frame). Inconsistency in the spatial reference of monitoring data will affect the performance of subsequent precise fusion of deformation patterns.

[0086] Therefore, the embodiment of the present application needs to unify the spatial references of Beidou and InSAR monitoring data through coordinate frame conversion, so that a deformation monitoring dataset containing Beidou and InSAR deformation values ​​under a unified coordinate reference can be constructed, namely:

[0087]

[0088] in, P BDS and( B BDS , L BDS ) represents the Beidou monitoring point name and longitude and latitude coordinates; P InSAR and( B InSAR , L InSAR ) represents the InSAR deformation pixel point name and longitude and latitude position; d BDS and d InSAR Represent the BeiDou and InSAR deformations respectively, t BDS and t InSAR They represent the corresponding moments of BeiDou and InSAR deformation respectively.

[0089] In step S103, based on the Beidou monitoring dataset and the InSAR monitoring dataset, a deformation law fusion model corresponding to the Beidou and InSAR at the same time is constructed to obtain the InSAR surface deformation of the target monitoring area according to the deformation law fusion model and the Beidou monitoring dataset.

[0090] Furthermore, the embodiment of the present application also needs to construct a deformation law fusion model driven by Beidou and InSAR monitoring data to obtain high-time-resolution InSAR surface deformation through the deformation law fusion model.

[0091] Optionally, in one embodiment of the present application, based on the BeiDou monitoring data set and the InSAR monitoring data set, a deformation law fusion model corresponding to the BeiDou and InSAR at the same time is constructed to obtain the InSAR surface deformation of the target monitoring area according to the deformation law fusion model and the BeiDou monitoring data set, including: obtaining the deformation law fusion model of the target monitoring area; m Beidou deformation monitoring points M Monitor deformation over time, and n InSAR monitoring pixels N Monitor the deformation periodically and determine m and n The first quantitative relationship between N and M The second quantitative relationship between n 、 m 、 N 、 M are all positive integers, andn > m , N < M Based on the first quantitative relationship, the second quantitative relationship, the BeiDou monitoring dataset and the InSAR monitoring dataset, a monitoring training dataset is constructed; the pre-constructed neural network model is trained with the monitoring training dataset to generate a deformation law fusion model; the deformation law fusion model is interpolated n InSAR monitoring pixels ( M - N ) period of monitoring deformation, and combined n InSAR monitoring pixels ( M - N ) period monitoring deformation and N The deformation is monitored regularly to generate the InSAR surface deformation of the target monitoring area.

[0092] It should be noted that based on the high-temporal-resolution Beidou monitoring dataset and the high-spatial-resolution InSAR monitoring dataset obtained above, a neural network is used to construct a conversion model between the Beidou monitoring multi-point deformation law and the InSAR monitoring surface deformation law at the same time.

[0093] As a possible implementation method, the embodiment of the present application may assume that the monitoring area includes m Beidou deformation monitoring points M Periodic monitoring of deformation and n InSAR deformed pixels N Monitor deformation regularly, and n > m , N < M , n 、 m 、 N 、 M If both are positive integers, the following monitoring training data set can be constructed:

[0094]

[0095] Afterwards, the embodiment of the present application can divide the monitoring training data set into a training set and a validation set according to the proportion, and then perform deep neural network model training, set parameters such as the number of layers, number of nodes, learning rate, loss function, etc., and use the training set data to train the model, and finally obtain a trained deep neural network model that integrates the deformation laws of Beidou and InSAR (i.e., a deformation law fusion model).

[0096] Optionally, in one embodiment of the present application, the mathematical expression of the deformation law fusion model is:

[0097]

[0098] in, f DNN Represents a pre-built neural network model; w Represents the weight parameter of the deformation law fusion model; b Represents the bias parameter of the deformation law fusion model; d BDS,m×1 Indicates LOS to one-dimensional data; d InSAR,n×1 Represents one-dimensional surface deformation information.

[0099] In the embodiment of the present application, the deformation law fusion model constructed above is shown as follows:

[0100] (7)

[0101] in, f DNN Represents a pre-built neural network model; w Represents the weight parameter of the deformation law fusion model; b Represents the bias parameter of the deformation law fusion model; d BDS,m×1 Indicates LOS to one-dimensional data; d InSAR,n×1 Represents one-dimensional surface deformation information.

[0102] In addition, the embodiment of the present application may use the mean value MEAN, the root mean square error (RMSE) and the determination coefficient R 2 The performance of the deformation law fusion model is evaluated by indicators such as , to verify the accuracy and reliability of the deformation law fusion model. The calculation formula is as follows:

[0103]

[0104] Therefore, based on the high temporal resolution BeiDou point deformation dataset obtained above, and combined with the BeiDou and InSAR deformation law fusion model, as shown in the following example: Figure 4 As shown in the figure, the InSAR surface deformation with high temporal resolution can be obtained to improve the problem of low temporal resolution of InSAR surface deformation monitoring.

[0105] In the embodiments of this application, it is assumed that the monitoring area includes m Beidou deformation monitoring points M Periodic monitoring of deformation and n InSAR monitoring pixels N Monitor deformation regularly, and n > m , N < M, then the embodiment of the present application can interpolate the deformation law fusion model established n InSAR monitoring pixels remaining ( MN ) period, that is:

[0106]

[0107] Afterwards, the embodiment of the present application can interpolate the n InSAR monitoring pixels remaining ( MN ) period and the deformation obtained by actual monitoring N The deformation combination is obtained as M A long-term high-temporal-resolution InSAR surface deformation monitoring sequence (i.e., InSAR surface deformation amount).

[0108] It is understood that in the embodiments of the present application, the fusion of InSAR data and GNSS data is mainly divided into the following three aspects:

[0109] (1) Using GNSS data to correct the atmospheric errors of InSAR data, thereby improving the accuracy of InSAR data;

[0110] (2) Using GNSS 3D deformation data to expand InSAR 3D deformation monitoring capabilities to obtain high spatial resolution 3D deformation data;

[0111] (3) Fuse high spatial resolution InSAR data and high temporal resolution GNSS data to generate high spatiotemporal resolution data.

[0112] Therefore, the embodiments of the present application can integrate and invert the three-dimensional surface deformation of areas such as large landslides through the complementary advantages of InSAR and Beidou / GNSS, and obtain three-dimensional deformation with high temporal and spatial resolution, thereby better assisting in the monitoring and early warning of geological disasters.

[0113] According to the method for optimizing the temporal resolution of BeiDou-assisted InSAR surface deformation series proposed in the embodiment of the present application, one-dimensional surface deformation information of the target monitoring area is extracted based on a preset SBAS-InSAR strategy; BeiDou / GNSS receivers are deployed in the target monitoring area to obtain BeiDou three-dimensional point deformation information of the target monitoring area, and the BeiDou three-dimensional point deformation information is converted into LOS-direction one-dimensional data, and a BeiDou monitoring data set and an InSAR monitoring data set of the target detection area are constructed based on the LOS-direction one-dimensional data and the one-dimensional surface deformation information; based on the BeiDou monitoring data set and the InSAR monitoring data set, a deformation law fusion model corresponding to BeiDou and InSAR at the same time is constructed to obtain the InSAR surface deformation of the target monitoring area based on the deformation law fusion model and the BeiDou monitoring data set. This application uses the respective advantages of InSAR and BeiDou / GNSS to integrate and invert the three-dimensional surface deformation of areas such as large landslides, thereby obtaining three-dimensional deformation with high temporal and spatial resolution to better assist in the monitoring and early warning of geological disasters.

[0114] Next, a device for optimizing the temporal resolution of a BeiDou-assisted InSAR surface deformation sequence proposed in an embodiment of the present application will be described with reference to the accompanying drawings.

[0115] Figure 5 It is a block diagram of a device for optimizing the time resolution of Beidou-assisted InSAR surface deformation series according to an embodiment of the present application.

[0116] like Figure 5 As shown, the device 10 for optimizing the time resolution of BeiDou-assisted InSAR surface deformation series includes: an extraction module 100 , a conversion module 200 and a fusion module 300 .

[0117] The extraction module 100 is used to extract one-dimensional surface deformation information of the target monitoring area based on a preset SBAS-InSAR strategy.

[0118] The conversion module 200 is used to deploy Beidou / GNSS receivers in the target monitoring area to obtain Beidou three-dimensional point deformation information of the target monitoring area, convert the Beidou three-dimensional point deformation information into LOS one-dimensional data, and construct a Beidou monitoring data set and an InSAR monitoring data set of the target detection area based on the LOS one-dimensional data and the one-dimensional surface deformation information.

[0119] The fusion module 300 is used to construct a deformation law fusion model corresponding to Beidou and InSAR at the same time based on the Beidou monitoring data set and the InSAR monitoring data set, so as to obtain the InSAR surface deformation of the target monitoring area according to the deformation law fusion model and the Beidou monitoring data set.

[0120] Optionally, in one embodiment of the present application, the conversion module 200 includes: a first acquisition unit, a determination unit, and a spatial reference unification unit.

[0121] Among them, the first acquisition unit is used to obtain the azimuth and altitude of the InSAR remote sensing satellite image data corresponding to the target detection area, and based on the azimuth and altitude, convert the Beidou three-dimensional point deformation information into LOS one-dimensional data.

[0122] The determination unit is used to respectively determine the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference corresponding to the LOS one-dimensional data and the one-dimensional planar deformation information.

[0123] The spatial reference unification unit is used to perform spatial reference unification operations on the LOS one-dimensional data and one-dimensional surface deformation information based on the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference, so as to obtain the LOS one-dimensional data and one-dimensional surface deformation information under the unified coordinate reference, and to establish the Beidou monitoring data set and the InSAR monitoring data set based on the LOS one-dimensional data and one-dimensional surface deformation information under the unified coordinate reference.

[0124] Optionally, in one embodiment of the present application, the fusion module 300 includes: a second acquisition unit, a construction unit, a training unit and a combination unit.

[0125] The second acquisition unit is used to obtain the target detection area m Beidou deformation monitoring points M Monitor deformation over time, and n InSAR monitoring pixels N Monitor the deformation periodically and determine m and n The first quantitative relationship between N and M The second quantitative relationship between n 、 m 、 N 、 M are all positive integers, and n > m , N < M .

[0126] The construction unit is used to construct a monitoring training data set based on the first quantitative relationship, the second quantitative relationship, the Beidou monitoring data set and the InSAR monitoring data set.

[0127] The training unit is used to train a pre-built neural network model by monitoring the training data set to generate a deformation law fusion model.

[0128] Combination unit, used to interpolate the fusion model according to the deformation law n InSAR monitoring pixels ( M - N ) period of monitoring deformation, and combined n InSAR monitoring pixels ( M - N ) period monitoring deformation and N The deformation is monitored regularly to generate the InSAR surface deformation of the target monitoring area.

[0129] Optionally, in one embodiment of the present application, the mathematical expression of the deformation law fusion model is:

[0130]

[0131] in, f DNN Represents a pre-built neural network model; w Represents the weight parameter of the deformation law fusion model; b Represents the bias parameter of the deformation law fusion model; d BDS,m×1 Indicates LOS to one-dimensional data; d InSAR,n×1 Represents one-dimensional surface deformation information.

[0132] It should be noted that the above explanation of the embodiment of the method for optimizing the time resolution of the Beidou-assisted InSAR surface deformation series is also applicable to the device for optimizing the time resolution of the Beidou-assisted InSAR surface deformation series in this embodiment, and will not be repeated here.

[0133] According to the embodiment of the present application, the device for optimizing the temporal resolution of the BeiDou-assisted InSAR surface deformation sequence proposed includes an extraction module 100 for extracting one-dimensional surface deformation information of the target monitoring area based on a preset SBAS-InSAR strategy; a conversion module 200 for deploying BeiDou / GNSS receivers in the target monitoring area to obtain BeiDou three-dimensional point deformation information of the target monitoring area, and converting the BeiDou three-dimensional point deformation information into LOS-direction one-dimensional data, and constructing a BeiDou monitoring data set and an InSAR monitoring data set of the target detection area based on the LOS-direction one-dimensional data and the one-dimensional surface deformation information; a fusion module 300 for constructing a deformation law fusion model corresponding to BeiDou and InSAR at the same time based on the BeiDou monitoring data set and the InSAR monitoring data set, so as to obtain the InSAR surface deformation of the target monitoring area based on the deformation law fusion model and the BeiDou monitoring data set. The present application utilizes the respective advantages of InSAR and BeiDou / GNSS to fuse and invert three-dimensional surface deformation of areas such as large landslides, thereby obtaining three-dimensional deformation with high temporal and spatial resolution to better assist in the monitoring and early warning of geological disasters.

[0134] Figure 6 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. The electronic device may include:

[0135] A memory 601 , a processor 602 , and a computer program stored in the memory 601 and executable on the processor 602 .

[0136] When the processor 602 executes the program, the method for optimizing the temporal resolution of the BeiDou-assisted InSAR surface deformation series provided in the above embodiment is implemented.

[0137] Furthermore, the electronic device further includes:

[0138] The communication interface 603 is used for communication between the memory 601 and the processor 602 .

[0139] The memory 601 is used to store computer programs that can be run on the processor 602 .

[0140] The memory 601 may include a high-speed RAM memory, and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.

[0141] If the memory 601, processor 602, and communication interface 603 are implemented independently, the communication interface 603, memory 601, and processor 602 can be connected to each other via a bus and communicate with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus. Buses can be divided into address buses, data buses, control buses, etc. For ease of representation, Figure 6 Only one thick line is used in the diagram, but this does not mean that there is only one bus or one type of bus.

[0142] Optionally, in a specific implementation, if the memory 601, the processor 602 and the communication interface 603 are integrated on a chip, the memory 601, the processor 602 and the communication interface 603 can communicate with each other through an internal interface.

[0143] The processor 602 may be a central processing unit (CPU), an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application.

[0144] An embodiment of the present application also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-mentioned method for optimizing the time resolution of the Beidou-assisted InSAR surface deformation series.

[0145] An embodiment of the present application further provides a computer program product, including a computer program, which, when executed, is used to implement the above-mentioned method for optimizing the time resolution of the Beidou-assisted InSAR surface deformation series.

[0146] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or N embodiments or examples in a suitable manner. In addition, those skilled in the art can combine and combine different embodiments or examples described in this specification and features of different embodiments or examples without contradiction.

[0147] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be understood to indicate or imply relative importance or implicitly specify the number of technical features indicated. Thus, a feature specified as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of this application, "N" means at least two, for example, two, three, etc., unless otherwise specifically defined.

[0148] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or N executable instructions for implementing a custom logical function or process step, and the scope of the preferred embodiments of the present application includes alternative implementations in which functions may be performed in a different order than shown or discussed, including performing functions in a substantially simultaneous manner or in a reverse order depending on the functions involved, which should be understood by those skilled in the art to which the embodiments of the present application pertain.

[0149] The logic and / or steps represented in a flowchart or otherwise described herein, for example, can be considered a sequenced list of executable instructions for implementing the logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (e.g., a computer-based system, a system including a processor, or other system that can fetch and execute instructions from an instruction execution system, apparatus, or device). For purposes of this specification, a "computer-readable medium" is any device that can contain, store, communicate, propagate, or transport a program for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (not exhaustive) of computer-readable media include: an electrical connection with one or more wires (electronic devices), a portable computer disk cartridge (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and programmable read-only memory (EPROM or flash memory), fiber optic devices, and a portable compact disc read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program can be obtained electronically by optically scanning the paper or other medium and then editing, interpreting or processing it in other suitable ways as necessary, and then storing it in a computer memory.

[0150] It should be understood that various parts of this application can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiment, the N steps or methods can be implemented using software or firmware stored in a memory and executed by a suitable instruction execution system. If implemented using hardware, as in another embodiment, any of the following technologies known in the art or a combination thereof can be used: a discrete logic circuit having logic gate circuits for implementing logical functions on data signals, an application-specific integrated circuit having suitable combinational logic gate circuits, a programmable gate array (PGA), a field-programmable gate array (FPGA), etc.

[0151] Those skilled in the art will understand that all or part of the steps in the method of the above embodiment can be completed by instructing related hardware through a program, and the program can be stored in a computer-readable storage medium. When the program is executed, it includes one or a combination of the steps of the method embodiment.

[0152] In addition, the functional units in the various embodiments of the present application may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into a module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0153] The storage medium mentioned above may be a read-only memory, a magnetic disk, or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present application. Persons skilled in the art may make changes, modifications, substitutions, and variations to the above embodiments within the scope of the present application.

Claims

1. A method for optimizing the temporal resolution of BeiDou-assisted InSAR surface deformation series, characterized in that: The following steps are involved: Based on the preset SBAS-InSAR strategy, the one-dimensional surface deformation information of the target monitoring area is extracted; Deploy a BeiDou / GNSS receiver in the target monitoring area to obtain BeiDou three-dimensional point deformation information of the target monitoring area, convert the BeiDou three-dimensional point deformation information into one-dimensional data in the LOS direction, and construct a BeiDou monitoring dataset and an InSAR monitoring dataset of the target monitoring area based on the one-dimensional data in the LOS direction and the one-dimensional surface deformation information; Based on the Beidou monitoring dataset and the InSAR monitoring dataset, a deformation law fusion model corresponding to the Beidou and InSAR at the same time is constructed, so as to obtain the InSAR surface deformation of the target monitoring area according to the deformation law fusion model and the Beidou monitoring dataset; The method includes: constructing a deformation law fusion model corresponding to Beidou and InSAR at the same time based on the Beidou monitoring dataset and the InSAR monitoring dataset, so as to obtain the InSAR surface deformation of the target monitoring area according to the deformation law fusion model and the Beidou monitoring dataset, including: Get the target monitoring area m Beidou deformation monitoring points M Monitor deformation over time, and n InSAR monitoring pixels N Monitor the deformation periodically and determine m and n The first quantitative relationship between N and M The second quantitative relationship between n 、 m 、 N 、 M are all positive integers, and n > m , N < M ; Constructing a monitoring training dataset based on the first quantitative relationship, the second quantitative relationship, the Beidou monitoring dataset, and the InSAR monitoring dataset; Training a pre-built neural network model using the monitoring training data set to generate the deformation law fusion model; According to the deformation law, the fusion model is interpolated n InSAR monitoring pixels ( M - N ) period of monitoring deformation, and combining the n InSAR monitoring pixels ( M - N ) period of monitoring deformation and the N The deformation is monitored regularly to generate the InSAR surface deformation of the target monitoring area.

2. The method for optimizing the temporal resolution of BeiDou-assisted InSAR surface deformation series according to claim 1, characterized in that: The converting of the Beidou three-dimensional point deformation information into one-dimensional data in the LOS direction, and constructing a Beidou monitoring data set and an InSAR monitoring data set of the target monitoring area according to the one-dimensional data in the LOS direction and the one-dimensional surface deformation information, comprises: Acquire the azimuth and elevation of the InSAR remote sensing satellite image data corresponding to the target monitoring area, and convert the Beidou three-dimensional point deformation information into the LOS direction one-dimensional data based on the azimuth and the elevation; Determine the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference corresponding to the LOS one-dimensional data and the one-dimensional planar deformation information respectively; Based on the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference, a spatial reference unification operation is performed on the LOS one-dimensional data and the one-dimensional surface deformation information to obtain the LOS one-dimensional data and the one-dimensional surface deformation information under a unified coordinate reference, and the Beidou monitoring dataset and the InSAR monitoring dataset are established based on the LOS one-dimensional data and the one-dimensional surface deformation information under the unified coordinate reference.

3. The method for optimizing the temporal resolution of BeiDou-assisted InSAR surface deformation series according to claim 1, characterized in that: The mathematical expression of the deformation law fusion model is: in, f DNN Representing the pre-built neural network model; w Represents the weight parameter of the deformation law fusion model; b represents the bias parameter of the deformation law fusion model; d BDS,m×1 Representing the LOS as one-dimensional data; d InSAR,n×1 Indicates the one-dimensional surface deformation information.

4. A device for optimizing the temporal resolution of BeiDou-assisted InSAR surface deformation series, characterized in that: include: The extraction module is used to extract the one-dimensional surface deformation information of the target monitoring area based on the preset SBAS-InSAR strategy; a conversion module, configured to deploy a BeiDou / GNSS receiver in the target monitoring area to obtain BeiDou three-dimensional point deformation information of the target monitoring area, convert the BeiDou three-dimensional point deformation information into one-dimensional data in the LOS direction, and construct a BeiDou monitoring dataset and an InSAR monitoring dataset of the target monitoring area based on the one-dimensional data in the LOS direction and the one-dimensional surface deformation information; a fusion module, configured to construct a deformation law fusion model corresponding to the Beidou and InSAR data sets at the same time based on the Beidou monitoring data set and the InSAR monitoring data set, so as to obtain the InSAR surface deformation of the target monitoring area according to the deformation law fusion model and the Beidou monitoring data set; The fusion module includes: The second acquisition unit is used to obtain the m Beidou deformation monitoring points M Monitor deformation over time, and n InSAR monitoring pixels N Monitor the deformation periodically and determine m and n The first quantitative relationship between N and M The second quantitative relationship between n 、 m 、 N 、 M are all positive integers, and n > m , N < M ; A construction unit, configured to construct a monitoring training dataset based on the first quantitative relationship, the second quantitative relationship, the Beidou monitoring dataset, and the InSAR monitoring dataset; A training unit, configured to train a pre-built neural network model using the monitoring training data set to generate the deformation law fusion model; Combining unit, used to interpolate the n InSAR monitoring pixels ( M - N ) period of monitoring deformation, and combining the n InSAR monitoring pixels ( M - N ) period of monitoring deformation and the N The deformation is monitored regularly to generate the InSAR surface deformation of the target monitoring area.

5. The device for optimizing the temporal resolution of BeiDou-assisted InSAR surface deformation series according to claim 4, characterized in that: The conversion module includes: A first acquisition unit is configured to acquire the azimuth and elevation of the InSAR remote sensing satellite image data corresponding to the target monitoring area, and convert the Beidou three-dimensional point deformation information into the LOS direction one-dimensional data based on the azimuth and the elevation; A determination unit, configured to respectively determine a Beidou deformation monitoring coordinate reference and an InSAR monitoring data coordinate reference corresponding to the LOS one-dimensional data and the one-dimensional planar deformation information; A spatial reference unification unit is configured to perform a spatial reference unification operation on the one-dimensional LOS data and the one-dimensional surface deformation information based on the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference to obtain the one-dimensional LOS data and the one-dimensional surface deformation information under a unified coordinate reference, and to establish the Beidou monitoring dataset and the InSAR monitoring dataset based on the one-dimensional LOS data and the one-dimensional surface deformation information under the unified coordinate reference.

6. An electronic device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method for optimizing the time resolution of a BeiDou-assisted InSAR surface deformation series according to any one of claims 1 to 3.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that: The program is executed by a processor to implement the method for optimizing the time resolution of the Beidou-assisted InSAR surface deformation series as described in any one of claims 1 to 3.

8. A computer program product comprising a computer program, characterized in that The computer program is executed to implement the method for optimizing the time resolution of the BeiDou-assisted InSAR surface deformation series according to any one of claims 1 to 3.

Images