Geological disaster hidden danger identification method and system based on InSAR technology

Abstract

The invention discloses a geological disaster hidden danger identification method and system based on an InSAR technology, and is applied to the technical field of intelligent monitoring, and the method comprises the steps: obtaining the SBAS deformation of a plurality of monitoring points of a target region, and dividing the target region into a plurality of monitoring grids according to the terrain of the target region; acquiring a monitoring point closest to the geometric center point of the monitoring grid in the monitoring grid as a reference monitoring point, and taking the SBAS deformation amount of the reference monitoring point as the reference deformation amount of the monitoring grid; and judging the geological disaster hidden danger risk of the monitoring grid according to the difference between the SBAS deformation of all the monitoring points in the monitoring grid and the reference deformation. According to the method, errors generated in the InSAR technology are synchronously eliminated in a data relative analysis mode, the actual deformation condition in a certain area can be effectively analyzed, the geological disaster hidden danger can be found out, the safety of the geological disaster risk body is improved, and extraction and early warning can be carried out on the risk body in time.

Description

TECHNICAL FIELD

[0001] The present application relates to intelligent monitoring technology, in particular to a geological disaster hidden danger identification method and system based on InSAR technology. BACKGROUND

[0002] SBAS (Small Baseline Subset) data is a surface deformation product obtained after InSAR interference-phase unwrapping-time inversion on multiple SAR images of the same area, which belongs to the post-processing output in the InSAR technical chain. Although InSAR technology can capture surface deformation at the millimeter level, atmospheric delay and DEM error can easily introduce centimeter-level pseudo deformation in geological disaster hidden danger identification, and long revisit period and large gradient deformation area often fail to phase unwrapping and miss reporting.

[0003] In the prior art, a time series InSAR method suitable for landslide multi-dimensional deformation monitoring is disclosed in Chinese patent No. CN202510168657.7. First, the deformation observation quantity of the line of sight of the measured landslide is obtained, and then it is input into the monitoring model containing the fractal composite structure deformation constraint model, thereby obtaining the deformation time series of the landslide body in the vertical direction, the east-west direction, the slope direction, and the vertical slope direction. This method has obvious advantages, can simultaneously obtain the deformation time series characteristics of the landslide in four dimensions of the vertical direction, the east-west direction, the slope direction, and the vertical slope direction, and the constructed model covers linear, acceleration, periodic, high-frequency, and thermal expansion deformation, fully considers the complexity of landslide deformation, and has strong adaptability.

[0004] From the prior art, it can be seen that the prior art focuses on the time series effect of the data generated by InSAR, which is affected by the above-mentioned atmospheric delay and DEM error, thereby generating a large error and affecting the final analysis result. SUMMARY

[0005] In order to at least overcome the above-mentioned deficiencies in the prior art, the purpose of the present application is to provide a geological disaster hidden danger identification method and system based on InSAR technology.

[0006] In a first aspect, the embodiments of the present application provide a geological disaster hidden danger identification method based on InSAR technology, comprising: obtaining the SBAS deformation of a plurality of monitoring points in a target area, and dividing the target area into a plurality of monitoring grids according to the topography of the target area; obtaining the monitoring point closest to the geometric center point of the monitoring grid in the monitoring grid as a reference monitoring point, and taking the SBAS deformation of the reference monitoring point as the reference deformation of the monitoring grid; judging the geological disaster hidden danger risk of the monitoring grid according to the difference between the SBAS deformation of all monitoring points in the monitoring grid and the reference deformation.

[0007] In a possible implementation, the dividing the target area into a plurality of monitoring grids according to the topography of the target area comprises: obtaining a moving direction of the target area when a disaster occurs, and constructing a plurality of longitudinal grid lines along the moving direction; constructing a plurality of transverse grid lines in a direction orthogonal to the longitudinal grid lines; surrounding the target area by the longitudinal grid lines, the transverse grid lines and a boundary line of the target area as the monitoring grid.

[0008] In a possible implementation, the judging the geological disaster hidden danger risk of the monitoring grid comprises: calculating physical distances from all the monitoring points in the monitoring grid to a reference monitoring point, and calculating absolute values of differences between each SBAS deformation and a reference deformation as deformation difference values; calculating ratios of the deformation difference values and corresponding physical distances as risk difference values of the corresponding monitoring points; calculating an average value of the risk difference values of all the monitoring points in the monitoring grid as a risk warning value of the monitoring grid; when the risk warning value exceeds a preset warning value, determining that the monitoring grid has a geological disaster hidden danger risk.

[0009] In a possible implementation, the risk warning value is calculated according to the following formula:

[0010] In the formula, R is the risk warning value of the monitoring grid, R i is the risk difference value of the i th monitoring point, x i is a horizontal coordinate value of the i th monitoring point, y i is a vertical coordinate value of the i th monitoring point, D i is the SBAS deformation of the i th monitoring point, x 0 is a horizontal coordinate value of the reference monitoring point, y 0 is a vertical coordinate value of the reference monitoring point, D 0 is the reference deformation of the reference monitoring point, and n is the number of monitoring points in the monitoring grid.

[0011] In a possible implementation, the obtaining of the reference monitoring point comprises: obtaining a coordinate of a geometric center point of the monitoring grid, and calculating physical distances from all the monitoring points in the monitoring grid to the geometric center point; taking the monitoring point with the smallest physical distance as the reference monitoring point.

[0012] In a second aspect, the application further provides a geological disaster hidden danger identification system based on InSAR technology, comprising: an acquisition unit configured to acquire SBAS deformation amounts of a plurality of monitoring points in a target area, and divide the target area into a plurality of monitoring grids according to a topography of the target area; a reference unit configured to acquire a monitoring point closest to a geometric center point of a monitoring grid as a reference monitoring point, and acquire the SBAS deformation amount of the reference monitoring point as a reference deformation amount of the monitoring grid; a judgment unit configured to judge a geological disaster risk of the monitoring grid according to a difference between the SBAS deformation amounts of all monitoring points in the monitoring grid and the reference deformation amount.

[0013] In a possible implementation, the acquisition unit is further configured to: acquire a motion direction of the target area when a disaster occurs, and construct a plurality of longitudinal grid lines along the motion direction; construct a plurality of transverse grid lines in a direction orthogonal to the longitudinal grid lines; surround an area surrounded by the longitudinal grid lines, the transverse grid lines and a boundary line of the target area as the monitoring grid.

[0014] In a possible implementation, the judgment unit is further configured to: calculate physical distances from the reference monitoring point to all monitoring points in the monitoring grid, and calculate absolute values of differences between each of the SBAS deformation amounts and the reference deformation amount as deformation difference values; calculate ratios of the deformation difference values to corresponding physical distances as risk difference values of the corresponding monitoring points; calculate an average value of the risk difference values of all monitoring points in the monitoring grid as a risk warning value of the monitoring grid; when the risk warning value exceeds a preset warning value, determine that the monitoring grid has a geological disaster risk.

[0015] In a possible implementation, the judgment unit is further configured to: calculate the risk warning value according to the following formula:

[0016] wherein R is the risk warning value of the monitoring grid, R i is the risk difference value of the i th monitoring point, x i is a horizontal coordinate value of the i th monitoring point, y i is a vertical coordinate value of the i th monitoring point, D iSBAS deformation of the i th monitoring point, x0 is the horizontal coordinate value of the reference monitoring point, y0 is the vertical coordinate value of the reference monitoring point, D0 is the reference deformation of the reference monitoring point, and n is the number of monitoring points in the monitoring grid.

[0017] In a possible implementation, the reference unit is further configured to: obtain the coordinates of the geometric center point of the monitoring grid, and calculate the physical distance from all monitoring points in the monitoring grid to the geometric center point; the monitoring point with the smallest physical distance is taken as the reference monitoring point.

[0018] Compared with the prior art, the present application has the following advantages and beneficial effects: The geological disaster hidden danger identification method and system based on InSAR technology can effectively analyze the actual deformation in a certain area to find out the geological disaster hidden danger, improve the safety of the geological disaster risk body, and can timely extract and warn the risk body. BRIEF DESCRIPTION OF DRAWINGS

[0019] The accompanying drawings, which are included to provide a further understanding of the embodiments of the present application and constitute a part of the application, do not constitute a limitation to the embodiments of the present application. In the drawings: Figure 1 The present application is a schematic diagram of the method steps of the embodiments; Figure 2 The present application is a schematic diagram of the embodiments. DETAILED DESCRIPTION

[0020] To make the purpose, technical scheme and advantages of the embodiments of the present application clearer, the technical scheme in the embodiments of the present application will be described clearly and completely below in combination with the drawings in the embodiments of the present application. It should be understood that the drawings in the present application only serve the purpose of illustration and description, and do not serve to limit the protection scope of the present application. In addition, it should be understood that the schematic drawings are not drawn to scale. The flowchart shows the operations implemented according to some embodiments of the present application. It should be understood that the operations of the flowchart can not be implemented in sequence, and the steps without logical context relationship can be reversed in sequence or implemented simultaneously. In addition, one or more other operations can be added to the flowchart or one or more operations can be removed from the flowchart under the guidance of the content of the present application.

[0021] In addition, the described embodiments are only some of the embodiments of the present application, not all the embodiments. The components of the embodiments of the present application described and shown in the drawings herein can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present application provided in the drawings is not intended to limit the scope of the claimed application, but only represents selected embodiments of the present application. Based on the embodiments of the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0022] Please refer to the figure, the flow chart of the geological disaster hidden danger identification method based on InSAR technology provided by the embodiment of the present application is shown, further, the geological disaster hidden danger identification method based on InSAR technology can specifically include the contents described in the following steps S1-S3.

[0023] S1: obtaining the SBAS deformation of a plurality of monitoring points in a target area, and dividing the target area into a plurality of monitoring grids according to the topography of the target area; S2: obtaining the monitoring point closest to the geometric center point of the monitoring grid in the monitoring grid as the reference monitoring point, and taking the SBAS deformation of the reference monitoring point as the reference deformation of the monitoring grid; S3: judging the geological disaster hidden danger risk of the monitoring grid according to the difference between the SBAS deformation of all monitoring points in the monitoring grid and the reference deformation.

[0024] When the embodiments of the present application are implemented, the SBAS deformation of a plurality of monitoring points in a target area needs to be obtained by InSAR technology, which can be realized by joint decomposition of ascending and descending track observation and other technical means, which belongs to mature existing technology. The present application does not make more limitations, and can refer to the specification T / CAGHP 013-2018 geological disaster InSAR monitoring technology guide for monitoring. At the same time, the target area in the embodiments of the present application is generally a complete geological disaster risk body, such as a complete landslide body. In order to reduce the error caused by InSAR technology, the target area is divided into a plurality of monitoring grids in the embodiments of the present application, and each monitoring grid is evaluated independently. For InSAR technology, in satellite photography, atmospheric delay and DEM error generally remain consistent within a monitoring grid, so independent evaluation of each monitoring grid can effectively reduce the evaluation error.

[0025] In the embodiment of the present application, the reference monitoring points are needed to be acquired first, each of which is the position closest to the geometric center of the monitoring grid in physical distance, and the overall deformation in the monitoring grid can be analyzed based on the reference deformation generated by the reference monitoring points and the deformations of other monitoring points in the grid, and then the geological disaster hidden danger risk of the monitoring grid can be judged. The error generated in the InSAR technology is synchronously eliminated by the relative data analysis manner, the actual deformation in a region can be effectively analyzed to find out the geological disaster hidden danger, the safety of the geological disaster risk body is improved, and the risk body can be extracted and warned in time.

[0026] In a possible implementation manner, the target region is divided into a plurality of monitoring grids according to the topography of the target region, including: acquiring a motion direction of the target region when a disaster occurs, and constructing a plurality of longitudinal grid lines along the motion direction; constructing a plurality of transverse grid lines in a direction orthogonal to the longitudinal grid lines; surrounding the target region by the longitudinal grid lines, the transverse grid lines and the boundary line as the monitoring grid.

[0027] In the embodiment of the present application, please refer to Figure 2 , which gives an example of dividing the monitoring grid, wherein the motion direction of the target region when a disaster occurs is needed to be acquired first, for the landslide body, the motion direction is generally the direction of the sliding surface, Figure 2 , the position of the landslide body is marked by green lines, and the entire landslide body is the target region. At this time, the possible sliding direction of the landslide body is acquired, in the example, the possible sliding direction of the landslide body is from the top of the slope to the slope foot of the river bank, so the transverse grid lines are constructed based on this. Meanwhile, the longitudinal grid lines orthogonal to the transverse grid lines are constructed, and a plurality of monitoring grids are segmented. As can be seen from Figure 2 , the monitoring points in different monitoring grids are represented by points from blue to red, and the color closer to red indicates that the deformation monitored is larger. It should be understood that the deformation monitored is the cumulative deformation in a time period, rather than the real-time quantity, because the satellite needs to fly over the region to generate a group of data. It should be understood that the specific interval division of the transverse grid lines and the longitudinal grid lines is generally based on experience, such as being divided in the width range of 50-100m to ensure that the atmospheric influence in the range is as consistent as possible, and the grid width shown in Figure 2 is 70m.

[0028] In a possible implementation manner, judging the geological disaster hidden danger risk of the monitoring grid includes: calculating the absolute value of the difference between the SBAS deformation of each of the monitoring points and the reference deformation as a deformation difference value; calculating the ratio of the deformation difference value and the corresponding physical distance as a risk difference value of the corresponding monitoring point; calculating the average value of the risk difference values of all the monitoring points in the monitoring grid as a risk warning value of the monitoring grid; determining that the monitoring grid has a geological disaster hidden danger risk when the risk warning value exceeds a preset warning value.

[0029] In the implementation of the embodiments of the present application, in the determination of the geological disaster hidden danger risk of the monitoring grid, the relative deformation calculation method is adopted. For the same monitoring grid, the error values of the SBAS deformation are similar, so the physical distances from different monitoring points to the reference monitoring point and the absolute difference values between the SBAS deformation of each monitoring point and the reference deformation need to be calculated first. The absolute difference value can effectively eliminate the errors in the monitoring grid, so the absolute difference value can better express the deformation difference of different monitoring points and the center point of the monitoring grid.

[0030] In the embodiments of the present application, since the deformation of a monitoring grid is continuous, the ratio of the above deformation difference and the specific physical distance can express the deformation of the landslide body in the unit distance, which expresses the stability of the region. At the same time, the ratio can also effectively reflect the factors such as the steepness of the slope and the influence of the vegetation in the monitoring grid, which have an impact on the deformation difference of different points. It is a comprehensive index and can well express the stability of the monitoring grid. It should be understood that for the InSAR technology, the satellite image generated is a two-dimensional graph, so the physical distance in the embodiments of the present application is a physical distance in the two-dimensional graph, that is, the physical distance calculated by the horizontal and vertical coordinates in the two-dimensional graph.

[0031] In the implementation of the embodiments of the present application, the risk warning value in the monitoring grid can be obtained by calculating the average value of the risk difference values of all the monitoring points in the monitoring grid. At this time, the risk warning value and the warning value obtained by statistics can realize the risk condition judgment in the monitoring grid. The warning value needs to be selected according to the period of the above SBAS deformation, and the period warning value in years is preferably 8.3 mm / km.

[0032] In a possible implementation, the risk warning value is calculated according to the following formula:

[0033] In the formula, R is a risk warning value of a monitoring grid, R i is a risk difference value of the ith monitoring point, x i is a horizontal coordinate value of the ith monitoring point, y i is a vertical coordinate value of the ith monitoring point, D i is an SBAS deformation amount of the ith monitoring point, x0 is a horizontal coordinate value of a reference monitoring point, y0 is a vertical coordinate value of the reference monitoring point, D0 is a reference deformation amount of the reference monitoring point, and n is a number of monitoring points in the monitoring grid.

[0034] In the implementation of the embodiments of the present application, a specific risk warning value calculation scheme is provided, in which the horizontal coordinate value is generally selected as a coordinate value along the latitude direction, and the vertical coordinate value is generally selected as a coordinate value along the longitude direction.

[0035] In a possible implementation, the obtaining of the reference monitoring point comprises: obtaining a coordinate of a geometric center point of the monitoring grid and calculating physical distances of all monitoring points in the monitoring grid to the geometric center point; taking the monitoring point with the smallest physical distance as the reference monitoring point.

[0036] Based on the same inventive concept, the present application also provides a geological disaster hidden danger identification system based on InSAR technology, comprising: an obtaining unit configured to obtain SBAS deformation amounts of multiple monitoring points in a target region and divide the target region into multiple monitoring grids according to a topography of the target region; a reference unit configured to obtain a monitoring point closest to a geometric center point of the monitoring grid as a reference monitoring point and take the SBAS deformation amount of the reference monitoring point as a reference deformation amount of the monitoring grid; a judging unit configured to judge a geological disaster hidden danger risk of the monitoring grid according to differences between SBAS deformation amounts of all monitoring points in the monitoring grid and the reference deformation amount.

[0037] In a possible implementation, the obtaining unit is further configured to: obtain a movement direction of the target region when a disaster occurs and construct multiple longitudinal grid lines along the movement direction; construct multiple transverse grid lines in a direction orthogonal to the longitudinal grid lines; enclose an area surrounded by the longitudinal grid lines, the transverse grid lines, and a boundary line of the target region as the monitoring grid.

[0038] In a possible implementation, the judging unit is further configured to: calculating a physical distance from the reference monitoring point to each monitoring point in the monitoring grid, and calculating a risk difference value of each monitoring point as a ratio of the deformation difference value and the corresponding physical distance; calculating a physical distance from the reference monitoring point to each monitoring point in the monitoring grid, and calculating a risk difference value of each monitoring point as a ratio of the deformation difference value and the corresponding physical distance; calculating a risk warning value of the monitoring grid as an average of risk difference values of all monitoring points in the monitoring grid; determining that the monitoring grid has a geological disaster risk when the risk warning value exceeds a preset warning value.

[0039] In a possible implementation, the determining unit is further configured to: calculate the risk warning value according to the following formula:

[0040] wherein, R is the risk warning value of the monitoring grid, R i is the risk difference value of the ith monitoring point, x i is the horizontal coordinate value of the ith monitoring point, y i is the vertical coordinate value of the ith monitoring point, D i is the SBAS deformation value of the ith monitoring point, x0 is the horizontal coordinate value of the reference monitoring point, y0 is the vertical coordinate value of the reference monitoring point, D0 is the reference deformation value of the reference monitoring point, and n is the number of monitoring points in the monitoring grid.

[0041] In a possible implementation, the reference unit is further configured to: obtain a coordinate of a geometric center point of the monitoring grid, and calculate a physical distance from the geometric center point to each monitoring point in the monitoring grid; determine the monitoring point with the minimum physical distance as the reference monitoring point.

[0042] Those skilled in the art can appreciate that the units and algorithm steps of the examples described in conjunction with the embodiments disclosed herein can be realized by electronic hardware, computer software or a combination of both. In order to clearly illustrate the interchangeability of hardware and software, the components and steps of the examples have been described in a general manner in the foregoing description. Whether the functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the present application.

[0043] In several embodiments provided in the present application, it should be understood that the disclosed apparatus and method can be implemented by other manners. For example, the apparatus embodiments described above are merely illustrative, for example, the division of the units is merely a logical function division, and actual implementation can have another division manner, for example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the units shown or discussed can be indirect coupling or communication connection through some interfaces, devices or units, and can also be electrical, mechanical or other form of connection.

[0044] The units described as separate components can or can not be physically separated, and it is obvious to those skilled in the art that the units and algorithm steps of each example described in combination with the embodiments disclosed herein can be realized by electronic hardware, computer software or a combination of both. In order to clearly illustrate the interchangeability of hardware and software, the components and steps of each example have been described in the above description in general terms. Whether the functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. The skilled person can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the present application.

[0045] In addition, each functional unit in each embodiment of the present application can be integrated into one processing unit, or each unit can exist physically, or two or more units can be integrated into one unit. The integrated unit can be realized in the form of hardware or software functional unit.

[0046] When the integrated unit is realized in the form of software functional unit and sold or used as an independent product, it can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the part that contributes to the prior art, or the whole or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a grid device, etc.) to execute all or part of the steps of the method described in each embodiment of the present application. The foregoing storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various program code storage media.

[0047] The above detailed description of the specific embodiments of the present application has been given to understand the purpose, technical solutions and beneficial effects of the present application. It should be understood that the above description is only a specific embodiment of the present application and is not used to limit the protection scope of the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for identifying geological disaster hidden dangers based on InSAR technology, characterized in that, The method comprises the following steps: obtaining SBAS deformation of multiple monitoring points in a target area, and dividing the target area into multiple monitoring grids according to the topography of the target area; obtaining a monitoring point closest to a geometric center point of the monitoring grid in the monitoring grid as a reference monitoring point, and obtaining the SBAS deformation of the reference monitoring point as a reference deformation of the monitoring grid; judging geological disaster hidden danger risk of the monitoring grid according to the difference between the SBAS deformation of all monitoring points in the monitoring grid and the reference deformation; judging geological disaster hidden danger risk of the monitoring grid comprises: calculating physical distances from all monitoring points in the monitoring grid to the reference monitoring point, and calculating absolute values of the difference between each SBAS deformation and the reference deformation as deformation difference values; calculating ratios of the deformation difference values and corresponding physical distances as risk difference values of the corresponding monitoring points; calculating an average value of the risk difference values of all monitoring points in the monitoring grid as a risk warning value of the monitoring grid; when the risk warning value exceeds a preset warning value, it is determined that the monitoring grid has geological disaster hidden danger risk. 2.The InSAR technology-based geological disaster hidden danger identification method according to claim 1, characterized in that, dividing the target area into multiple monitoring grids according to the topography of the target area comprises: obtaining a movement direction of the target area when a disaster occurs, and constructing multiple longitudinal grid lines along the movement direction; constructing multiple transverse grid lines in a direction perpendicular to the longitudinal grid lines; the area surrounded by the longitudinal grid lines, the transverse grid lines and the boundary line of the target area is taken as the monitoring grid. 3.The InSAR technology-based geological disaster hidden danger identification method according to claim 1, characterized in that, The risk warning value is calculated according to the following formula: ； In the formula, R is a risk early warning value of a monitoring grid, R i is a risk difference value of the ith monitoring point, x i is a horizontal coordinate value of the ith monitoring point, y i is a vertical coordinate value of the ith monitoring point, D i is an SBAS deformation amount of the ith monitoring point, x0 is a horizontal coordinate value of a reference monitoring point, y0 is a vertical coordinate value of the reference monitoring point, D0 is a reference deformation amount of the reference monitoring point, and n is a number of monitoring points in the monitoring grid. 4.The InSAR technology-based geological disaster hidden danger identification method according to claim 1, characterized in that, the reference monitoring point comprises: obtaining coordinates of the geometric center point of the monitoring grid, and calculating physical distances from all monitoring points in the monitoring grid to the geometric center point; the monitoring point with the smallest physical distance is taken as the reference monitoring point.

5. A geological disaster hazard identification system based on InSAR technology, characterized in that, The method comprises the following steps: an obtaining unit is configured to obtain SBAS deformation of multiple monitoring points in a target area, and divide the target area into multiple monitoring grids according to the topography of the target area; a reference unit is configured to obtain a monitoring point closest to a geometric center point of the monitoring grid in the monitoring grid as a reference monitoring point, and obtain the SBAS deformation of the reference monitoring point as a reference deformation of the monitoring grid; a judging unit is configured to judge geological disaster hidden danger risk of the monitoring grid according to the difference between the SBAS deformation of all monitoring points in the monitoring grid and the reference deformation; the judging unit is further configured to: calculate physical distances from all monitoring points in the monitoring grid to the reference monitoring point, and calculate absolute values of the difference between each SBAS deformation and the reference deformation as deformation difference values; calculate ratios of the deformation difference values and corresponding physical distances as risk difference values of the corresponding monitoring points; calculate an average value of the risk difference values of all monitoring points in the monitoring grid as a risk warning value of the monitoring grid; when the risk warning value exceeds a preset warning value, it is determined that the monitoring grid has geological disaster hidden danger risk. 6.The InSAR technology-based geological disaster hidden danger identification system according to claim 5, characterized in that, the obtaining unit is further configured to: acquire a moving direction of the target area when a disaster occurs, and construct a plurality of longitudinal grid lines along the moving direction; construct a plurality of transverse grid lines in a direction orthogonal to the longitudinal grid lines; surround an area formed by the longitudinal grid lines, the transverse grid lines and a boundary line of the target area as the monitoring grid. 7.The InSAR technology-based geological disaster hidden danger identification system according to claim 5, characterized in that, The judging unit is further configured to: calculate the risk warning value according to the following formula: ； In the formula, R is a risk early warning value of a monitoring grid, R i is a risk difference value of the ith monitoring point, x i is a horizontal coordinate value of the ith monitoring point, y i is a vertical coordinate value of the ith monitoring point, D i is an SBAS deformation amount of the ith monitoring point, x0 is a horizontal coordinate value of a reference monitoring point, y0 is a vertical coordinate value of the reference monitoring point, D0 is a reference deformation amount of the reference monitoring point, and n is a number of monitoring points in the monitoring grid. 8.The InSAR technology based geological disaster hidden danger identification system according to claim 5, characterized in that, The reference unit is further configured to: acquire a coordinate of a geometric center point of the monitoring grid, and calculate a physical distance from all monitoring points in the monitoring grid to the geometric center point; take the monitoring point with the minimum physical distance as the reference monitoring point.

Images