Microscale deformation monitoring system based on interferometric synthetic aperture radar

Abstract

The present invention relates to the technical field of interferometric synthetic aperture radar (ISA), and discloses a micro-scale deformation monitoring system based on ISA, comprising a data acquisition module, a deformation analysis module, and a monitoring and early warning module. The data acquisition module comprises a data acquisition unit, a data preprocessing unit, and a data error early warning unit, and the deformation analysis module comprises a first analysis unit, a second analysis unit, and a comprehensive analysis unit. The system uses ISA to collect deformation variables at several different points of a target location, simultaneously obtains physical environmental factors during the data collection process, and analyzes and obtains an error influence coefficient. Simultaneously, the system calculates an overall deformation coefficient based on the deformation variables at different points of the target location, calculates an overall deformation trend based on the deformation variables at different points of the target location, and comprehensively calculates a comprehensive risk coefficient to determine the micro-scale deformation risk of the target location, reduce interference with the data caused by environmental and physical factors, and ensure the accuracy of the analysis.

Description

Technical Field

[0001] The present invention relates to the technical field of interferometric synthetic aperture radar, and in particular to a micro-scale deformation monitoring system based on interferometric synthetic aperture radar. Background Art

[0002] Interferometric synthetic aperture radar (ISAR) can obtain three-dimensional topographic maps of the Earth's surface with high spatial resolution and elevation accuracy around the clock, in all weather conditions. Using this information, changes in the Earth's land surface and ice and snow surfaces can be monitored, disasters such as earthquakes, volcanic eruptions, landslides, and floods can be predicted, and this can aid production in agriculture, forestry, fisheries, and other sectors, providing information support for various activities. In short, three-dimensional topographic maps obtained using ISAR have important application value and a wide range of applications.

[0003] Although there are a variety of actual data resources available from abroad, their main purpose is to promote the signal processing research of interferometric synthetic aperture radar. At present, my country is in the research and development stage of interferometric synthetic aperture radar systems, and is more concerned about system-level issues, such as: overall parameter selection, system performance evaluation, system error analysis and compensation, etc. A simulation system with high fidelity and high efficiency will provide a powerful tool for the overall design of the system. The flexibility and controllability of the simulation can facilitate comparative tests of different system schemes, provide a basis for correcting the theoretical analysis model of the system, and provide a rich and targeted data source for signal processing research. For interferometric synthetic aperture radar systems, the existing technology has not taken into account the errors caused by subtle meteorological environments, which may lead to inaccurate deformation detection data of dams or bridges. Summary of the Invention

[0004] In view of the shortcomings of the existing technology, the present invention provides a micro-scale deformation monitoring system based on interferometric synthetic aperture radar, which has the advantages of accurate monitoring and solves the above technical problems.

[0005] To achieve the above-mentioned object, the present invention provides the following technical solution: a micro-scale deformation monitoring system based on interferometric synthetic aperture radar, comprising a data acquisition module, a deformation analysis module and a monitoring and early warning module;

[0006] The data acquisition module includes a data acquisition unit, a data preprocessing unit, and a data error warning unit. The data acquisition unit collects deformation variables of several different points of the target location through interferometric synthetic aperture radar. The data preprocessing unit is used to preprocess the data collected by the data acquisition unit and filter out interference. The data error warning unit obtains an error influence coefficient based on a comprehensive analysis of the screening process in the data preprocessing unit and sends the error influence coefficient to the monitoring and warning module.

[0007] The monitoring and early warning module includes an error early warning unit, which determines whether to issue an early warning based on the error influence coefficient;

[0008] The deformation analysis module is data-connected to the data acquisition module. The deformation analysis module includes a first analysis unit, a second analysis unit, and a comprehensive analysis unit. The first analysis unit calculates the overall deformation coefficient based on the deformation variables at different points of the target location. The second analysis unit calculates the overall deformation trend based on the deformation variables at different points of the target location. The comprehensive analysis unit calculates the comprehensive risk coefficient based on the overall deformation coefficient and the overall deformation trend.

[0009] The monitoring and early warning module also includes a comprehensive risk early warning unit and a scheduling unit. The comprehensive risk early warning unit determines whether to issue an early warning based on the comprehensive risk coefficient. The scheduling unit is called when the comprehensive risk unit issues an early warning and executes personnel scheduling.

[0010] As a preferred technical solution of the present invention, the specific expression for the data acquisition unit to collect the deformation variables of several different points at the target location is as follows: ;

[0011] in, Represents a point dataset, Respectively represent the deformation of the first target point collected, , No. The deformation of the target point, , No. The deformation of the target point, .

[0012] As a preferred technical solution of the present invention, the data preprocessing unit is used to preprocess the data collected by the data acquisition unit and filter out interference in the following specific steps:

[0013] Step A1: Obtain atmospheric error similarity coefficient ;

[0014] Step A2: Obtaining the Interferometric Synthetic Aperture Radar Noise Error Coefficient ;

[0015] Step A3: Based on the atmospheric error similarity coefficient in step A1 and the interferometric synthetic aperture radar noise error coefficient in step A2 Construction error influence coefficient , the specific expression is as follows: ;

[0016] in, and They represent weight coefficients that sum to 1, represents the error influence coefficient, represents the atmospheric error similarity coefficient, represents the noise error coefficient of interferometric synthetic aperture radar.

[0017] As a preferred technical solution of the present invention, the step A1 obtains the atmospheric error similarity coefficient The specific steps are as follows:

[0018] Step A1.1: Obtain real-time meteorological monitoring data from ground meteorological stations, including water vapor content and temperature;

[0019] Step A1.2: Calculate the atmospheric error similarity coefficient based on real-time meteorological monitoring data obtained from ground meteorological stations The specific expression is as follows:

[0020] ;

[0021] in, represents the atmospheric error similarity coefficient, It represents the maximum water vapor content that does not affect the detection accuracy of interferometric synthetic aperture radar. express The water vapor content sampled at each moment, represents the water vapor content sampled at time t, express The absolute value of represents the temperature sampled at time t, Indicates the temperature sampled at time t-1.

[0022] As a preferred technical solution of the present invention, the step A2 obtains the noise error coefficient of the interferometric synthetic aperture radar The specific steps are as follows:

[0023] Step A2.1: Obtain the electrical interference coefficient during the interferometric synthetic aperture radar imaging process, including the signal-to-noise ratio;

[0024] Step A2.2: Calculate the noise error coefficient of the interferometric synthetic aperture radar based on the signal-to-noise ratio obtained in step A2.1 The specific expression is as follows: ;

[0025] in, represents the signal-to-noise ratio, represents the noise error coefficient of interferometric synthetic aperture radar.

[0026] As a preferred technical solution of the present invention, the specific steps of the error warning unit judging whether to issue a warning based on the error influence coefficient are as follows: Exceeding the error impact threshold When the warning is issued, the point data set The deformation variables and error influence coefficients of all points in Perform corresponding storage and terminate subsequent work. When the error influence coefficient Does not exceed the error impact threshold When , the error warning unit is terminated.

[0027] As a preferred technical solution of the present invention, the first analysis unit calculates the overall deformation coefficient based on the deformation variables at different points of the target location. The specific expression is as follows:

[0028] ;

[0029] in, represents the overall deformation coefficient, Represents a point dataset The maximum value of Indicates total The deformation variables of the target points are summed.

[0030] As a preferred technical solution of the present invention, the second analysis unit calculates the overall deformation trend based on the deformation variables at different points of the target location, and the specific expression is as follows:

[0031] ;

[0032] in, Indicates the overall deformation trend, Indicates total The deformation variables of the target points are summed up. Indicates the last The deformation of the target point, Indicates the last sampling The deformation variables of the target points are summed.

[0033] As a preferred technical solution of the present invention, the comprehensive analysis unit calculates the comprehensive risk coefficient based on the overall deformation coefficient and the overall deformation trend, and the specific expression is as follows:

[0034] ;

[0035] in, represents the overall deformation coefficient, Indicates the overall deformation trend, represents the comprehensive risk factor, represents a natural constant, Represents a logarithmic function whose base is a natural constant.

[0036] As a preferred technical solution of the present invention, the specific steps of the comprehensive risk warning unit to determine whether to issue a warning based on the comprehensive risk coefficient are as follows: Exceeding the comprehensive risk threshold When the comprehensive risk factor Does not exceed the comprehensive risk threshold No warning is given when

[0037] The scheduling unit is called when the comprehensive risk unit issues an early warning, and the specific strategy for executing personnel scheduling is to schedule personnel based on distance as the priority.

[0038] Compared with the existing technology, the present invention provides a micro-scale deformation monitoring system based on interferometric synthetic aperture radar, which has the following beneficial effects:

[0039] The present invention collects deformation variables of several different points of a target location through interferometric synthetic aperture radar, and simultaneously obtains physical environmental factors in the data collection process, and performs comprehensive analysis to obtain an error influence coefficient. At the same time, the overall deformation coefficient is calculated based on the deformation variables of different points of the target location, the overall deformation trend is calculated based on the deformation variables of different points of the target location, and a comprehensive risk coefficient is obtained by comprehensive calculation, thereby judging the micro-scale deformation risk of the target point, reducing the interference of environmental and physical factors on the data, and ensuring the accuracy of the analysis. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 Schematic diagram of the system framework of the present invention. DETAILED DESCRIPTION

[0041] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0042] See also Figure 1 , a micro-scale deformation monitoring system based on interferometric synthetic aperture radar, including data acquisition module, deformation analysis module and monitoring and early warning module;

[0043] The data acquisition module includes a data acquisition unit, a data preprocessing unit, and a data error warning unit. The data acquisition unit collects deformation variables of several different points at the target location through interferometric synthetic aperture radar. The data preprocessing unit is used to preprocess the data collected by the data acquisition unit and filter out interference. The data error warning unit obtains the error influence coefficient based on the comprehensive analysis of the screening process in the data preprocessing unit and sends the error influence coefficient to the monitoring and warning module.

[0044] The specific expression of the deformation variables of several different points at the target location collected by the data acquisition unit is as follows: ;

[0045] in, Represents a point dataset, Respectively represent the deformation of the first target point collected, , No. The deformation of the target point, , No. The deformation of the target point, .

[0046] The data preprocessing unit is used to preprocess the data collected by the data acquisition unit and filter out interference. The specific steps are as follows:

[0047] Step A1: Obtain atmospheric error similarity coefficient ;

[0048] Step A1: Obtain atmospheric error similarity coefficient The specific steps are as follows:

[0049] Step A1.1: Obtain real-time meteorological monitoring data from ground meteorological stations, including water vapor content and temperature;

[0050] Step A1.2: Calculate the atmospheric error similarity coefficient based on real-time meteorological monitoring data obtained from ground meteorological stations The specific expression is as follows:

[0051] ;

[0052] in, represents the atmospheric error similarity coefficient, It represents the maximum water vapor content that does not affect the detection accuracy of interferometric synthetic aperture radar. express The water vapor content sampled at each moment, represents the water vapor content sampled at time t, express The absolute value of represents the temperature sampled at time t, Indicates the temperature sampled at time t-1;

[0053] Step A2: Obtaining the Interferometric Synthetic Aperture Radar Noise Error Coefficient ;

[0054] Step A2: Obtaining the Interferometric Synthetic Aperture Radar Noise Error Coefficient The specific steps are as follows:

[0055] Step A2.1: Obtain the electrical interference coefficient during the interferometric synthetic aperture radar imaging process, including the signal-to-noise ratio;

[0056] Step A2.2: Calculate the noise error coefficient of the interferometric synthetic aperture radar based on the signal-to-noise ratio obtained in step A2.1 The specific expression is as follows: ;

[0057] in, represents the signal-to-noise ratio, represents the noise error coefficient of interferometric synthetic aperture radar;

[0058] Step A3: Based on the atmospheric error similarity coefficient in step A1 and the interferometric synthetic aperture radar noise error coefficient in step A2 Construction error influence coefficient , the specific expression is as follows: ;

[0059] in, and They represent weight coefficients that sum to 1, represents the error influence coefficient, represents the atmospheric error similarity coefficient, represents the noise error coefficient of interferometric synthetic aperture radar;

[0060] The monitoring and early warning module includes an error early warning unit. The error early warning unit determines whether to issue an early warning based on the error influence coefficient. The specific steps of the error early warning unit are as follows: Exceeding the error impact threshold When the warning is issued, the point data set The deformation variables and error influence coefficients of all points in Perform corresponding storage and terminate subsequent work. When the error influence coefficient Does not exceed the error impact threshold When , the error warning unit is terminated;

[0061] The deformation analysis module is connected to the data acquisition module. The deformation analysis module includes a first analysis unit, a second analysis unit, and a comprehensive analysis unit. The first analysis unit calculates the overall deformation coefficient based on the deformation variables of different points at the target location. The second analysis unit calculates the overall deformation trend based on the deformation variables of different points at the target location. The comprehensive analysis unit calculates the comprehensive risk coefficient based on the overall deformation coefficient and the overall deformation trend, thereby determining the micro-scale deformation risk of the target location and reducing the interference of environmental and physical factors on the data, ensuring the accuracy of the analysis.

[0062] The specific expression of the overall deformation coefficient calculated by the first analysis unit based on the deformation variables at different points of the target location is as follows:

[0063] ;

[0064] in, represents the overall deformation coefficient, Represents a point dataset The maximum value of Indicates total The deformation variables of each target point are summed up, and the second analysis unit calculates the overall deformation trend based on the deformation variables of different points in the target location as follows:

[0065] ;

[0066] in, Indicates the overall deformation trend, Indicates total The deformation variables of the target points are summed up. Indicates the last The deformation of the target point, Indicates the last sampling The deformation variables of each target point are summed up, and the comprehensive analysis unit calculates the comprehensive risk coefficient based on the overall deformation coefficient and the overall deformation trend. The specific expression is as follows: ;

[0067] in, represents the overall deformation coefficient, Indicates the overall deformation trend, represents the comprehensive risk factor, represents a natural constant, Represents a logarithmic function whose base is a natural constant;

[0068] The monitoring and early warning module also includes a comprehensive risk warning unit and a scheduling unit. The comprehensive risk warning unit determines whether to issue an early warning based on the comprehensive risk coefficient. The scheduling unit is called when the comprehensive risk unit issues an early warning and executes personnel scheduling.

[0069] The specific steps for the comprehensive risk warning unit to determine whether to issue a warning based on the comprehensive risk coefficient are as follows: Exceeding the comprehensive risk threshold When the comprehensive risk factor Does not exceed the comprehensive risk threshold No warning is given when

[0070] The dispatch unit is called when the comprehensive risk unit issues an early warning, and the specific strategy for executing personnel dispatch is to dispatch personnel based on distance as the priority.

[0071] Example 1: Obtaining the atmospheric error similarity coefficient during this implementation and noise error coefficient of interferometric synthetic aperture radar The data involved are shown in Table 1 below:

[0072] Table 1

[0073] ;

[0074] ; Now calculate , now exceeds , the point dataset The deformation variables and error influence coefficients of all points in Perform corresponding storage, terminate subsequent work, issue an early warning, and wait for the next cycle day for sampling;

[0075] Example 2: During this implementation In line with the scope, no display, At this time, the deformation analysis module obtains the deformation variables of several different points. In this application, all deformation variables are the results after absolute value processing, specifically 20 corresponding points, see Table 2 below:

[0076] Table 2

[0077] ;

[0078] At this time, the calculation , , , Exceeding the comprehensive risk threshold At this time, personnel are dispatched based on distance as the priority to conduct a secondary survey of the on-site environment and take corresponding measures, so as to accurately estimate the risk situation in the corresponding area;

[0079] This embodiment only provides a feasible solution and does not represent the optimal solution. The setting of the threshold value is for the convenience of comparison. The value of the threshold value depends on the amount of sample data and the number of cardinality set by those skilled in the art for each set of sample data. As long as it does not affect the proportional relationship between the parameter and the quantized value, the value of the weight can be determined by those skilled in the art based on each sample data and multiple rounds of experiments.

[0080] While embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions, and variations may be made to these embodiments without departing from the principles and spirit of the invention, and that the scope of the invention is defined by the appended claims and their equivalents.

Claims

1. A micro-scale deformation monitoring system based on interferometric synthetic aperture radar, characterized by: Including data acquisition module, deformation analysis module and monitoring and early warning module; The data acquisition module includes a data acquisition unit, a data preprocessing unit, and a data error warning unit. The data acquisition unit collects deformation variables of several different points of the target location through interferometric synthetic aperture radar. The data preprocessing unit is used to preprocess the data collected by the data acquisition unit and filter out interference. The data error warning unit obtains an error influence coefficient based on a comprehensive analysis of the screening process in the data preprocessing unit and sends the error influence coefficient to the monitoring and warning module. The data preprocessing unit is used to preprocess the data collected by the data collection unit and filter out interference. The specific steps are as follows: Step A1: Obtain atmospheric error similarity coefficient , the specific steps are as follows: Step A1.1: Obtain real-time meteorological monitoring data from ground meteorological stations, including water vapor content and temperature; Step A1.2: Calculate the atmospheric error similarity coefficient based on real-time meteorological monitoring data obtained from ground meteorological stations The specific expression is as follows: ； in, represents the atmospheric error similarity coefficient, It represents the maximum water vapor content that does not affect the detection accuracy of interferometric synthetic aperture radar. express The water vapor content sampled at each moment, represents the water vapor content sampled at time t, express The absolute value of represents the temperature sampled at time t, Indicates the temperature sampled at time t-1; Step A2: Obtaining the Interferometric Synthetic Aperture Radar Noise Error Coefficient , the specific steps are as follows: Step A2.1: Obtain the electrical interference coefficient during the interferometric synthetic aperture radar imaging process, including the signal-to-noise ratio; Step A2.2: Calculate the noise error coefficient of the interferometric synthetic aperture radar based on the signal-to-noise ratio obtained in step A2.1 The specific expression is as follows: ;in, represents the signal-to-noise ratio, represents the noise error coefficient of interferometric synthetic aperture radar; Step A3: Based on the atmospheric error similarity coefficient in step A1 and the interferometric synthetic aperture radar noise error coefficient in step A2 Construction error influence coefficient , the specific expression is as follows: ; in, and They represent weight coefficients that sum to 1, represents the error influence coefficient, represents the atmospheric error similarity coefficient, represents the noise error coefficient of interferometric synthetic aperture radar; The monitoring and early warning module includes an error early warning unit, which determines whether to issue an early warning based on the error influence coefficient; The deformation analysis module is data-connected to the data acquisition module. The deformation analysis module includes a first analysis unit, a second analysis unit, and a comprehensive analysis unit. The first analysis unit calculates the overall deformation coefficient based on the deformation variables at different points of the target location. The second analysis unit calculates the overall deformation trend based on the deformation variables at different points of the target location. The comprehensive analysis unit calculates the comprehensive risk coefficient based on the overall deformation coefficient and the overall deformation trend. The monitoring and early warning module also includes a comprehensive risk early warning unit and a scheduling unit. The comprehensive risk early warning unit determines whether to issue an early warning based on the comprehensive risk coefficient. The scheduling unit is called when the comprehensive risk unit issues an early warning and executes personnel scheduling.

2. The micro-scale deformation monitoring system based on interferometric synthetic aperture radar according to claim 1, characterized in that: The specific expression for the data acquisition unit to collect the deformation variables of several different points at the target location is as follows: ; in, Represents a point dataset, Respectively represent the deformation of the first target point collected, , No. The deformation of the target point, , No. The deformation of the target point, .

3. The micro-scale deformation monitoring system based on interferometric synthetic aperture radar according to claim 2, characterized in that: The specific steps of the error warning unit judging whether to issue a warning based on the error influence coefficient are as follows: Exceeding the error impact threshold When the warning is issued, the point data set The deformation variables and error influence coefficients of all points in Perform corresponding storage and terminate subsequent work. When the error influence coefficient Does not exceed the error impact threshold When , the error warning unit is terminated.

4. The micro-scale deformation monitoring system based on interferometric synthetic aperture radar according to claim 2, characterized in that: The specific expression of the overall deformation coefficient calculated by the first analysis unit based on the deformation variables at different points of the target location is as follows: ； in, represents the overall deformation coefficient, Represents a point dataset The maximum value of Indicates total The deformation variables of the target points are summed.

5. The micro-scale deformation monitoring system based on interferometric synthetic aperture radar according to claim 4 is characterized in that: The second analysis unit calculates the overall deformation trend based on the deformation variables at different points of the target location and obtains the following specific expression: ； in, Indicates the overall deformation trend, Indicates total The deformation variables of the target points are summed up. Indicates the last The deformation of the target point, Indicates the last sampling The deformation variables of the target points are summed.

6. The micro-scale deformation monitoring system based on interferometric synthetic aperture radar according to claim 5, characterized in that: The comprehensive analysis unit calculates the comprehensive risk coefficient based on the overall deformation coefficient and the overall deformation trend, and obtains the following specific expression: ; in, represents the overall deformation coefficient, Indicates the overall deformation trend, represents the comprehensive risk factor, represents a natural constant, Represents a logarithmic function whose base is a natural constant.

7. The micro-scale deformation monitoring system based on interferometric synthetic aperture radar according to claim 6, characterized in that: The specific steps of the comprehensive risk warning unit to determine whether to issue a warning based on the comprehensive risk coefficient are as follows: Exceeding the comprehensive risk threshold When the comprehensive risk factor Does not exceed the comprehensive risk threshold No warning is given when The scheduling unit is called when the comprehensive risk unit issues an early warning, and the specific strategy for executing personnel scheduling is to schedule personnel based on distance as the priority.

Images