CN112540369A - Landslide three-dimensional deformation resolving method and system integrating GNSS and lifting rail time sequence InSAR - Google Patents
Landslide three-dimensional deformation resolving method and system integrating GNSS and lifting rail time sequence InSAR Download PDFInfo
- Publication number
- CN112540369A CN112540369A CN202011355738.1A CN202011355738A CN112540369A CN 112540369 A CN112540369 A CN 112540369A CN 202011355738 A CN202011355738 A CN 202011355738A CN 112540369 A CN112540369 A CN 112540369A
- Authority
- CN
- China
- Prior art keywords
- deformation
- dimensional
- insar
- gnss
- landslide
- 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
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000012544 monitoring process Methods 0.000 claims abstract description 117
- 230000001174 ascending effect Effects 0.000 claims abstract description 38
- 238000004364 calculation method Methods 0.000 claims abstract description 37
- 238000001514 detection method Methods 0.000 claims abstract description 21
- 230000000007 visual effect Effects 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000012806 monitoring device Methods 0.000 claims description 4
- 230000007774 longterm Effects 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 15
- 238000000354 decomposition reaction Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Classifications
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/86—Combinations of radar systems with non-radar systems, e.g. sonar, direction finder
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
The invention discloses a landslide three-dimensional deformation calculation method fusing GNSS and lifting rail time sequence InSAR, belongs to the technical field of risk identification, and is used for solving the technical problem that the existing detection mode is incomplete in result. The method comprises the following steps: acquiring a plurality of original images of the ascending and descending radar satellite in a GNSS monitoring period of a target monitoring area; acquiring the surface deformation of a target monitoring area under the view angles of an ascending and descending rail radar through a time sequence InSAR method; decomposing the ground surface deformation of the target monitoring area under the visual angles of the rail ascending and descending radar in the three-dimensional direction; and calculating the three-dimensional deformation of the target monitoring area by fusing an InSAR-GNSS three-dimensional joint solution model, and identifying the landslide of the target monitoring area according to the three-dimensional deformation. The method effectively solves the problem that the InSAR is difficult to obtain the three-dimensional deformation, effectively unifies the high time resolution and the high plane position precision of the GNSS and the high space resolution and the elevation deformation precision of the InSAR technology, and has high practicability for obtaining the three-dimensional deformation field of the landslide.
Description
Technical Field
The invention relates to the technical field of risk identification, in particular to a landslide three-dimensional deformation calculation method and system fusing GNSS and lifting rail time sequence InSAR.
Background
The method is used for accurately and effectively monitoring the landslide, is the basis for early forecasting of landslide disasters, and is also an important means for landslide disaster risk assessment and disaster prevention and reduction. The traditional side slope contact type safety monitoring mode has high manufacturing cost, difficult equipment maintenance and low reliability, can only acquire 'point-shaped' monitoring information, is difficult to acquire 'surface-shaped' dynamic change information of the side slope, and is not suitable for large-range deformation monitoring of the bank slope in a reservoir area. The Synthetic Aperture Radar interferometry (InSAR) technology can stably and continuously observe the earth surface in a fixed revisit period for a long time in a large range, is convenient to reveal the space-time evolution rule of the geoscience phenomenon, and has great potential in the landslide monitoring field in recent years.
However, the InSAR result of a single platform or track can only provide the deformation of the ground surface on a one-dimensional radar Line of Sight (LOS), and the deformation of a landslide often slides down along a slope surface under the influence of gravity, so that the LOS of the InSAR reflects only one projection of the true deformation of the landslide to the deformation, the actual deformation condition of the bank slope is difficult to reflect, and if the true three-dimensional deformation field of the ground surface is to be reconstructed, at least three measurement results in different directions or equivalent prior information are required. However, in some current experiments and researches, only one type of SAR data of a monitoring area for rail ascending or rail descending is often used, and in many cases, only monitoring deformation in one direction cannot completely reflect the real deformation condition of a monitoring area, which is easy to cause misjudgment or missing judgment of ground surface deformation monitoring information, and even may cause deviation. If the deformation happens to occur in the sensor heading, it is likely that the deformation is not monitored at all.
To remedy such a drawback, in recent years, researchers at home and abroad have conducted related studies. Herrera and the like find that on a local scale, the fusion of multi-platform InSAR results is helpful for distinguishing different landslide deformation directions, reflecting different deformation rate modes and distinguishing slow landslides caused by natural factors and fast landslides caused by artificial factors; on the basis of the landslide parallel displacement hypothesis, the MT-InSAR deformation monitoring results of ALOS rail ascending and ENVISAT rail descending are fused, and the slope deformation rate of the navicular drainage slope is obtained; HE and the like fuse multi-Aperture InSAR (Multiple-Aperture InSAR, MAI) and SBAS technologies, extract azimuth time sequence deformation of the landslide of the compliant mining area, fuse the azimuth time sequence deformation with LOS (LoS) direction time sequence deformation, and obtain a two-dimensional time sequence deformation field of the landslide; because the Offset-Tracking method (Offset-Tracking) can simultaneously acquire LOS direction deformation and azimuth direction deformation, Raucoules and the like process TerrasAR-X lifting rail data in an Alps region in France by utilizing Offset-Tracking and acquire a three-dimensional deformation result of the La Valette huge landslide in the region. Due to the side-view imaging mode of the SAR satellite and the particularity of landslide deformation, the research on monitoring the three-dimensional deformation of the landslide by utilizing the time series InSAR technology is very limited at present.
The global navigation satellite system (GNSS, GPS/Beidou system) can obtain continuous surface three-dimensional deformation monitoring results, but is limited by ground receiving equipment conditions, often cannot achieve high net distribution density, belongs to discrete point observation, has low spatial resolution and is not enough to meet the requirement of high spatial resolution deformation monitoring, and InSAR can provide planar continuous information of the whole research area; the InSAR is insensitive to deformation information in the horizontal direction, the horizontal deformation precision of a global navigation satellite system (GNSS, GPS/Beidou system) is high, data fusion of the InSAR and the global navigation satellite system can correct errors which are difficult to eliminate of the InSAR data, and effective unification of high time resolution and high plane position precision of the GNSS and high space resolution and elevation deformation precision of an InSAR technology can be achieved. Therefore, how to effectively fuse two complementary earth observation technologies to obtain a three-dimensional deformation field with high spatial-temporal resolution and high precision in a research area, perform rapid and effective deformation monitoring and early warning on the bank slope of the hydropower engineering, provide technical support for guaranteeing the safety of the hydropower engineering and the life safety of people, and have important significance for promoting the progress and application popularization of the bank slope deformation monitoring technology.
Disclosure of Invention
The invention aims to provide a landslide three-dimensional deformation calculation method and system integrating GNSS and lifting rail time sequence InSAR and having high spatial and temporal resolution and high precision.
In order to achieve the above purpose, the invention provides the following technical scheme:
the landslide three-dimensional deformation calculation method fusing the GNSS and the lifting rail time sequence InSAR comprises the following steps:
acquiring a plurality of original images of the ascending and descending radar satellite in a GNSS monitoring period of a target monitoring area;
acquiring the surface deformation of a target monitoring area under the view angles of an ascending and descending rail radar through a time sequence InSAR method;
decomposing the ground surface deformation of the target monitoring area under the visual angles of the rail ascending and descending radar in the three-dimensional direction;
and calculating the three-dimensional deformation of the target monitoring area by fusing an InSAR-GNSS three-dimensional joint solution model, and identifying the landslide of the target monitoring area according to the three-dimensional deformation.
According to the method, the original images of the orbit ascending and descending radar satellites are processed through a time sequence InSAR method, landmark deformation under the view angles of the orbit ascending and descending radar in a target detection area can be obtained, and then a three-dimensional joint solution model can be better constructed in a combined manner with GNSS monitoring through decomposition of the three-dimensional directions of the earth surface deformation under the view angles of the orbit ascending and descending radar, so that the three-dimensional deformation of the target detection area is determined, and landslide identification is carried out on the target area. According to the method, the GNSS monitoring is added while the target monitoring area is monitored in the rail ascending and descending mode, errors which are difficult to eliminate in InSAR data are effectively corrected, the high time resolution and high plane position accuracy of the GNSS and the high space resolution and elevation deformation accuracy of the InSAR technology can be effectively unified, the progress and application and popularization of the bank slope deformation monitoring technology are effectively promoted, and the method has high practicability.
Preferably, in the above technical solution, before the obtaining of the plurality of orbit ascending and descending radar satellite original images within the target monitoring area GNSS monitoring period, the method further includes the following steps:
and deploying a GNSS monitoring device in the target monitoring area to acquire continuous three-dimensional deformation monitoring data in the monitoring period of the target detection area.
Preferably, in the above technical solution, the method includes the steps of calculating the three-dimensional deformation of the target monitoring area by fusing an InSAR-GNSS three-dimensional joint solution model, and identifying the landslide of the target monitoring area according to the three-dimensional deformation, including:
determining the relation between the deformation rate of an interferogram of the lifting rail InSAR to the earth surface and the deformation rate of the interferogram of the lifting rail InSAR in the three-dimensional orthogonal direction according to the SAR image imaging geometric principle;
constructing a three-dimensional joint resolving model based on the deformation rate relation and the GNSS monitoring data;
and calculating three-dimensional deformation based on the three-dimensional joint calculation model, and carrying out early identification on the landslide.
Preferably, in the above technical solution, a relation between a deformation rate of the interferogram of the lifting rail InSAR to the ground surface and a deformation rate of the interferogram of the lifting rail InSAR in the three-dimensional orthogonal direction is as follows:
wherein:
and:
the deformation rates to the ground surface of the interferogram of the ascending orbit InSAR and the deformation rates to the ground surface of the interferogram of the descending orbit InSAR are respectively;
and
unit projection vectors of the rail ascending and descending realizing directions respectively; [ v ] ofe vn vu]TThe deformation rates of the observation points in the three-dimensional orthogonal direction are respectively.
Preferably, in the above technical solution, the calculation formula of the three-dimensional joint solution model is as follows:
written in matrix form as:
L＝AX
wherein:
for the observation value of GPS in three-dimensional direction, A is matrix parameter, X ═ ve vnvu]TThe deformation rates of the observation points in the three-dimensional orthogonal direction are respectively, and L is three-dimensional deformation.
The invention also provides a landslide three-dimensional deformation calculation system integrating the GNSS and the lifting rail time sequence InSAR, which comprises a lifting rail image acquisition device, a lifting rail InSAR data processing device and a three-dimensional combined calculation model device:
the GNSS variable acquisition and lifting orbit image acquisition device is used for acquiring long-term three-dimensional deformation of a GNSS in a target monitoring area and a plurality of lifting orbit radar satellite original images in a monitoring period;
the lifting rail InSAR data processing device is used for acquiring the surface deformation of a target monitoring area under the lifting and lowering rail radar visual angle by a time sequence InSAR method;
the three-dimensional joint solution model device is used for decomposing the ground surface deformation of the target monitoring area under the view angles of the rail ascending and descending radar in the three-dimensional direction;
and calculating the three-dimensional deformation of the target monitoring area by fusing an InSAR-GNSS three-dimensional joint solution model, and identifying the landslide of the target monitoring area according to the three-dimensional deformation.
The system provided by the invention can acquire the deformation of the landmarks in the elevation and the descent radar in the target detection area under the elevation and the descent radar visual angle through the processing of the original images of the elevation and the descent radar satellites by the time sequence InSAR method, and can better jointly construct a three-dimensional joint solution model with GNSS monitoring through the decomposition of the deformation of the ground surface in the elevation and the descent radar visual angle in the three-dimensional direction, so as to determine the three-dimensional deformation of the target detection area and recognize the landslide of the target area. The system adds GNSS monitoring while performing rail ascending and descending monitoring on a target monitoring area, effectively corrects errors which are difficult to eliminate of InSAR data, can effectively unify high time resolution and high plane position precision of the GNSS and high space resolution and elevation deformation precision of InSAR technology, effectively promotes progress and application and popularization of the bank slope deformation monitoring technology, and has high practicability.
Preferably, in the above technical solution, the GNSS monitoring apparatus further includes:
the system is deployed in a target monitoring area and used for acquiring continuous three-dimensional deformation monitoring data in a target detection area monitoring period.
Preferably, in the above technical solution, the three-dimensional joint solution model device is further configured to determine a relationship between a deformation rate of the interferogram of the lifting rail InSAR to the earth surface and a deformation rate of the interferogram of the lifting rail InSAR in a three-dimensional orthogonal direction according to an imaging geometry principle of the SAR image;
constructing a three-dimensional joint resolving model based on the deformation rate relation and the GNSS monitoring data;
and calculating three-dimensional deformation based on the three-dimensional joint calculation model, and carrying out early identification on the landslide.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention and not to limit the invention. In the drawings:
fig. 1 is a schematic flow chart of a landslide three-dimensional deformation calculation method incorporating GNSS and lifting rail timing InSAR according to an embodiment of the present invention;
FIG. 2 is a schematic logic flow diagram of a landslide three-dimensional deformation calculation method fusing a GNSS and a lifting rail time sequence InSAR in FIG. 1;
fig. 3 is a schematic block diagram of a landslide three-dimensional deformation calculation method fusing a GNSS and a lifting rail timing InSAR according to another embodiment of the present invention.
Detailed Description
In order to make the technical problems, technical solutions and advantageous effects to be solved by the present invention more clearly apparent, the present invention is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
It will be understood that when an element is referred to as being "secured to" or "disposed on" another element, it can be directly on the other element or be indirectly on the other element. When an element is referred to as being "connected to" another element, it can be directly connected to the other element or be indirectly connected to the other element.
Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of the present invention, "a plurality" means two or more unless specifically defined otherwise. The meaning of "a number" is one or more unless specifically limited otherwise.
In the description of the present invention, it is to be understood that the terms "upper", "lower", "front", "rear", "left", "right", and the like indicate orientations or positional relationships based on those shown in the drawings, and are only for convenience of description and simplicity of description, but do not indicate or imply that the referred device or element must have a specific orientation, be constructed in a specific orientation, and be operated, and thus, should not be construed as limiting the present invention.
In the description of the present invention, it should be noted that, unless otherwise explicitly specified or limited, the terms "mounted," "connected," and "connected" are to be construed broadly, e.g., as meaning either a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; either directly or indirectly through intervening media, either internally or in any other relationship. The specific meanings of the above terms in the present invention can be understood by those skilled in the art according to specific situations.
As shown in fig. 1 and fig. 2, the landslide three-dimensional deformation calculation method fusing GNSS and lifting rail timing InSAR provided by the invention includes the following steps:
acquiring a plurality of original images of the ascending and descending radar satellite in a GNSS monitoring period of a target monitoring area;
acquiring the surface deformation of a target monitoring area under the view angles of an ascending and descending rail radar through a time sequence InSAR method;
decomposing the ground surface deformation of the target monitoring area under the visual angles of the rail ascending and descending radar in the three-dimensional direction;
and calculating the three-dimensional deformation of the target monitoring area by fusing an InSAR-GNSS three-dimensional joint solution model, and identifying the landslide of the target monitoring area according to the three-dimensional deformation.
According to the method, the original images of the orbit ascending and descending radar satellites are processed through a time sequence InSAR method, landmark deformation under the view angles of the orbit ascending and descending radar in a target detection area can be obtained, and then a three-dimensional joint solution model can be better constructed in a combined manner with GNSS monitoring through decomposition of the three-dimensional directions of the earth surface deformation under the view angles of the orbit ascending and descending radar, so that the three-dimensional deformation of the target detection area is determined, and landslide identification is carried out on the target area. According to the method, the GNSS monitoring is added while the target monitoring area is monitored in the rail ascending and descending mode, errors which are difficult to eliminate in InSAR data are effectively corrected, the high time resolution and high plane position accuracy of the GNSS and the high space resolution and elevation deformation accuracy of the InSAR technology can be effectively unified, the progress and application and popularization of the bank slope deformation monitoring technology are effectively promoted, and the method has high practicability.
As an implementation manner, before the acquiring a plurality of raw images of the orbiting radar satellite within a GNSS monitoring period of a target monitoring area, the method further includes the following steps:
step S05: and deploying a GNSS monitoring device in the target monitoring area to acquire continuous three-dimensional deformation monitoring data in the monitoring period of the target detection area.
Through the acquired continuous three-dimensional deformation monitoring data in the target detection area detection period, data support is provided for the combination of the rail ascending and descending data and the GNSS monitoring data, and the effective execution of the landslide identification method for the target detection area is ensured.
As an implementation mode, the method comprises the following steps of calculating the three-dimensional deformation of a target monitoring area by fusing an InSAR-GNSS three-dimensional joint solution model, and identifying landslide of the target monitoring area according to the three-dimensional deformation:
determining the relation between the deformation rate of an interferogram of the lifting rail InSAR to the earth surface and the deformation rate of the interferogram of the lifting rail InSAR in the three-dimensional orthogonal direction according to the SAR image imaging geometric principle;
constructing a three-dimensional joint resolving model based on the deformation rate relation and the GNSS monitoring data;
and calculating three-dimensional deformation based on the three-dimensional joint calculation model, and carrying out early identification on the landslide.
By determining the relation between the deformation rate of the interferogram of the lifting rail InSAR to the earth surface and the deformation rate of the interferogram of the lifting rail InSAR in the three-dimensional orthogonal direction, data and parameter support is provided for the components of the three-dimensional joint calculation model, the components of the three-dimensional joint calculation model are ensured, and early landslide identification is realized.
As an implementation manner, the deformation rate of the interferogram of the lifting rail InSAR to the ground surface and the deformation rate of the interferogram in the three-dimensional orthogonal direction have the following relationship:
wherein:
and:
the deformation rates to the ground surface of the interferogram of the ascending orbit InSAR and the deformation rates to the ground surface of the interferogram of the descending orbit InSAR are respectively;
and
unit projection vectors of the rail ascending and descending realizing directions respectively; [ v ] ofe vn vu]TThe deformation rates of the observation points in the three-dimensional orthogonal direction are respectively.
As an implementation manner, the calculation formula of the three-dimensional joint solution model is as follows:
written in matrix form as:
L＝AX
wherein:
for the observation value of GPS in three-dimensional direction, A is matrix parameter, X ═ ve vnvu]TThe deformation rates of the observation points in the three-dimensional orthogonal direction are respectively, and L is three-dimensional deformation.
As shown in fig. 3, the invention further provides a landslide three-dimensional deformation calculation system integrating GNSS and lifting rail timing InSAR, a lifting rail image acquisition device, a lifting rail InSAR data processing device, and a three-dimensional joint calculation model device:
the lifting orbit image acquisition device is used for acquiring long-term three-dimensional deformation of a GNSS (global navigation satellite system) in a target monitoring area and a plurality of lifting orbit radar satellite original images in a monitoring period;
the lifting rail InSAR data processing device is used for acquiring the surface deformation of a target monitoring area under the lifting and lowering rail radar visual angle by a time sequence InSAR method;
the three-dimensional joint solution model device is used for decomposing the ground surface deformation of the target monitoring area under the view angles of the rail ascending and descending radar in the three-dimensional direction;
and calculating the three-dimensional deformation of the target monitoring area by fusing an InSAR-GNSS three-dimensional joint solution model, and identifying the landslide of the target monitoring area according to the three-dimensional deformation.
The system provided by the invention can acquire the deformation of the landmarks in the elevation and the descent radar in the target detection area under the elevation and the descent radar visual angle through the processing of the original images of the elevation and the descent radar satellites by the time sequence InSAR method, and can better jointly construct a three-dimensional joint solution model with GNSS monitoring through the decomposition of the deformation of the ground surface in the elevation and the descent radar visual angle in the three-dimensional direction, so as to determine the three-dimensional deformation of the target detection area and recognize the landslide of the target area. The system adds GNSS monitoring while performing rail ascending and descending monitoring on a target monitoring area, effectively corrects errors which are difficult to eliminate of InSAR data, can effectively unify high time resolution and high plane position precision of the GNSS and high space resolution and elevation deformation precision of InSAR technology, effectively promotes progress and application and popularization of the bank slope deformation monitoring technology, and has high practicability.
As an implementation manner, the GNSS monitoring apparatus is further included:
the system is deployed in a target monitoring area and used for acquiring continuous three-dimensional deformation monitoring data in a target detection area monitoring period.
Through the acquired continuous three-dimensional deformation monitoring data in the target detection area detection period, data support is provided for the combination of the rail ascending and descending data and the GNSS monitoring data, and the effective execution of the landslide identification method for the target detection area is ensured.
As an implementable embodiment, the three-dimensional joint solution model device is further configured to determine a relation between a deformation rate to the earth surface of an interferogram of the lifting rail InSAR and a deformation rate in a three-dimensional orthogonal direction according to an SAR image imaging geometry principle;
constructing a three-dimensional joint resolving model based on the deformation rate relation and the GNSS monitoring data;
and calculating three-dimensional deformation based on the three-dimensional joint calculation model, and carrying out early identification on the landslide.
By determining the relation between the deformation rate of the interferogram of the lifting rail InSAR to the earth surface and the deformation rate of the interferogram of the lifting rail InSAR in the three-dimensional orthogonal direction, data and parameter support is provided for the components of the three-dimensional joint calculation model, the components of the three-dimensional joint calculation model are ensured, and early landslide identification is realized.
In the foregoing description of embodiments, the particular features, structures, materials, or characteristics may be combined in any suitable manner in any one or more embodiments or examples.
The above description is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present invention, and all the changes or substitutions should be covered within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the appended claims.
Claims (8)
1. A landslide three-dimensional deformation calculation method fusing GNSS and lifting rail time sequence InSAR is characterized by comprising the following steps:
acquiring a plurality of original images of the ascending and descending radar satellite in a GNSS monitoring period of a target monitoring area;
acquiring the surface deformation of a target monitoring area under the view angles of an ascending and descending rail radar through a time sequence InSAR method;
decomposing the ground surface deformation of the target monitoring area under the visual angles of the rail ascending and descending radar in the three-dimensional direction;
and calculating the three-dimensional deformation of the target monitoring area by fusing an InSAR-GNSS three-dimensional joint solution model, and identifying the landslide of the target monitoring area according to the three-dimensional deformation.
2. The landslide three-dimensional deformation calculation method integrating the GNSS and the ascending and descending orbit timing InSAR according to claim 1, further comprising, before the obtaining of the plurality of ascending and descending orbit radar satellite raw images within the GNSS monitoring period of the target monitoring area, the steps of:
and deploying a GNSS monitoring device in the target monitoring area to acquire continuous three-dimensional deformation monitoring data in the monitoring period of the target detection area.
3. The landslide three-dimensional deformation calculation method fusing the GNSS and the lifting rail time sequence InSAR as claimed in claim 1, wherein the three-dimensional deformation of the target monitoring area is calculated by fusing an InSAR-GNSS three-dimensional joint calculation model, and landslide identification is performed on the target monitoring area according to the three-dimensional deformation, comprising the following steps:
determining the relation between the deformation rate of an interferogram of the lifting rail InSAR to the earth surface and the deformation rate of the interferogram of the lifting rail InSAR in the three-dimensional orthogonal direction according to the SAR image imaging geometric principle;
constructing a three-dimensional joint resolving model based on the deformation rate relation and the GNSS monitoring data;
and calculating three-dimensional deformation based on the three-dimensional joint calculation model, and carrying out early identification on the landslide.
4. The landslide three-dimensional deformation calculation method integrating the GNSS and the lifting rail time sequence InSAR is characterized in that the relation between the deformation rate of the interferogram of the lifting rail InSAR to the earth surface and the deformation rate of the interferogram of the lifting rail InSAR in the three-dimensional orthogonal direction is as follows:
wherein:
and:
the deformation rates to the ground surface of the interferogram of the ascending orbit InSAR and the deformation rates to the ground surface of the interferogram of the descending orbit InSAR are respectively;
and
unit projection vectors of the rail ascending and descending realizing directions respectively; [ v ] ofevnvu]TThe deformation rates of the observation points in the three-dimensional orthogonal direction are respectively.
5. The landslide three-dimensional deformation calculation method integrating GNSS and lifting rail time sequence InSAR as claimed in claim 4, wherein the calculation formula of the three-dimensional joint calculation model is as follows:
written in matrix form as:
L＝AX
wherein:
for the observation value of GPS in three-dimensional direction, A is matrix parameter, X ═ vevnvu]TThe deformation rates of the observation points in the three-dimensional orthogonal direction are respectively, and L is three-dimensional deformation.
6. A landslide three-dimensional deformation calculation system integrating GNSS and lifting rail time sequence InSAR is characterized by comprising a lifting rail image acquisition device, a lifting rail InSAR data processing device and a three-dimensional joint calculation model device;
the GNSS variable acquisition and lifting orbit image acquisition device is used for acquiring long-term three-dimensional deformation of a GNSS in a target monitoring area and a plurality of lifting orbit radar satellite original images in a monitoring period;
the lifting rail InSAR data processing device is used for acquiring the surface deformation of a target monitoring area under the lifting and lowering rail radar visual angle by a time sequence InSAR method;
the three-dimensional joint solution model device is used for decomposing the ground surface deformation of the target monitoring area under the view angles of the rail ascending and descending radar in the three-dimensional direction;
and calculating the three-dimensional deformation of the target monitoring area by fusing an InSAR-GNSS three-dimensional joint solution model, and identifying the landslide of the target monitoring area according to the three-dimensional deformation.
7. The landslide three-dimensional deformation calculation system integrating the GNSS and the lifting rail time sequence InSAR as claimed in claim 6, further comprising a GNSS monitoring device:
the system is deployed in a target monitoring area and used for acquiring continuous three-dimensional deformation monitoring data in a target detection area monitoring period.
8. The landslide three-dimensional deformation calculation system integrating the GNSS and the lifting rail time sequence InSAR as claimed in claim 6, wherein the three-dimensional joint calculation model device is further used for determining a relation between a deformation rate to the earth surface and a deformation rate in a three-dimensional orthogonal direction of an interferogram of the lifting rail InSAR according to an SAR image imaging geometry principle;
constructing a three-dimensional joint resolving model based on the deformation rate relation and the GNSS monitoring data;
and calculating three-dimensional deformation based on the three-dimensional joint calculation model, and carrying out early identification on the landslide.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011355738.1A CN112540369A (en) | 2020-11-27 | 2020-11-27 | Landslide three-dimensional deformation resolving method and system integrating GNSS and lifting rail time sequence InSAR |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011355738.1A CN112540369A (en) | 2020-11-27 | 2020-11-27 | Landslide three-dimensional deformation resolving method and system integrating GNSS and lifting rail time sequence InSAR |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112540369A true CN112540369A (en) | 2021-03-23 |
Family
ID=75015166
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011355738.1A Pending CN112540369A (en) | 2020-11-27 | 2020-11-27 | Landslide three-dimensional deformation resolving method and system integrating GNSS and lifting rail time sequence InSAR |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112540369A (en) |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113091598A (en) * | 2021-04-06 | 2021-07-09 | 中国矿业大学 | Method for defining stability grade range of goaf building site by InSAR |
| CN113280750A (en) * | 2021-06-09 | 2021-08-20 | 武汉大学 | Three-dimensional deformation monitoring method and device |
| CN113281742A (en) * | 2021-06-02 | 2021-08-20 | 西南交通大学 | SAR landslide early warning method based on landslide deformation information and meteorological data |
| CN113532374A (en) * | 2021-06-30 | 2021-10-22 | 中铁第四勘察设计院集团有限公司 | Ground settlement detection method, device, equipment and storage medium |
| CN113624122A (en) * | 2021-08-10 | 2021-11-09 | 中咨数据有限公司 | Bridge deformation monitoring method fusing GNSS data and InSAR technology |
| CN113740844A (en) * | 2021-09-09 | 2021-12-03 | 甘肃中星鸿图科技有限公司 | Dam body three-dimensional deformation monitoring-oriented two-foundation radar combined observation method |
| CN113762623A (en) * | 2021-09-09 | 2021-12-07 | 甘肃中星鸿图科技有限公司 | Landslide direction and trend identification and prediction method and device and electronic equipment |
| CN113790667A (en) * | 2021-11-18 | 2021-12-14 | 中大检测(湖南)股份有限公司 | Dam deformation detection method based on radar |
| CN113848551A (en) * | 2021-09-24 | 2021-12-28 | 成都理工大学 | Landslide depth inversion method using InSAR lifting rail deformation data |
| CN113960596A (en) * | 2021-10-20 | 2022-01-21 | 苏州深蓝空间遥感技术有限公司 | Landslide three-dimensional deformation monitoring method based on Beidou and PS-InSAR |
| CN114089335A (en) * | 2021-11-16 | 2022-02-25 | 安徽理工大学 | Mountain area mining subsidence three-dimensional deformation extraction method based on monorail InSAR |
| CN114863643A (en) * | 2022-04-25 | 2022-08-05 | 浙江科技学院 | Early warning method for rainfall type landslide based on GIS area |
| CN115469308A (en) * | 2022-08-25 | 2022-12-13 | 中国地震局地质研究所 | Multi-track InSAR inter-seismic deformation rate field splicing method, device, equipment and medium |
| CN116659429A (en) * | 2023-08-01 | 2023-08-29 | 齐鲁空天信息研究院 | Multi-source data high-precision time sequence earth surface three-dimensional deformation resolving method and system |
| CN116819525A (en) * | 2023-05-19 | 2023-09-29 | 东北大学 | InSAR landslide monitoring and evaluating method and system for high steep side slope |
| CN117109426A (en) * | 2023-08-28 | 2023-11-24 | 兰州交通大学 | Three-dimensional deformation field modeling method fusing GNSS/InSAR observation data |
| CN117168373A (en) * | 2023-07-20 | 2023-12-05 | 中国卫通集团股份有限公司 | Reservoir dam body deformation monitoring system based on satellite leads to remote integration |
| CN117553737A (en) * | 2023-11-13 | 2024-02-13 | 内蒙古高新科技控股有限责任公司 | Slope safety monitoring method and device based on lead telemetry information |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101770027A (en) * | 2010-02-05 | 2010-07-07 | 河海大学 | Ground surface three-dimensional deformation monitoring method based on InSAR and GPS data fusion |
| US20160047940A1 (en) * | 2014-08-15 | 2016-02-18 | California Institute Of Technology | Systems and methods for advanced rapid imaging and analysis for earthquakes |
| CN110058236A (en) * | 2019-05-21 | 2019-07-26 | 中南大学 | It is a kind of towards three-dimensional Ground Deformation estimation InSAR and GNSS determine Quan Fangfa |
| CN111458709A (en) * | 2020-06-08 | 2020-07-28 | 河南大学 | Satellite-borne radar wide-area earth surface two-dimensional deformation field monitoring method and device |
| GB202012738D0 (en) * | 2019-08-15 | 2020-09-30 | China Institute Of Water Resources And Res | System and method for monitoring deformation of dam slope |
-
2020
- 2020-11-27 CN CN202011355738.1A patent/CN112540369A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101770027A (en) * | 2010-02-05 | 2010-07-07 | 河海大学 | Ground surface three-dimensional deformation monitoring method based on InSAR and GPS data fusion |
| US20160047940A1 (en) * | 2014-08-15 | 2016-02-18 | California Institute Of Technology | Systems and methods for advanced rapid imaging and analysis for earthquakes |
| CN110058236A (en) * | 2019-05-21 | 2019-07-26 | 中南大学 | It is a kind of towards three-dimensional Ground Deformation estimation InSAR and GNSS determine Quan Fangfa |
| GB202012738D0 (en) * | 2019-08-15 | 2020-09-30 | China Institute Of Water Resources And Res | System and method for monitoring deformation of dam slope |
| CN111458709A (en) * | 2020-06-08 | 2020-07-28 | 河南大学 | Satellite-borne radar wide-area earth surface two-dimensional deformation field monitoring method and device |
Non-Patent Citations (4)
| Title |
|---|
| 曹海坤: "GPS、InSAR数据联合解算地表三维形变场", 中国优秀硕士学位论文全文数据库基础科学辑, no. 02, pages 10 - 12 * |
| 曹海坤;赵丽华;张勤;瞿伟;聂建亮;: "利用附加系统误差参数的升降轨InSAR-GPS数据融合方法建立三维形变场", 武汉大学学报(信息科学版), no. 09 * |
| 杨晓杰等: "矿山工程地质学", 中国矿业大学出版社, pages: 269 - 271 * |
| 范青松;汤翠莲;陈于;张晓东;: "GPS与InSAR技术在滑坡监测中的应用研究", 测绘科学, no. 05 * |
Cited By (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113091598A (en) * | 2021-04-06 | 2021-07-09 | 中国矿业大学 | Method for defining stability grade range of goaf building site by InSAR |
| CN113281742A (en) * | 2021-06-02 | 2021-08-20 | 西南交通大学 | SAR landslide early warning method based on landslide deformation information and meteorological data |
| CN113280750A (en) * | 2021-06-09 | 2021-08-20 | 武汉大学 | Three-dimensional deformation monitoring method and device |
| CN113532374B (en) * | 2021-06-30 | 2022-09-09 | 中铁第四勘察设计院集团有限公司 | Ground settlement detection method, device, equipment and storage medium |
| CN113532374A (en) * | 2021-06-30 | 2021-10-22 | 中铁第四勘察设计院集团有限公司 | Ground settlement detection method, device, equipment and storage medium |
| CN113624122A (en) * | 2021-08-10 | 2021-11-09 | 中咨数据有限公司 | Bridge deformation monitoring method fusing GNSS data and InSAR technology |
| CN113762623A (en) * | 2021-09-09 | 2021-12-07 | 甘肃中星鸿图科技有限公司 | Landslide direction and trend identification and prediction method and device and electronic equipment |
| CN113762623B (en) * | 2021-09-09 | 2024-04-05 | 雷添杰 | Landslide direction and trend identification and prediction method and device and electronic equipment |
| CN113740844B (en) * | 2021-09-09 | 2024-04-02 | 雷添杰 | Dam three-dimensional deformation monitoring-oriented two-foundation radar combined observation method |
| CN113740844A (en) * | 2021-09-09 | 2021-12-03 | 甘肃中星鸿图科技有限公司 | Dam body three-dimensional deformation monitoring-oriented two-foundation radar combined observation method |
| CN113848551A (en) * | 2021-09-24 | 2021-12-28 | 成都理工大学 | Landslide depth inversion method using InSAR lifting rail deformation data |
| CN113848551B (en) * | 2021-09-24 | 2023-03-21 | 成都理工大学 | Landslide depth inversion method using InSAR lifting rail deformation data |
| CN113960596A (en) * | 2021-10-20 | 2022-01-21 | 苏州深蓝空间遥感技术有限公司 | Landslide three-dimensional deformation monitoring method based on Beidou and PS-InSAR |
| CN114089335B (en) * | 2021-11-16 | 2022-09-06 | 安徽理工大学 | Mountain area mining subsidence three-dimensional deformation extraction method based on monorail InSAR |
| CN114089335A (en) * | 2021-11-16 | 2022-02-25 | 安徽理工大学 | Mountain area mining subsidence three-dimensional deformation extraction method based on monorail InSAR |
| CN113790667A (en) * | 2021-11-18 | 2021-12-14 | 中大检测(湖南)股份有限公司 | Dam deformation detection method based on radar |
| CN114863643A (en) * | 2022-04-25 | 2022-08-05 | 浙江科技学院 | Early warning method for rainfall type landslide based on GIS area |
| CN114863643B (en) * | 2022-04-25 | 2023-09-22 | 浙江科技学院 | Rainfall landslide early warning method based on GIS region |
| CN115469308A (en) * | 2022-08-25 | 2022-12-13 | 中国地震局地质研究所 | Multi-track InSAR inter-seismic deformation rate field splicing method, device, equipment and medium |
| CN116819525B (en) * | 2023-05-19 | 2024-02-09 | 东北大学 | InSAR landslide monitoring and evaluating method and system for high steep side slope |
| CN116819525A (en) * | 2023-05-19 | 2023-09-29 | 东北大学 | InSAR landslide monitoring and evaluating method and system for high steep side slope |
| CN117168373A (en) * | 2023-07-20 | 2023-12-05 | 中国卫通集团股份有限公司 | Reservoir dam body deformation monitoring system based on satellite leads to remote integration |
| CN117168373B (en) * | 2023-07-20 | 2024-07-09 | 中国卫通集团股份有限公司 | Reservoir dam body deformation monitoring system based on satellite leads to remote integration |
| CN116659429A (en) * | 2023-08-01 | 2023-08-29 | 齐鲁空天信息研究院 | Multi-source data high-precision time sequence earth surface three-dimensional deformation resolving method and system |
| CN117109426B (en) * | 2023-08-28 | 2024-03-22 | 兰州交通大学 | Three-dimensional deformation field modeling method fusing GNSS/InSAR observation data |
| CN117109426A (en) * | 2023-08-28 | 2023-11-24 | 兰州交通大学 | Three-dimensional deformation field modeling method fusing GNSS/InSAR observation data |
| CN117553737A (en) * | 2023-11-13 | 2024-02-13 | 内蒙古高新科技控股有限责任公司 | Slope safety monitoring method and device based on lead telemetry information |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112540369A (en) | Landslide three-dimensional deformation resolving method and system integrating GNSS and lifting rail time sequence InSAR | |
| WO2022213673A1 (en) | Method for extracting three-dimensional surface deformation by combining unmanned aerial vehicle doms and satellite-borne sar images | |
| Brasington et al. | Monitoring and modelling morphological change in a braided gravel‐bed river using high resolution GPS‐based survey | |
| Rosenau et al. | A processing system to monitor Greenland outlet glacier velocity variations at decadal and seasonal time scales utilizing the Landsat imagery | |
| Brenner | Extraction of features from mobile laser scanning data for future driver assistance systems | |
| Corsini et al. | Integrating airborne and multi‐temporal long‐range terrestrial laser scanning with total station measurements for mapping and monitoring a compound slow moving rock slide | |
| Liu et al. | Detection of crustal movement from TerraSAR-X intensity images for the 2011 Tohoku, Japan earthquake | |
| CN110232683A (en) | A kind of landslide detection method based on unmanned plane point cloud | |
| Dematteis et al. | Integration of robotic total station and digital image correlation to assess the three-dimensional surface kinematics of a landslide | |
| Aati et al. | A new approach for 2-D and 3-D precise measurements of ground deformation from optimized registration and correlation of optical images and ICA-based filtering of image geometry artifacts | |
| Roncella et al. | A fixed terrestrial photogrammetric system for landslide monitoring | |
| Bertin et al. | Geomorphic change detection using cost-effective structure-from-motion photogrammetry: Evaluation of direct georeferencing from consumer-grade UAS at Orewa beach (New Zealand) | |
| CN117968631A (en) | Pavement subsidence detection method based on unmanned aerial vehicle DOM and satellite-borne SAR image | |
| Dreier et al. | The potential of UAV-based laser scanning for deformation monitoring. Case study on a water dam | |
| Qiao et al. | Assessment of geo-positioning capability of high resolution satellite imagery for densely populated high buildings in metropolitan areas | |
| Hormes et al. | Innovative methods to monitor rock and mountain slope deformation | |
| Previtali et al. | Accurate 3D surface measurement of mountain slopes through a fully automated image-based technique | |
| Meisina et al. | Choice of surveying methods for landslides monitoring | |
| Sefercik et al. | Country-scale discontinuity analysis of AW3D30 and SRTM Global DEMS: case study in Turkey | |
| CN114185048A (en) | Method, system and storage medium for extracting landslide displacement vector by foundation InSAR | |
| Tamimi et al. | Performance Assessment of a Mini Mobile Mapping System: Iphone 14 pro Installed on a e-Scooter | |
| Tampubolon et al. | Investigations on InSAR data processing standard for volcano island monitoring in Indonesia | |
| Lee et al. | Automatic building reconstruction with satellite images and digital maps | |
| Li et al. | Terrestrial mobile mapping towards real-time geospatial data collection | |
| Mancini et al. | Contribution of geomatics engineering and VGI within the landslide risk assessment procedures |
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 |
Application publication date: 20210323 |