Landslide susceptibility partition map making method, system and device and storage medium

Abstract

A landslide susceptibility partition map making method, system and device and a storage medium, the method comprises the steps of firstly obtaining geographic data of a geological disaster susceptibility area, then preprocessing radar data to obtain an initial surface deformation rate map, then based on historical landslide data and geological factor data, adopting an evidence weight method to carry out assignment weighting to obtain an initial surface deformation rate map; making an initial landslide susceptibility partition map, then performing assignment updating on the initial earth surface deformation rate map and the initial landslide susceptibility partition map based on the contribution matrix, and coupling the updated two to obtain a final landslide susceptibility partition map; in the application, based on the static geological data of which the detection area is unchanged and common, and in combination with the recent dynamic performance of the earth surface deformation rate map, the landslide hidden danger screening of the detection area is comprehensively realized, and meanwhile, the influence of each factor on the final result is displayed through the landslide susceptibility partition map generated by the evidence weight method, so that the human subjective influence is reduced. The method is not only objective, but also accurate.

Description

Technical Field

[0001] The invention relates to a means for making a zoning map, belongs to the field of landslide hazard identification, and in particular to a method, system, equipment and storage medium for making a landslide susceptibility zoning map. Background Art

[0002] With the development of social economy and the acceleration of urbanization, the relationship between man and nature is becoming closer and closer, and the frequency of landslide disasters is also increasing. There is an urgent need for effective landslide identification and evaluation technology to identify potential landslide and geological disaster risk areas at an early stage on a regional scale, so as to avoid geological risks.

[0003] Surface deformation is the most intuitive and comprehensive manifestation of landslide occurrence and is the preferred indicator for landslide identification and evaluation. However, landslides identified by InSAR technology may have errors, and the excessive surface deformation caused by many landslides will cause InSAR decoherence, making InSAR impossible to identify. In addition, landslide susceptibility assessment is also an effective means to identify landslide hazards. It aims to partition the possibility of landslide occurrence in the survey area by analyzing a series of disaster-causing factors in the survey area, but it has certain subjectivity and limitations of a small evaluation range.

[0004] In the prior art, there is a method of using InSAR surface deformation results to assist in verifying the landslide susceptibility zoning map, but the InSAR results are not integrated into the production process and are only used as a verification method, which does not improve the accuracy of the landslide susceptibility zoning map. There is also a method of using MT-InSAR results as one of the factors to participate in the landslide susceptibility zoning map extraction using the hierarchical analysis method, but this method is based on expert experience and has strong subjectivity. There is also a method of using InSAR settlement results as one of the inputs to train a neural network and extract potential landslides, but the training set production is relatively cumbersome and the engineering workload is large. Therefore, an objective and accurate method is urgently needed to solve the above problems existing in the prior art. Summary of the invention

[0005] The purpose of the present invention is to overcome the above-mentioned defects and problems existing in the prior art and to provide an objective and accurate method, system, device and storage medium for making a landslide susceptibility zoning map.

[0006] To achieve the above objectives, the technical solution of the present invention is: a method for making a landslide susceptibility zoning map, comprising:

[0007] S1. Obtaining geographic data of areas prone to geological disasters; the geographic data includes historical landslide data, three-land survey data, geological factor data, road data, digital elevation model and radar data;

[0008] S2, preprocessing based on radar data to obtain the initial surface deformation rate map of the geological disaster prone area;

[0009] S3. Based on historical landslide data and geological factor data, the weight of evidence method is used to assign weights and produce an initial landslide susceptibility zoning map for geological disaster-prone areas;

[0010] S4. Based on the contribution matrix, the initial surface deformation rate map and the initial landslide susceptibility zoning map are assigned and updated to obtain an updated surface deformation rate map and an updated landslide susceptibility zoning map, and the two are coupled to obtain a final landslide susceptibility zoning map for geological disaster-prone areas.

[0011] The step S2 specifically includes:

[0012] S21, performing synthetic aperture radar interferometry processing on the radar data, and then using the least square method to calculate and obtain the surface deformation time series;

[0013] S22. Based on the singular value decomposition method, the surface deformation time series is temporally and spatially fitted to obtain the initial surface deformation rate map.

[0014] The step S3 specifically includes:

[0015] S31, presetting grid units as the spatial scale limit of geological disaster prone areas, and selecting evidence factors from historical landslide data and geological factor data and assigning corresponding layers; the evidence factors include elevation, slope, slope aspect, fault distance, engineering rock group, water system distance, precipitation, road distance and land classification;

[0016] S32, calculating the contribution weight of each evidence factor to the occurrence of geological hazards through conditional probability, assigning a comprehensive weight to each evidence factor layer, and then superimposing and summing the comprehensive weights of each evidence factor layer to obtain the susceptibility index of the grid unit;

[0017] S33, assigning the susceptibility index to each evidence factor layer for weighted superposition to obtain the initial landslide susceptibility zoning map; the calculation expression is as follows:

[0018]

[0019]

[0020]

[0021] in: is the weight value of the evidence factor existence area; is the weight value of the area where the evidence factor does not exist; W iis the weight value of the influence of the evidence factor on the occurrence of landslide; P is the conditional probability; F is the evidence factor; L is the number of landslide hazard units; is the number of non-landslide hazard units; (F|L) and are the number of units where the evidence factor exists in the landslide hazard unit and the non-landslide hazard unit; and are the numbers of units where the evidence factor does not exist in the landslide hazard unit and the non-landslide hazard unit, respectively.

[0022] In step S4, obtaining and updating the surface deformation rate map specifically includes:

[0023] S41, based on the nearest neighbor method, downsampling the initial surface deformation rate map to the initial landslide susceptibility partition Figure 1 consistent grid cells;

[0024] S42, dividing the initial surface deformation rate map into a plurality of deformation rate intervals according to the surface deformation rate, and classifying and counting the pixel points in different deformation rate intervals to obtain the deformation severity of the different deformation rate intervals;

[0025] S43. Assigning a contribution degree to the initial surface deformation rate map based on the deformation severity to obtain an updated surface deformation rate map.

[0026] In step S4, obtaining and updating the landslide susceptibility zoning map specifically includes:

[0027] S44. Assign different contribution degrees to the initial landslide susceptibility zoning map according to the hazard level to obtain an updated landslide susceptibility zoning map.

[0028] In step S4, obtaining the final landslide susceptibility zoning map of the geological disaster prone area specifically includes:

[0029] S45, calculating the contribution of a certain area on the updated landslide susceptibility zoning map and the updated surface deformation rate map to obtain an updated landslide susceptibility index LSI; the calculation formula is as follows:

[0030] LSI=C LSM *C GDM ;

[0031] Where: C LSM Contribution value for updating the landslide susceptibility zoning map; C GDM Contribution value for updating the surface deformation rate map;

[0032] S46. Zoning the updated landslide susceptibility index and assigning hazard levels to different areas respectively, thereby obtaining a final landslide susceptibility zoning map for geological disaster-prone areas.

[0033] A landslide susceptibility zoning map making system, which is applied to the landslide susceptibility zoning map making method, and the system comprises:

[0034] A geographic data acquisition module is used to acquire geographic data of areas prone to geological disasters; the geographic data includes historical landslide data, three-land survey data, geological factor data, road data, digital elevation model and radar data;

[0035] An initial surface deformation rate map construction module is used to perform preprocessing based on radar data to obtain an initial surface deformation rate map of geological disaster-prone areas;

[0036] The initial landslide susceptibility zoning map construction module is used to produce the initial landslide susceptibility zoning map of the geological disaster prone area based on the historical landslide data and geological factor data and the weight of evidence method for weighting;

[0037] The coupling update module is used to assign values ​​to and update the initial surface deformation rate map and the initial landslide susceptibility zoning map based on the contribution matrix, obtain the updated surface deformation rate map and the updated landslide susceptibility zoning map, and couple the two to obtain the final landslide susceptibility zoning map of the geological disaster-prone area.

[0038] A device for making a landslide susceptibility zoning map, the device comprising a processor and a memory;

[0039] The memory is used to store computer program code and transmit the computer program code to the processor;

[0040] The processor is used to execute the method for preparing a landslide susceptibility zoning map according to the instructions in the computer program code.

[0041] A computer-readable storage medium stores computer-executable instructions. When the computer-executable instructions are executed on a computer, the method for making a landslide susceptibility zoning map is implemented.

[0042] Compared with the prior art, the present invention has the following beneficial effects:

[0043] In the method, system, device and storage medium for making a landslide susceptibility zoning map of the present invention, the method first obtains geographic data of a geological disaster prone area, then pre-processes radar data to obtain an initial surface deformation rate map, then uses the evidence weight method to assign weights based on historical landslide data and geological factor data to make an initial landslide susceptibility zoning map, then assigns and updates the initial surface deformation rate map and the initial landslide susceptibility zoning map based on a contribution matrix, and couples the updated two to obtain a final landslide susceptibility zoning map of a geological disaster prone area; in the application, the design is based on the static geological data of the detection area that is unchanged and common, and combined with the recent dynamic performance of the surface deformation rate map, to comprehensively realize the screening of landslide hazards in the detection area, and the landslide susceptibility zoning map generated by the evidence weight method can show the influence of each factor on the final extraction result, reduce the influence of human subjectivity, and improve objectivity. Therefore, the present invention is not only more objective, but also more accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] Figure 1 It is a schematic diagram of the process of the present invention.

[0045] Figure 2 It is a schematic diagram of the system structure of the present invention.

[0046] Figure 3 It is a schematic diagram of the device structure of the present invention.

[0047] In the figure: geographic data acquisition module 1, initial surface deformation rate map construction module 2, initial landslide susceptibility zoning map construction module 3, coupling update module 4, processor 5, memory 6, computer program code 61. DETAILED DESCRIPTION

[0048] The present invention is further described in detail below in conjunction with the accompanying drawings and specific implementation methods.

[0049] Embodiment 1:

[0050] See also Figure 1 , a method for preparing a landslide susceptibility zoning map, comprising:

[0051] S1. Obtaining geographic data of areas prone to geological disasters; the geographic data includes historical landslide data, three-land survey data, geological factor data, road data, digital elevation model and radar data;

[0052] S2, preprocessing based on radar data to obtain the initial surface deformation rate map of the geological disaster prone area;

[0053] Furthermore, the step S2 specifically includes:

[0054] S21. Small baseline aggregate synthetic aperture radar interferometry (SBAS-InSAR) processing is performed on the radar data to overcome the decoherence phenomenon caused by the long temporal and spatial baseline, and to alleviate the influence of atmospheric effect and terrain on the difference. Then, the surface deformation time series is obtained by the least squares method.

[0055] S22. Based on the singular value decomposition (SVD) method, the surface deformation time series is spatiotemporally fitted to obtain the initial surface deformation rate map with an accuracy of up to millimeter level.

[0056] S3. Based on historical landslide data and geological factor data, the weight of evidence method is used to assign weights and produce an initial landslide susceptibility zoning map for geological disaster-prone areas;

[0057] Furthermore, the step S3 specifically includes:

[0058] S31. Preset a grid unit of 100 m × 100 m as the spatial scale limit of the geological disaster-prone area, and select elevation, slope, slope direction, fault distance, engineering rock group, water system distance, precipitation, road distance and land classification from the historical landslide data and geological factor data as evidence factors and assign corresponding layers;

[0059] S32, calculating the contribution weight of each evidence factor to the occurrence of geological hazards through conditional probability, assigning a comprehensive weight to each evidence factor layer, and then superimposing and summing the comprehensive weights of each evidence factor layer to obtain the susceptibility index of the grid unit;

[0060] S33, assigning the susceptibility index to each evidence factor layer for weighted superposition to obtain an initial landslide susceptibility zoning map; the calculation expression of the evidence weight method is as follows:

[0061]

[0062]

[0063]

[0064] in: is the weight value of the evidence factor existence area; is the weight value of the area where the evidence factor does not exist; W i is the weight value of the influence of the evidence factor on the occurrence of landslide; P is the conditional probability; F is the evidence factor; L is the number of landslide hazard units; is the number of non-landslide hazard units; (F|L) and are the number of units where the evidence factor exists in the landslide hazard unit and the non-landslide hazard unit; and are the numbers of units where the evidence factor does not exist in the landslide hazard unit and the non-landslide hazard unit, respectively.

[0065] S4. Based on the contribution matrix, the initial surface deformation rate map and the initial landslide susceptibility zoning map are assigned and updated to obtain an updated surface deformation rate map and an updated landslide susceptibility zoning map, and the two are coupled to obtain a final landslide susceptibility zoning map for geological disaster-prone areas.

[0066] Further, step S4 specifically includes:

[0067] S41, based on the nearest neighbor method, downsampling the initial surface deformation rate map to the initial landslide susceptibility partition Figure 1 consistent grid cells;

[0068] S42, dividing the initial surface deformation rate map into a plurality of deformation rate intervals according to the surface deformation rate, and classifying and counting the pixel points in different deformation rate intervals to obtain the deformation severity of the different deformation rate intervals;

[0069] S43, assigning a contribution to the initial surface deformation rate map based on the severity of the deformation, wherein the greater the deformation severity, the higher the contribution, and finally obtaining an updated surface deformation rate map;

[0070] S44, assigning different contribution degrees to the initial landslide susceptibility zoning map according to the hazard level to obtain an updated landslide susceptibility zoning map; the hazard level includes very high, high, medium, and low;

[0071] S45, calculating the contribution of a certain area on the updated landslide susceptibility zoning map and the updated surface deformation rate map to obtain an updated landslide susceptibility index LSI; the calculation formula is as follows:

[0072] LSI=C LSM *C GDM ;

[0073] Where: C LSM Contribution value for updating the landslide susceptibility zoning map; C GDM Contribution value for updating the surface deformation rate map;

[0074] S46. Zoning the updated landslide susceptibility index and assigning hazard levels to different areas respectively, to obtain a final landslide susceptibility zoning map for geological disaster prone areas; the hazard levels include very high, very high, high, medium, and low.

[0075] Embodiment 2:

[0076] See also Figure 2 , a landslide susceptibility zoning map making system, the system is applied to the landslide susceptibility zoning map making method, the system comprises:

[0077] Geographic data acquisition module 1, used to acquire geographic data of areas prone to geological disasters; the geographic data includes historical landslide data, three-land survey data, geological factor data, road data, digital elevation model and radar data;

[0078] Initial surface deformation rate map construction module 2, used for preprocessing based on radar data to obtain an initial surface deformation rate map of a geological disaster prone area;

[0079] Further, the initial surface deformation rate map construction module 2 obtains the initial surface deformation rate map according to the following steps;

[0080] S21, performing synthetic aperture radar interferometry processing on radar data and using the least square method to obtain the surface deformation time series;

[0081] S22, performing spatiotemporal fitting on the surface deformation time series based on the singular value decomposition method to obtain an initial surface deformation rate map;

[0082] The initial landslide susceptibility zoning map construction module 3 is used to make the initial landslide susceptibility zoning map of the geological disaster prone area based on the historical landslide data and geological factor data and weighted by using the weight of evidence method;

[0083] Further, the initial landslide susceptibility zoning map construction module 3 obtains the initial landslide susceptibility zoning map according to the following steps;

[0084] S31, presetting grid units as the spatial scale limit of geological disaster prone areas, and selecting evidence factors from historical landslide data and geological factor data and assigning corresponding layers; the evidence factors include elevation, slope, slope aspect, fault distance, engineering rock group, water system distance, precipitation, road distance and land classification;

[0085] S32, calculating the contribution weight of each evidence factor to the occurrence of geological hazards through conditional probability, assigning a comprehensive weight to each evidence factor layer, and then superimposing and summing the comprehensive weights of each evidence factor layer to obtain the susceptibility index of the grid unit;

[0086] S33, assigning the susceptibility index to each evidence factor layer for weighted superposition to obtain the initial landslide susceptibility zoning map; the calculation expression is as follows:

[0087]

[0088]

[0089]

[0090] in: is the weight value of the evidence factor existence area; is the weight value of the area where the evidence factor does not exist; W i is the weight value of the influence of the evidence factor on the occurrence of landslide; P is the conditional probability; F is the evidence factor; L is the number of landslide hazard units; is the number of non-landslide hazard units; (F|L) and are the number of units where the evidence factor exists in the landslide hazard unit and the non-landslide hazard unit; and are the numbers of units where the evidence factor does not exist in the landslide hazard unit and the non-landslide hazard unit, respectively.

[0091] A coupling update module 4 is used to assign and update the initial surface deformation rate map and the initial landslide susceptibility zoning map based on the contribution matrix, obtain an updated surface deformation rate map and an updated landslide susceptibility zoning map, and couple the two to obtain a final landslide susceptibility zoning map of the geological disaster prone area;

[0092] Further, the coupling update module 4 obtains the final landslide susceptibility zoning map of the geological disaster prone area according to the following steps;

[0093] S41, based on the nearest neighbor method, downsampling the initial surface deformation rate map to the initial landslide susceptibility partition Figure 1 consistent grid cells;

[0094] S42, dividing the initial surface deformation rate map into a plurality of deformation rate intervals according to the surface deformation rate, and classifying and counting the pixel points in different deformation rate intervals to obtain the deformation severity of the different deformation rate intervals;

[0095] S43, assigning a contribution to the initial surface deformation rate map based on the deformation severity, to obtain an updated surface deformation rate map;

[0096] S44, assigning different contribution degrees to the initial landslide susceptibility zoning map according to the hazard level, to obtain an updated landslide susceptibility zoning map;

[0097] S45, calculating the contribution of a certain area on the updated landslide susceptibility zoning map and the updated surface deformation rate map to obtain an updated landslide susceptibility index LSI; the calculation formula is as follows:

[0098] LSI=C LSM *C GDM ;

[0099] Where: C LSM Contribution value for updating the landslide susceptibility zoning map; C GDM Contribution value for updating the surface deformation rate map;

[0100] S46. Zoning the updated landslide susceptibility index and assigning hazard levels to different areas respectively, thereby obtaining a final landslide susceptibility zoning map for geological disaster-prone areas.

[0101] Embodiment 3:

[0102] See also Figure 3 , a device for making a landslide susceptibility zoning map, the device comprising a processor 5 and a memory 6;

[0103] The memory 6 is used to store the computer program code 61 and transmit the computer program code 61 to the processor 5;

[0104] The processor 5 is used to execute the method for making a landslide susceptibility zoning map according to the instructions in the computer program code 61.

[0105] Embodiment 4:

[0106] A computer-readable storage medium stores computer-executable instructions. When the computer-executable instructions are executed on a computer, the method for producing a landslide susceptibility zoning map described in Example 1 is implemented.

[0107] Generally speaking, the computer instructions for implementing the method of the present invention may be carried in any combination of one or more computer-readable storage media. Non-transitory computer-readable storage media may include any computer-readable media, except for the signal itself that is temporarily propagating.

[0108] The computer-readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or device, or any combination thereof. More specific examples of computer-readable storage media (a non-exhaustive list) include: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EKROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present invention, a computer-readable storage medium may be any tangible medium containing or storing a program that may be used by or in conjunction with an instruction execution system, device, or device.

[0109] Computer program code for performing the operation of the present invention can be written in one or more programming languages ​​or a combination thereof, including object-oriented programming languages ​​such as Java, SMalltalk, C++, and conventional procedural programming languages ​​such as "C" or similar programming languages, in particular, Python suitable for neural network computing and platform frameworks based on TensorFlow, PyTorch, etc. can be used. The program code can be executed entirely on the user's computer, partially on the user's computer, as an independent software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In the case of a remote computer, the remote computer can be connected to the user's computer or to an external computer (for example, using an Internet service provider to connect via the Internet) through any type of network, including a local area network (LAN) or a wide area network (WAN).

[0110] The above-mentioned device and non-temporary computer-readable storage medium can be found in the detailed description of a method for preparing a landslide susceptibility zoning map and its beneficial effects, which will not be repeated here.

[0111] Although the embodiments of the present invention have been shown and described above, it should be understood that the above embodiments are exemplary and should not be construed as limitations on the present invention. A person skilled in the art may change, modify, replace and vary the above embodiments within the scope of the present invention.

Claims

1. A method for preparing a landslide susceptibility zoning map, characterized in that: include: S1. Obtaining geographic data of areas prone to geological disasters; the geographic data includes historical landslide data, three-land survey data, geological factor data, road data, digital elevation model and radar data; S2, preprocessing based on radar data to obtain the initial surface deformation rate map of the geological disaster prone area; S3. Based on historical landslide data and geological factor data, the weight of evidence method is used to assign weights and produce an initial landslide susceptibility zoning map for geological disaster-prone areas; S4. Based on the contribution matrix, the initial surface deformation rate map and the initial landslide susceptibility zoning map are assigned and updated to obtain an updated surface deformation rate map and an updated landslide susceptibility zoning map, and the two are coupled to obtain a final landslide susceptibility zoning map for geological disaster-prone areas.

2. The method for preparing a landslide susceptibility zoning map according to claim 1, characterized in that: The step S2 specifically includes: S21, performing synthetic aperture radar interferometry processing on the radar data, and then using the least square method to calculate and obtain the surface deformation time series; S22. Based on the singular value decomposition method, the surface deformation time series is temporally and spatially fitted to obtain the initial surface deformation rate map.

3. The method for preparing a landslide susceptibility zoning map according to claim 1, characterized in that: The step S3 specifically includes: S31, presetting grid units as the spatial scale limit of geological disaster prone areas, and selecting evidence factors from historical landslide data and geological factor data and assigning corresponding layers; the evidence factors include elevation, slope, slope aspect, fault distance, engineering rock group, water system distance, precipitation, road distance and land classification; S32, calculating the contribution weight of each evidence factor to the occurrence of geological hazards through conditional probability, assigning a comprehensive weight to each evidence factor layer, and then superimposing and summing the comprehensive weights of each evidence factor layer to obtain the susceptibility index of the grid unit; S33, assigning the susceptibility index to each evidence factor layer for weighted superposition to obtain an initial landslide susceptibility zoning map; the calculation expression is as follows: IN i =In i + -IN i - ; Where: W i + is the weight value of the evidence factor existence area; W i - is the weight value of the area where the evidence factor does not exist; W i is the weight value of the influence of the evidence factor on the occurrence of landslide; p is the conditional probability; F is the evidence factor; L is the number of landslide hazard units; is the number of non-landslide hazard units; (F|L) and are the number of units where the evidence factor exists in the landslide hazard unit and the non-landslide hazard unit; and are the numbers of units where the evidence factor does not exist in the landslide hazard unit and the non-landslide hazard unit, respectively.

4. The method for preparing a landslide susceptibility zoning map according to claim 1, characterized in that: In step S4, obtaining and updating the surface deformation rate map specifically includes: S41, downsampling the initial surface deformation rate map to a grid unit consistent with the initial landslide susceptibility zoning map based on the nearest neighbor method; S42, dividing the initial surface deformation rate map into a plurality of deformation rate intervals according to the surface deformation rate, and classifying and counting the pixel points in different deformation rate intervals to obtain the deformation severity of the different deformation rate intervals; S43. Assigning a contribution to the initial surface deformation rate map based on the deformation severity to obtain an updated surface deformation rate map.

5. The method for preparing a landslide susceptibility zoning map according to claim 4, characterized in that: In step S4, obtaining and updating the landslide susceptibility zoning map specifically includes: S44. Assign different contribution degrees to the initial landslide susceptibility zoning map according to the hazard level to obtain an updated landslide susceptibility zoning map.

6. The method for preparing a landslide susceptibility zoning map according to claim 5, characterized in that: In step S4, obtaining the final landslide susceptibility zoning map of the geological disaster prone area specifically includes: S45, calculating the contribution of a certain area on the updated landslide susceptibility zoning map and the updated surface deformation rate map to obtain an updated landslide susceptibility index LSI; the calculation formula is as follows: LSI=C LSM *C GDM ; Where: C LSM Contribution value for updating the landslide susceptibility zoning map; C GDM Contribution value for updating the surface deformation rate map; S46. Zoning the updated landslide susceptibility index and assigning hazard levels to different areas respectively, thereby obtaining a final landslide susceptibility zoning map for geological disaster-prone areas.

7. A landslide susceptibility zoning map making system, characterized in that: The system is applied to the method described in any one of claims 1 to 6, and the system comprises: A geographic data acquisition module (1) is used to acquire geographic data of areas prone to geological disasters; the geographic data includes historical landslide data, three-land survey data, geological factor data, road data, digital elevation model and radar data; An initial surface deformation rate map construction module (2) is used to perform preprocessing based on radar data to obtain an initial surface deformation rate map of a geological disaster prone area; An initial landslide susceptibility zoning map construction module (3) is used to generate an initial landslide susceptibility zoning map of a geological disaster prone area based on historical landslide data and geological factor data and using an evidence weight method for weighting; A coupling update module (4) is used to assign values ​​to and update the initial surface deformation rate map and the initial landslide susceptibility zoning map based on the contribution matrix, obtain an updated surface deformation rate map and an updated landslide susceptibility zoning map, and couple the two to obtain a final landslide susceptibility zoning map of the geological disaster prone area.

8. A device for making a landslide susceptibility zoning map, characterized in that: The device comprises a processor (5) and a memory (6); The memory (6) is used to store computer program code (61) and transmit the computer program code (61) to the processor (5); The processor (5) is used to execute the method for preparing a landslide susceptibility zoning map according to any one of claims 1 to 6 according to the instructions in the computer program code (61).

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer-executable instructions, and when the computer-executable instructions are executed on a computer, the method for preparing a landslide susceptibility zoning map according to any one of claims 1 to 6 is implemented.

Images