CN111854596A - Ground surface deformation real-time monitoring method and system based on ground-based interferometric radar data - Google Patents
Ground surface deformation real-time monitoring method and system based on ground-based interferometric radar data Download PDFInfo
- Publication number
- CN111854596A CN111854596A CN202010868109.2A CN202010868109A CN111854596A CN 111854596 A CN111854596 A CN 111854596A CN 202010868109 A CN202010868109 A CN 202010868109A CN 111854596 A CN111854596 A CN 111854596A
- Authority
- CN
- China
- Prior art keywords
- deformation
- pixel
- accumulated
- time
- rate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000001133 acceleration Effects 0.000 claims abstract description 48
- 238000012545 processing Methods 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000012937 correction Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 7
- 238000009825 accumulation Methods 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000012300 Sequence Analysis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000005305 interferometry Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 241001252483 Kalimeris Species 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/16—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a ground surface deformation real-time monitoring method and system based on ground interference radar data. The method comprises the following steps: solving a differential interference pattern of radar single-vision complex images of target monitoring areas at two adjacent moments; determining the deformation quantity of each pixel at adjacent moments according to the differential interference pattern; according to the deformation quantity of each pixel at each adjacent moment, the accumulated deformation quantity of each pixel relative to the initial moment is obtained; solving the deformation rate of each pixel according to the accumulated deformation quantity of each pixel; according to the deformation rate of each pixel, solving the deformation acceleration of each pixel at the current moment; and determining whether to send out an early warning or not according to the accumulated deformation amount of each pixel at the current moment relative to the initial moment, the deformation rate of each pixel at the current moment and the deformation acceleration of each pixel at the current moment by combining the accumulated deformation threshold, the deformation rate threshold and the deformation acceleration threshold. The method and the system for monitoring the deformation of the earth surface in real time can realize accurate monitoring and early warning of the deformation of the earth surface.
Description
Technical Field
The invention relates to the field of ground surface deformation monitoring, in particular to a ground surface deformation real-time monitoring method and system based on ground interference radar data.
Background
Ground-based synthetic aperture radar interferometry (GB-InSAR) is a deformation monitoring technology based on a microwave detection active imaging mode, micro-deformation monitoring superior to millimeter level can be realized through the radar interferometry technology, earth surface deformation information can be accurately acquired without contacting a target area, and the method has the technical advantages of high space-time resolution, high flexibility, high precision and the like. In recent years, landslide, glacier movement and mining area settlement often occur, and it is very important to provide a method capable of accurately monitoring and early warning ground surface deformation.
Disclosure of Invention
The invention aims to provide a ground surface deformation real-time monitoring method and system based on ground-based interferometric radar data, which can realize accurate monitoring and early warning of ground surface deformation.
In order to achieve the purpose, the invention provides the following scheme:
a ground surface deformation real-time monitoring method based on ground-based interferometric radar data comprises the following steps:
solving a differential interference pattern of radar single-vision complex images of target monitoring areas at two adjacent moments;
determining the deformation quantity of each pixel at adjacent moments according to the differential interference pattern;
according to the deformation quantity of each pixel at each adjacent moment, the accumulated deformation quantity of each pixel relative to the initial moment is obtained;
solving the deformation rate of each pixel according to the accumulated deformation quantity of each pixel;
according to the deformation rate of each pixel, solving the deformation acceleration of each pixel at the current moment;
and determining whether to send out an early warning or not according to the accumulated deformation amount of each pixel at the current moment relative to the initial moment, the deformation rate of each pixel at the current moment and the deformation acceleration of each pixel at the current moment by combining the accumulated deformation threshold, the deformation rate threshold and the deformation acceleration threshold.
Optionally, the obtaining of the differential interferogram of the radar single-view complex image of the target monitoring area at two adjacent moments specifically includes:
acquiring radar single-view complex images of two adjacent moments of a target monitoring area;
and carrying out differential interference processing on the radar single-view complex images at two adjacent moments to obtain differential interference images at two adjacent moments.
Optionally, the determining the deformation amount of each pixel at adjacent time according to the differential interference pattern specifically includes:
and carrying out filtering, phase unwrapping and atmospheric correction processing on the differential interference pattern, and determining the deformation amount of each pixel at adjacent time according to the processed differential interference pattern.
Optionally, the calculating of the accumulated deformation amount of each pixel relative to the initial time according to the deformation amount of each pixel at each adjacent time specifically includes:
and the accumulated deformation quantity of each pixel is obtained by accumulating the deformation quantities of the pixels at each adjacent moment.
Optionally, the obtaining of the deformation rate of each pixel according to the accumulated deformation amount of each pixel specifically includes:
and in a set time period, performing least square linear fitting on the accumulated deformation quantity of the pixel, and determining the deformation rate of the pixel in the set time period according to the slope of a fitting straight line.
Optionally, the obtaining of the deformation rate of each pixel according to the accumulated deformation amount of each pixel specifically includes:
according to the accumulated deformation amount of the pixel at each moment, establishing a curve of the accumulated deformation amount of the pixel along with the change of time, and recording the curve as an accumulated deformation-time curve;
deriving the accumulated deformation-time curve to obtain a deformation rate-time curve of the pixel;
and determining the deformation rate of the pixel based on the deformation rate-time curve.
Optionally, the obtaining of the deformation acceleration of each pixel at the current moment according to the deformation rate of each pixel specifically includes:
the deformation rate-time curve of the pixel is derived to obtain a deformation acceleration-time curve of the pixel;
and determining the deformation acceleration of the pixel based on the deformation acceleration-time curve.
The invention also provides a ground surface deformation real-time monitoring system based on the ground interference radar data, which comprises the following components:
the differential interference image solving module is used for solving the differential interference image of the radar single-vision complex image of the target monitoring area at two adjacent moments;
the deformation quantity solving module is used for determining the deformation quantity of each pixel at adjacent moments according to the differential interference pattern;
the accumulated deformation quantity calculation module is used for calculating the accumulated deformation quantity of each pixel relative to the initial time according to the deformation quantity of each pixel at each adjacent time;
the deformation rate calculation module is used for calculating the deformation rate of each pixel according to the accumulated deformation quantity of each pixel;
the deformation acceleration calculation module is used for calculating the deformation acceleration of each pixel at the current moment according to the deformation rate of each pixel;
and the early warning module is used for determining whether to send out early warning or not according to the accumulated deformation amount of each pixel at the current moment relative to the initial moment, the deformation rate of each pixel at the current moment and the deformation acceleration of each pixel at the current moment by combining the accumulated deformation threshold, the deformation rate threshold and the deformation acceleration threshold.
Alternatively to this, the first and second parts may,
the differential interferogram solving module specifically comprises:
the radar image acquisition unit is used for acquiring radar single-view complex images of two adjacent moments in a target monitoring area;
the differential interference processing unit is used for carrying out differential interference processing on the radar single-view complex images at two adjacent moments to obtain differential interference images at two adjacent moments;
the deformation quantity calculating module specifically comprises:
and the deformation quantity calculating unit is used for carrying out filtering, phase unwrapping and atmospheric correction processing on the differential interference pattern and determining the deformation quantity of each pixel at adjacent time according to the processed differential interference pattern.
Alternatively to this, the first and second parts may,
the accumulated deformation quantity calculating module specifically comprises:
the accumulated deformation quantity calculation unit is used for calculating the accumulated deformation quantity of each pixel through the accumulation of the deformation quantity of the pixel at each adjacent moment;
the deformation rate solving module specifically comprises:
the first deformation rate solving unit is used for performing least square linear fitting on the accumulated deformation quantity of the pixel in a set time period and determining the deformation rate of the pixel in the set time period according to the slope of a fitting straight line;
the second deformation rate obtaining unit is used for establishing an accumulated deformation-time curve of the accumulated deformation quantity of the pixel along with the change of time according to the accumulated deformation quantity of the pixel at each moment, deriving the accumulated deformation-time curve to obtain a deformation rate-time curve of the pixel, and determining the deformation rate of the pixel based on the deformation rate-time curve;
the deformation acceleration calculation module specifically includes:
and the deformation acceleration calculation unit is used for calculating the derivative of the deformation rate-time curve of the pixel element to obtain the deformation acceleration-time curve of the pixel element, and determining the deformation acceleration of the pixel element based on the deformation acceleration-time curve.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects: according to the ground surface deformation real-time monitoring method and system based on the ground interference radar data, the obtained radar observation data are processed in real time, the accumulated deformation quantity, the deformation rate and the deformation acceleration of each pixel in the target monitoring area are obtained by solving the deformation of an interference pattern and using a time sequence analysis processing method, and ground surface deformation early warning is carried out by integrating three threshold parameters of an accumulated deformation quantity threshold, a deformation rate threshold and a deformation acceleration threshold, so that the monitoring precision and the early warning precision of ground surface deformation are improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
Fig. 1 is a flowchart of a ground surface deformation real-time monitoring method based on ground-based interferometric radar data according to embodiment 1 of the present invention;
fig. 2 is a structural diagram of a ground surface deformation real-time monitoring system based on ground-based interferometric radar data according to embodiment 2 of the present invention;
FIG. 3 is a flow chart of the data processing of the ground-based interferometric radar of the present invention;
FIG. 4 is a graph of cumulative deformation in the present invention;
FIG. 5 is a graph of average deformation rate in the present invention;
FIG. 6 is a line graph of three curves in the present invention, FIG. 6(a) is a graph of cumulative deformation versus time, FIG. 6(b) is a graph of deformation rate versus time, and FIG. 6(c) is a graph of deformation acceleration versus time;
FIG. 7 is a graph of data processing accuracy of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Example 1
Referring to fig. 1, the present embodiment provides a ground surface deformation real-time monitoring method based on ground-based interferometric radar data, including:
step 101: solving a differential interference pattern of radar single-vision complex images of target monitoring areas at two adjacent moments;
step 102: determining the deformation quantity of each pixel at adjacent moments according to the differential interference pattern;
step 103: according to the deformation quantity of each pixel at each adjacent moment, the accumulated deformation quantity of each pixel relative to the initial moment is obtained;
step 104: solving the deformation rate of each pixel according to the accumulated deformation quantity of each pixel;
step 105: according to the deformation rate of each pixel, solving the deformation acceleration of each pixel at the current moment;
step 106: and determining whether to send out an early warning or not according to the accumulated deformation amount of each pixel at the current moment relative to the initial moment, the deformation rate of each pixel at the current moment and the deformation acceleration of each pixel at the current moment by combining the accumulated deformation threshold, the deformation rate threshold and the deformation acceleration threshold.
All the images in the steps are stored in a numerical matrix form, a real numerical value or a complex value is stored at the position of each pixel point, and the image can be formed if all the pixel points have numerical values, so that the image is different from a color image storing indexes of various colors.
In step 101, performing differential interferogram processing on radar single-view complex images of target monitoring areas at two adjacent moments, which may specifically be: and carrying out conjugate multiplication on complex values of all corresponding pixels of two adjacent scenes of SLC data (radar single-view complex data) in the time sequence to obtain an interference phase. Suppose y1、y2Respectively the complex value, y, of a pixel point of two images2 *Is y2The conjugated complex number of (1) then has y1＝a1+b1i、y2 *＝a2-b2i, conjugate multiplication formula is shown as the following formula:
y1y2 *＝(a1+b1i)(a2-b2i)＝a0+b0i。
in step 102, filtering, phase unwrapping and atmospheric correction processing are firstly carried out on the differential interference pattern, and then the deformation amount of each pixel at adjacent time is determined according to the processed differential interference pattern. Filtering is used to attenuate the effect of noise on the interference phase. Since the phase obtained by interference does not have information of the whole number of cycles, i.e. is in a phase winding state, phase unwrapping is required to restore to an absolute interference phase. And the atmospheric states of the two scenes of SLC data acquisition are different, so that the interferogram is influenced by the atmospheric phase, therefore, based on the unwrapping phase, a high-quality point (a point which keeps high phase quality on a time sequence and can be a high-coherence point generally) is selected, a function model y of distance and phase is established as a r to solve the atmospheric phase, in the formula, y is the unwrapped phase, a is an unknown coefficient, r is the slant range of the radar to a high-quality target point, the unknown coefficient a can be solved by a least square method, and the atmospheric phase can be solved according to the radar slant range of the rest each point. And subtracting the atmospheric phase corresponding to the pixel point from the unwrapped phase of each pixel point on the unwrapped image to finish atmospheric correction, thus obtaining the deformation quantity of the interference pattern, wherein the deformation quantity of all the pixels is the deformation image.
In step 103, the method for calculating the accumulated deformation amount may specifically be: selecting a plurality of interferograms with continuous time sequence, sequentially adding the deformation quantities of the pixel points at the same position in different interferograms to obtain the accumulated deformation quantity of the pixel points, and calculating the accumulated deformation quantity of all the pixel points to obtain the accumulated deformation image between the corresponding moments of the first interferogram and the last interferogram. In the loop process of the method of the present invention, each time new SLC data is obtained, the deformation graph of the interference graph at the latest time can be obtained through step 101, and the total accumulated deformation graph (accumulated deformation graph between the previous time and the initial time) at the previous time is calculated in step 103 of the last iteration, if there is no previous iteration, the accumulated deformation graph with all pixel values being 0 is used, and the accumulated deformation graph and the value of the deformation image element at the latest time are correspondingly added to obtain the total accumulated deformation graph at the latest time.
As an embodiment, in step 104, a specific method for calculating the deformation rate may be: on the premise of assuming that the accumulated deformation amount of the target monitoring area in a certain period of time is changed linearly, the least square linear fitting is carried out on the accumulated deformation graph contained in the period of time to obtain the average deformation rate graph of the period of time. If the spatial-temporal incoherence between the data can reduce the precision of the result, in order to avoid the spatial-temporal incoherence problem caused by overlong time interval, the length of the time interval should be controlled, namely the quantity of the accumulated deformation graphs in the processing process.
The formula of the least square linear fitting is shown as the formula, and the formula is a solving formula of unknown quantity x in the formula:
Ax＝b
wherein,
a is a time coefficient matrix, b is a deformation value matrix obtained by accumulating the deformation graphs in the second step, and x is a matrix of the quantity to be solved. n denotes the nth SLC data and m denotes the mth data point in the first SLC data. t is tnTime interval between nth scene data and 1 st scene data, dmnIs a difference value of an m-th point between the n-th scene data and the cumulative distortion map of the 1-th scene data,
is the mean deformation rate (first order term), c, of the mth data point fittedmIs the constant term of the linear equation fitted to the mth data point. The matrix x of the quantity to be solved can be obtained by a deformation formula of the formula:
x＝(ATA)-1ATb
and obtaining the deformation rate of all points in the graph by calculating x, and obtaining a deformation rate graph.
As another embodiment, in step 104, the method for obtaining the image element shape shift may further include:
according to the accumulated deformation amount of the pixel at each moment, establishing a curve of the accumulated deformation amount of the pixel along with the change of time, and recording the curve as an accumulated deformation-time curve; deriving the accumulated deformation-time curve to obtain a deformation rate-time curve of the pixel; and determining the deformation rate of the pixel based on the deformation rate-time curve. On this basis, the method for solving the deformation acceleration of the pixel can be as follows: the deformation rate-time curve of the pixel is derived to obtain a deformation acceleration-time curve of the pixel; and determining the deformation acceleration of the pixel based on the deformation acceleration-time curve.
The cumulative deformation-time curve is longitudinal with the cumulative deformation amountAnd the axis takes the radar image acquisition time as a horizontal axis, arranges the accumulated deformation graphs according to the corresponding image acquisition time sequence, selects a certain same pixel point position to extract the accumulated deformation according to the time sequence, so that each accumulated deformation graph can obtain the coordinate of a curve point, the radar image acquisition time of the image is the abscissa of the curve point, the accumulated deformation quantity extracted from the image is the ordinate, and a plurality of curve points are connected into a line according to the image time sequence to obtain the accumulated deformation-time curve. If the curve formed by a certain pixel point on the graph is not accurate enough, a plurality of numerical values can be converted into one numerical value by averaging or taking a median value and the like for the numerical values of the points in a certain same area on the graph, the processed accumulated deformation values are extracted according to the time sequence to form different points on the curve, and then the curve points are connected into a line to obtain the accumulated deformation-time curve. The deformation rate-time curve can be obtained from the deformation rate map in the same manner as the cumulative deformation-time curve is obtained from the cumulative deformation map. The deformation acceleration-time curve is obtained on the basis of the deformation rate-time curve, and the calculation formula is shown as the formula, wherein t isnShowing the acquisition time of the radar image corresponding to the nth scene cumulative deformation diagram,
representing the average deformation rate of the latest period of time obtained in step three at that moment, anIndicating the acceleration calculated at that time. The method comprises the steps of taking a certain radar image acquisition moment as an abscissa, taking an acceleration calculated at the moment as an ordinate, forming points on a deformation acceleration-time curve, calculating accelerations at other moments to obtain more curve points, and connecting the curve points into a line according to a time sequence, namely the deformation acceleration-time curve.
In addition, the curve can more intuitively represent the trend of the region deformation in the form of an image, and provides convenience for adjusting the early warning threshold required in the step 106.
As an implementation manner, the early warning method in step 106 may specifically be: and comparing the values on the three curves with a preset early warning threshold respectively, and if the value corresponding to the curve point at the latest moment on one curve exceeds the threshold range, giving out early warning to send out warning information, or giving out early warning when the values of two or more curves exceed the threshold range simultaneously. And if the values of the three curves do not exceed the threshold range, continuing the observation.
Example 2
Referring to fig. 2, the present embodiment provides a ground surface deformation real-time monitoring system based on ground-based interferometric radar data, which includes:
the differential interference pattern solving module 201 is used for solving a differential interference pattern of radar single-vision complex images in target monitoring areas at two adjacent moments;
the deformation quantity solving module 202 is used for determining the deformation quantity of each pixel at adjacent time according to the differential interference pattern;
the accumulated deformation quantity calculation module 203 is used for calculating the accumulated deformation quantity of each pixel relative to the initial time according to the deformation quantity of each pixel at each adjacent time;
the deformation rate calculation module 204 is used for calculating the deformation rate of each pixel according to the accumulated deformation amount of each pixel;
the deformation acceleration calculation module 205 is configured to calculate a deformation acceleration of each pixel at the current time according to the deformation rate of each pixel;
and the early warning module 206 is configured to determine whether to send an early warning according to the accumulated deformation amount of each pixel at the current time relative to the initial time, the deformation rate of each pixel at the current time, and the deformation acceleration of each pixel at the current time, in combination with the accumulated deformation threshold, the deformation rate threshold, and the deformation acceleration threshold.
The differential interferogram obtaining module 201 specifically includes:
the radar image acquisition unit is used for acquiring radar single-view complex images of two adjacent moments in a target monitoring area;
the differential interference processing unit is used for carrying out differential interference processing on the radar single-view complex images at two adjacent moments to obtain differential interference images at two adjacent moments;
the deformation quantity calculating module 202 specifically includes:
and the deformation quantity calculating unit is used for carrying out filtering, phase unwrapping and atmospheric correction processing on the differential interference pattern and determining the deformation quantity of each pixel at adjacent time according to the processed differential interference pattern.
The accumulated deformation quantity calculating module 203 specifically includes:
the accumulated deformation quantity calculation unit is used for calculating the accumulated deformation quantity of each pixel through the accumulation of the deformation quantity of the pixel at each adjacent moment;
the deformation rate obtaining module 204 specifically includes:
the first deformation rate solving unit is used for performing least square linear fitting on the accumulated deformation quantity of the pixel in a set time period and determining the deformation rate of the pixel in the set time period according to the slope of a fitting straight line;
the second deformation rate obtaining unit is used for establishing an accumulated deformation-time curve of the accumulated deformation quantity of the pixel along with the change of time according to the accumulated deformation quantity of the pixel at each moment, deriving the accumulated deformation-time curve to obtain a deformation rate-time curve of the pixel, and determining the deformation rate of the pixel based on the deformation rate-time curve;
the deformation acceleration calculation module 205 specifically includes:
and the deformation acceleration calculation unit is used for calculating the derivative of the deformation rate-time curve of the pixel element to obtain the deformation acceleration-time curve of the pixel element, and determining the deformation acceleration of the pixel element based on the deformation acceleration-time curve.
The invention is explained below by way of experiments
GB-InSAR images with the time interval of 1 hour collected in the mining area of the Hetian Tangshan kalimeris are selected in the experiment, the collection equipment is MIMO, and the wave band is Ku wave band. The data processing flow is shown in fig. 3.
The method comprises the following steps: and carrying out differential interference processing on the latest acquired SLC and the SLC acquired at the previous time to calculate the deformation of the interference pattern.
Step two: and accumulating all the generated interference pattern deformation to obtain a new accumulated deformation pattern.
Step three: and (4) taking the deformation of the interference pattern generated by the latest 20 scenes, and carrying out time sequence processing to obtain a deformation rate diagram of the time period.
Step four: and generating three statistical curves according to all the accumulated deformation graphs and the deformation rate graphs.
Step five: and judging to alarm or continue observation according to the comparison of the set threshold value and the curve.
The specific software running time is influenced by the configuration of the computer where the software is located, 382 SLC data with the specification of 2096 × 246 are processed on an i7-7700HQ CPU computer provided with Intel, the time consumed by processing single data is within 25 seconds, the speed is high, and when the performance of the computer is enough, the real-time processing of the software can meet the requirement of quick drawing.
The calculated accumulated deformation graph is shown in fig. 4, and it is obvious by combining the legend that the deformation value of the blue region in the red box is large, and the blue region is the actual deformation region. The deformation rate map obtained from the previous 20 scenes is shown in fig. 5, and it can be seen that the region with the larger rate substantially coincides with the region with the larger deformation in fig. 4. Taking the average value of a small area of a scene image in the fourth step as the value of a certain point on the curve, and taking the small area in fig. 4 as the processed area, the three curves obtained are shown in fig. 6, and it can be seen that the accumulated deformation quantity is accumulated continuously with time, and the deformation rate and the acceleration are larger in the first 35 scene values and then fluctuate around the 0 value. The three curves can reflect the trend of the deformation of the region from unstable to gradually stable.
Ideally, the accumulated deformation amount calculated by the region stabilization point in the monitoring period should be 0. And when the accuracy of the processing result is analyzed, because the third party synchronously monitors the data to be lacked, whether the accumulated deformation curve of the stable point of the monitoring area is stable or not and whether the standard deviation of the residual error obtained by subtracting the model value from the interference value is smaller or not are used for judging the accuracy of the processing result. The model value is composed of two parts, an atmospheric model generated based on a distance function and a deformation model for time fitting the deformation. As the monitoring precision of the ground SAR system is mm level, the point precision with the standard deviation of data smaller than 1mm in real-time processing can be considered. The standard deviation calculation of the accumulated data residuals is performed by selecting 100 points with the size of 5 × 5 at four corners of the stable point region, and the statistical result of the standard deviation is shown in fig. 7. The total number of high-precision points with standard deviation of residual errors smaller than 1mm is 77, and the percentage is 77%, so that most point positions can meet the monitoring precision of submillimeter, and the result meets the requirements for judging the trend of the accumulated deformation curve and the landslide early warning.
According to the method and the system for monitoring the surface deformation in real time, the acquired radar observation data are processed in real time, the accumulated deformation quantity, the deformation rate and the deformation acceleration of each pixel in a target monitoring area are acquired by solving the deformation of an interference pattern and using a time sequence analysis processing method, and the surface deformation is early warned by integrating three threshold parameters of an accumulated deformation quantity threshold, a deformation rate threshold and a deformation acceleration threshold, so that the monitoring precision and the early warning precision of the surface deformation are improved.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (10)
1. A ground surface deformation real-time monitoring method based on ground-based interferometric radar data is characterized by comprising the following steps:
solving a differential interference pattern of radar single-vision complex images of target monitoring areas at two adjacent moments;
determining the deformation quantity of each pixel at adjacent moments according to the differential interference pattern;
according to the deformation quantity of each pixel at each adjacent moment, the accumulated deformation quantity of each pixel relative to the initial moment is obtained;
solving the deformation rate of each pixel according to the accumulated deformation quantity of each pixel;
according to the deformation rate of each pixel, solving the deformation acceleration of each pixel at the current moment;
and determining whether to send out an early warning or not according to the accumulated deformation amount of each pixel at the current moment relative to the initial moment, the deformation rate of each pixel at the current moment and the deformation acceleration of each pixel at the current moment by combining the accumulated deformation threshold, the deformation rate threshold and the deformation acceleration threshold.
2. The ground-based interferometric radar data-based ground surface deformation real-time monitoring method according to claim 1, wherein the obtaining of the differential interferogram of the radar single vision complex image of the target monitoring area at two adjacent moments specifically comprises:
acquiring radar single-view complex images of two adjacent moments of a target monitoring area;
and carrying out differential interference processing on the radar single-view complex images at two adjacent moments to obtain differential interference images at two adjacent moments.
3. The ground-based interferometric radar data-based surface deformation real-time monitoring method according to claim 1, wherein the determining the deformation amount of each pixel at adjacent time according to the differential interferogram specifically comprises:
and carrying out filtering, phase unwrapping and atmospheric correction processing on the differential interference pattern, and determining the deformation amount of each pixel at adjacent time according to the processed differential interference pattern.
4. The ground-based interferometric radar data-based ground surface deformation real-time monitoring method according to claim 1, wherein the accumulated deformation quantity of each pixel relative to the initial time is obtained according to the deformation quantity of each pixel at each adjacent time, and specifically comprises:
and the accumulated deformation quantity of each pixel is obtained by accumulating the deformation quantities of the pixels at each adjacent moment.
5. The ground-based interferometric radar data-based surface deformation real-time monitoring method according to claim 1, wherein the obtaining of the deformation rate of each pixel according to the accumulated deformation quantity of each pixel specifically comprises:
and in a set time period, performing least square linear fitting on the accumulated deformation quantity of the pixel, and determining the deformation rate of the pixel in the set time period according to the slope of a fitting straight line.
6. The ground-based interferometric radar data-based surface deformation real-time monitoring method according to claim 1, wherein the obtaining of the deformation rate of each pixel according to the accumulated deformation quantity of each pixel specifically comprises:
according to the accumulated deformation amount of the pixel at each moment, establishing a curve of the accumulated deformation amount of the pixel along with the change of time, and recording the curve as an accumulated deformation-time curve;
deriving the accumulated deformation-time curve to obtain a deformation rate-time curve of the pixel;
and determining the deformation rate of the pixel based on the deformation rate-time curve.
7. The ground-based interferometric radar data-based ground surface deformation real-time monitoring method according to claim 6, wherein the obtaining of the deformation acceleration of each pixel at the current time according to the deformation rate of each pixel specifically comprises:
the deformation rate-time curve of the pixel is derived to obtain a deformation acceleration-time curve of the pixel;
and determining the deformation acceleration of the pixel based on the deformation acceleration-time curve.
8. The utility model provides a ground surface deformation real-time monitoring system based on ground interference radar data which characterized in that includes:
the differential interference image solving module is used for solving the differential interference image of the radar single-vision complex image of the target monitoring area at two adjacent moments;
the deformation quantity solving module is used for determining the deformation quantity of each pixel at adjacent moments according to the differential interference pattern;
the accumulated deformation quantity calculation module is used for calculating the accumulated deformation quantity of each pixel relative to the initial time according to the deformation quantity of each pixel at each adjacent time;
the deformation rate calculation module is used for calculating the deformation rate of each pixel according to the accumulated deformation quantity of each pixel;
the deformation acceleration calculation module is used for calculating the deformation acceleration of each pixel at the current moment according to the deformation rate of each pixel;
and the early warning module is used for determining whether to send out early warning or not according to the accumulated deformation amount of each pixel at the current moment relative to the initial moment, the deformation rate of each pixel at the current moment and the deformation acceleration of each pixel at the current moment by combining the accumulated deformation threshold, the deformation rate threshold and the deformation acceleration threshold.
9. The ground-based interferometric radar data based surface deformation real-time monitoring system of claim 8,
the differential interferogram solving module specifically comprises:
the radar image acquisition unit is used for acquiring radar single-view complex images of two adjacent moments in a target monitoring area;
the differential interference processing unit is used for carrying out differential interference processing on the radar single-view complex images at two adjacent moments to obtain differential interference images at two adjacent moments;
the deformation quantity calculating module specifically comprises:
and the deformation quantity calculating unit is used for carrying out filtering, phase unwrapping and atmospheric correction processing on the differential interference pattern and determining the deformation quantity of each pixel at adjacent time according to the processed differential interference pattern.
10. The ground-based interferometric radar data based surface deformation real-time monitoring system of claim 8,
the accumulated deformation quantity calculating module specifically comprises:
the accumulated deformation quantity calculation unit is used for calculating the accumulated deformation quantity of each pixel through the accumulation of the deformation quantity of the pixel at each adjacent moment;
the deformation rate solving module specifically comprises:
the first deformation rate solving unit is used for performing least square linear fitting on the accumulated deformation quantity of the pixel in a set time period and determining the deformation rate of the pixel in the set time period according to the slope of a fitting straight line;
the second deformation rate obtaining unit is used for establishing an accumulated deformation-time curve of the accumulated deformation quantity of the pixel along with the change of time according to the accumulated deformation quantity of the pixel at each moment, deriving the accumulated deformation-time curve to obtain a deformation rate-time curve of the pixel, and determining the deformation rate of the pixel based on the deformation rate-time curve;
the deformation acceleration calculation module specifically includes:
and the deformation acceleration calculation unit is used for calculating the derivative of the deformation rate-time curve of the pixel element to obtain the deformation acceleration-time curve of the pixel element, and determining the deformation acceleration of the pixel element based on the deformation acceleration-time curve.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010868109.2A CN111854596A (en) | 2020-08-26 | 2020-08-26 | Ground surface deformation real-time monitoring method and system based on ground-based interferometric radar data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010868109.2A CN111854596A (en) | 2020-08-26 | 2020-08-26 | Ground surface deformation real-time monitoring method and system based on ground-based interferometric radar data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111854596A true CN111854596A (en) | 2020-10-30 |
Family
ID=72968177
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010868109.2A Pending CN111854596A (en) | 2020-08-26 | 2020-08-26 | Ground surface deformation real-time monitoring method and system based on ground-based interferometric radar data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111854596A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113466846A (en) * | 2021-06-10 | 2021-10-01 | 上海铁路北斗测量工程技术有限公司 | Migration monitoring method and system for high-speed railway bridge |
| CN114894082A (en) * | 2022-04-21 | 2022-08-12 | 北方雷科(安徽)科技有限公司 | Mine-based slope deformation landslide early warning method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001083243A (en) * | 1999-09-13 | 2001-03-30 | Mitsubishi Electric Corp | Extraction apparatus for three-dimensional information on landform by interference-type synthetic aperture radar |
| CN104515988A (en) * | 2014-12-16 | 2015-04-15 | 中国安全生产科学研究院 | Side slope safety monitoring and prewarning method based on ground-based synthetic aperture radar |
| CN105974411A (en) * | 2016-04-28 | 2016-09-28 | 电子科技大学 | DINSAR (differential interferometric synthetic aperture radar)-based high-voltage transmission tower top tilt displacement monitoring method |
| CN108548497A (en) * | 2018-06-12 | 2018-09-18 | 西南交通大学 | It is a kind of that system and monitoring method are monitored based on the high gradient slope of radar and infrared ray |
| CN109709550A (en) * | 2019-01-17 | 2019-05-03 | 武汉大学 | A kind of reservoir stability deformation monitoring processing method based on InSAR image data |
-
2020
- 2020-08-26 CN CN202010868109.2A patent/CN111854596A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001083243A (en) * | 1999-09-13 | 2001-03-30 | Mitsubishi Electric Corp | Extraction apparatus for three-dimensional information on landform by interference-type synthetic aperture radar |
| CN104515988A (en) * | 2014-12-16 | 2015-04-15 | 中国安全生产科学研究院 | Side slope safety monitoring and prewarning method based on ground-based synthetic aperture radar |
| CN105974411A (en) * | 2016-04-28 | 2016-09-28 | 电子科技大学 | DINSAR (differential interferometric synthetic aperture radar)-based high-voltage transmission tower top tilt displacement monitoring method |
| CN108548497A (en) * | 2018-06-12 | 2018-09-18 | 西南交通大学 | It is a kind of that system and monitoring method are monitored based on the high gradient slope of radar and infrared ray |
| CN109709550A (en) * | 2019-01-17 | 2019-05-03 | 武汉大学 | A kind of reservoir stability deformation monitoring processing method based on InSAR image data |
Non-Patent Citations (3)
| Title |
|---|
| 周吕 等: "GBInSAR 河流滑坡变形监测应用", 《地理空间信息》 * |
| 杨红磊 等: "GB-InSAR监测大型露天矿边坡形变", 《地球物理学进展》 * |
| 黄胜军: "《复杂岩体边坡变形与失稳预测研究》", 31 August 2016, 冶金工业出版社 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113466846A (en) * | 2021-06-10 | 2021-10-01 | 上海铁路北斗测量工程技术有限公司 | Migration monitoring method and system for high-speed railway bridge |
| CN113466846B (en) * | 2021-06-10 | 2024-03-26 | 上海铁路北斗测量工程技术有限公司 | High-speed railway bridge deviation monitoring method and system |
| CN114894082A (en) * | 2022-04-21 | 2022-08-12 | 北方雷科(安徽)科技有限公司 | Mine-based slope deformation landslide early warning method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2021227423A1 (en) | Insar digital elevation model construction method and system based on dynamic baseline | |
| CN115856889B (en) | InSAR time sequence deformation monitoring method capable of automatically correcting errors | |
| CN111059998B (en) | High-resolution-based time sequence InSAR deformation monitoring method and system | |
| CN107102333B (en) | Satellite-borne InSAR long and short baseline fusion unwrapping method | |
| Reigber et al. | Refined estimation of time-varying baseline errors in airborne SAR interferometry | |
| CN106249236B (en) | A kind of spaceborne InSAR long-short baselines image joint method for registering | |
| CN107014399B (en) | Combined calibration method for satellite-borne optical camera-laser range finder combined system | |
| CN111854596A (en) | Ground surface deformation real-time monitoring method and system based on ground-based interferometric radar data | |
| CN109633648A (en) | A kind of more baseline phase estimation devices and method based on possibility predication | |
| CN108196267B (en) | GNSS CP technology-based uninterrupted time transfer method | |
| CN110109106A (en) | A kind of InSAR interferometric phase unwrapping method in region with a varied topography | |
| CN109752715B (en) | SAR data total-dispersion body detection method and device | |
| DE102020211472A1 (en) | SYSTEM AND METHOD FOR IMPROVING A NON-INERTIAL TRACKING SYSTEM WITH INERTIAL RESTRICTIONS | |
| CN112050725A (en) | Surface deformation monitoring method fusing InSAR and GPS | |
| KR20190043478A (en) | System and method for generating high resolution digital elevation model | |
| CN110109105A (en) | A method of the InSAR technical monitoring Ground Deformation based on timing | |
| CN107869968A (en) | A kind of quick three-dimensional scan method and system suitable for complex object surface | |
| CN111398959A (en) | InSAR time sequence earth surface deformation monitoring method based on earth surface stress strain model | |
| CN103454636A (en) | Differential interferometric phase estimation method based on multi-pixel covariance matrixes | |
| CN113204023B (en) | Dual-polarization phase optimization earth surface deformation monitoring method combining PS target and DS target | |
| CN113848550A (en) | Ground radar adaptive threshold permanent scatterer identification method, device and storage medium | |
| CN103323840A (en) | Method for time alignment between interference SAR echo data and platform motion and gesture data | |
| CN103325092A (en) | Method and device for generating two-dimensional phase disentanglement quality picture | |
| CN109752698A (en) | A kind of inertial navigation method for estimating error of airborne synthetic aperture radar | |
| CN114332191A (en) | Three-dimensional point cloud error compensation method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201030 |