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

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

A post-earthquake landslide disaster monitoring method and device based on Insar data

technical field

The invention relates to the technical field of landslide disaster monitoring, in particular to a post-earthquake landslide disaster monitoring method and device based on Insar data.

Background technique

Many studies have shown that strong earthquake shaking not only triggers co-seismic landslides but also amplifies post-earthquake landslide activity, which may be due to the decrease in shear strength of slope materials or disturbed slope morphology. Thus, the observed increased landslide susceptibility in the post-earthquake period actually means that rainfall events that would not have triggered any landslides before the earthquake may have triggered landslides after the earthquake, which can be expressed as the overall landslide susceptibility of a given landscape The level of increase, this concept is defined as the earthquake legacy effect. However, the current technology and patents have not yet involved the overall scheme of post-earthquake landslide activity detection. The relevant research only involves insar's terrain deformation detection, and has not systematically and completely given the technical solution for post-earthquake landslide activity detection.

For the above problems, no effective solution has been proposed yet.

Contents of the invention

In view of this, the object of the present invention is to provide a post-earthquake landslide disaster monitoring method and device based on Insar data, to alleviate the technical problem that the existing earthquake landslide disaster monitoring method is difficult to monitor post-earthquake landslide activities.

In the first aspect, an embodiment of the present invention provides a post-earthquake landslide disaster monitoring method based on Insar data, including: acquiring remote sensing data of the area to be monitored before and after the earthquake, performing interference processing and permanent scatterer interferometry on the remote sensing data processing to obtain the target data of the permanent scatterer corresponding to the remote sensing data, wherein the target data includes: deformation time series and average annual line-of-sight velocity map; based on the target data, determine the target data of the permanent scatterer corresponding to the remote sensing data The annual LOS deformation velocity, wherein, the annual LOS deformation velocity includes: the annual LOS deformation velocity before the earthquake and the annual LOS deformation velocity after the earthquake; based on the DEM data of the area to be monitored, the slope unit in the area to be monitored is determined, Wherein, the area of the slope unit is within a preset area range and the circular variance of the slope direction of the grid unit in the slope unit is less than a preset threshold; based on the annual LOS deformation rate of the permanent scatterer corresponding to the remote sensing data , to determine the state of the slope unit before and after the earthquake, wherein the state includes: a stable state and an active state.

Further, performing interference processing and permanent scatterer interferometry processing on the remote sensing data to obtain the target data of the permanent scatterer corresponding to the remote sensing data includes: using a preset open source tool to perform interference processing on the remote sensing data, Obtaining an interferogram; performing permanent scatterer interferometry processing on the interferogram to obtain target data of the permanent scatterer corresponding to the remote sensing data.

Further, determining the slope units in the area to be monitored based on the DEM data of the area to be monitored includes: determining the slope units of each grid unit in the area to be monitored based on the DEM data of the area to be monitored. Slope and aspect: Determine the slope unit in the area to be monitored based on the slope and aspect of each grid unit in the area to be monitored.

Further, based on the annual LOS deformation rate of the permanent scatterer corresponding to the remote sensing data, determining the state of the slope unit before and after the earthquake includes: determining the state of the slope unit before and after the earthquake based on the LOS deformation rate of the year before the earthquake. The state before the earthquake; based on the annual LOS deformation velocity after the earthquake, the state of the slope unit after the earthquake is determined.

Further, determining the state of the slope unit before the earthquake based on the LOS deformation rate in the year before the earthquake includes: determining the number of first permanent scatterers contained in the slope unit based on the LOS deformation rate in the year before the earthquake , wherein, the first permanent scatterer is a permanent scatterer whose LOS deformation rate in the year before the earthquake is greater than or equal to a first preset threshold; if the number of first permanent scatterers included in the slope unit is greater than the first preset number, Then the state of the slope unit before the earthquake is an active state; if the number of first permanent scatterers contained in the slope unit is less than or equal to the first preset number, the state of the slope unit before the earthquake is steady state.

Further, determining the state of the slope unit after the earthquake based on the annual LOS deformation rate after the earthquake includes: determining the second permanent scattering contained in the slope unit based on the annual LOS deformation rate after the earthquake. The number of scatterers, wherein, the second permanent scatterer is a permanent scatterer whose annual LOS deformation rate after the earthquake is greater than or equal to a second preset threshold; if the number of second permanent scatterers contained in the slope unit is greater than the second preset threshold If the number is set, the state of the slope unit after the earthquake is an active state; if the number of second permanent scatterers contained in the slope unit is less than or equal to the second preset number, the slope unit will be active after the earthquake. state is a stable state.

In the second aspect, the embodiment of the present invention also provides a post-earthquake landslide disaster monitoring device based on Insar data, including: an acquisition unit, configured to acquire remote sensing data of the area to be monitored before and after the earthquake, and perform interference processing on the remote sensing data and permanent scatterer interferometry processing to obtain the target data of the permanent scatterer corresponding to the remote sensing data, wherein the target data includes: a deformation time series and an average annual line-of-sight velocity map; a first determination unit is configured to The target data is used to determine the annual LOS deformation velocity of the permanent scatterer corresponding to the remote sensing data, wherein the annual LOS deformation velocity includes: the annual LOS deformation velocity before the earthquake and the annual LOS deformation velocity after the earthquake; The DEM data of the area to be monitored determines the slope unit in the area to be monitored, wherein the area of the slope unit is within the preset area range and the circular variance of the slope direction of the grid unit in the slope unit less than the preset threshold; the third determination unit is configured to determine the state of the slope unit before and after the earthquake based on the annual LOS deformation rate of the permanent scatterer corresponding to the remote sensing data, wherein the state includes: a stable state and active state.

Further, the acquisition unit is configured to: use a preset open source tool to perform interference processing on the remote sensing data to obtain an interferogram; perform permanent scatterer interferometry processing on the interferogram to obtain the remote sensing data corresponding to Target data for permanent scatterers.

In a third aspect, an embodiment of the present invention also provides an electronic device, including a memory and a processor, the memory is used to store a program that supports the processor to execute the method described in the first aspect above, and the processor is configured to for executing programs stored in the memory.

In a fourth aspect, the embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored.

In the embodiment of the present invention, by acquiring the remote sensing data of the area to be monitored before and after the earthquake, performing interference processing and permanent scatterer interferometry processing on the remote sensing data, the target data of the permanent scatterers corresponding to the remote sensing data are obtained, wherein , the target data includes: deformation time series and average annual line-of-sight velocity map; based on the target data, the annual LOS deformation velocity of the permanent scatterer corresponding to the remote sensing data is determined, wherein the annual LOS deformation velocity includes: the year before the earthquake LOS deformation velocity and annual LOS deformation velocity after the earthquake; based on the DEM data of the area to be monitored, determine the slope unit in the area to be monitored, wherein the area of the slope unit is within a preset area range and the The circular variance of the slope aspect of the grid unit in the slope unit is less than a preset threshold; based on the annual LOS deformation rate of the permanent scatterer corresponding to the remote sensing data, the state of the slope unit before and after the earthquake is determined, wherein the The state includes: stable state and active state, which achieves the purpose of monitoring landslide activity after the earthquake, and then solves the technical problem that the existing earthquake landslide hazard monitoring method is difficult to monitor landslide activity after the earthquake, thus realizing the post-earthquake disaster and landslide monitoring. The assessment of risk provides support for the technical effects.

Additional features and advantages of the invention will be set forth in the description which follows, and in part will be apparent 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-mentioned objects, features and advantages of the present invention more comprehensible, preferred embodiments will be described in detail below together with the accompanying drawings.

Description of drawings

In order to more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the following will briefly introduce the drawings that need to be used in the description of the specific embodiments or the prior art. Obviously, the accompanying drawings in the following description The drawings show some implementations of the present invention, and those skilled in the art can also obtain other drawings based on these drawings without creative work.

Fig. 1 is the flow chart of a kind of post-earthquake landslide disaster monitoring method based on Insar data that the embodiment of the present invention provides;

Fig. 2 is a schematic diagram of a post-earthquake landslide disaster monitoring device based on Insar data provided by an embodiment of the present invention;

Fig. 3 is a schematic diagram of an electronic device provided by an embodiment of the present invention.

Detailed ways

In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below in conjunction with the accompanying drawings. Obviously, the described embodiments are part of the embodiments of the present invention, not all of them. the embodiment. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts belong to the protection scope of the present invention.

Embodiment one:

According to an embodiment of the present invention, an embodiment of a post-earthquake landslide disaster monitoring method based on Insar data is provided. system, and, although a logical order is shown in the flowcharts, in some cases the steps shown or described may be performed in an order different from that shown or described herein.

Fig. 1 is the flow chart of a kind of post-earthquake landslide hazard monitoring method based on Insar data according to an embodiment of the present invention, as shown in Fig. 1, the method comprises the following steps:

Step S102, obtaining the remote sensing data of the area to be monitored before and after the earthquake, performing interference processing and permanent scatterer interferometry processing on the remote sensing data, and obtaining the target data of the permanent scatterer corresponding to the remote sensing data, wherein the target data Includes: deformation time series and mean annual line-of-sight velocity maps;

It should be noted that the remote sensing data in the embodiment of the present invention is Sentinel-1 SLC data.

Step S104, based on the target data, determine the annual LOS deformation velocity of the permanent scatterer corresponding to the remote sensing data, wherein the annual LOS deformation velocity includes: the annual LOS deformation velocity before the earthquake and the annual LOS deformation velocity after the earthquake;

It should be noted that the deformation results in the radar line-of-sight direction (LOS direction) are processed by InSAR technology, and the deformation results are calculated on an annual basis to obtain the annual LOS deformation speed.

Step S106, based on the DEM data of the area to be monitored, determine the slope unit in the area to be monitored, wherein the area of the slope unit is within the preset area range and the grid unit in the slope unit The circular variance of aspect is less than the preset threshold;

It should be noted that the DEM data of the area to be monitored is the DEM data before the earthquake.

Step S108, based on the annual LOS deformation velocity of the permanent scatterer corresponding to the remote sensing data, determine the state of the slope unit before and after the earthquake, wherein the state includes: a stable state and an active state.

In the embodiment of the present invention, by acquiring the remote sensing data of the area to be monitored before and after the earthquake, performing interference processing and permanent scatterer interferometry processing on the remote sensing data, the target data of the permanent scatterers corresponding to the remote sensing data are obtained, wherein , the target data includes: deformation time series and average annual line-of-sight velocity map; based on the target data, the annual LOS deformation velocity of the permanent scatterer corresponding to the remote sensing data is determined, wherein the annual LOS deformation velocity includes: the year before the earthquake LOS deformation velocity and annual LOS deformation velocity after the earthquake; based on the DEM data of the area to be monitored, determine the slope unit in the area to be monitored, wherein the area of the slope unit is within a preset area range and the The circular variance of the slope aspect of the grid unit in the slope unit is less than a preset threshold; based on the annual LOS deformation rate of the permanent scatterer corresponding to the remote sensing data, the state of the slope unit before and after the earthquake is determined, wherein the The state includes: stable state and active state, which achieves the purpose of monitoring landslide activity after the earthquake, and then solves the technical problem that the existing earthquake landslide hazard monitoring method is difficult to monitor landslide activity after the earthquake, thus realizing the post-earthquake disaster and landslide monitoring. The assessment of risk provides support for the technical effects.

In the embodiment of the present invention, step S102 includes the following steps:

Using preset open source tools to perform interference processing on the remote sensing data to obtain an interferogram;

Perform permanent scatterer interferometry processing on the interferogram to obtain target data of permanent scatterers corresponding to the remote sensing data.

In the embodiment of the present invention, preferably, the open source tool adopts SNAP2StaMPS.

In the embodiment of the present invention, step S106 includes the following steps:

Based on the DEM data of the area to be monitored, determine the slope and aspect of each grid unit in the area to be monitored;

A slope unit in the area to be monitored is determined based on the slope and aspect of each grid unit in the area to be monitored.

In the embodiment of the present invention, after obtaining the DEM data of the area to be monitored, the maximum area threshold m, the minimum area threshold a, and the circular variance c of the slope unit are set as the control parameters of the extraction results. The smaller the circular variance of the slope , indicating that the slope aspects of all grid cells inside the slope cell are closer. When the circular variance of the grid unit in the small area meets the set threshold c, it means that the small area has a uniform slope direction, which is regarded as a slope unit.

In the embodiment of the present invention, step S108 includes the following steps:

Determine the state of the slope unit before the earthquake based on the deformation velocity of the LOS in the year before the earthquake;

Based on the annual LOS deformation rate after the earthquake, the state of the slope unit after the earthquake is determined.

Specifically, based on the LOS deformation rate in the year before the earthquake, the number of first permanent scatterers contained in the slope unit is determined, wherein the first permanent scatterer is that the LOS deformation rate in the year before the earthquake is greater than or equal to the first preset Threshold permanent scatterers;

If the number of first permanent scatterers contained in the slope unit is greater than the first preset number, the state of the slope unit before the earthquake is an active state;

If the number of first permanent scatterers included in the slope unit is less than or equal to the first preset number, the state of the slope unit before the earthquake is a stable state.

Based on the post-earthquake annual LOS deformation rate, determine the number of second permanent scatterers contained in the slope unit, wherein the second permanent scatterer is that the post-earthquake annual LOS deformation rate is greater than or equal to a second preset threshold permanent scatterers;

If the number of second permanent scatterers included in the slope unit is greater than a second preset number, the state of the slope unit after the earthquake is an active state;

If the number of second permanent scatterers included in the slope unit is less 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 above-mentioned first preset threshold and the second preset threshold may be the same or different, and the above-mentioned first preset number and the second preset number may be the same or different. The situation sets itself.

In the embodiment of the present invention, if both the pre-earthquake state and the post-earthquake state of the slope unit are stable, it means that the slope corresponding to the slope unit is stable before and after the earthquake and has not been affected by the earthquake.

If both the pre-earthquake state and the post-earthquake state of the slope unit are active, it means that the slope corresponding to the slope unit is stable before and after the earthquake and is not affected by dynamic changes.

If the pre-earthquake state of the slope unit is stable and the post-earthquake state is active, it means that the state of the slope corresponding to the slope unit changes from stable to active due to the impact of the earthquake.

If the pre-earthquake state of the slope unit is active and the post-earthquake state is stable, it means that the state of the slope corresponding to the slope unit has changed from active to stable due to the impact of the earthquake.

In the embodiment of the present invention, the permanent scatterers are aggregated by the slope unit instead of the standard pixel density aggregation method, so that the whole area to be monitored can be analyzed holistically.

In the embodiment of the present invention, the above method can be used to systematically and consistently examine the hillslope evolution process in the post-earthquake period, and it can be applied to other areas affected by seismic images. Understanding the evolution of mountain slopes affected by strong earthquakes helps to better assess post-earthquake hazards and risks, and to plan management and mitigation measures.

Embodiment two:

The embodiment of the present invention also provides a post-earthquake landslide disaster monitoring device based on Insar data, which is used to implement the post-earthquake landslide disaster monitoring method based on Insar data provided by the above-mentioned content of the embodiment of the present invention. The following is the implementation of the present invention The specific introduction of the post-earthquake landslide hazard monitoring device based on Insar data is provided in the example.

As shown in Figure 2, Figure 2 is a schematic diagram of the above-mentioned post-earthquake landslide disaster monitoring device based on Insar data. The post-earthquake landslide disaster monitoring device based on Insar data includes:

The acquiring unit 10 is configured to acquire the remote sensing data of the area to be monitored before and after the earthquake, perform interference processing and permanent scatterer interferometry processing on the remote sensing data, and obtain the target data of the permanent scatterer corresponding to the remote sensing data, wherein the The target data include: deformation time series and mean annual line-of-sight velocity map;

The first determination unit 20 is configured to determine the annual LOS deformation rate of the permanent scatterer corresponding to the remote sensing data based on the target data, wherein the annual LOS deformation rate includes: the annual LOS deformation rate before the earthquake and the annual LOS deformation rate after the earthquake speed;

The second determination unit 30 is configured to determine the 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 within a preset area range and within the slope unit The circular variance of the aspect of the raster cells is less than the preset threshold;

The third determination unit 40 is configured to determine the state of the slope unit before and after the earthquake based on the annual LOS deformation rate of the permanent scatterer corresponding to the remote sensing data, wherein the state includes: a stable state and an active state.

In the embodiment of the present invention, by acquiring the remote sensing data of the area to be monitored before and after the earthquake, performing interference processing and permanent scatterer interferometry processing on the remote sensing data, the target data of the permanent scatterers corresponding to the remote sensing data are obtained, wherein , the target data includes: deformation time series and average annual line-of-sight velocity map; based on the target data, the annual LOS deformation velocity of the permanent scatterer corresponding to the remote sensing data is determined, wherein the annual LOS deformation velocity includes: the year before the earthquake LOS deformation velocity and annual LOS deformation velocity after the earthquake; based on the DEM data of the area to be monitored, determine the slope unit in the area to be monitored, wherein the area of the slope unit is within a preset area range and the The circular variance of the slope aspect of the grid unit in the slope unit is less than a preset threshold; based on the annual LOS deformation rate of the permanent scatterer corresponding to the remote sensing data, the state of the slope unit before and after the earthquake is determined, wherein the The state includes: stable state and active state, which achieves the purpose of monitoring landslide activity after the earthquake, and then solves the technical problem that the existing earthquake landslide hazard monitoring method is difficult to monitor landslide activity after the earthquake, thus realizing the post-earthquake disaster and landslide monitoring. The assessment of risk provides support for the technical effects.

Embodiment three:

An embodiment of the present invention also provides an electronic device, including a memory and a processor, the memory is used to store a program that supports the processor to execute the method described in the first embodiment above, and the processor is configured to execute the program stored in memory.

Referring to Fig. 3, the embodiment of the present invention also provides an electronic device 100, comprising: 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 are connected through the bus 52; processing The processor 50 is used to execute executable modules stored in the memory 51, such as computer programs.

Wherein, the memory 51 may include a high-speed random access memory (RAM, Random Access Memory), and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory. The communication connection between the system network element and at least one other network element is realized through at least one communication interface 53 (which may be wired or wireless), and the Internet, wide area network, local network, metropolitan area network, etc. can be used.

The bus 52 can be an ISA bus, a PCI bus or an EISA bus, etc. The bus can be divided into address bus, data bus, control bus and so on. For ease of representation, only one double-headed arrow is used in FIG. 3 , but it does not mean that there is only one bus or one type of bus.

Wherein, the memory 51 is used to store the program, and the processor 50 executes the program after receiving the execution instruction, and the method performed by the flow process definition device disclosed in any of the embodiments of the present invention described above can be applied to processing In the device 50, or implemented by the processor 50.

The processor 50 may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by an integrated logic circuit of hardware in the processor 50 or instructions in the form of software. Above-mentioned processor 50 can be general-purpose processor, comprises central processing unit (Central Processing Unit, be called for short CPU), network processor (Network Processor, be called for short NP) etc.; Circuit (Application Specific Integrated Circuit, referred to as ASIC), off-the-shelf programmable gate array (Field-Programmable Gate Array, referred to as FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components. Various methods, steps and logic block diagrams disclosed in the embodiments of the present invention may be implemented or executed. A general-purpose processor may be a microprocessor, or the processor may be any conventional processor, or the like. The steps of the methods disclosed in the embodiments of the present invention may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module can be located in a mature storage medium in the field such as random access memory, flash memory, read-only memory, programmable read-only memory or electrically erasable programmable memory, register. The storage medium is located in the memory 51, and the processor 50 reads the information in the memory 51, and completes the steps of the above method in combination with its hardware.

Embodiment four:

An embodiment of the present invention also provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is run by a processor, the steps of the method described in the first embodiment above are executed.

In addition, in the description of the embodiments of the present invention, unless otherwise specified and limited, the terms "installation", "connection" and "connection" should be interpreted in a broad sense, for example, it can be a fixed connection or a detachable connection , or integrally connected; it may be mechanically connected or electrically connected; it may be directly connected or indirectly connected through an intermediary, and it may be the internal communication of two components. Those of ordinary skill in the art can understand the specific meanings of the above terms in the present invention in specific situations.

In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer" etc. The indicated orientation or positional relationship is based on the orientation or positional relationship shown in the drawings, and is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the referred device or element must have a specific orientation, or in a specific orientation. construction and operation, therefore, should not be construed as limiting the invention. In addition, the terms "first", "second", and "third" are used for descriptive purposes only, and should not 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 ways. The device embodiments described above are only illustrative. For example, the division of the units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components can be combined or May be integrated into another system, or some features may be ignored, or not implemented. In another point, the mutual coupling or direct coupling or communication connection shown or discussed may be through some communication interfaces, and the indirect coupling or communication connection of devices or units may be in electrical, mechanical or other forms.

The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place, or may be distributed to multiple network units. Part or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, each unit may exist separately physically, or two or more units may be integrated into one unit.

Finally, it should be noted that: the above-described embodiments are only specific implementations of the present invention, used to illustrate the technical solutions of the present invention, rather than limiting them, and the scope of protection of the present invention is not limited thereto, although referring to the foregoing The embodiment has described the present invention in detail, and those of ordinary skill in the art should understand that any person familiar with the technical field can still modify the technical solutions described in the foregoing embodiments within the technical scope disclosed in the present invention Changes can be easily thought of, or equivalent replacements are made to some of the technical features; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should be included in the scope of the present invention within the scope of protection. Therefore, the protection scope of the present invention should be based on 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.

Images