CN116243269A - Post-earthquake landslide hazard monitoring method and device based on Insar data - Google Patents

Post-earthquake landslide hazard monitoring method and device based on Insar data Download PDF

Info

Publication number
CN116243269A
CN116243269A CN202310499274.9A CN202310499274A CN116243269A CN 116243269 A CN116243269 A CN 116243269A CN 202310499274 A CN202310499274 A CN 202310499274A CN 116243269 A CN116243269 A CN 116243269A
Authority
CN
China
Prior art keywords
earthquake
unit
slope
remote sensing
permanent
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.)
Granted
Application number
CN202310499274.9A
Other languages
Chinese (zh)
Other versions
CN116243269B (en
Inventor
向杰
胡辉
夏鹏
颜秋宇
邹谢华
朱清
严煦
奚晓谦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Mining News
Nanjing Aerospace Hongtu Information Technology Co ltd
China University of Geosciences
Original Assignee
China Mining News
Nanjing Aerospace Hongtu Information Technology Co ltd
China University of Geosciences
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Mining News, Nanjing Aerospace Hongtu Information Technology Co ltd, China University of Geosciences filed Critical China Mining News
Priority to CN202310499274.9A priority Critical patent/CN116243269B/en
Publication of CN116243269A publication Critical patent/CN116243269A/en
Application granted granted Critical
Publication of CN116243269B publication Critical patent/CN116243269B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/50—Context or environment of the image
    • G06V20/52—Surveillance or monitoring of activities, e.g. for recognising suspicious objects
    • 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
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/40—Extraction of image or video features
    • G06V10/62—Extraction of image or video features relating to a temporal dimension, e.g. time-based feature extraction; Pattern tracking
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/10—Terrestrial scenes
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Electromagnetism (AREA)
  • Alarm Systems (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a post-earthquake landslide hazard monitoring method and device based on Insar data, which relate to the technical field of landslide hazard monitoring and comprise the following steps: remote sensing data of the area to be monitored before and after the earthquake is obtained, interference processing and permanent scatterer interference measurement processing are carried out on the remote sensing data, and target data of the permanent scatterer corresponding to the remote sensing data are obtained; determining the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data based on the target data; determining a slope unit in the area to be monitored based on DEM data of the area to be monitored, wherein the area of the slope unit is in a preset area range, and the circular variance of the slope direction of the grid unit in the slope unit is smaller than a preset threshold value; based on the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data, the state of the slope unit before and after the earthquake is determined, and the technical problem that the landslide activity is difficult to monitor after the earthquake in the existing earthquake landslide disaster monitoring method is solved.

Description

Post-earthquake landslide hazard monitoring method and device based on Insar data
Technical Field
The invention relates to the technical field of landslide hazard monitoring, in particular to a post-earthquake landslide hazard monitoring method and device based on Insar data.
Background
Many studies have shown that intense seismic vibrations not only induce co-seismic landslides, but also amplify post-seismic landslide activity, possibly due to reduced shear strength of the slope material or disturbing the morphology of the slope. Thus, an increase in landslide sensitivity is observed during post-earthquake periods, which in effect means that rainfall events that do not induce any landslide before the earthquake may induce landslide after the earthquake, which may be expressed as an increase in the overall landslide sensitivity level for a given landscape, a concept defined as the earthquake carry-over effect. However, the current technology and patent do not relate to the whole scheme of post-earthquake landslide activity detection, and related researches only relate to the terrain deformation detection of insar, and the technical scheme of post-earthquake landslide activity detection cannot be systematically and completely given.
An effective solution to the above-mentioned problems has not been proposed yet.
Disclosure of Invention
In view of the above, the invention aims to provide a post-earthquake landslide hazard monitoring method and device based on Insar data, so as to solve the technical problem that the existing earthquake landslide hazard monitoring method is difficult to monitor post-earthquake landslide activity.
In a first aspect, an embodiment of the present invention provides a post-earthquake landslide hazard monitoring method based on Insar data, including: remote sensing data of a region to be monitored before and after an earthquake are obtained, interference processing and permanent scatterer interference measurement processing are carried out on the remote sensing data, and target data of a permanent scatterer corresponding to the remote sensing data are obtained, wherein the target data comprise: a deformation time series and an average annual line-of-sight velocity map; and determining the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data based on the target data, wherein the annual LOS deformation speed comprises: the LOS deformation speed before earthquake and the LOS deformation speed after earthquake; determining a slope unit in the area to be monitored based on DEM data of the area to be monitored, wherein the area of the slope unit is in a preset area range, and the circular variance of the slope direction of the grid unit in the slope unit is smaller than a preset threshold value; and determining the state of the slope unit before and after the earthquake based on the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data, wherein the state comprises the following steps: steady state and active state.
Further, performing interference processing and permanent scatterer interference measurement processing on the remote sensing data to obtain target data of a permanent scatterer corresponding to the remote sensing data, including: performing interference processing on the remote sensing data by using a preset open source tool to obtain an interference pattern; and carrying out permanent scatterer interferometry processing on the interferogram to obtain target data of the permanent scatterer corresponding to the remote sensing data.
Further, determining a ramp unit in the area to be monitored based on DEM data of the area to be monitored, including: determining the gradient and the slope direction of each grid unit in the area to be monitored based on the DEM data of the area to be monitored; and determining the slope units in the area to be monitored based on the gradient and the slope direction of each grid unit in the area to be monitored.
Further, determining the state of the slope unit before and after the earthquake based on the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data includes: determining the state of the slope unit before the earthquake based on the LOS deformation speed of the previous earthquake; and determining the state of the slope unit after the earthquake based on the LOS deformation speed of the later year of the earthquake.
Further, determining a state of the ramp unit before the earthquake based on the annual LOS deformation speed before the earthquake includes: determining the number of first permanent scatterers contained in the slope unit based on the previous-year LOS deformation speed of the earthquake, wherein the first permanent scatterers are permanent scatterers with the previous-year LOS deformation speed of the earthquake being greater than or equal to a first preset threshold value; if the number of the first permanent scatterers contained in the slope unit is larger than a first preset number, the state of the slope unit before an earthquake is an active state; and if the number of the first permanent scatterers contained in the slope unit is smaller than or equal to the first preset number, the state of the slope unit before earthquake is a stable state.
Further, determining a state of the ramp unit after the earthquake based on the post-earthquake annual LOS deformation speed includes: determining the number of second permanent scatterers contained in the slope unit based on the post-earthquake annual LOS deformation speed, wherein the second permanent scatterers are permanent scatterers with the post-earthquake annual LOS deformation speed greater than or equal to a second preset threshold value; if the number of the second permanent scatterers contained in the slope unit is larger than a second preset number, the state of the slope unit after the earthquake is an active state; and if the number of the second permanent scatterers contained in the slope unit is smaller than or equal to the second preset number, the state of the slope unit after the earthquake is a stable state.
In a second aspect, an embodiment of the present invention further provides a post-earthquake landslide hazard monitoring device based on the instr data, including: the device comprises an acquisition unit, a processing unit and a processing unit, wherein the acquisition unit is used for acquiring remote sensing data of an area to be monitored before and after an earthquake, carrying out interference processing and permanent scatterer interference measurement processing on the remote sensing data to obtain target data of a permanent scatterer corresponding to the remote sensing data, wherein the target data comprises: a deformation time series and an average annual line-of-sight velocity map; the first determining unit is configured to determine an annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data based on the target data, where the annual LOS deformation speed includes: the LOS deformation speed before earthquake and the LOS deformation speed after earthquake; the second determining unit is used for determining a slope unit in the area to be monitored based on the DEM data of the area to be monitored, wherein the area of the slope unit is in a preset area range, and the circular variance of the slope direction of the grid unit in the slope unit is smaller than a preset threshold value; the third determining unit is configured to determine, based on an annual LOS deformation speed of a permanent scatterer corresponding to the remote sensing data, a state of the ramp unit before and after an earthquake, where the state includes: steady state and active state.
Further, the acquiring unit is configured to: performing interference processing on the remote sensing data by using a preset open source tool to obtain an interference pattern; and carrying out permanent scatterer interferometry processing on the interferogram to obtain target data of the permanent scatterer corresponding to the remote sensing data.
In a third aspect, an embodiment of the present invention further provides an electronic device, including a memory and a processor, where the memory is configured to store a program for supporting the processor to execute the method described in the first aspect, and the processor is configured to execute the program stored in the memory.
In a fourth aspect, embodiments of the present invention also provide a computer-readable storage medium having a computer program stored thereon.
In the embodiment of the invention, remote sensing data of an area to be monitored before and after an earthquake are acquired, interference processing and permanent scatterer interference measurement processing are carried out on the remote sensing data, and target data of a permanent scatterer corresponding to the remote sensing data are obtained, wherein the target data comprise: a deformation time series and an average annual line-of-sight velocity map; and determining the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data based on the target data, wherein the annual LOS deformation speed comprises: the LOS deformation speed before earthquake and the LOS deformation speed after earthquake; determining a slope unit in the area to be monitored based on DEM data of the area to be monitored, wherein the area of the slope unit is in a preset area range, and the circular variance of the slope direction of the grid unit in the slope unit is smaller than a preset threshold value; and determining the state of the slope unit before and after the earthquake based on the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data, wherein the state comprises the following steps: the method has the advantages that the aim of monitoring the post-earthquake landslide activity is fulfilled due to the stable state and the active state, the technical problem that the post-earthquake landslide activity is difficult to monitor by the existing earthquake landslide disaster monitoring method is solved, and the technical effect of providing support for the post-earthquake disaster and risk assessment is achieved.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
In order to make the above objects, features and advantages of the present invention more comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the description of the embodiments or the prior art will be briefly described, and it is obvious that the drawings in the description below are some embodiments of the present invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
Fig. 1 is a flowchart of a post-earthquake landslide hazard monitoring method based on Insar data provided by an embodiment of the invention;
fig. 2 is a schematic diagram of a post-earthquake landslide hazard monitoring device based on Insar data according to an embodiment of the present invention;
fig. 3 is a schematic diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Embodiment one:
in accordance with an embodiment of the present invention, there is provided an embodiment of a post-earthquake landslide hazard monitoring method based on Insar data, it being noted that the steps shown in the flowchart of the figures may be performed in a computer system such as a set of computer-executable instructions, and, although a logical order is shown in the flowchart, in some cases, the steps shown or described may be performed in an order other than that shown herein.
Fig. 1 is a flowchart of a post-earthquake landslide hazard monitoring method based on Insar data according to an embodiment of the invention, and as shown in fig. 1, the method includes the steps of:
step S102, remote sensing data of a region to be monitored before and after an earthquake is obtained, interference processing and permanent scatterer interference measurement processing are carried out on the remote sensing data, and target data of a permanent scatterer corresponding to the remote sensing data are obtained, wherein the target data comprise: a deformation time series and an average annual line-of-sight velocity map;
it should be noted that the remote sensing data in the embodiment of the present invention is Sentinel-1 SLC data.
Step S104, determining the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data based on the target data, wherein the annual LOS deformation speed comprises the following steps: the LOS deformation speed before earthquake and the LOS deformation speed after earthquake;
the deformation result of the radar line-of-sight direction (LOS direction) obtained by the InSAR technique processing is calculated year by year to obtain the annual LOS deformation speed.
Step S106, determining a slope unit in the area to be monitored based on DEM data of the area to be monitored, wherein the area of the slope unit is in a preset area range, and the circular variance of the slope direction of the grid unit in the slope unit is smaller than a preset threshold;
it should be noted that DEM data of the area to be monitored is DEM data before an earthquake.
Step S108, determining the state of the slope unit before and after an earthquake based on the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data, wherein the state comprises: steady state and active state.
In the embodiment of the invention, remote sensing data of an area to be monitored before and after an earthquake are acquired, interference processing and permanent scatterer interference measurement processing are carried out on the remote sensing data, and target data of a permanent scatterer corresponding to the remote sensing data are obtained, wherein the target data comprise: a deformation time series and an average annual line-of-sight velocity map; and determining the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data based on the target data, wherein the annual LOS deformation speed comprises: the LOS deformation speed before earthquake and the LOS deformation speed after earthquake; determining a slope unit in the area to be monitored based on DEM data of the area to be monitored, wherein the area of the slope unit is in a preset area range, and the circular variance of the slope direction of the grid unit in the slope unit is smaller than a preset threshold value; and determining the state of the slope unit before and after the earthquake based on the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data, wherein the state comprises the following steps: the method has the advantages that the aim of monitoring the post-earthquake landslide activity is fulfilled due to the stable state and the active state, the technical problem that the post-earthquake landslide activity is difficult to monitor by the existing earthquake landslide disaster monitoring method is solved, and the technical effect of providing support for the post-earthquake disaster and risk assessment is achieved.
In the embodiment of the present invention, step S102 includes the following steps:
performing interference processing on the remote sensing data by using a preset open source tool to obtain an interference pattern;
and carrying out permanent scatterer interferometry processing on the interferogram to obtain target data of the permanent scatterer corresponding to the remote sensing data.
In the embodiment of the invention, the SNAP2StaMPS is preferably adopted as the open source tool.
In the embodiment of the present invention, step S106 includes the following steps:
determining the gradient and the slope direction of each grid unit in the area to be monitored based on the DEM data of the area to be monitored;
and determining the slope units in the area to be monitored based on the gradient and the slope direction of each grid unit in the area to be monitored.
In the embodiment of the invention, after DEM data of an area to be monitored is acquired, a maximum area threshold value m, a minimum area threshold value a and a circular variance c of a slope unit are set as control parameters of an extraction result, and the smaller the circular variance of the slope direction is, the closer the slope directions of all grid units in the slope unit are. When the round variance of the grid unit in the small area meets the set threshold value c, the small area has uniform slope direction as a slope unit.
In the embodiment of the present invention, step S108 includes the steps of:
determining the state of the slope unit before the earthquake based on the LOS deformation speed of the previous earthquake;
and determining the state of the slope unit after the earthquake based on the LOS deformation speed of the later year of the earthquake.
Specifically, based on the previous-year LOS deformation speed of the earthquake, determining the number of first permanent scatterers contained in the slope unit, wherein the first permanent scatterers are permanent scatterers with the previous-year LOS deformation speed of the earthquake being greater than or equal to a first preset threshold value;
if the number of the first permanent scatterers contained in the slope unit is larger than a first preset number, the state of the slope unit before an earthquake is an active state;
and if the number of the first permanent scatterers contained in the slope unit is smaller than or equal to the first preset number, the state of the slope unit before earthquake is a stable state.
Determining the number of second permanent scatterers contained in the slope unit based on the post-earthquake annual LOS deformation speed, wherein the second permanent scatterers are permanent scatterers with the post-earthquake annual LOS deformation speed greater than or equal to a second preset threshold value;
if the number of the second permanent scatterers contained in the slope unit is larger than a second preset number, the state of the slope unit after the earthquake is an active state;
and if the number of the second permanent scatterers contained in the slope unit is smaller than or equal to the second preset number, the state of the slope unit after the earthquake is a stable state.
It should be noted that, the first preset threshold value and the second preset threshold value may be the same or different, and the first preset number and the second preset number may be the same or different, and specific numerical values are set by the staff according to actual situations.
In the embodiment of the invention, if the pre-earthquake state and the post-earthquake state of the slope unit are both stable states, the slope corresponding to the slope unit is stable before and after an earthquake, and the slope is not affected by the earthquake.
If the pre-earthquake state and the post-earthquake state of the slope unit are both active states, the slope corresponding to the slope unit is stable before and after the earthquake and is not influenced by dynamic changes.
If the pre-earthquake state and the post-earthquake state of the slope unit are both active states, the state of the slope corresponding to the slope unit is changed from stable to active due to the influence of earthquake.
If the pre-earthquake state and the post-earthquake state of the slope unit are both the active state and the steady state, the state of the slope corresponding to the slope unit is changed from active to steady due to the influence of the earthquake.
In the embodiment of the invention, the whole area to be monitored can be analyzed integrally by the slope unit to aggregate the permanent scatterers instead of using the method of standard pixel density collection.
In the embodiment of the invention, the hillside evolution process in the post-earthquake period can be systematically and consistently checked by the method, and the method is applied to other areas subjected to earthquake images. Knowing the evolution of hills affected by strong earthquakes helps to better assess post-earthquake disasters and risks, as well as planning management and mitigation measures.
Embodiment two:
the embodiment of the invention also provides a post-earthquake landslide hazard monitoring device based on the Insar data, which is used for executing the post-earthquake landslide hazard monitoring method based on the Insar data provided by the embodiment of the invention, and the following is a specific introduction of the post-earthquake landslide hazard monitoring device based on the Insar data provided by the embodiment of the invention.
As shown in fig. 2, fig. 2 is a schematic diagram of the post-earthquake landslide hazard monitoring device based on the instr data, where the post-earthquake landslide hazard monitoring device based on the instr data includes:
the acquiring unit 10 is configured to acquire remote sensing data of an area to be monitored before and after an earthquake, perform interference processing and permanent scatterer interferometry processing on the remote sensing data, and obtain target data of a permanent scatterer corresponding to the remote sensing data, where the target data includes: a deformation time series and an average annual line-of-sight velocity map;
a first determining unit 20, configured to determine an annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data based on the target data, where the annual LOS deformation speed includes: the LOS deformation speed before earthquake and the LOS deformation speed after earthquake;
a second determining unit 30, configured to determine a ramp unit in the area to be monitored based on DEM data of the area to be monitored, where an area of the ramp unit is within a preset area range and a circular variance of a slope direction of a grid unit in the ramp unit is smaller than a preset threshold;
a third determining unit 40, configured to determine, based on the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data, a state of the ramp unit before and after the earthquake, where the state includes: steady state and active state.
In the embodiment of the invention, remote sensing data of an area to be monitored before and after an earthquake are acquired, interference processing and permanent scatterer interference measurement processing are carried out on the remote sensing data, and target data of a permanent scatterer corresponding to the remote sensing data are obtained, wherein the target data comprise: a deformation time series and an average annual line-of-sight velocity map; and determining the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data based on the target data, wherein the annual LOS deformation speed comprises: the LOS deformation speed before earthquake and the LOS deformation speed after earthquake; determining a slope unit in the area to be monitored based on DEM data of the area to be monitored, wherein the area of the slope unit is in a preset area range, and the circular variance of the slope direction of the grid unit in the slope unit is smaller than a preset threshold value; and determining the state of the slope unit before and after the earthquake based on the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data, wherein the state comprises the following steps: the method has the advantages that the aim of monitoring the post-earthquake landslide activity is fulfilled due to the stable state and the active state, the technical problem that the post-earthquake landslide activity is difficult to monitor by the existing earthquake landslide disaster monitoring method is solved, and the technical effect of providing support for the post-earthquake disaster and risk assessment is achieved.
Embodiment III:
an embodiment of the present invention further provides an electronic device, including a memory and a processor, where the memory is configured to store a program that supports the processor to execute the method described in the first embodiment, and the processor is configured to execute the program stored in the memory.
Referring to fig. 3, an embodiment of the present invention further provides an electronic device 100, including: a processor 50, a memory 51, a bus 52 and a communication interface 53, the processor 50, the communication interface 53 and the memory 51 being connected by the bus 52; the processor 50 is arranged to execute executable modules, such as computer programs, stored in the memory 51.
The memory 51 may include a high-speed random access memory (RAM, random Access Memory), and may further include a non-volatile memory (non-volatile memory), such as at least one magnetic disk memory. The communication connection between the system network element and at least one other network element is achieved via at least one communication interface 53 (which may be wired or wireless), and the internet, wide area network, local network, metropolitan area network, etc. may be used.
Bus 52 may be an ISA bus, a PCI bus, an EISA bus, or the like. The buses may be classified as address buses, data buses, control buses, etc. For ease of illustration, only one bi-directional arrow is shown in FIG. 3, but not only one bus or type of bus.
The memory 51 is configured to store a program, and the processor 50 executes the program after receiving an execution instruction, and the method executed by the apparatus for flow defining disclosed in any of the foregoing embodiments of the present invention may be applied to the processor 50 or implemented by the processor 50.
The processor 50 may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuitry in hardware in the processor 50 or by instructions in the form of software. The processor 50 may be a general-purpose processor, including a Central Processing Unit (CPU), a Network Processor (NP), and the like; but may also be a digital signal processor (Digital Signal Processing, DSP for short), application specific integrated circuit (Application Specific Integrated Circuit, ASIC for short), off-the-shelf programmable gate array (Field-Programmable Gate Array, FPGA for short), or other programmable logic device, discrete gate or transistor logic device, discrete hardware components. The disclosed methods, steps, and logic blocks in the embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present invention may be embodied directly in the execution of a hardware decoding processor, or in the execution of a combination of hardware and software modules in a decoding processor. The software modules may be located in a random access memory, flash memory, read only memory, programmable read only memory, or electrically erasable programmable memory, registers, etc. as well known in the art. The storage medium is located in a memory 51 and the processor 50 reads the information in the memory 51 and in combination with its hardware performs the steps of the above method.
Embodiment four:
the embodiment of the invention also provides a computer readable storage medium, and a computer program is stored on the computer readable storage medium, and when the computer program is executed by a processor, the steps of the method in the first embodiment are executed.
In addition, in the description of embodiments of the present invention, unless explicitly stated and limited otherwise, the terms "mounted," "connected," and "connected" are to be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the above terms in the present invention will be understood in specific cases by those of ordinary skill in the art.
In the description of the present invention, it should be noted that the directions or positional relationships indicated by the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc. are based on the directions or positional relationships shown in the drawings, are merely for convenience of describing the present invention and simplifying the description, and do not indicate or imply that the devices or elements referred to must have a specific orientation, be configured and operated in a specific orientation, and thus should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
In the several embodiments provided in this application, it should be understood that the disclosed systems, devices, and methods may be implemented in other manners. The above-described apparatus embodiments are merely illustrative, for example, the division of the units is merely a logical function division, and there may be other manners of division in actual implementation, and for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be through some communication interface, device or unit indirect coupling or communication connection, which may be in electrical, mechanical or other form.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit.
Finally, it should be noted that: the above examples are only specific embodiments of the present invention, and are not intended to limit the scope of the present invention, but it should be understood by those skilled in the art that the present invention is not limited thereto, and that the present invention is described in detail with reference to the foregoing examples: any person skilled in the art may modify or easily conceive of the technical solution described in the foregoing embodiments, or perform equivalent substitution of some of the technical features, while remaining within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention, and are intended to be included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (10)

1. The post-earthquake landslide hazard monitoring method based on Insar data is characterized by comprising the following steps of:
remote sensing data of a region to be monitored before and after an earthquake are obtained, interference processing and permanent scatterer interference measurement processing are carried out on the remote sensing data, and target data of a permanent scatterer corresponding to the remote sensing data are obtained, wherein the target data comprise: a deformation time series and an average annual line-of-sight velocity map;
and determining the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data based on the target data, wherein the annual LOS deformation speed comprises: the LOS deformation speed before earthquake and the LOS deformation speed after earthquake;
determining a slope unit in the area to be monitored based on DEM data of the area to be monitored, wherein the area of the slope unit is in a preset area range, and the circular variance of the slope direction of the grid unit in the slope unit is smaller than a preset threshold value;
and determining the state of the slope unit before and after the earthquake based on the annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data, wherein the state comprises the following steps: steady state and active state.
2. The method of claim 1, wherein performing interferometry and permanent scatterer interferometry on the remote sensing data to obtain target data for a permanent scatterer corresponding to the remote sensing data comprises:
performing interference processing on the remote sensing data by using a preset open source tool to obtain an interference pattern;
and carrying out permanent scatterer interferometry processing on the interferogram to obtain target data of the permanent scatterer corresponding to the remote sensing data.
3. The method of claim 1, wherein determining a ramp unit in the area to be monitored based on DEM data of the area to be monitored comprises:
determining the gradient and the slope direction of each grid unit in the area to be monitored based on the DEM data of the area to be monitored;
and determining the slope units in the area to be monitored based on the gradient and the slope direction of each grid unit in the area to be monitored.
4. The method of claim 1, wherein determining the state of the ramp unit before and after the earthquake based on the annual LOS deformation velocity of the permanent scatterer corresponding to the remote sensing data comprises:
determining the state of the slope unit before the earthquake based on the LOS deformation speed of the previous earthquake;
and determining the state of the slope unit after the earthquake based on the LOS deformation speed of the later year of the earthquake.
5. The method of claim 4, wherein determining the pre-seismic state of the ramp unit based on the pre-seismic LOS deformation velocity comprises:
determining the number of first permanent scatterers contained in the slope unit based on the previous-year LOS deformation speed of the earthquake, wherein the first permanent scatterers are permanent scatterers with the previous-year LOS deformation speed of the earthquake being greater than or equal to a first preset threshold value;
if the number of the first permanent scatterers contained in the slope unit is larger than a first preset number, the state of the slope unit before an earthquake is an active state;
and if the number of the first permanent scatterers contained in the slope unit is smaller than or equal to the first preset number, the state of the slope unit before earthquake is a stable state.
6. The method of claim 4, wherein determining the post-earthquake state of the ramp unit based on the post-earthquake LOS deformation velocity comprises:
determining the number of second permanent scatterers contained in the slope unit based on the post-earthquake annual LOS deformation speed, wherein the second permanent scatterers are permanent scatterers with the post-earthquake annual LOS deformation speed greater than or equal to a second preset threshold value;
if the number of the second permanent scatterers contained in the slope unit is larger than a second preset number, the state of the slope unit after the earthquake is an active state;
and if the number of the second permanent scatterers contained in the slope unit is smaller than or equal to the second preset number, the state of the slope unit after the earthquake is a stable state.
7. Post-earthquake landslide hazard monitoring device based on Insar data, characterized by comprising:
the device comprises an acquisition unit, a processing unit and a processing unit, wherein the acquisition unit is used for acquiring remote sensing data of an area to be monitored before and after an earthquake, carrying out interference processing and permanent scatterer interference measurement processing on the remote sensing data to obtain target data of a permanent scatterer corresponding to the remote sensing data, wherein the target data comprises: a deformation time series and an average annual line-of-sight velocity map;
the first determining unit is configured to determine an annual LOS deformation speed of the permanent scatterer corresponding to the remote sensing data based on the target data, where the annual LOS deformation speed includes: the LOS deformation speed before earthquake and the LOS deformation speed after earthquake;
the second determining unit is used for determining a slope unit in the area to be monitored based on the DEM data of the area to be monitored, wherein the area of the slope unit is in a preset area range, and the circular variance of the slope direction of the grid unit in the slope unit is smaller than a preset threshold value;
the third determining unit is configured to determine, based on an annual LOS deformation speed of a permanent scatterer corresponding to the remote sensing data, a state of the ramp unit before and after an earthquake, where the state includes: steady state and active state.
8. The apparatus of claim 7, wherein the acquisition unit is configured to:
performing interference processing on the remote sensing data by using a preset open source tool to obtain an interference pattern;
and carrying out permanent scatterer interferometry processing on the interferogram to obtain target data of the permanent scatterer corresponding to the remote sensing data.
9. An electronic device comprising a memory for storing a program supporting the processor to perform the method of any one of claims 1 to 6, and a processor configured to execute the program stored in the memory.
10. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, performs the steps of the method according to any of the preceding claims 1 to 6.
CN202310499274.9A 2023-05-06 2023-05-06 Post-earthquake landslide hazard monitoring method and device based on Insar data Active CN116243269B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310499274.9A CN116243269B (en) 2023-05-06 2023-05-06 Post-earthquake landslide hazard monitoring method and device based on Insar data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310499274.9A CN116243269B (en) 2023-05-06 2023-05-06 Post-earthquake landslide hazard monitoring method and device based on Insar data

Publications (2)

Publication Number Publication Date
CN116243269A true CN116243269A (en) 2023-06-09
CN116243269B CN116243269B (en) 2023-07-28

Family

ID=86628093

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310499274.9A Active CN116243269B (en) 2023-05-06 2023-05-06 Post-earthquake landslide hazard monitoring method and device based on Insar data

Country Status (1)

Country Link
CN (1) CN116243269B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021008282A1 (en) * 2019-07-12 2021-01-21 清华大学 Seismic landslide quick report analysis method and apparatus based on actually-measured seismic motion
CN113064188A (en) * 2020-08-07 2021-07-02 国网浙江省电力有限公司 Transformer substation geological deformation monitoring method based on SAR satellite and Beidou satellite
CN113096005A (en) * 2021-04-06 2021-07-09 中国科学院生态环境研究中心 Radar time sequence differential interferometry method for monitoring mountain body lifting speed at present
CN113192086A (en) * 2021-05-11 2021-07-30 中国自然资源航空物探遥感中心 Generation method of geological disaster hidden danger deformation intensity distribution diagram and storage medium
CN113343563A (en) * 2021-05-27 2021-09-03 中交第二公路勘察设计研究院有限公司 Landslide susceptibility evaluation method based on automatic sample selection and surface deformation rate
CN114036794A (en) * 2021-11-12 2022-02-11 中国海洋大学 Method for analyzing stability of seabed slope after earthquake
CN114493319A (en) * 2022-02-11 2022-05-13 中国科学院沈阳应用生态研究所 Cross-time-scale combined ancient landslide resurgence risk grading evaluation method and device
CN114694030A (en) * 2022-04-21 2022-07-01 中煤航测遥感集团有限公司 Landslide detection method, device, equipment and storage medium
CN114791273A (en) * 2022-06-24 2022-07-26 中国铁道科学研究院集团有限公司铁道建筑研究所 InSAR deformation monitoring result interpretation method for landslide

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021008282A1 (en) * 2019-07-12 2021-01-21 清华大学 Seismic landslide quick report analysis method and apparatus based on actually-measured seismic motion
CN113064188A (en) * 2020-08-07 2021-07-02 国网浙江省电力有限公司 Transformer substation geological deformation monitoring method based on SAR satellite and Beidou satellite
CN113096005A (en) * 2021-04-06 2021-07-09 中国科学院生态环境研究中心 Radar time sequence differential interferometry method for monitoring mountain body lifting speed at present
CN113192086A (en) * 2021-05-11 2021-07-30 中国自然资源航空物探遥感中心 Generation method of geological disaster hidden danger deformation intensity distribution diagram and storage medium
CN113343563A (en) * 2021-05-27 2021-09-03 中交第二公路勘察设计研究院有限公司 Landslide susceptibility evaluation method based on automatic sample selection and surface deformation rate
CN114036794A (en) * 2021-11-12 2022-02-11 中国海洋大学 Method for analyzing stability of seabed slope after earthquake
CN114493319A (en) * 2022-02-11 2022-05-13 中国科学院沈阳应用生态研究所 Cross-time-scale combined ancient landslide resurgence risk grading evaluation method and device
CN114694030A (en) * 2022-04-21 2022-07-01 中煤航测遥感集团有限公司 Landslide detection method, device, equipment and storage medium
CN114791273A (en) * 2022-06-24 2022-07-26 中国铁道科学研究院集团有限公司铁道建筑研究所 InSAR deformation monitoring result interpretation method for landslide

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
周振凯: "２ ０ １ ５ 年１ ２ 月７ 日帕米尔高原Ｍ ｗ ７ ．１ 地震同震变形场 Ｄ － Ｉ ｎ Ｓ Ａ Ｒ 观测及构造稳定性分析", 工程地质学报 *
张诗茄;蒋建军;缪亚敏;白世彪;: "基于SBAS技术的岷江流域潜在滑坡识别", 山地学报, no. 01 *
彭令;牛瑞卿;: "三峡库区白家包滑坡变形特征与影响因素分析", 中国地质灾害与防治学报, no. 04 *
王凯 等: "斜坡单元提取方法研究进展和展望", 长江科学院院报, vol. 37, no. 6 *

Also Published As

Publication number Publication date
CN116243269B (en) 2023-07-28

Similar Documents

Publication Publication Date Title
CN109413565B (en) Method and device for determining service cell of coverage scene and electronic equipment
Corsini et al. Use of ROC curves for early warning of landslide displacement rates in response to precipitation (Piagneto landslide, Northern Apennines, Italy)
Sedan et al. Armagedom—a tool for seismic risk assessment illustrated with applications
CN109326087B (en) Urban waterlogging early warning method and device based on drainage pipe network monitoring
CN109669798B (en) Crash analysis method, crash analysis device, electronic equipment and storage medium
CN110149654B (en) Method and device for determining faults of base station antenna feeder system
CN114611941A (en) Landslide disaster evaluation method and system considering space dangerousness
CN116243269B (en) Post-earthquake landslide hazard monitoring method and device based on Insar data
CN112751910A (en) Information collection method and device
US20100174490A1 (en) Seismic Detection in Electricity Meters
CN115953034A (en) Method and device for evaluating and zoning risk of gale disaster, electronic equipment and medium
Boken et al. Agricultural water use estimation using geospatial modeling and a geographic information system
CN110661529A (en) Method and device for generating step amplitude sequence
CN116362624A (en) Building earthquake resistance evaluation method, device, terminal and storage medium
CN111143500A (en) Visualized area calculation method, terminal, control device and storage medium
RU2464594C2 (en) Method of estimating main characteristics of anticipated strong tsunamigenic earthquake and system for realising said method
CN116909907A (en) Memory leakage detection method and device, electronic equipment and readable storage medium
CN115331215B (en) Three-dimensional identification and matching method and device for drought event
CN108287890B (en) Data management method and device
CN109612670B (en) Protective net monitoring method, system and terminal equipment
JP6453784B2 (en) Reliability evaluation apparatus, reliability evaluation method, and program
CN112949442B (en) Abnormal event pre-recognition method and device, electronic equipment and monitoring system
CN114332766A (en) Dangerous behavior early warning method and device based on action rule and application thereof
CN112665730A (en) Method, device, equipment and storage medium for detecting pre-earthquake temperature anomaly
Oynakov et al. Spatial variation of precursory seismic quiescence observed before Earthquake from 01.04. 2010 in the region of crete

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
GR01 Patent grant
GR01 Patent grant