CN112200355B - Landslide prediction method, landslide prediction device, landslide prediction equipment and storage medium - Google Patents
Landslide prediction method, landslide prediction device, landslide prediction equipment and storage medium Download PDFInfo
- Publication number
- CN112200355B CN112200355B CN202011063199.4A CN202011063199A CN112200355B CN 112200355 B CN112200355 B CN 112200355B CN 202011063199 A CN202011063199 A CN 202011063199A CN 112200355 B CN112200355 B CN 112200355B
- Authority
- CN
- China
- Prior art keywords
- landslide
- grid
- size
- rainfall
- grids
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 74
- 238000003860 storage Methods 0.000 title claims abstract description 14
- 230000015654 memory Effects 0.000 claims description 13
- 230000003203 everyday effect Effects 0.000 claims description 11
- 230000003068 static effect Effects 0.000 claims description 11
- 238000000605 extraction Methods 0.000 claims description 10
- 230000002354 daily effect Effects 0.000 claims description 7
- 239000002689 soil Substances 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 239000013049 sediment Substances 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 description 22
- 230000006870 function Effects 0.000 description 17
- 238000012549 training Methods 0.000 description 14
- 238000004422 calculation algorithm Methods 0.000 description 11
- 238000009826 distribution Methods 0.000 description 10
- 238000012544 monitoring process Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 10
- 238000013145 classification model Methods 0.000 description 8
- 238000005070 sampling Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000012937 correction Methods 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 6
- 210000002569 neuron Anatomy 0.000 description 6
- 238000012952 Resampling Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000010276 construction Methods 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 4
- 238000003066 decision tree Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000010606 normalization Methods 0.000 description 4
- 238000012876 topography Methods 0.000 description 4
- 238000012795 verification Methods 0.000 description 4
- 238000003062 neural network model Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 108010014173 Factor X Proteins 0.000 description 2
- 230000004913 activation Effects 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 210000004027 cell Anatomy 0.000 description 2
- 238000003708 edge detection Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- NIWYLNFSJJLYHU-UHFFFAOYSA-N n-ethyl-1-(2-fluorophenyl)propan-2-amine Chemical compound CCNC(C)CC1=CC=CC=C1F NIWYLNFSJJLYHU-UHFFFAOYSA-N 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000000306 recurrent effect Effects 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000013527 convolutional neural network Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000007477 logistic regression Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000007637 random forest analysis Methods 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000006403 short-term memory Effects 0.000 description 1
- 238000004513 sizing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Computational Linguistics (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the invention discloses a landslide prediction method, device, equipment and storage medium, which are used for acquiring a first grid size of a target area and a second grid size resampled by the first grid size, dividing the second grid size into a plurality of grids, extracting characteristic information of each divided grid, and improving the accuracy of landslide prediction by taking the landslide occurrence probability as a landslide prediction result of the target area.
Description
Technical Field
The embodiment of the invention relates to the technical field of landslide monitoring, in particular to a landslide prediction method, a landslide prediction device, landslide prediction equipment and a storage medium.
Background
Landslide is one of the most common disastrous natural disasters, has the characteristics of wide distribution range, high occurrence frequency, multiple occurrence, regional property, severity and the like, and can cause a large amount of casualties and serious environmental and infrastructure losses every year. The method has important significance in evaluating the liability of landslide.
The existing landslide susceptibility prediction can be divided into a deterministic method and a non-deterministic method according to the difference of theoretical basis on which the landslide susceptibility prediction is based. The deterministic method is mainly a directional analysis based on expert experience and knowledge and a landslide process or physical model analysis method, and the prediction accuracy is poor. With the rapid development of computer technology and 3S technology in recent years, non-deterministic methods have been widely used, mainly including fuzzy logic methods, analytic hierarchy processes, decision trees, and the like. However, the time precision of the landslide factor processed by the method is poor, especially the rainfall factor, only one year of rainfall is usually considered, so that the prediction precision is not ideal, meanwhile, due to the fact that the influence factors of landslide are numerous, the prediction efficiency of prediction by only adopting one non-determining method is poor, and the accuracy of landslide prediction by directly utilizing landslide data resampled from one grid size to another grid size is poor due to the fact that some characteristic differences of landslide data under grids of different sizes are large.
Disclosure of Invention
The invention provides a landslide prediction method, device, equipment and storage medium, which realize feature re-extraction of landslide data of resampled grids, and landslide prediction is carried out according to the extracted feature information and landslide data of each grid every day, so that the prediction precision is improved, and meanwhile, the landslide prediction is carried out based on a landslide prediction model, so that the prediction efficiency is improved.
In a first aspect, an embodiment of the present invention provides a method for predicting a landslide, where the method for predicting a landslide includes:
acquiring a first grid size of a target area, and determining a second grid size resampled by the first grid size, wherein the second grid size is larger than the first grid size;
dividing an initial grid corresponding to the second grid size into a plurality of grids according to the first grid size, and extracting characteristic information of landslide influence factors in the grids corresponding to the second grid size based on the first grid size, wherein the landslide influence factors are daily landslide data in a set time period of each grid under the grids corresponding to the second grid size;
and determining the occurrence probability of each grid landslide according to the characteristic information and the landslide influence factor based on a landslide prediction model, and taking the occurrence probability as a landslide prediction result of the target area.
In a second aspect, an embodiment of the present invention further provides a device for predicting a landslide, where the device for predicting a landslide includes:
a grid size determining module, configured to obtain a first grid size of a target area, and determine a second grid size resampled by the first grid size, where the second grid size is greater than the first grid size;
The characteristic information extraction module is used for dividing the initial grid corresponding to the second grid size into a plurality of grids according to the first grid size, and extracting characteristic information of landslide influence factors in the grids corresponding to the second grid size based on the first grid size, wherein the landslide influence factors are daily landslide data in a set time period of each grid under the grids corresponding to the second grid size;
and the landslide prediction module is used for determining the occurrence probability of landslide of each grid according to the characteristic information and the landslide influence factor based on a landslide prediction model, and taking the occurrence probability as a landslide prediction result of the target area.
In a third aspect, an embodiment of the present invention further provides a landslide prediction apparatus, including:
one or more processors;
a memory for storing one or more programs;
when the one or more programs are executed by the one or more processors, the one or more processors implement the landslide prediction method provided by any embodiment of the present invention.
In a fourth aspect, embodiments of the present invention also provide a storage medium containing computer-executable instructions, which when executed by a computer processor, are used to perform the method of predicting landslide provided by any of the embodiments of the present invention.
According to the technical scheme, the first grid size of the target area and the second grid size resampled by the first grid size are obtained, the second grid size is divided into a plurality of grids, and the characteristic information of each divided grid is extracted, so that the balance of the grids, the richness of landslide influence factors and the advantage of small calculated amount can be simultaneously considered, and in the process of predicting landslide probability in a follow-up mode, the landslide influence factors are input into a landslide prediction model based on the characteristic information and the landslide influence factors and in a unit of day, the landslide occurrence probability is predicted, and the landslide occurrence probability is used as a landslide prediction result of the target area, so that the landslide prediction accuracy and the landslide prediction efficiency are improved.
Drawings
FIG. 1 is a flow chart of a method for predicting landslide in accordance with a first embodiment of the invention;
FIG. 2 is a flow chart of a method for predicting landslide in accordance with a second embodiment of the present invention;
FIG. 3 is a flowchart of a landslide prediction method in accordance with a third embodiment of the present invention;
fig. 4 is a schematic structural diagram of a landslide prediction device according to a fourth embodiment of the present invention;
fig. 5 is a schematic structural diagram of a landslide prediction apparatus in a fifth embodiment of the invention.
Detailed Description
The invention is described in further detail below with reference to the drawings and examples. It is to be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention. It should be further noted that, for convenience of description, only some, but not all of the structures related to the present invention are shown in the drawings.
Example 1
Fig. 1 is a flowchart of a landslide prediction method according to an embodiment of the present invention, where the method may be performed by a landslide prediction device, and the method includes the following steps:
s110, acquiring a first grid size of the target area, and determining a second grid size resampled by the first grid size.
The target area is usually an area where landslide occurs, and may be any designated area. Raster data is a form of data that divides space into regular grids, each grid being a grid or cell, and assigns corresponding attribute values to each cell to represent an entity. The second grid size is greater than the first grid size.
Optionally, the method for determining the first grid size includes: determining an initial size and an initial grid number of each grid of the target area, and determining a landslide grid number and a non-landslide grid number, wherein the landslide grid number is the number of grids on which landslide occurs, and the non-landslide grid number is the number of grids on which landslide does not occur; if the ratio of the landslide grid number to the non-landslide grid number is in a set balance interval, taking the initial size as the first grid size; if the ratio of the landslide grid number to the non-landslide grid number is not in the set balance interval, adjusting the size of each grid of the target area, and re-determining the landslide grid number, the non-landslide grid number and the current grid number; and if the ratio of the redetermined landslide grid number to the non-landslide grid number is in the set balance interval and the ratio of the current grid number to the initial grid number is smaller than a set threshold value, taking the size of the grid after adjustment as the first grid size.
The initial dimension may be any dimension, for example 1km by 1km. The initial number of grids is the number of grids corresponding to the initial size. It can be understood that when the target area is rasterized, landslide of the target area falls into the grids, the grids where landslide occurs and the grids where landslide does not occur are determined according to the landslide position and the landslide area, and the number of landslide grids and the number of non-landslide grids are determined. Further, calculating the ratio of the number of landslide grids to the number of non-landslide grids, if the ratio of the number of landslide grids to the number of non-landslide grids is in a set balance interval, indicating that the number of landslide grids in the grids corresponding to the initial size is balanced with the number of non-landslide grids, and taking the initial size as the first grid size; if the ratio of the number of landslide grids to the number of non-landslide grids is not in a set balance interval, the size of each grid is adjusted, the number of landslide grids, the number of non-landslide grids and the current grid number are redetermined, if the ratio of the number of redetermined landslide grids to the number of non-landslide grids is in the set balance interval, the number of landslide grids and the number of grids are in a balance state, if the ratio of the current grid number to the number of initial grids is smaller than a set threshold, the total number of landslide grids in a target area is in a balance state, and the size of each grid after adjustment can be used as a first grid size.
Optionally, the present embodiment may further determine the first grid size by: acquiring rainfall data of each rainfall station of the target area, an initial size and an initial grid number of each grid of the target area, and determining a landslide grid number and a non-landslide grid number, wherein the landslide grid number is the number of grids with landslide, and the non-landslide grid number is the number of grids without landslide; calculating rainfall predicted values of the rainfall stations corresponding to any current grid according to rainfall data of the rainfall stations with a specific number, and calculating rainfall errors based on the rainfall predicted values and actually measured rainfall values of the rainfall stations corresponding to the current grid; if the rainfall error is smaller than a set error threshold, and the ratio of the landslide grid number to the non-landslide grid number is in a set balance interval, taking the initial size as the first grid size; if the rainfall error is larger than a set error threshold, adjusting the size of each grid of the target area, and recalculating the rainfall error according to the rainfall predicted value and the actually measured rainfall value of the rainfall station corresponding to any grid after adjustment; and if the recalculated rainfall error is smaller than the set error threshold, the ratio of the redetermined landslide grid number to the redetermined non-landslide grid number is in a set balance interval, the ratio of the current grid number to the initial grid number is smaller than the set threshold, and the size of the grid after adjustment is taken as the first grid size.
Specifically, a plurality of monitoring points can be set in the target area, and each monitoring point is provided with a rainfall station so as to acquire rainfall data of each monitoring point in real time. Specifically, the target area is rasterized, a plurality of rainfall stations are distributed in different grids, when rainfall of a certain monitoring point is predicted, rainfall of a plurality of other monitoring points can be obtained, the rainfall of the monitoring point is determined based on a collaborative Kriging interpolation method and by taking elevation information as auxiliary information, the rainfall is used as predicted rainfall of the monitoring point, the predicted rainfall and the actually measured rainfall obtained by the rainfall stations of the monitoring point are subjected to rainfall error, and if the rainfall error is smaller than a set error threshold, the predicted rainfall of the rainfall station in the grids under the initial size is more accurate; meanwhile, calculating the ratio of the number of landslide grids to the number of non-landslide grids, judging whether the ratio is in a set balance interval, if so, indicating that the number of landslide grids corresponding to the initial size is balanced with the number of non-landslide grids, and taking the initial size as a first grid size; if the rainfall error is not less than the set error threshold, the predicted rainfall of the rainfall station in the grids with the initial size is inaccurate, the size of each grid needs to be adjusted, the rainfall error is recalculated according to the rainfall predicted value and the actually measured rainfall value of the rainfall station in the grids with the adjusted grid size until the recalculated rainfall error is less than the set error threshold, if the ratio of the redetermined landslide grid number to the redetermined non-landslide grid number is also in the set balance interval and the ratio of the current grid number to the initial grid number is less than the set threshold, the actually measured rainfall of the rainfall station of the adjusted grid is accurate, the landslide grid number corresponding to the adjusted grid size is balanced with the non-landslide grid number, and the adjusted total grid number is balanced, and the size of the adjusted grid is taken as the first grid size.
The second grid size in this embodiment may also be determined in the above manner, or the initialization grid size is taken as the second grid size, and the specific determination manner of the second grid size is not described in detail.
S120, dividing the initial grid corresponding to the second grid size into a plurality of grids according to the first grid size, and extracting characteristic information of landslide impact factors in the grids corresponding to the second grid size based on the first grid size.
And the landslide influence factor is the landslide data of each day in a set time period of each grid under the grid corresponding to the second grid size. The landslide impact factors mainly include various landslide impact factors of landslide, such as elevation factors and environmental factors. The set period of time may be one day, three days, one week, one month, or other period of time. The landslide impact factor comprises a dynamic factor and a static factor, wherein the dynamic factor comprises at least one of rainfall and soil humidity, and the static factor comprises at least one of elevation, slope direction, planar curvature, profile curvature, topography humidity index, water current intensity index, sediment transport index, topography roughness index, distance to fault, distance to river, distance to road, lithology, land utilization, vegetation coverage. The characteristic information includes: at least one of a weighted value, a maximum value, a minimum value, a differential mean value and a standard deviation characteristic of the landslide impact factor.
Optionally, the feature information extracting method includes: acquiring characteristic values of landslide impact factors in eight neighbor grids of the current grid; and determining the characteristic information of the current grid according to the characteristic values in the eight-neighborhood grid. For example, if the characteristic information of the current grid is gradient, the gradient and the gradient height in the eight-neighborhood grid of the current grid can be obtained, the acute angle formed by the current grid and any eight-neighborhood grid is calculated based on the gradient and the gradient length of the eight-neighborhood grid of the current grid by adopting the Pythagorean theorem, and the gradient of the current grid is determined based on the acute angle formed by the current grid and each eight-neighborhood grid. Alternatively, the gradients of other eight neighbor grids can be calculated in the same manner as described above until the gradient calculation of all grids is completed, so as to obtain the characteristic information of the landslide impact factor under the grids of the first grid size.
It can be understood that under different grid sizes, the characterization degree of landslide impact factors is different, the grids with smaller sizes can keep the richness of the characteristic information, and the grids with larger sizes can reduce the redundancy of data. After resampling the landslide impact factor of the target area from the first grid size to the second grid size, some characteristics of the landslide impact factor at the second grid size differ significantly from the landslide impact factor at the first grid size (e.g., slope and slope direction in the static factor). In order to simultaneously consider the balance of grids, the richness of landslide influence factors and the advantage of small calculated amount, in the embodiment, after resampling the landslide influence factors, grids corresponding to the second grid size after resampling are divided into a plurality of grids according to the first grid size before resampling, and feature information of the landslide influence factors in the grids corresponding to the second grid size is extracted based on the first grid size, so that feature re-extraction of the grids corresponding to the second grid size is realized. Before extracting the characteristic information of the landslide impact factor in the grids corresponding to the second grid size based on the first grid size, the difference between the landslide impact factors before and after resampling needs to be judged, if the difference is large, the characteristic is extracted again for the grids corresponding to the second grid size by adopting the mode of S120, and if the difference is small, the landslide impact factor under the first grid size or the second grid size is directly obtained.
Illustratively, the first grid size is 30m×30m, the second grid size is 1km×1km, the grids corresponding to the first grid size are resampled to the grids corresponding to the second grid size, the gradient and the slope direction are greatly different under the two grid sizes, the initial grids corresponding to the second grid size are divided into 1111 grids according to the first grid size, and at least one of the weighted value, the minimum value, the differential mean value and the standard deviation of the gradient and the slope direction is extracted under 1111 grids, so that the feature re-extraction of the grids corresponding to the second size is completed.
S130, determining occurrence probability of landslide of each grid based on the landslide prediction model according to the characteristic information and the landslide influence factor.
The landslide prediction model can be a single neural network model, or can be a complex neural network model formed by sequentially connecting an input end and an output end. Specifically, the training process of the landslide prediction model is as follows: basic information of landslide of a target area or all areas is extracted from files such as landslide field investigation reports, typical landslide monitoring reports and the like, the basic information comprises landslide information such as landslide occurrence time, longitude and latitude, disaster scale and the like and actual measurement rainfall of the landslide area, landslide influence factors of historical time periods of various landslide occurrence lands are determined according to the landslide information, the landslide influence factors of the historical time periods can comprise landslide influence factors and landslide characteristic information under different grid sizes, and the landslide influence factors of the historical time periods form a training set and a verification set according to a set proportion, such as 8:2 or 7:3; and carrying out data correction and registration on each landslide influence factor, carrying out coordinate unified processing and grid unified processing, thus obtaining landslide influence factors with consistent grid sizes, and training a landslide prediction model by using the landslide influence factors processed by the steps. The landslide influence factors are extracted through characteristic engineering to form an input characteristic matrix of a landslide prediction model; initializing parameters of a landslide prediction model, inputting the input feature matrix into the landslide prediction model, performing model training to obtain historical landslide information in a historical time period, adjusting the parameters of the landslide prediction model according to an evaluation result based on an F1-value (F1-Score) and an ROC (Receiver Operating Characteristic) as evaluation indexes, and performing model verification through a verification set when the parameters are satisfied, wherein the trained landslide prediction model is obtained after verification. Wherein the historical landslide information includes a probability of each grid landslide occurrence.
Optionally, before inputting the landslide prediction model according to the feature information and the landslide impact factor, the landslide impact factor and the feature information may be further preprocessed, where the preprocessing includes at least one of coordinate unification processing, correction processing, data discretization processing, and grid unification processing. According to the embodiment, the occurrence probability of each grid landslide is determined by combining the characteristic information and the landslide influence factors, the data size of the landslide influence factors can be increased, the occurrence probability of each grid landslide can be determined more accurately according to the characteristic information, and the accuracy of landslide prediction is improved.
According to the technical scheme provided by the embodiment, the first grid size of the target area and the second grid size resampled by the first grid size are obtained, the second grid size is divided into a plurality of grids, and the characteristic information of each divided grid is extracted, so that the balance of the grids, the richness of landslide influence factors and the advantage of small calculated amount can be simultaneously considered, and in the process of predicting landslide probability subsequently, the landslide influence factors are input into a landslide prediction model based on the characteristic information and the landslide influence factors and in a unit of day to predict the landslide occurrence probability, and the landslide occurrence probability is used as a landslide prediction result of the target area, so that the landslide prediction accuracy and the landslide prediction efficiency are improved.
Example two
Fig. 2 is a flowchart of a landslide prediction method according to a second embodiment of the present invention, where the previous embodiment is further refined. Optionally, the determining, based on the landslide prediction model, the occurrence probability of each grid landslide according to the feature information and the landslide influence factor includes: determining the occurrence probability of landslide every day according to the characteristic information and the landslide influence factor based on a first preset model; and if the occurrence probability of the current sky landslide is larger than a preset probability threshold value, determining the occurrence probability of each grid landslide based on the characteristic information and the landslide influence factor on the basis of a second preset model. For parts which are not described in detail in this method embodiment, reference is made to the above-described embodiments. Referring specifically to fig. 2, the method may include the steps of:
s210, acquiring a first grid size of the target area, and determining a second grid size resampled by the first grid size.
S220, dividing the initial grid corresponding to the second grid size into a plurality of grids according to the first grid size, and extracting characteristic information of landslide impact factors in the grids corresponding to the second grid size based on the first grid size.
And the landslide influence factor is the landslide data of each day in a set time period of each grid under the grid corresponding to the second grid size.
S230, determining the occurrence probability of landslide every day according to the characteristic information and the landslide influence factor based on the first preset model.
The first preset model may be a neural network model, or other learning algorithm. By way of example, the first preset model may be a support vector machine algorithm (Support Vector Machine, SVM), a Long Short Term Memory Network (LSTM), a logistic regression model (Logistics Regression, LR), an XGBoost (Extreme Gradient Boosting, extreme gradient boost decision tree) algorithm, a GBDT (Gradient Boosting Decision Tree, gradient boost decision tree) algorithm, a full convolution Network (Fully Convolutional Networks, FCN), a cyclic convolution Network (Recurrent Neural Network, RNN), a Residual Network (ResNet), a gated loop unit (Gate Recurrent Unit, GRU), and the like. The training method of the first preset model is the same as that of the landslide prediction model in the foregoing embodiment, and will not be described in detail herein.
Optionally, the method for determining the occurrence probability of the landslide every day is as follows: calculating construction characteristics of the dynamic factors of each grid every day, wherein the construction characteristics comprise the sum, average value, maximum value, minimum value, range, quartile and rainfall time of each dynamic factor of each grid, and obtaining a first characteristic matrix of the first preset model according to the construction characteristics and the static factors every day; and inputting the first feature matrix into the first preset model to determine the occurrence probability of landslide every day according to the first preset model.
As in the previous embodiments, the landslide impact factors include dynamic factors including rainfall (precipitation) and soil humidity, and static factors including elevation, slope direction, plane curvature, section curvature, topography humidity index, water current intensity index, sediment transport index, topography roughness index, distance to fault, distance to river, distance to road, lithology, land utilization, vegetation coverage, and the like.
Illustratively, the landslide impact factor X of each grid acquired daily is: x= { X (1) ,X (2) ,…,X (d) And, wherein the superscript indicates the number of days, X (c) C=1, 2, …, d, representing the landslide impact factor of each grid in the target area on day c, where matrix X (c) The rows of (1) represent a grid, the columns represent landslide factors, and the matrix X (c) Is m x n, i.e. comprises m grids, n landslide factors, where i=1, 2, … n 1 Represents a dynamic factor (rainfall, soil humidity, etc.), i=n 1 +1,n 1 +2, … n, represent the statics factor. Constructing feature set X in days 1 (first feature matrix) constructed as follows: for each X (c) Calculating n by grid 1 The sum, average, maximum, minimum, range, upper quartile, lower quartile and other structural features of each dynamic factor are used to obtain feature set X 1 . For training data, landslide tag vector Y may also be constructed 1 To indicate whether or not a landslide has occurred every day, wherein 1 indicates that a landslide has occurred and 0 indicates that no landslide has occurred.
S240, if the occurrence probability of the current sky landslide is larger than a preset probability threshold value, determining the occurrence probability of each grid landslide based on the second preset model according to the characteristic information and the landslide influence factor.
Optionally, the method for determining the occurrence probability of each grid landslide is as follows: feature extraction is carried out on landslide influence factors of grids corresponding to the current day so as to generate a second feature matrix of the second preset model; normalizing the second feature matrix; and inputting the normalized second feature matrix into the second preset model to determine the occurrence probability of each grid landslide in the current day according to the second preset model.
Specifically, assuming that the landslide occurrence probability on the c-th day is greater than the prediction probability threshold, the landslide influence factor X on the c-th day (c) Send to a second preset model to construct a second presetThe specific construction mode of the feature set of the model is as follows: and calculating the construction characteristics of the dynamic factors of each grid on the day, such as the corresponding sum, the mean, the variance, the median, the mean and the variance of the difference, the skewness, the kurtosis and the like, wherein the total sum, the mean, the variance, the median, the mean and the variance of the difference are respectively calculated on 3 days, 7 days, 15 days and 30 days, and the kurtosis is the early rainfall index calculated based on the rainfall attenuation index. Meanwhile, taking the characteristics of whether the current grid is the maximum value or the minimum value in the set range, whether the current grid exceeds the average value corresponding to the set range, whether the grid in the set range has landslide, the number of times of landslide of the current grid and the like into consideration in the set range taking the current grid as the center, such as the 3*3 area range, and finally obtaining a second characteristic matrix X of the second preset model 2 fea 。
Specifically, the second feature matrix X can be normalized based on max-min 2 fea Is normalized. Of course, other normalization algorithms can be selected for normalization processing. The embodiment of the invention does not limit the normalization algorithm of the first feature matrix and the second feature matrix. Further, after the normalization process, the method further includes: feature selection is carried out on the second feature matrix based on a multiple collinearity method so as to carry out feature selection according to the collinearity degree, and a screened second feature matrix is obtained
Specifically, the second preset model can be a one-dimensional convolutional neural network model, and comprises a convolutional layer, a batch standardization layer, an activation function and an optimization layer, wherein the value range of the number of convolutional kernels of the convolutional layer is 32-512, and the step length is 16; the activation function may include any one of a ReLU function (Rectified Linear Unit, linear rectification function), a Linear function (Linear function), a Sigmoid function, and a Tanh function (hyperbolic function); the optimization method comprises any one of SGD (Stochastic Gradient Descent, gradient update rule), adam (Adaptive Moment Estimation ), nadam (Nesterov Adaptive Moment Estimation, nesterov acceleration adaptive moment estimation), adagard (Adaptive Gradient Algorithm ), RMSprop (Root Mean Square Prop, forward root mean square gradient descent algorithm) and other optimization algorithms; the initial learning rate may take a value of 0.0001, 0.001, 0.01, or 0.1; the value range of the number of the hidden layer neurons is 4-256, the step length is 4, the value range of the number of the hidden layer neurons is 3-8, and the step length is 1; the random discarding rate is 0-0.8, and the step length is 0.05. The training method of the first preset model is the same as that of the landslide prediction model in the foregoing embodiment, and will not be described in detail herein.
According to the technical scheme provided by the embodiment, the first grid size of the target area and the second grid size resampled by the first grid size are obtained, the second grid size is divided into a plurality of grids, the characteristic information of each divided grid is extracted, landslide prediction is performed based on the characteristic information and landslide influence factors, the accuracy of data is improved, especially, the dynamic factors in landslide are predicted in units of days, the time sequence information of the dynamic factors is fully considered, and the prediction accuracy is improved; by setting two preset models, the landslide prediction is carried out in two stages, whether the probability of landslide occurrence in the current day is larger than a set value or not is predicted in the first stage, if so, the landslide prediction in the second stage is carried out, the landslide occurrence probability of a specific grid is carried out by taking the data in the current day as a reference, the landslide occurrence probability is taken as a landslide prediction result of a target area, the waste of calculation resources is greatly reduced, and the landslide prediction efficiency is improved.
Example III
Fig. 3 is a flowchart of a landslide prediction method according to a third embodiment of the present invention. The embodiment is used for analyzing the whole landslide prediction flow. For parts which are not described in detail in this method embodiment, reference is made to the above-described embodiments. Referring specifically to fig. 3, the method may include the steps of:
S310, acquiring a first grid size of the target area, and determining a second grid size resampled by the first grid size.
Optionally, before the first grid size of the target area is obtained, remote sensing image information and terrain information of the target area can be obtained, each functional area of the target area is determined based on the functional area classification model according to the remote sensing image information and the terrain information, and the target functional area is screened according to labels corresponding to each functional area.
The remote sensing image information refers to satellite images and can be obtained through land reflectivity products of Landsat (terrestrial satellite) series, and can reflect the land feature types of grids of a target area, such as grassland types, woodland types, rice types and building types; the terrain information may be obtained from a digital elevation model (Digital Elevation Model, DEM) of the target area, reflecting the maximum degree of variation of the grid information of a certain grid and surrounding grids of the target area. The terrain information may include elevation information, gradient information, curvature information, etc., which may reflect a terrain feature Gao Chengdian of a certain grid of the target area, elevation points near an important geographic target, or key elevation points on a distribution range, the gradient information may reflect a degree of steepness of a certain grid of the target area, and the curvature information may reflect a concave-convex condition of a certain grid of the target area.
Optionally, the method for determining each functional area of the target area includes: preprocessing the remote sensing image information, wherein the preprocessing comprises at least one of atmospheric correction, radiation correction and edge detection processing; and respectively extracting features of the preprocessed remote sensing image information and the preprocessed topographic information, inputting the features of the remote sensing image and the features of the topographic information into the functional area classification model, and determining each functional area of the target area, wherein the features of the remote sensing image comprise the mean value, the maximum value and the minimum value of the remote sensing image information, the topographic information comprises the elevation information, the gradient information and the curvature information, and the features of the topographic information comprise the mean value, the maximum value and the minimum value of the elevation information, the gradient information and the curvature information.
The functional area classification model is a random forest model, and the training method of the functional area classification model comprises the following steps: acquiring an initial forest model; sample remote sensing image information and sample terrain information of each functional area in a historical time period are input into the initial forest model, and a prediction label and the probability of the prediction label of each functional area are determined; and adjusting parameters of the initial forest model based on the prediction label, the probability of the prediction label and the sample label of each functional area until the prediction label is consistent with the sample label and the probability of the prediction label reaches a set threshold value, so as to obtain the functional area classification model. The probability refers to the probability of a predicted label of the output of the initial forest model, parameters of the initial forest model are adjusted through iteration, if the predicted label is consistent with the sample label and the predicted label probability reaches a set threshold, the probability that the predicted label is the sample label is high, and the initial forest model under the iteration number is used as a functional class classification model. The set threshold may be a larger value, for example, the set threshold takes a value of 0.9. Optionally, before the sample remote sensing image information is input into the initial forest model, the sample remote sensing image information may be preprocessed, for example, atmospheric correction, radiation correction, edge detection, false color synthesis processing, etc. are performed on the sample remote sensing image information, so as to improve the training accuracy of the functional area classification model of the sample remote sensing information.
It can be understood that each functional area output by the functional area classification model comprises all functional areas of a target area, including functional areas of town buildings, urban green lands, water bodies, farmlands, bare soil, mountain forests and the like, and landslide occurrence probabilities of the functional areas are different. For example, urban buildings and urban green lands have extremely low probability of landslide, and water bodies and mountain forests have relatively high probability of landslide. By comparing the historical prediction probability corresponding to each functional area with the historical probability threshold, the target functional area is screened out, and the target functional area is subjected to targeted landslide prediction, so that the calculated amount can be reduced, and the landslide prediction efficiency of the target area is improved.
S320, dividing the initial grid corresponding to the second grid size into a plurality of grids according to the first grid size, and extracting characteristic information of landslide impact factors in the grids corresponding to the second grid size based on the first grid size.
If each functional area of the target area is determined in S310 and the target functional area is screened according to the label corresponding to each functional area, S320 may be replaced by: dividing an initial grid of a target functional area corresponding to the second grid size into a plurality of grids according to the first grid size, and extracting characteristic information of landslide influence factors in the grids corresponding to the second grid size based on the first grid size.
The landslide impact factor includes a dynamic factor and a static factor, wherein the dynamic factor includes at least one of a rainfall and a soil humidity, and the rainfall can be determined by inputting the obtained geographical environment data of the target area into a rainfall interpolation model. The geographic environment data includes geographic location data, atmospheric data, terrain data, and underlying data. The rainfall interpolation model can be a back propagation model (BP, back propagation neural network), and is a multi-layer feedforward network consisting of nonlinear transformation units based on an error back propagation algorithm, wherein the BP generally consists of an input layer, an implicit layer and an output layer, each layer further comprises N neurons, the neurons at the same layer are mutually independent, and the output of the neurons between each layer only affects the input of the neurons at the lower layer after a specific excitation function.
Specifically, the training method of the rainfall interpolation model comprises the following steps: acquiring an initial model, determining an initial weight matrix and an initial threshold value of the initial model, and calculating the fitness of the initial model; sample geographic environment information in a historical time period is input into the initial model, predicted rainfall is determined, and according to the predicted rainfall and the actually measured rainfall in the historical time period, an initial weight matrix and an initial threshold value of the initial model are iteratively adjusted based on a genetic algorithm; and adjusting the initial model and calculating the fitness of the adjusted model based on the iteratively adjusted weight matrix threshold value until the recalculated fitness reaches an expected value, and taking the model corresponding to the fitness reaching the expected value as the rainfall interpolation model.
Specifically, the initial weight matrix may include a connection weight between an input layer and an implicit layer and a connection weight between the implicit layer and an output layer;the initial threshold may include a threshold of an hidden layer and a threshold of an output layer. The fitness of the initial model is calculated by the following formula:where n is the number of samples, y (i) is the actual output of the ith sample of the initial model, and t (i) is the expected output of the ith sample of the initial model. Specifically, sample geographic environment information in a historical time period is input into an initial model to obtain a predicted rainfall, a root mean square error between the predicted rainfall and an actually measured rainfall is calculated, if the root mean square error is larger than a preset error threshold, weight coding and threshold coding of the initial model are iteratively adjusted, an adjusted weight matrix and a threshold are determined based on the adjusted weight coding and threshold coding until the root mean square error is not smaller than the preset error threshold, the initial model is adjusted to be in a stable state, the fitness of the adjusted model is calculated by adopting the fitness calculation formula, and if the recalculated fitness reaches a desired value, a model corresponding to the fitness reaching the desired value is used as the rainfall interpolation model.
Optionally, the weight and the threshold are coded by adopting a binary coding, real number coding or Gray code coding mode, and the adjusted weight matrix and the adjusted threshold are determined according to the weight coding and the threshold coding. Optionally, the initial model may be adjusted to a stable state by calculating an average error, an average absolute error or a linear correlation coefficient between the predicted rainfall and the actually measured rainfall, and iteratively adjusting the weight code and the threshold code of the initial model according to the average error, the average absolute error or the linear correlation coefficient until the average error, the average absolute error or the linear correlation coefficient is not less than a preset error threshold, and further calculating the fitness of the adjusted model and determining the rainfall interpolation model.
It should be noted that, the sample geographic environment data in the historical time period is multidimensional, and when the rainfall interpolation model is trained, the precision and the robustness of the rainfall interpolation model can be improved by the multidimensional sample geographic environment data, so that after the daily geographic environment data is obtained, the prediction precision of the rainfall probability can be improved based on the rainfall interpolation model prediction, and the prediction precision of landslide prediction can be improved.
S330, determining the occurrence probability of landslide every day according to the characteristic information and the landslide influence factor based on the first preset model.
Unlike the previous embodiments, the first preset model may be a self-classifying learning model that iteratively adjusts a self-step factor and a downsampling ratio determination based on historical landslide data, historical non-landslide data. The self-step factors are determined according to the number of sub-boxes and the iteration times of the historical non-landslide data, and the downsampling proportion is determined according to the self-step factors of each sub-box and the self-step factors of all sub-boxes. Optionally, the method for determining the historical landslide data and the historical non-landslide data includes: determining labels corresponding to grids of a history area, and determining historical rainfall and specific correlation factors corresponding to landslide grids and non-landslide grids at each history time point respectively; and generating a landslide data set comprising the historical rainfall and the specific correlation factor respectively corresponding to the landslide grids and the non-landslide grids based on the labels, and determining the historical landslide data and the historical non-landslide data according to the labels corresponding to the grids in the landslide data set.
Optionally, the method for determining the historical landslide data and the historical non-landslide data includes: determining labels corresponding to grids of a history area, and determining historical rainfall and specific correlation factors corresponding to landslide grids and non-landslide grids at each history time point respectively; and generating a landslide data set comprising the historical rainfall and the specific correlation factor respectively corresponding to the landslide grids and the non-landslide grids based on the labels, and determining the historical landslide data and the historical non-landslide data according to the labels corresponding to the grids in the landslide data set.
Optionally, the training method of the first preset model includes: acquiring an initial prediction model, and acquiring the rainfall of a landslide grid and the rainfall of a non-landslide grid at each historical landslide time point, and specific correlation factors under each label; and inputting the rainfall of the landslide grids at the historical landslide time points, the rainfall of the non-landslide grids and the specific correlation factors under the labels into the initial prediction model according to the days, and adjusting a loss function of the initial prediction model based on the landslide probability output by the initial prediction model and the landslide probability corresponding to the historical landslide time points until the loss function reaches a set threshold value to obtain the first preset model.
The step of adjusting the loss function of the initial prediction model based on the landslide probability output by the initial prediction model and the landslide probability corresponding to the historical landslide time point until the loss function reaches a set threshold value, and obtaining the first preset model comprises the following steps: determining an initial classification hardness of the initial prediction model based on a rainfall of the landslide grid at the historical landslide time point, a rainfall of the non-landslide grid, each of the labels and a specific correlation factor under each of the labels; determining the rainfall capacity of the non-landslide grids and the number of the sub-boxes of the specific correlation factors of the non-landslide grids according to the initial classification hardness, and determining the self-step factors of the sub-boxes of the initial prediction model based on the number of the sub-boxes; determining a downsampling proportion of each sub-box based on the self-step factors, determining the rainfall capacity of downsampled non-landslide grids in each sub-box based on the downsampling proportion, and determining specific correlation factors of the downsampled non-landslide grids; inputting the rainfall of the downsampled non-landslide grids, the specific correlation factors of the downsampled non-landslide grids and the rainfall of the landslide grids into the initial prediction model, determining the loss function based on the landslide probability output by the initial prediction model and the landslide probability corresponding to the historical landslide time point, and iteratively adjusting the self-walking factors and the downsampling proportion of each sub-bin based on the training classification hardness of the loss function on a single sample; and adjusting the initial prediction model according to the self-step factors adjusted in an iteration mode and the downsampling proportion of each sub-bin until the loss function reaches a set threshold value, and obtaining the first preset model.
The downsampling proportion of each sub-box is the ratio of the self-step factor of each sub-box to the self-step factors of all sub-boxes, and the self-step factors are obtained by adding 1 to the sum of the number of sub-boxes, the opposite number of sub-box labels and the iteration number.
The landslide prediction method based on the self-step classification learning model has the advantages that: the self-step factors of the self-step classification learning model are determined according to the number of sub-boxes of the historical non-landslide data and the iteration times, and the down-sampling proportion is determined according to the self-step factors of each sub-box and the self-step factors of all sub-boxes, so that when the self-step factors and the down-sampling proportion are adjusted based on the historical landslide data and the historical non-landslide data in an iteration way, the sampling proportion of each sub-box gradually changes uniformly from the lower down-sampling number of the low hardness to the last lower down-sampling number of each sub-box along with the increase of the iteration times, and the lower down-sampling proportion of each box is uniformly changed, and meanwhile, the lower down-sampling number of the sub-boxes with small hardness can be always higher than the lower down-sampling number of the sub-boxes with large hardness, so that the first preset model with high diversity, robustness and strong inclusion is obtained.
S340, if the occurrence probability of the current sky landslide is larger than a preset probability threshold value, determining the occurrence probability of each grid landslide based on the second preset model according to the characteristic information and the landslide influence factor.
Alternatively, the second preset model may be a self-organizing learning model. The determination manner of the historical landslide data for training the second preset model and the training method of the second preset model can refer to S330, and this step will not be described in detail.
Optionally, after determining the occurrence probability of each grid landslide in the current day, the occurrence probability of each grid in the current day may be compared with an early warning threshold of the target area, and the occurrence level of the landslide in the target area is determined based on the obtained early warning level corresponding to the occurrence probability of each grid, where the early warning threshold is determined according to the first coefficient and the second coefficient of the grid landslide in each area.
The method for determining the first coefficient and the second coefficient comprises the following steps: acquiring historical probability of landslide occurrence in a set time period of grid landslide of each area; determining the first determined probability interval and the second determined probability interval based on the distribution characteristics of the historical probabilities; respectively calculating the sum of products of probability threshold values and corresponding coefficients of the same risk level in the first determined probability interval and the second determined probability interval to obtain intermediate determined probabilities; and if the intermediate determined probability does not reach the set evaluation index, iteratively adjusting the corresponding coefficient of the probability threshold value of each risk level until the intermediate determined probability reaches the set evaluation index, and determining a first coefficient and a second coefficient according to the corresponding coefficient reaching the set evaluation index.
The method for determining the early warning threshold value comprises the following steps: determining a first prediction probability interval and a second prediction probability interval based on the distribution characteristics of occurrence probability of each grid landslide; calculating a first product of each risk level in the first prediction probability interval and the first coefficient, calculating a second product of each risk level in the second prediction probability interval and the second coefficient, and taking the sum of the first product and the second product as an early warning threshold value of each risk level.
Wherein the determining the first determined probability interval and the second determined probability interval based on the distribution characteristics of the historical probabilities comprises: determining the historical probability and the times corresponding to the historical probability; determining the density distribution characteristics and breakpoint distribution characteristics of the historical probabilities according to the historical probabilities and the landslide occurrence times corresponding to the historical probabilities; the first determined probability interval is determined based on the intensity distribution characteristics and the second determined probability interval is determined based on the breakpoint distribution characteristics.
It can be understood that the above manner is based on the historical probability of the landslide occurring in the set time period of the grid landslide of each region, the first determined probability interval and the second determined probability interval are determined according to the distribution characteristics of the historical probability, the first coefficient and the second coefficient of each region are determined according to the probability threshold value of each risk level in the first determined probability interval and the second determined probability interval and the set evaluation index, and after the occurrence probability of different regions is obtained, the early warning threshold value of the different regions can be flexibly determined according to the first coefficient, the second coefficient and the occurrence probability corresponding to the different regions, so that the risk level of the region can be accurately determined according to the early warning threshold value corresponding to each region.
The determining the occurrence level of the landslide of the target area based on the obtained early warning level corresponding to the occurrence probability of each grid comprises the following steps: screening target grades larger than a first grade in the early warning threshold; calculating the average grade of the grids corresponding to the target grade; and determining the occurrence level of the landslide of the target area based on the average level, the set coefficient and the ratio of the grid number corresponding to the target level to all the grid numbers in the target area.
The calculation formula of the occurrence level of the landslide of the target area is as follows:
wherein alpha is a set coefficient, and is obtained by performing Bayesian calculation on the historical occurrence level of the landslide of the target area,and p is the ratio of the number of grids corresponding to the target grade to the number of all grids in the target zone.
Example IV
Fig. 4 is a schematic diagram of the result of a landslide prediction device according to a fourth embodiment of the present invention, where, as shown in fig. 4, the landslide prediction device includes: a grid sizing module 410, a feature information extraction module 420, and a landslide prediction module 430.
Wherein, the grid size determining module 410 is configured to obtain a first grid size of the target area, and determine a second grid size resampled by the first grid size, where the second grid size is greater than the first grid size;
The feature information extraction module 420 is configured to divide an initial grid corresponding to the second grid size into a plurality of grids according to the first grid size, and extract feature information of a landslide impact factor in the grid corresponding to the second grid size based on the first grid size, where the landslide impact factor is landslide data of each day in a set time period of each grid under the grid corresponding to the second grid size;
and the landslide prediction module 430 is configured to determine occurrence probability of landslide of each grid according to the feature information and the landslide influence factor based on a landslide prediction model, and take the occurrence probability as a landslide prediction result of the target area.
According to the technical scheme, the first grid size of the target area and the second grid size resampled by the first grid size are obtained, the second grid size is divided into a plurality of grids, and the characteristic information of each divided grid is extracted, so that the balance of the grids, the richness of landslide influence factors and the advantage of small calculated amount can be simultaneously considered, and in the process of predicting landslide probability in a follow-up mode, the landslide influence factors are input into a landslide prediction model based on the characteristic information and the landslide influence factors and in a unit of day, the landslide occurrence probability is predicted, and the landslide occurrence probability is used as a landslide prediction result of the target area, so that the landslide prediction accuracy and the landslide prediction efficiency are improved.
Optionally, the feature information includes: at least one of a weighted value, a maximum value, a minimum value, a differential mean value and a standard deviation characteristic of the landslide impact factor.
Optionally, the feature information extraction module 420 is further configured to obtain a feature value of a landslide impact factor in an eight-neighborhood grid of the current grid;
and determining the characteristic information of the current grid according to the characteristic values in the eight-neighborhood grid.
Optionally, the grid size determining module 410 is further configured to determine an initial size and an initial number of grids of each grid of the target area, and determine a landslide grid number and a non-landslide grid number, where the landslide grid number is a number of grids that have undergone landslide, and the non-landslide grid number is a number of grids that have not undergone landslide;
if the ratio of the landslide grid number to the non-landslide grid number is in a set balance interval, taking the initial size as the first grid size;
if the ratio of the landslide grid number to the non-landslide grid number is not in the set balance interval, adjusting the size of each grid of the target area, and re-determining the landslide grid number, the non-landslide grid number and the current grid number;
And if the ratio of the redetermined landslide grid number to the non-landslide grid number is in the set balance interval and the ratio of the current grid number to the initial grid number is smaller than a set threshold value, taking the size of the grid after adjustment as the first grid size.
Optionally, the grid size determining module 410 is further configured to obtain rainfall data of each rainfall station in the target area, an initial size and an initial grid number of each grid in the target area, and determine a landslide grid number and a non-landslide grid number, where the landslide grid number is a grid number where landslide occurs, and the non-landslide grid number is a grid number where no landslide occurs;
calculating rainfall predicted values of the rainfall stations corresponding to any current grid according to rainfall data of the rainfall stations with a specific number, and calculating rainfall errors based on the rainfall predicted values and actually measured rainfall values of the rainfall stations corresponding to the current grid;
if the rainfall error is smaller than a set error threshold, and the ratio of the landslide grid number to the non-landslide grid number is in a set balance interval, taking the initial size as the first grid size;
if the rainfall error is larger than a set error threshold, adjusting the size of each grid of the target area, and recalculating the rainfall error according to the rainfall predicted value and the actually measured rainfall value of the rainfall station corresponding to any grid after adjustment;
And if the recalculated rainfall error is smaller than the set error threshold, the ratio of the redetermined landslide grid number to the redetermined non-landslide grid number is in a set balance interval, the ratio of the current grid number to the initial grid number is smaller than the set threshold, and the size of the grid after adjustment is taken as the first grid size.
Optionally, the landslide prediction module 430 is further configured to determine, based on a first preset model, an occurrence probability of a landslide every day according to the feature information and the landslide impact factor;
and if the occurrence probability of the current sky landslide is larger than a preset probability threshold value, determining the occurrence probability of each grid landslide based on the characteristic information and the landslide influence factor on the basis of a second preset model.
Optionally, the landslide impact factor comprises a dynamic factor and a static factor, wherein the dynamic factor comprises at least one of rainfall and soil moisture, and the static factor comprises at least one of elevation, slope direction, planar curvature, profile curvature, terrain moisture index, water current intensity index, sediment transport index, terrain roughness index, distance to fault, distance to river, distance to road, lithology, land utilization, and vegetation coverage.
The landslide prediction device provided by the embodiment of the invention can execute the landslide prediction method provided by any embodiment of the invention, and has the corresponding functional modules and beneficial effects of the execution method.
Example five
Fig. 5 is a schematic structural diagram of a landslide prediction apparatus according to a fifth embodiment of the present invention, where, as shown in fig. 5, the apparatus includes a processor 510, a memory 520, an input device 530, and an output device 540; the number of device processors 510 may be one or more, one processor 510 being illustrated in fig. 5; the processor 510, memory 520, input means 530 and output means 540 in the device may be connected by a bus or other means, for example in fig. 5.
The memory 520 is a computer readable storage medium, and may be used to store a software program, a computer executable program, and a module, such as program instructions/modules corresponding to a landslide prediction method in an embodiment of the present invention (for example, the grid size determining module 410, the feature information extracting module 420, and the landslide prediction module 430 in a landslide prediction device). The processor 510 executes various functional applications of the apparatus and data processing, i.e., implements the landslide prediction method described above, by running software programs, instructions, and modules stored in the memory 520.
Memory 520 may include primarily a program storage area and a data storage area, wherein the program storage area may store an operating system, at least one application program required for functionality; the storage data area may store data created according to the use of the terminal, etc. In addition, memory 520 may include high-speed random access memory, and may also include non-volatile memory, such as at least one magnetic disk storage device, flash memory device, or other non-volatile solid-state storage device. In some examples, memory 520 may further include memory located remotely from processor 510, which may be connected to the device/terminal/server via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The input means 530 may be used to receive input numeric or character information and to generate key signal inputs related to user settings and function control of the device. The output 540 may include a display device such as a display screen.
Example six
A sixth embodiment of the invention also provides a storage medium containing computer-executable instructions, which when executed by a computer processor, are for performing a method of predicting a landslide, the method comprising:
Acquiring a first grid size of a target area, and determining a second grid size resampled by the first grid size, wherein the second grid size is larger than the first grid size;
dividing an initial grid corresponding to the second grid size into a plurality of grids according to the first grid size, and extracting characteristic information of landslide influence factors in the grids corresponding to the second grid size based on the first grid size, wherein the landslide influence factors are daily landslide data in a set time period of each grid under the grids corresponding to the second grid size;
and determining the occurrence probability of each grid landslide according to the characteristic information and the landslide influence factor based on a landslide prediction model.
Of course, the storage medium containing the computer executable instructions provided in the embodiments of the present invention is not limited to the method operations described above, and may also perform the related operations in the landslide prediction method provided in any embodiment of the present invention.
From the above description of embodiments, it will be clear to a person skilled in the art that the present invention may be implemented by means of software and necessary general purpose hardware, but of course also by means of hardware, although in many cases the former is a preferred embodiment. Based on such understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art in the form of a software product, which may be stored in a computer readable storage medium, such as a floppy disk, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), a FLASH Memory (FLASH), a hard disk or an optical disk of a computer, etc., and include several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the method according to the embodiments of the present invention.
It should be noted that, in the embodiment of the landslide prediction apparatus described above, each unit and module included are only divided according to the functional logic, but not limited to the above-described division, so long as the corresponding functions can be implemented; in addition, the specific names of the functional units are also only for distinguishing from each other, and are not used to limit the protection scope of the present invention.
Note that the above is only a preferred embodiment of the present invention and the technical principle applied. It will be understood by those skilled in the art that the present invention is not limited to the particular embodiments described herein, but is capable of various obvious changes, rearrangements and substitutions as will now become apparent to those skilled in the art without departing from the scope of the invention. Therefore, while the invention has been described in connection with the above embodiments, the invention is not limited to the embodiments, but may be embodied in many other equivalent forms without departing from the spirit or scope of the invention, which is set forth in the following claims.
Claims (8)
1. A method of predicting landslide, comprising:
acquiring a first grid size of a target area, and determining a second grid size resampled by the first grid size, wherein the second grid size is larger than the first grid size;
Dividing an initial grid corresponding to the second grid size into a plurality of grids according to the first grid size, and extracting characteristic information of landslide influence factors in the grids corresponding to the second grid size based on the first grid size, wherein the landslide influence factors are daily landslide data in a set time period of each grid under the grids corresponding to the second grid size;
determining occurrence probability of each grid landslide according to the characteristic information and the landslide influence factor based on a landslide prediction model;
the first grid size determining method comprises the following steps:
determining an initial size and an initial grid number of each grid of the target area, and determining a landslide grid number and a non-landslide grid number, wherein the landslide grid number is the number of grids on which landslide occurs, and the non-landslide grid number is the number of grids on which landslide does not occur;
if the ratio of the landslide grid number to the non-landslide grid number is in a set balance interval, taking the initial size as the first grid size;
if the ratio of the landslide grid number to the non-landslide grid number is not in the set balance interval, adjusting the size of each grid of the target area, and re-determining the landslide grid number, the non-landslide grid number and the current grid number;
If the ratio of the redetermined landslide grid number to the non-landslide grid number is in the set balance interval and the ratio of the current grid number to the initial grid number is smaller than a set threshold value, taking the size of the adjusted grid as the first grid size;
the extracting the characteristic information of the landslide impact factor in the grids corresponding to the second grid size based on the first grid size comprises the following steps:
acquiring characteristic values of landslide impact factors in eight neighbor grids of the current grid;
and determining the characteristic information of the current grid according to the characteristic values in the eight-neighborhood grid.
2. The prediction method according to claim 1, wherein the characteristic information includes: at least one of a weighted value, a maximum value, a minimum value, a differential mean value and a standard deviation characteristic of the landslide impact factor.
3. The prediction method according to claim 1, wherein the determining method of the first grid size includes:
acquiring rainfall data of each rainfall station of the target area, an initial size and an initial grid number of each grid of the target area, and determining a landslide grid number and a non-landslide grid number, wherein the landslide grid number is the number of grids with landslide, and the non-landslide grid number is the number of grids without landslide;
Calculating rainfall predicted values of the rainfall stations corresponding to any current grid according to rainfall data of the rainfall stations with a specific number, and calculating rainfall errors based on the rainfall predicted values and actually measured rainfall values of the rainfall stations corresponding to the current grid;
if the rainfall error is smaller than a set error threshold, and the ratio of the landslide grid number to the non-landslide grid number is in a set balance interval, taking the initial size as the first grid size;
if the rainfall error is larger than a set error threshold, adjusting the size of each grid of the target area, and recalculating the rainfall error according to the rainfall predicted value and the actually measured rainfall value of the rainfall station corresponding to any grid after adjustment;
and if the recalculated rainfall error is smaller than the set error threshold, the ratio of the redetermined landslide grid number to the redetermined non-landslide grid number is in a set balance interval, the ratio of the current grid number to the initial grid number is smaller than the set threshold, and the size of the grid after adjustment is taken as the first grid size.
4. The prediction method according to claim 1, wherein determining occurrence probability of each grid landslide based on the feature information and the landslide influence factor based on a landslide prediction model includes:
Determining the occurrence probability of landslide every day according to the characteristic information and the landslide influence factor based on a first preset model;
and if the occurrence probability of the current sky landslide is larger than a preset probability threshold value, determining the occurrence probability of each grid landslide based on the characteristic information and the landslide influence factor on the basis of a second preset model.
5. The prediction method according to claim 1, wherein the landslide impact factor comprises a dynamic factor and a static factor, wherein the dynamic factor comprises at least one of rainfall and soil humidity, and the static factor comprises at least one of elevation, slope direction, planar curvature, profile curvature, terrain humidity index, water current intensity index, sediment transport index, terrain roughness index, fault distance, river distance, road distance, lithology, land utilization, and vegetation coverage.
6. A landslide prediction apparatus comprising:
a grid size determining module, configured to obtain a first grid size of a target area, and determine a second grid size resampled by the first grid size, where the second grid size is greater than the first grid size;
The grid size determining module is further configured to determine an initial size and an initial grid number of each grid of the target area, and determine a landslide grid number and a non-landslide grid number, where the landslide grid number is a grid number where landslide occurs, and the non-landslide grid number is a grid number where landslide does not occur;
if the ratio of the landslide grid number to the non-landslide grid number is in a set balance interval, taking the initial size as the first grid size;
if the ratio of the landslide grid number to the non-landslide grid number is not in the set balance interval, adjusting the size of each grid of the target area, and re-determining the landslide grid number, the non-landslide grid number and the current grid number;
if the ratio of the redetermined landslide grid number to the non-landslide grid number is in the set balance interval and the ratio of the current grid number to the initial grid number is smaller than a set threshold value, taking the size of the adjusted grid as the first grid size;
the characteristic information extraction module is used for dividing the initial grid corresponding to the second grid size into a plurality of grids according to the first grid size, and extracting characteristic information of landslide influence factors in the grids corresponding to the second grid size based on the first grid size, wherein the landslide influence factors are daily landslide data in a set time period of each grid under the grids corresponding to the second grid size;
The characteristic information extraction module is also used for obtaining characteristic values of landslide influence factors in eight neighbor grids of the current grid;
determining the characteristic information of the current grid according to the characteristic values in the eight neighborhood grids;
and the landslide prediction module is used for determining the occurrence probability of landslide of each grid according to the characteristic information and the landslide influence factor based on a landslide prediction model, and taking the occurrence probability as a landslide prediction result of the target area.
7. A prediction apparatus for landslide, characterized in that the apparatus comprises:
one or more processors;
a memory for storing one or more programs;
when executed by the one or more processors, causes the one or more processors to implement the landslide prediction method of any one of claims 1-5.
8. A storage medium containing computer executable instructions which, when executed by a computer processor, are for performing a method of landslide prediction as claimed in any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011063199.4A CN112200355B (en) | 2020-09-30 | 2020-09-30 | Landslide prediction method, landslide prediction device, landslide prediction equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011063199.4A CN112200355B (en) | 2020-09-30 | 2020-09-30 | Landslide prediction method, landslide prediction device, landslide prediction equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112200355A CN112200355A (en) | 2021-01-08 |
| CN112200355B true CN112200355B (en) | 2024-01-05 |
Family
ID=74013561
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011063199.4A Active CN112200355B (en) | 2020-09-30 | 2020-09-30 | Landslide prediction method, landslide prediction device, landslide prediction equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112200355B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112785078A (en) * | 2021-02-02 | 2021-05-11 | 武汉中地云申科技有限公司 | Landslide prediction method and terminal combining soil moisture information and effective rainfall |
| CN117974927B (en) * | 2024-03-28 | 2024-06-25 | 煤炭科学研究总院有限公司 | Real-time early warning method for coal mine disasters of four-dimensional attribute geological model |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102915616A (en) * | 2012-05-11 | 2013-02-06 | 新疆大学 | Grid-accumulation-type regional land desertification early-warning method |
| KR101721114B1 (en) * | 2016-06-27 | 2017-03-30 | 서울대학교산학협력단 | Method for Determining the Size of Grid for Clustering on Multi-Scale Web Map Services using Location-Based Point Data |
| CN110362949A (en) * | 2019-07-23 | 2019-10-22 | 电子科技大学 | A kind of landslide sensitivity assessment method neural network based |
| CN111144651A (en) * | 2019-12-26 | 2020-05-12 | 杭州鲁尔物联科技有限公司 | Geological disaster prediction method, device and equipment |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7143098B2 (en) * | 2002-05-10 | 2006-11-28 | International Business Machines Corporation | Systems, methods, and computer program products to reduce computer processing in grid cell size determination for indexing of multidimensional databases |
-
2020
- 2020-09-30 CN CN202011063199.4A patent/CN112200355B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102915616A (en) * | 2012-05-11 | 2013-02-06 | 新疆大学 | Grid-accumulation-type regional land desertification early-warning method |
| KR101721114B1 (en) * | 2016-06-27 | 2017-03-30 | 서울대학교산학협력단 | Method for Determining the Size of Grid for Clustering on Multi-Scale Web Map Services using Location-Based Point Data |
| CN110362949A (en) * | 2019-07-23 | 2019-10-22 | 电子科技大学 | A kind of landslide sensitivity assessment method neural network based |
| CN111144651A (en) * | 2019-12-26 | 2020-05-12 | 杭州鲁尔物联科技有限公司 | Geological disaster prediction method, device and equipment |
Non-Patent Citations (1)
| Title |
|---|
| The influence of DEM spatial resolution on landslide susceptibility mapping in the Baxie River basin, NW China;Zhuo Chen et al.;Natural Hazards;853-877 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112200355A (en) | 2021-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112200358B (en) | Landslide level prediction method, landslide level prediction device, landslide level prediction equipment and storage medium | |
| CN112200362B (en) | Landslide prediction method, landslide prediction device, landslide prediction equipment and storage medium | |
| Soma et al. | Optimization of causative factors using logistic regression and artificial neural network models for landslide susceptibility assessment in Ujung Loe Watershed, South Sulawesi Indonesia | |
| CN112200354A (en) | Landslide prediction method, device, equipment and storage medium | |
| Zhao et al. | A Comparative Study of Landslide Susceptibility Mapping Using SVM and PSO‐SVM Models Based on Grid and Slope Units | |
| CN111784044A (en) | Landslide prediction method, device, equipment and storage medium | |
| Chu et al. | Comparison of landslide susceptibility maps using random forest and multivariate adaptive regression spline models in combination with catchment map units | |
| Collini et al. | Predicting and understanding landslide events with explainable AI | |
| CN111310898A (en) | Landslide hazard susceptibility prediction method based on RNN | |
| CN115688404B (en) | Rainfall landslide early warning method based on SVM-RF model | |
| CN116151483B (en) | Regional rainfall landslide probabilistic prediction method and prediction terminal | |
| CN112200356A (en) | Landslide prediction method, device, equipment and storage medium | |
| Biard et al. | Automated detection of weather fronts using a deep learning neural network | |
| CN112200355B (en) | Landslide prediction method, landslide prediction device, landslide prediction equipment and storage medium | |
| Meena et al. | Assessing the importance of conditioning factor selection in landslide susceptibility for the province of Belluno (region of Veneto, northeastern Italy) | |
| CN117556197B (en) | Typhoon vortex initialization method based on artificial intelligence | |
| KR20200052398A (en) | Method and apparatus for landslide susceptibility mapping using machine-learning architecture | |
| Ngo et al. | Landslide susceptibility mapping using single machine learning models: a case study from Pithoragarh District, India | |
| Bui et al. | Landslide susceptibility prediction mapping with advanced ensemble models: Son La province, Vietnam | |
| Akinci et al. | Evaluation of landslide susceptibility of Şavşat District of Artvin Province (Turkey) using machine learning techniques | |
| Viet Du et al. | Deep learning to assess the effects of land use/land cover and climate change on landslide susceptibility in the Tra Khuc river basin of Vietnam | |
| CN112200363A (en) | Landslide prediction method, device, equipment and storage medium | |
| CN117233869B (en) | Site short-term wind speed prediction method based on GRU-BiTCN | |
| Momeni et al. | Pattern‐based calibration of cellular automata by genetic algorithm and Shannon relative entropy | |
| CN112200357B (en) | Landslide prediction method, landslide prediction device, landslide prediction equipment and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |