Beidou-assisted InSAR earth surface deformation sequence time resolution optimization method

Abstract

The invention relates to the technical field of high-precision deformation monitoring, in particular to a Beidou-assisted InSAR (Interferometric Synthetic Aperture Radar) earth surface deformation sequence time resolution optimization method, which comprises the following steps: extracting one-dimensional planar deformation based on SBAS-InSAR; three-dimensional deformation based on the Beidou / GNSS monitoring point is extracted; converting the Beidou three-dimensional deformation sequence into InSAR radar sight line one-dimensional deformation; a Beidou and InSAR high-quality monitoring data set is constructed; constructing a deformation rule fusion model driven by Beidou and InSAR monitoring data; and the InSAR planar deformation with high time resolution is obtained. Therefore, the problems of relatively high monitoring cost, relatively long station distribution period, relatively low temporal-spatial resolution and the like of an existing earth surface deformation monitoring technology are solved.

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] The frequent occurrence of geological disasters in China has caused serious harm to people's lives and property. In order to monitor surface deformation and prevent the occurrence of geological disasters, monitoring technologies such as Interferometric Synthetic Aperture Radar (InSAR) and Global Navigation Satellite System (GNSS) have emerged.

[0003] InSAR technology has the advantages of working around the clock, all-weather and high spatial resolution, but 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 the highest accuracy reaching millimeter level. However, compared with InSAR technology, GNSS ground monitoring equipment is more expensive and has a longer station 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 need to be urgently addressed. Summary of the invention

[0005] The present application provides a method for optimizing the temporal resolution of a BeiDou-assisted InSAR surface deformation sequence, in order to solve the problems of the existing surface deformation monitoring technology, 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 data set and an InSAR monitoring data set 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 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.

[0007] Optionally, in an embodiment of the present application, the conversion of the Beidou three-dimensional point deformation information into one-dimensional LOS data, and the construction of the Beidou monitoring dataset and the InSAR monitoring dataset for the target detection area based on the one-dimensional LOS data and the one-dimensional planar deformation information include: obtaining the azimuth angle and altitude angle of the InSAR remote sensing satellite image data corresponding to the target detection area, and based on the azimuth angle and the altitude angle, converting the Beidou three-dimensional point deformation information into the one-dimensional LOS data; respectively determining the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference corresponding to the one-dimensional LOS data and the one-dimensional planar 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 one-dimensional LOS data and the one-dimensional planar deformation information to obtain the one-dimensional LOS data and the one-dimensional planar deformation information under the unified coordinate reference, and establishing the Beidou monitoring dataset and the InSAR monitoring dataset according to the one-dimensional LOS data and the one-dimensional planar deformation information under the unified coordinate reference.

[0008] Optionally, in an embodiment of the present application, the construction of the deformation law fusion model corresponding to Beidou and InSAR at the same moment based on the Beidou monitoring dataset and the InSAR monitoring dataset to obtain the InSAR planar deformation amount of the target monitoring area according to the deformation law fusion model and the Beidou monitoring dataset includes: obtaining the M-phase monitoring deformation amounts of m Beidou deformation monitoring points in the target detection area, and the N-phase monitoring deformation amounts of n InSAR monitoring pixel points, and respectively determining the first quantitative relationship between m and n and the second quantitative relationship between N and M, where n, m, N, and 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-constructed neural network model through the monitoring training dataset to generate the deformation law fusion model; interpolating the (M-N)-phase monitoring deformation amounts of the n InSAR monitoring pixel points according to the deformation law fusion model, and combining the (M-N)-phase monitoring deformation amounts of the n InSAR monitoring pixel points and the N-phase monitoring deformation amounts to generate the InSAR planar deformation amount of the target monitoring area.

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

[0010] d InSAR,n×1 =f DNN (w, b, d BDS,m×1 )

[0011] Among them, f DNN represents 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 Represents the LOS as one-dimensional data; d InSAR,n×1 Represents 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, used to extract one-dimensional surface deformation information of a target monitoring area based on a preset SBAS-InSAR strategy; a conversion module, used to 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-direction one-dimensional data, and construct a Beidou monitoring data set and an InSAR monitoring data set of the target detection area according to the LOS one-dimensional data and the one-dimensional surface deformation information; a fusion module, 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 amount of the target monitoring area according to 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 the 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 an embodiment of the present application, the fusion module includes: a second acquisition unit configured to acquire the M-phase monitoring deformation amounts of m Beidou deformation monitoring points and the N-phase monitoring deformation amounts of n InSAR monitoring pixel points within the target detection area, and respectively determine a first quantitative relationship between m and n and a second quantitative relationship between N and M, where n, m, N, and 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-constructed neural network model through the monitoring training dataset to generate the deformation law fusion model; and a combination unit configured to interpolate the (M - N)-phase monitoring deformation amounts of the n InSAR monitoring pixel points according to the deformation law fusion model, and combine the (M - N)-phase monitoring deformation amounts of the n InSAR monitoring pixel points and the N-phase monitoring deformation amounts to generate the InSAR areal deformation amount of the target monitoring area.

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

[0016] d InSAR，n×1 = f DNN (w, b, d BDS，m×1 )

[0017] where f DNN represents the pre-constructed 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 represents the LOS-direction one-dimensional data; d InSAR,n×1 represents the one-dimensional areal deformation information.

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

[0019] An embodiment of the fourth aspect of the present application provides a computer-readable storage medium storing a computer program, and when the program is executed by a processor, it implements the method for optimizing the time resolution of the Beidou-assisted InSAR surface deformation sequence as described above.

[0020] An embodiment of the fifth aspect of the present application provides a computer program product including a computer program, and when the computer program is executed, it is used to implement the method for optimizing the time resolution of the Beidou-assisted InSAR surface deformation sequence as described above.

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

[0022] The embodiment of the present application can extract the one-dimensional surface deformation information of the target monitoring area based on the 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 Beidou three-dimensional point deformation information into LOS one-dimensional data, and construct Beidou monitoring data set and InSAR monitoring data set of the target detection area according to LOS one-dimensional data and one-dimensional surface deformation information; based on Beidou monitoring data set and InSAR monitoring data set, construct the deformation law fusion model corresponding to Beidou and InSAR at the same time, so as to obtain the InSAR surface deformation amount of the target monitoring area according to the deformation law fusion model and 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 large landslides and other areas, so as to obtain three-dimensional deformation with high spatiotemporal resolution, so as to better help the monitoring and early warning of geological disasters. Thus, the existing surface deformation monitoring technology solves the problems of high monitoring cost, long station deployment cycle, and low spatiotemporal resolution.

[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 the 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 provided according to an embodiment of the present application;

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

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

[0028] Figure 4 A schematic diagram of a deformation law fusion model construction process provided for an 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 sequence according to an embodiment of the present application;

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

[0031] Among them, 10-Beidou-assisted InSAR surface deformation sequence time resolution optimization device; 100-extraction module, 200-conversion module, 300-fusion module; 601-memory, 602-processor, 603-communication interface. DETAILED DESCRIPTION

[0032] The embodiments of the present application are described in detail below, and examples of the embodiments 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 the optimization method of the time resolution of the Beidou-assisted InSAR surface deformation sequence of the embodiment of the present application with reference to the accompanying drawings. In view of the problems mentioned in the above background technology, the present application provides a method for optimizing the time resolution of the Beidou-assisted InSAR surface deformation sequence, in which 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 one-dimensional data in the LOS (Line Of Sight, line-of-sight propagation) direction, and a Beidou monitoring data set and an InSAR monitoring data set of the target detection area are constructed according to the one-dimensional data in the LOS direction 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 the Beidou and InSAR at the same time is constructed to obtain the InSAR surface deformation amount of the target monitoring area according to 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 large landslides and other areas, thereby obtaining three-dimensional deformation with high temporal and spatial resolution to better assist in the monitoring and early warning of geological disasters. This solves the problems of existing surface deformation monitoring technologies, such as high monitoring costs, long station deployment cycles, and low temporal and spatial resolution.

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

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

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

[0037] It should be noted that, according to the surface features of natural scenes, the embodiment of the present application may adopt SBAS-InSAR (Small Baseline Subset InSAR, differential interferometry short baseline set timing analysis) technology to obtain the timing 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, is the 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 backscattering coefficient points using spectral phase stability, thereby effectively supplementing the shortcomings of coherence point selection in low coherence areas.

[0041] After selecting the target point, the target point is untangled 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 model estimation parameters, and improve the elevation error, deformation rate, baseline error, and atmospheric phase through multiple iterations of regression analysis, so as to obtain terrain elevation error and line-of-sight deformation sequence information (i.e., one-dimensional surface deformation information) d InSAR .

[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 data set and 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.

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

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

[0046]

[0047] Among them, P, are pseudorange and carrier phase observations respectively; dt r 、dt s are the receiver clock error and satellite clock error respectively; T and I are the tropospheric delay error and ionospheric delay error respectively; λ is the wavelength of the carrier phase; N is the integer ambiguity; t r ,t s are the time when the station receiver receives the satellite signal and the time when the satellite transmits the signal; M P , are the multipath effects corresponding to pseudorange and carrier phase respectively; ε P , are the measurement accuracy of pseudorange and carrier phase respectively; c is the constant speed of light; ρ(t s ,t r ) 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. Obtaining satellite observation values

[0050] Deploy Beidou / GNSS receivers in areas with large landslide disaster risks 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 / international GNSS research institutions.

[0053] 3. Unified data format

[0054] The data format is unified by decoding and converting the Beidou / GNSS satellite signals, receiver observation data and ephemeris data to obtain standard satellite signals, standard observation data and standard ephemeris data.

[0055] 4. Data quality analysis

[0056] Use teqc and other data quality analysis software 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 solution software to perform baseline solution on the observation data of the relevant area to 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 general adjustment software package COSAGPS to perform three-dimensional vector network adjustment to 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 (d n d e d u ).

[0063] The following is an explanation of 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 the three-dimensional deformation of the monitoring point based on Beidou / GNSS. Figure 3 As shown, in an embodiment of the present application, the extraction process of the three-dimensional deformation of the monitoring point based on Beidou / GNSS is as follows:

[0065] 1. Beidou observation data collation:

[0066] (1) Obtaining the original observation data of Beidou / GNSS stations in the large landslide area;

[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 raw observation data based on synchronous observation duration, multipath error, signal-to-noise ratio, and satellite visibility;

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

[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 is established;

[0074] 3. Extraction of 3D deformation information:

[0075] (1) Obtain “clean” Beidou observation data from the base station and the monitoring point;

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

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

[0078] (4) Performing 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 the monitoring points) is extracted.

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

[0081] Optionally, in one embodiment of the present application, Beidou three-dimensional point deformation information is converted into one-dimensional data in the LOS direction, and a Beidou monitoring data set and an InSAR monitoring data set of the target detection area are constructed based on the one-dimensional data in the LOS direction 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 one-dimensional data in the LOS direction; respectively determining the Beidou deformation monitoring coordinate reference and the InSAR monitoring data coordinate reference corresponding to the one-dimensional data in the LOS direction 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 one-dimensional data in the LOS direction and the one-dimensional surface deformation information to obtain the one-dimensional data in the LOS direction 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 one-dimensional data in the LOS direction 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 angle α and altitude angle θ of the InSAR remote sensing satellite image data; secondly, obtain the three-dimensional data (d n d e d u ); Finally, the three-dimensional data of Beidou monitoring points in large landslide areas are converted into one-dimensional data in the LOS direction.

[0083] InSAR radar line of sight (LOS) observation vector d BDS and the three-dimensional deformation component (d n d e d u ) 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 of Beidou deformation monitoring is the CGCS2000 coordinate system; the coordinate reference of InSAR monitoring data is determined by the reference point coordinate system, which can be CGCS2000, WGS84 and ITRF (International Terrestrial Reference Frame); the inconsistency of the spatial reference of monitoring data will affect the performance of the subsequent precise fusion of deformation laws.

[0086] Therefore, in the embodiments of the present application, it is necessary to unify the spatial reference of Beidou and InSAR monitoring data through coordinate frame transformation, and then a deformation monitoring data set including Beidou and InSAR deformation amounts under a unified coordinate reference can be constructed, that is:

[0087]

[0088] Among them, P BDS and (B BDS , L BDS ) represent the Beidou monitoring point name and longitude and latitude coordinates; P InSAR and (B InSAR , L InSAR ) represent the InSAR deformation pixel point name and longitude and latitude position; d BDS and d InSAR respectively represent the Beidou and InSAR deformation amounts, and t BDS and t InSAR respectively represent the corresponding moments of the Beidou and InSAR deformation amounts.

[0089] In step S103, 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 moment is constructed to obtain the InSAR areal deformation amount of the target monitoring area according to the deformation law fusion model and the Beidou monitoring data set.

[0090] Furthermore, in the embodiments of the present application, it is also necessary to construct a deformation law fusion model driven by Beidou and InSAR monitoring data to obtain the InSAR areal deformation amount with high time resolution through the deformation law fusion model.

[0091] Optionally, in an 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 Beidou and InSAR at the same moment is constructed to obtain the InSAR areal deformation amount of the target monitoring area according to the deformation law fusion model and the Beidou monitoring data set, including: obtaining the M-phase monitoring deformation amounts of m Beidou deformation monitoring points in the target detection area, and the N-phase monitoring deformation amounts of n InSAR monitoring pixel points, and respectively determining the first quantitative relationship between m and n and the second quantitative relationship between N and M, where n, m, N, and M are all positive integers, and n>m, N<M; constructing 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; training a pre-constructed neural network model through the monitoring training data set to generate a deformation law fusion model; interpolating the (M-N)-phase monitoring deformation amounts of n InSAR monitoring pixel points according to the deformation law fusion model, and combining the (M-N)-phase monitoring deformation amounts of n InSAR monitoring pixel points and the N-phase monitoring deformation amounts to generate the InSAR areal deformation amount of the target monitoring area.

[0092] It should be noted that based on the above-mentioned obtained Beidou monitoring dataset with high time resolution and InSAR monitoring dataset with high spatial resolution, a conversion model for the deformation law of multiple Beidou monitoring points and the deformation law of InSAR monitoring planar shape at the same moment is constructed through a neural network.

[0093] As a feasible implementation method, the embodiment of the present application may assume that the monitoring area contains the deformation amounts of M periods of m Beidou deformation monitoring points and the deformation amounts of N periods of n InSAR deformation pixel points, and n > m, N < M, and n, m, N, and M are all positive integers. Then, the following monitoring training dataset can be constructed:

[0094]

[0095] After that, the embodiment of the present application can divide the monitoring training dataset into a training set and a validation set according to a ratio, and then perform deep neural network model training. Set parameters such as the number of layers, the number of nodes, the learning rate, and the loss function, and use the training set data for model training. Finally, a trained deep neural network model that fuses the deformation laws of Beidou and InSAR (i.e., the deformation law fusion model) is obtained.

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

[0097] d InSAR,n×1 = f DNN (w, b, d BDS,m×1 ) where f DNN represents a pre-constructed 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 represents the LOS-direction one-dimensional data; d InSAR,n×1 represents the one-dimensional planar deformation information.

[0098] In the embodiment of the present application, the above-mentioned constructed deformation law fusion model is shown as the following formula:

[0099] d InSAR,n×1 = f DNN (w, b, d BDS,m×1 ) (7)

[0100] where f DNN represents a pre-constructed 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 represents the LOS-direction one-dimensional data; d InSAR,n×1 represents the one-dimensional planar deformation information.

[0101] In addition, the embodiments of the present application can use indicators such as the mean value MEAN, the root mean square error (RMSE), and the coefficient of determination R 2 to evaluate the performance of the deformation law fusion model, so as to verify the accuracy and reliability of the deformation law fusion model, etc. The calculation formulas are as follows:

[0102]

[0103] Therefore, based on the obtained Beidou point deformation dataset with high temporal resolution, and combined with the Beidou and InSAR deformation law fusion model, as Figure 4 shown, high-temporal-resolution InSAR areal deformation amounts can be obtained, so as to improve the problem of the low temporal resolution of InSAR surface deformation monitoring.

[0104] In the embodiments of the present application, it is assumed that the monitoring area contains the deformation amounts of M periods of m Beidou deformation monitoring points and the deformation amounts of N periods of n InSAR monitoring pixel points, and n > m, N < M. Then, the embodiments of the present application can interpolate the deformation amounts of the remaining (M - N) periods of the n InSAR monitoring pixel points through the established deformation law fusion model, that is:

[0105]

[0106] After that, the embodiments of the present application can combine the interpolated deformation amounts of the remaining (M - N) periods of the n InSAR monitoring pixel points with the actually monitored N-period deformation amounts, so as to obtain an M-period high-temporal-resolution InSAR areal deformation monitoring sequence (i.e., the InSAR areal deformation amount).

[0107] It can be 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:

[0108] (1) Using GNSS data to correct the atmospheric error of InSAR data, so as to improve the accuracy of InSAR data;

[0109] (2) Using GNSS three-dimensional deformation data to expand the three-dimensional deformation monitoring ability of InSAR, so as to obtain three-dimensional deformation data with high spatial resolution;

[0110] (3) Fusing high-spatial-resolution InSAR data and high-temporal-resolution GNSS data to generate high spatio-temporal resolution data.

[0111] Thus, the embodiments of the present application can fuse and invert the three-dimensional surface deformation in areas such as large landslides through the complementary advantages of InSAR and Beidou / GNSS, obtain three-dimensional deformation with high spatio-temporal resolution, and thus better assist in the monitoring and early warning of geological disasters.

[0112] According to the optimization method of the time resolution of the Beidou-assisted InSAR surface deformation sequence proposed in the embodiment of the present application, the one-dimensional surface deformation information of the target monitoring area is extracted based on the preset SBAS-InSAR strategy; Beidou / GNSS receivers are deployed in the target monitoring area to obtain the 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 data set and InSAR monitoring data set of the target detection area are constructed according to 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, a deformation law fusion model corresponding to Beidou and InSAR at the same time is constructed to obtain the InSAR surface deformation amount of the target monitoring area according to 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.

[0113] Secondly, a device for optimizing the time resolution of a Beidou-assisted InSAR surface deformation sequence proposed in an embodiment of the present application is described with reference to the accompanying drawings.

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

[0115] like Figure 5 As shown, the Beidou-assisted InSAR surface deformation sequence time resolution optimization device 10 includes: an extraction module 100, a conversion module 200 and a fusion module 300.

[0116] 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.

[0117] 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, 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.

[0118] 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.

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

[0120] 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.

[0121] 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 surface deformation information.

[0122] The spatial reference unification unit is used to perform spatial reference unification operation 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 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.

[0123] 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.

[0124] The second acquisition unit is used to obtain the M-period monitoring deformation of m Beidou deformation monitoring points and the N-period monitoring deformation of n InSAR monitoring pixel points in the target detection area, and respectively determine the first quantitative relationship between m and n and the second quantitative relationship between N and M, wherein n, m, N, and M are all positive integers, and n>m, N <M。

[0125] 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.

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

[0127] The combination unit is used to interpolate the monitoring deformation of n InSAR monitoring pixels (MN) period according to the deformation law fusion model, and combine the monitoring deformation of n InSAR monitoring pixels (MN) period and the monitoring deformation of N period to generate the InSAR surface deformation of the target monitoring area.

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

[0129] d InSAR,n×1 =fDNN (w, b, d BDS,m×1 )

[0130] Among them, f DNN represents 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 Indicates LOS one-dimensional data; d InSAR,n×1 Represents one-dimensional surface deformation information.

[0131] 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 sequence is also applicable to the device for optimizing the time resolution of the Beidou-assisted InSAR surface deformation sequence of this embodiment, which will not be repeated here.

[0132] According to the optimization device of the time resolution of the Beidou-assisted InSAR surface deformation sequence proposed in the embodiment of the present application, it includes an extraction module 100, which 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 200, which is used to deploy Beidou / GNSS receivers in the target monitoring area to obtain the 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 the Beidou monitoring data set and InSAR monitoring data set of the target detection area according to the LOS one-dimensional data and the one-dimensional surface deformation information; a fusion module 300, which is used to construct the 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 according to 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, so as to better assist in the monitoring and early warning of geological disasters.

[0133] Figure 6 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application. The electronic device may include:

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

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

[0136] Furthermore, the electronic device further comprises:

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

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

[0139] 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.

[0140] If the memory 601, the processor 602 and the communication interface 603 are implemented independently, the communication interface 603, the memory 601 and the processor 602 can be connected to each other through a bus and communicate with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component (PCI) bus or an Extended Industry Standard Architecture (EISA) bus. The bus can be divided into an address bus, a data bus, a control bus, 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 only one type of bus.

[0141] 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.

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

[0143] 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 method for optimizing the time resolution of a Beidou-assisted InSAR surface deformation sequence.

[0144] An embodiment of the present application also 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.

[0145] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. 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 may be combined in any one or N embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

[0146] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Therefore, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the features. In the description of this application, "N" means at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0147] 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 the steps of a custom logical function or process, and the scope of the preferred embodiments of the present application includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present application belong.

[0148] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purpose of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or N wirings (electronic devices), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), a fiber optic device, and a portable compact disk 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 may be obtained electronically by optically scanning the paper or other medium and then editing, interpreting or processing in other suitable ways as necessary and then storing it in a computer memory.

[0149] It should be understood that the various parts of the present application can be implemented by hardware, software, firmware or a combination thereof. In the above embodiment, the N steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. If implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0150] A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

[0151] In addition, each functional unit in each embodiment 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 one 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.

[0152] The storage medium mentioned above may be a read-only memory, a disk or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limiting the present application. A person of ordinary skill in the art may change, modify, replace and modify 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 data set and an InSAR monitoring data set of the target detection area according to the one-dimensional data in the LOS direction 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 according to the deformation law fusion model and the Beidou monitoring data set.

2. The method for optimizing the temporal resolution of the BeiDou-assisted InSAR surface deformation series according to claim 1, characterized in that: The method of converting 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 detection area according to the one-dimensional data in the LOS direction and the one-dimensional surface deformation information, comprises: Acquire the azimuth and altitude of the InSAR remote sensing satellite image data corresponding to the target detection area, and convert the Beidou three-dimensional point deformation information into the LOS one-dimensional data based on the azimuth and the altitude; 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 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 one-dimensional LOS data and the one-dimensional surface deformation information to obtain the one-dimensional LOS data and the one-dimensional surface deformation information under a unified coordinate reference, and the Beidou monitoring data set and the InSAR monitoring data set are established based on the one-dimensional LOS data and the one-dimensional surface deformation information under a unified coordinate reference.

3. The method for optimizing the temporal resolution of the BeiDou-assisted InSAR surface deformation series according to claim 2, characterized in that: The method of 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 according to the deformation law fusion model and the Beidou monitoring data set, includes: Obtain the M-period monitoring deformation of m Beidou deformation monitoring points and the N-period monitoring deformation of n InSAR monitoring pixel points in the target detection area, and determine the first quantitative relationship between m and n and the second quantitative relationship between N and M, respectively, wherein n, m, N, and 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, construct a monitoring training data set; Training a pre-built neural network model using the monitoring training data set to generate the deformation law fusion model; The monitoring deformation of the n InSAR monitoring pixel points (MN) period is interpolated according to the deformation law fusion model, and the monitoring deformation of the n InSAR monitoring pixel points (MN) period and the monitoring deformation of the N period are combined to generate the InSAR surface deformation of the target monitoring area.

4. The method for optimizing the temporal resolution of BeiDou-assisted InSAR surface deformation series according to claim 3, characterized in that: The mathematical expression of the deformation law fusion model is: d InSAR，n×1 =f DNN (w,b,d BDS，m×1 ) Among them, f DNN represents 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 Indicates that the LOS is one-dimensional data; d InSAR,n×1 Represents the one-dimensional surface deformation information.

5. A device for optimizing the temporal resolution of BeiDou-assisted InSAR surface deformation series, characterized in that: include: An extraction module is used to extract one-dimensional surface deformation information of the target monitoring area based on a preset SBAS-InSAR strategy; A conversion module is used to deploy a Beidou / GNSS receiver 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 according to the LOS one-dimensional data and the one-dimensional surface deformation information; A fusion module 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.

6. The device for optimizing the time resolution of BeiDou-assisted InSAR surface deformation series according to claim 5, characterized in that: The conversion module comprises: A first acquisition unit is used to acquire the azimuth and altitude of the InSAR remote sensing satellite image data corresponding to the target detection area, and convert the Beidou three-dimensional point deformation information into the LOS one-dimensional data based on the azimuth and the altitude; A determination unit, used 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; The spatial reference unification unit is used 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 data set and the InSAR monitoring data set according to the one-dimensional LOS data and the one-dimensional surface deformation information under the unified coordinate reference.

7. The device for optimizing the time resolution of BeiDou-assisted InSAR surface deformation series according to claim 6, characterized in that: The fusion module includes: The second acquisition unit is used to obtain the M-period monitoring deformation of m Beidou deformation monitoring points and the N-period monitoring deformation of n InSAR monitoring pixel points in the target detection area, and respectively determine the first quantitative relationship between m and n and the second quantitative relationship between N and M, wherein n, m, N, and M are all positive integers, and n>m, N <M; A construction unit, configured 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, used for training a pre-built neural network model through the monitoring training data set to generate the deformation law fusion model; A combining unit is used to interpolate the monitoring deformation of the n InSAR monitoring pixel points (MN) period according to the deformation law fusion model, and combine the monitoring deformation of the n InSAR monitoring pixel points (MN) period and the monitoring deformation of the N period to generate the InSAR surface deformation of the target monitoring area.

8. 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 sequence as described in any one of claims 1 to 4.

9. 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 a Beidou-assisted InSAR surface deformation sequence as described in any one of claims 1 to 4.

10. 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 as described in any one of claims 1 to 4.

Images