CN112200355B - 一种滑坡的预测方法、装置、设备和存储介质 - Google Patents
一种滑坡的预测方法、装置、设备和存储介质 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
本发明实施例公开了一种滑坡的预测方法、装置、设备和存储介质,获取目标区域的第一栅格尺寸和由第一栅格尺寸重采样后的第二栅格尺寸,将第二栅格尺寸划分成多个栅格后并提取划分后的每个栅格的特征信息,可以同时兼顾栅格的平衡性、滑坡影响因子的丰富性以及计算量较小的优势,并便于在后续预测滑坡概率时,基于特征信息和滑坡影响因子,并以天为单位将滑坡影响因子输入至滑坡预测模型中,预测滑坡发生概率,将滑坡发生概率作为目标区域的滑坡预测结果,提高了滑坡预测的准确度。
Description
技术领域
本发明实施例涉及滑坡监测技术领域,尤其涉及一种滑坡的预测方法、装置、设备和存储介质。
背景技术
山体滑坡是最常见的灾难性自然灾害之一,其具有分布范围广、发生频次高、多发性、区域性和严重性等特点,因滑坡每年会造成大量的人员伤亡和重大的环境和基础设施损失。对滑坡的易发性进行评估具有重要意义。
现有的滑坡易发性预测根据其所依据的理论基础的差异,可以分为确定性方法和非确定性方法。确定性方法主要是基于专家经验和知识的定向分析和基于滑坡过程或者物理学模型的进行分析的方法,预测准确度较差。近年来随着计算机技术和3S技术的高速发展,非确定性方法得到了广泛的应用,主要包括模糊逻辑法、层次分析法、决策树等。但是上述方法所处理的滑坡因子的时间精度较差,尤其是降雨量因子,通常仅考虑一年的降雨量,因此预测精度也不太理想,同时,由于滑坡的影响因子众多,仅采用一种非确定方法进行预测的预测效率较差,并且,不同尺寸的栅格下的滑坡数据的某些特征差异较大,直接利用由一个栅格尺寸重采样到另一个栅格尺寸的滑坡数据进行滑坡预测的准确度较差。
发明内容
本发明提供了一种滑坡的预测方法、装置、设备和存储介质,实现了对重采样后的栅格的滑坡数据进行特征再提取,根据提取的特征信息和每天的各个栅格的滑坡数据进行滑坡预测,提高了预测精度,同时基于滑坡预测模型进行滑坡预测,提高了预测效率。
第一方面,本发明实施例提供了一种滑坡的预测方法,该滑坡的预测方法包括:
获取目标区域的第一栅格尺寸,并确定由所述第一栅格尺寸重采样后的第二栅格尺寸,其中,所述第二栅格尺寸大于所述第一栅格尺寸;
根据所述第一栅格尺寸将所述第二栅格尺寸对应的初始栅格划分成多个栅格,并基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的特征信息,其中,所述滑坡影响因子为所述第二栅格尺寸对应的栅格下每个栅格的设定时间段内的每天的滑坡数据;
基于滑坡预测模型,根据所述特征信息和所述滑坡影响因子确定各栅格滑坡的发生概率,将所述发生概率作为所述目标区域的滑坡预测结果。
第二方面,本发明实施例还提供了一种滑坡的预测装置,该滑坡的预测装置包括:
栅格尺寸确定模块,用于获取目标区域的第一栅格尺寸,并确定由所述第一栅格尺寸重采样后的第二栅格尺寸,其中,所述第二栅格尺寸大于所述第一栅格尺寸;
特征信息提取模块,用于根据所述第一栅格尺寸将所述第二栅格尺寸对应的初始栅格划分成多个栅格,并基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的特征信息,其中,所述滑坡影响因子为所述第二栅格尺寸对应的栅格下每个栅格的设定时间段内的每天的滑坡数据;
滑坡预测模块,用于基于滑坡预测模型,根据所述特征信息和所述滑坡影响因子确定各栅格滑坡的发生概率,将所述发生概率作为所述目标区域的滑坡预测结果。
第三方面,本发明实施例还提供了一种滑坡的预测设备,该设备包括:
一个或多个处理器;
存储器,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现本发明任意实施例提供的滑坡的预测方法。
第四方面,本发明实施例还提供了一种包含计算机可执行指令的存储介质,所述计算机可执行指令在由计算机处理器执行时用于执行本发明任意实施例提供的滑坡的预测方法。
本发明实施例的技术方案,获取目标区域的第一栅格尺寸和由第一栅格尺寸重采样后的第二栅格尺寸,将第二栅格尺寸划分成多个栅格后并提取划分后的每个栅格的特征信息,可以同时兼顾栅格的平衡性、滑坡影响因子的丰富性以及计算量较小的优势,并便于在后续预测滑坡概率时,基于特征信息和滑坡影响因子,并以天为单位将滑坡影响因子输入至滑坡预测模型中,预测滑坡发生概率,将滑坡发生概率作为目标区域的滑坡预测结果,提高了滑坡预测的准确度和滑坡预测效率。
附图说明
图1是本发明实施例一中的一种滑坡的预测方法的流程图;
图2是本发明实施例二中的一种滑坡的预测方法的流程图;
图3是本发明实施例三中的一种滑坡的预测方法的流程图;
图4是本发明实施例四中的一种滑坡的预测装置的结构示意图;
图5是本发明实施例五中的一种滑坡的预测设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
图1是本发明实施例一提供的一种滑坡的预测方法的流程图,本实施例可适用于对滑坡易发性进行评估的情况,该方法可以由滑坡的预测装置来执行,如图1所示,该方法包括如下步骤:
S110、获取目标区域的第一栅格尺寸,并确定由第一栅格尺寸重采样后的第二栅格尺寸。
其中,目标区域通常是发生过滑坡的区域,也可以是指定的任意一个区域。栅格数据就是将空间分割成有规律的网格,每一个网格成为一个栅格或单元,并在各个单元上赋予相应的属性值来表示实体的一种数据形式。所述第二栅格尺寸大于所述第一栅格尺寸。
可选地,所述第一栅格尺寸的确定方法,包括:确定所述目标区域的每个栅格的初始尺寸和初始栅格数,以及确定滑坡栅格数和非滑坡栅格数,其中,所述滑坡栅格数是发生滑坡的栅格数量,所述非滑坡栅格数是未发生滑坡的栅格数量;如果所述滑坡栅格数与非滑坡栅格数的比值处于设定平衡区间内,将所述初始尺寸作为所述第一栅格尺寸;如果所述滑坡栅格数与所述非滑坡栅格数的比值未处于所述设定平衡区间内,调节所述目标区域的每个栅格的尺寸,重新确滑坡栅格数、非滑坡栅格数和当前栅格数;如果重新确定的滑坡栅格数与非滑坡栅格数的比值处于所述设定平衡区间内,且所述当前栅格数与所述初始栅格数的比值小于设定阈值,将调节后的栅格的尺寸作为所述第一栅格尺寸。
其中,初始尺寸可以是任意尺寸,例如1km*1km。所述初始栅格数是初始尺寸对应的栅格数量。可以理解的是,将目标区域进行栅格化时,目标区域的滑坡落入栅格中,根据滑坡位置和滑坡面积确定发生滑坡的栅格和未发生滑坡的栅格,并确定滑坡栅格数和非滑坡栅格数。进一步地,计算滑坡栅格数和非滑坡栅格数的比值,如果滑坡栅格数和非滑坡栅格数的比值处于设定平衡区间内,则说明初始尺寸对应的栅格内的滑坡栅格数和非滑坡栅格数平衡,可以将初始尺寸作为第一栅格尺寸;如果滑坡栅格数和非滑坡栅格数的比值未处于设定平衡区间内,调节每个栅格的尺寸,重新确定滑坡栅格数和非滑坡栅格数以及当前栅格数,如果重新确定的滑坡栅格数与非滑坡栅格数的比值处于所述设定平衡区间内,则说明栅格尺寸变化前后,滑坡栅格的数量以及栅格数量处于平衡状态,如果当前栅格数与初始栅格数的比值小于设定阈值,说明栅格尺寸变化前后,目标区域的总的滑坡栅格数处于平衡状态,可以将调节后的每个栅格的尺寸作为第一栅格尺寸。
可选地,本实施例还可以通过如下步骤确定第一栅格尺寸:获取所述目标区域的各雨量站的降雨数据、所述目标区域的每个栅格的初始尺寸和初始栅格数,以及确定滑坡栅格数和非滑坡栅格数,其中,所述滑坡栅格数是发生滑坡的栅格数量,所述非滑坡栅格数是未发生滑坡的栅格数量;根据特定数量的雨量站的降雨数据计算任一当前栅格对应的雨量站的降雨预测值,基于所述降雨预测值与当前栅格对应的雨量站的实测雨量值计算降雨误差;如果所述降雨误差小于设定误差阈值,且所述滑坡栅格数与所述非滑坡栅格数的比值处于设定平衡区间内,将所述初始尺寸作为所述第一栅格尺寸;如果所述降雨误差大于设定误差阈值,调节所述目标区域的每个栅格的尺寸,根据调节后的所述任一栅格对应的雨量站的降雨预测值与实测雨量值重新计算降雨误差;如果重新计算的降雨误差小于所述设定误差阈值,且重新确定的滑坡栅格数与重新确定的非滑坡栅格数的比值处于设定平衡区间内,且当前栅格数与所述初始栅格数的比值小于设定阈值,将调节后的栅格的尺寸作为所述第一栅格尺寸。
具体地,可以在目标区域设置多个监测点,每个监测点设置一个雨量站,以实时获取各个监测点的降雨数据。具体地,将目标区域栅格化,设置的多个雨量站分布在不同的栅格中,预测某个监测点的降雨量时,可以获取多个其他监测点的降雨量,基于协同克里金插值法并以高程信息作为辅助信息确定该监测点的降雨量,将所述降雨量作为该监测点的预测降雨量,将预测降雨量与通过该监测点的雨量站获取的实测降雨量求降雨误差,如果降雨误差小于设定误差阈值,则说明初始尺寸下的栅格中该雨量站的预测降雨量较准确;同时,计算滑坡栅格数与非滑坡栅格数的比值,判断该比值是否处于设定平衡区间内,如果是,则说明初始尺寸对应的滑坡栅格数和非滑坡栅格数平衡,将初始尺寸作为第一栅格尺寸;如果降雨误差不小于设定误差阈值,则说明初始尺寸下的栅格中该雨量站的预测降雨量不准确,需要调节每个栅格的尺寸,根据栅格尺寸调节后的栅格中该雨量站的降雨预测值与实测雨量值重新计算降雨误差,直至重新计算的降雨误差小于设定误差阈值,如果重新确定的滑坡栅格数与重新确定的非滑坡栅格数的比值也处于设定平衡区间内,且当前栅格数与初始栅格数的比值小于设定阈值,则说明调节后的栅格的雨量站的实测降雨量较准确,且调节后的栅格尺寸对应的滑坡栅格数与非滑坡栅格数平衡,并且调节后的总的栅格数平衡,则将调节后的栅格的尺寸作为第一栅格尺寸。
本实施例中的第二栅格尺寸也可以采用上述方式确定,或者将初始化栅格尺寸作为第二栅格尺寸,第二栅格尺寸的具体确定方式不再详细赘述。
S120、根据第一栅格尺寸将第二栅格尺寸对应的初始栅格划分成多个栅格,并基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的特征信息。
其中,所述滑坡影响因子为所述第二栅格尺寸对应的栅格下的每个栅格的设定时间段内的每天的滑坡数据。滑坡影响因子主要包括滑坡的各个滑坡影响因子,如高程因子和环境因子。设定时间段可以是一天、三天、一周、一个月或者其他时间段。所述滑坡影响因子包括动态因子和静态因子,其中,所述动态因子包括降雨量和土壤湿度中的至少一项,所述静态因子包括高程、坡度、坡向、平面曲率、剖面曲率、地形湿度指数、水流强度指数、沉积运输指数、地形粗糙度指数、与断层距离、与河流距离、与道路距离、岩性、土地利用、植被覆盖率中的至少一项。所述特征信息包括:所述滑坡影响因子的加权值、最大值、最小值、差分的均值和标准差特征中的至少一种。
可选地,特征信息的提取方法为:获取当前栅格的八邻域栅格内的滑坡影响因子的特征值;根据所述八邻域栅格内的特征值,确定所述当前栅格的特征信息。示例性地,如果当前栅格的特征信息为坡度,可以获取当前栅格的八邻域栅格内的坡长和坡高,基于当前栅格的坡高以及八邻域栅格的坡长,并采用勾股定理,计算当前栅格与任一八邻域栅格构成的锐角角度,基于当前栅格与各八邻域栅格构成的锐角角度,确定当前栅格的坡度。可选地,其他八邻域栅格也可以采用上述同样方式计算坡度,直至所有栅格的坡度计算完成,得到第一栅格尺寸的栅格下的滑坡影响因子的特征信息。
可以理解的是,在不同栅格尺寸下,滑坡影响因子的表征程度不同,尺寸较小的栅格能够保留特征信息的丰富性,尺寸较大的栅格能够精简数据的冗余度。将目标区域的滑坡影响因子由第一栅格尺寸重采样到第二栅格尺寸后,第二栅格尺寸下的滑坡影响因子与第一栅格尺寸下的滑坡影响因子的某些特征差异较大(例如静态因子中的坡度和坡向)。为了同时兼顾栅格的平衡性、滑坡影响因子的丰富性以及计算量较小的优势,本实施例将滑坡影响因子重采样后,根据重采样前的第一栅格尺寸将重采样后的第二栅格尺寸对应的栅格划分成多个栅格,并基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的滑坡影响因子的特征信息,实现对第二栅格尺寸对应的栅格进行特征再提取。需要说明的是,在基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的滑坡影响因子的特征信息之前,需要判断重采样前后的滑坡影响因子的差异,如果差异较大,采用S120的方式对第二栅格尺寸对应的栅格进行特征再提取,如果差异较小,直接获取第一栅格尺寸或第二栅格尺寸下的滑坡影响因子。
示例性地,第一栅格尺寸是30m*30m,第二栅格尺寸是1km*1km,将第一栅格尺寸对应的栅格重采样至第二栅格尺寸对应的栅格,坡度与坡向在两种栅格尺寸下的差异性较大,根据第一栅格尺寸将第二栅格尺寸对应的初始栅格划分成1111个栅格,并在1111个栅格下提取坡度和坡向的加权值、最小值、差分均值以及标准差中的至少一种,完成对第二尺寸对应的栅格的特征再提取。
S130、基于滑坡预测模型,根据特征信息和滑坡影响因子确定各栅格滑坡的发生概率。
其中,所述滑坡预测模型可以是单一神经网络模型,也可以是输入端和输出端依次连接形成的复杂神经网络模型。具体地,滑坡预测模型的训练过程为:从滑坡野外勘察报告、典型滑坡监测报告等文件提取目标区域或全部区域的滑坡的基本信息包括滑坡发生的时间、经纬度、受灾规模等滑坡信息以及滑坡区域的实测降雨量,并根据该滑坡信息确定各个滑坡发生地的历史时间段的滑坡影响因子,历史时间段的滑坡影响因子可以包括不同栅格大小下的滑坡影响因子以及滑坡特征信息,将历史时间段的滑坡影响因子以设定比例组成训练集和验证集,如8:2或7:3;将各个滑坡影响因子进行数据校正与配准,并进行坐标统一处理和栅格统一处理,从而得到栅格大小一致的滑坡影响因子,并将通过上述步骤处理后的滑坡影响因子训练滑坡预测模型。其中,该滑坡影响因子通过特征工程提取各个滑坡影响因子的特征,组成滑坡预测模型的输入特征矩阵;初始化滑坡预测模型的参数,并将该输入特征矩阵输入滑坡预测模型中,进行模型训练,得到历史时间段的历史滑坡信息,并基于F1-值(F1-Score)与ROC(Receiver Operating Characteristic,接收者操作特征)作为评价指标,根据评价结果进行滑坡预测模型参数的调整,当满足时,训练初步完成,通过验证集进行模型验证,验证通过则得到训练好的滑坡预测模型。其中,所述历史滑坡信息包括各栅格发生滑坡的概率。
可选地,在根据特征信息和滑坡影响因子输入至滑坡预测模型之前,还可以对滑坡影响因子和特征信息进行预处理,其中,所述预处理包括坐标统一处理、校正处理、数据离散化处理和栅格统一处理中的至少一项。本实施例通过将特征信息和滑坡影响因子结合确定各栅格滑坡的发生概率,可以增加滑坡影响因子的数据量,根据特征信息可以更准确的确定各栅格滑坡的发生概率,提高了滑坡预测的准确性。
本实施例提供的技术方案,获取目标区域的第一栅格尺寸和由第一栅格尺寸重采样后的第二栅格尺寸,将第二栅格尺寸划分成多个栅格后并提取划分后的每个栅格的特征信息,可以同时兼顾栅格的平衡性、滑坡影响因子的丰富性以及计算量较小的优势,并便于在后续预测滑坡概率时,基于特征信息和滑坡影响因子,并以天为单位将滑坡影响因子输入至滑坡预测模型中,预测滑坡发生概率,将滑坡发生概率作为目标区域的滑坡预测结果,提高了滑坡预测的准确度和滑坡预测效率。
实施例二
图2为本发明实施例二提供的一种滑坡的预测方法的流程图,本实施例是对上一实施例的进一步细化。可选地,所述基于滑坡预测模型,根据所述特征信息和所述滑坡影响因子确定各栅格滑坡的发生概率,包括:基于第一预设模型,根据所述特征信息和所述滑坡影响因子确定每天滑坡的发生概率;如果当前天滑坡的发生概率大于预设概率阈值,则基于第二预设模型,根据在所述特征信息和所述滑坡影响因子,确定各个栅格滑坡的发生概率。在该方法实施例中未详尽描述的部分请参考上述实施例。具体参见图2所示,该方法可以包括如下步骤:
S210、获取目标区域的第一栅格尺寸,并确定由第一栅格尺寸重采样后的第二栅格尺寸。
S220、根据第一栅格尺寸将第二栅格尺寸对应的初始栅格划分成多个栅格,并基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的特征信息。
其中,所述滑坡影响因子为所述第二栅格尺寸对应的栅格下每个栅格的设定时间段内的每天的滑坡数据。
S230、基于第一预设模型,根据特征信息和滑坡影响因子确定每天滑坡的发生概率。
其中,第一预设模型可以是神经网络模型,或者其他学习算法。示例性的,第一预设模型可以是支持向量机算法(Support Vector Machine,SVM)、长短期记忆网络(LongShort-Term Memory,LSTM)、逻辑回归模型(Logistics Regression,LR)、XGBoost(ExtremeGradient Boosting,极端梯度提升决策树)算法、GBDT(Gradient Boosting DecisionTree,梯度提升决策树)算法、全卷积网络(Fully Convolutional Networks,FCN)、循环卷积网络(Recurrent Neural Network,RNN)、残差网络(Residual Network,ResNet)、门控循环单元(Gate Recurrent Unit,GRU)等。所述第一预设模型的训练方法与前述实施例中的滑坡预测模型的训练方法一样,此处不再详细描述。
可选地,每天滑坡的发生概率的确定方法为:计算每天的各个栅格的所述动态因子的构造特征,其中,所述构造特征包括各个栅格的每个动态因子的总和、平均值、最大值、最小值、极差、四分位数和雨量时间,并根据每天的所述构造特征和静态因子得到所述第一预设模型的第一特征矩阵;将所述第一特征矩阵输入所述第一预设模型中,以根据所述第一预设模型确定每天滑坡的发生概率。
如前述实施例所述,所述滑坡影响因子包括动态因子和静态因子,动态因子包括降雨量(降水量)和土壤湿度,静态因子包括高程、坡度、坡向、平面曲率、剖面曲率、地形湿度指数、水流强度指数、沉积运输指数、地形粗糙度指数、与断层距离、与河流距离、与道路距离、岩性、土地利用、植被覆盖率等。
示例性的,所获取的每天的各个栅格的滑坡影响因子X为:X={X(1),X(2),…,X(d)},其中,上角标表示天数,X(c),c=1,2,…,d,表示第c天目标区域内各个栅格的滑坡影响因子,其中矩阵X(c)的行表示栅格,列表示滑坡因子,矩阵X(c)的大小为m×n,即包括m个栅格,n个滑坡因子,其中,i=1,2,…n1,表示动态因子(降雨量和土壤湿度等),i=n1+1,n1+2,…n,表示静态因子。以天为单位构造特征集X1(第一特征矩阵),构造方式如下:对于每一个X(c),按栅格计算n1个动态因子的总和、平均值、最大值、最小值、极差、上四分位数和下四位数等构造特征,由各个构造特征和静态因子,得到特征集X1。对于训练数据,还可以构造滑坡标签向量Y1,以表示每天是否发生滑坡,其中,1表示发生滑坡,0表示未发生滑坡。
S240、如果当前天滑坡的发生概率大于预设概率阈值,则基于第二预设模型,根据在特征信息和所述滑坡影响因子,确定各个栅格滑坡的发生概率。
可选地,每个栅格滑坡的发生概率的确定方法为:对所述当前天对应的各个栅格的滑坡影响因子进行特征提取,以生成所述第二预设模型的第二特征矩阵;对所述第二特征矩阵进行归一化处理;将归一化处理后的所述第二特征矩阵输入所述第二预设模型中,以根据所述第二预设模型确定当前天的各个栅格滑坡的发生概率。
具体地,假设第c天的滑坡发生概率大于预测概率阈值,则将第c天的滑坡影响因子X(c)发送至第二预设模型,构造第二预设模型的特征集,具体构造方式为:对该天的各个栅格的动态因子计算其3天、7天,15天和30天分别对应的总和、均值、方差、中位数以及差分的均值和方差、偏度、峰度等构造特征,所述峰度为基于降雨衰减指数计算的前期降雨指数。同时,考虑以当前栅格为中心的设定范围内,如3*3区域范围内,各个栅格的情况,包括当前栅格是否为设定范围内的最大值或最小值,是否超过设定范围对应的均值,以及设定范围的栅格是否发生过滑坡,当前栅格发生滑坡的次数等特征,最终得到第二预设模型的第二特征矩阵X2 fea。
具体地,可以基于max-min归一化算法对第二特征矩阵X2 fea的各个特征进行归一化处理。当然也可以选取其他归一化算法进行归一化处理。本发明实施例对第一特征矩阵和第二特征矩阵归一化算法不进行限定。进一步地,在归一化处理之后,还包括:基于多重共线性方法对第二特征矩阵进行特征选择,以根据共线性程度进行特征筛选,得到筛选的后的第二特征矩阵
具体地,第二预设模型可以是一维卷积神经网络模型,包括卷积层、批标准化层、激活函数和优化层,其中,卷积层的卷积核个数的取值范围为32~512,步长为16;激活函数可以包括ReLU函数(Rectified Linear Unit,线性整流函数)、Linear函数(线性函数)、Sigmoid函数和Tanh函数(双曲函数)中的任意一种;优化方法包括SGD(StochasticGradient Descent,梯度更新规则)、Adam(Adaptive Moment Estimation,自适应矩估计)、Nadam(Nesterov Adaptive Moment Estimation,Nesterov加速自适应矩估计)、Adagrad(Adaptive Gradient Algorithm,自适应梯度算法)和RMSprop(Root Mean Square Prop,前向均方根梯度下降演算法)等优化算法中任意一种;初始学习率的取值可以是0.0001、0.001、0.01或0.1;隐藏层神经元个数的取值范围为4~256,步长为4,隐藏层层数的取值范围为3~8,步长为1;随机丢弃率为0~0.8,步长为0.05。所述第一预设模型的训练方法与前述实施例中的滑坡预测模型的训练方法一样,此处不再详细描述。
本实施例提供的技术方案,获取目标区域的第一栅格尺寸和由第一栅格尺寸重采样后的第二栅格尺寸,将第二栅格尺寸划分成多个栅格后并提取划分后的每个栅格的特征信息,基于特征信息和滑坡影响因子进行滑坡预测,提高了数据的精度,尤其是对于滑坡中的动态因子,以天为单位进行预测,充分考虑了动态因子的时序信息,提高了预测精度;通过设置两个预设模型,分两阶段进行滑坡预测,第一阶段预测当前天发生滑坡的概率是否大于设定值,若是,则进行第二阶段滑坡预测,以当前天的数据为基准进行具体栅格的滑坡发生概率,将滑坡发生概率作为目标区域的滑坡预测结果,大幅度降低了计算资源的浪费,提高了滑坡预测的效率。
实施例三
图3为本发明实施例三提供的一种滑坡的预测方法的流程图。本实施例用于对滑坡预测的整个流程进行解析。在该方法实施例中未详尽描述的部分请参考上述实施例。具体参见图3所示,该方法可以包括如下步骤:
S310、获取目标区域的第一栅格尺寸,并确定由第一栅格尺寸重采样后的第二栅格尺寸。
可选地,在获取目标区域的第一栅格尺寸之前,还可以获取目标区域的遥感影像信息和地形信息,基于功能区分类模型,根据所述遥感影像信息和所述地形信息,确定目标区域的各功能区,并根据各功能区对应的标签筛选目标功能区。
其中,所述遥感影像信息指的是卫星影像,可以通过Landsat(陆地卫星)系列的地表反射率产品获取,所述遥感影像信息可以反映目标区域的各栅格的地物类型,例如,草地类型、林地类型、水稻类型以及建筑物类型;所述地形信息可以从目标区域的数字高程模型(Digital Elevation Model,简称DEM)获取,所述地形信息反映了目标区域的某个栅格和周围栅格的栅格信息的最大变化程度。所述地形信息可以包括高程信息、坡度信息以及曲率信息等,所述高程信息可以反映目标区域的某个栅格的地形特征高程点、重要地理目标附近的高程点或者分布范围上的关键高程点,所述坡度信息可以反映目标区域的某个栅格的陡缓的程度,所述曲率信息可以反映目标区域的某个栅格的凹凸情况。
可选地,确定所述目标区的各功能区的方法包括:对所述遥感影像信息进行预处理,其中,所述预处理包括大气校正、辐射校正以及边缘检测处理中的至少一种;分别对预处理后的遥感影像信息和地形信息进行特征提取,并将遥感影像特征和地形特征输入至所述功能区分类模型,确定所述目标区域的各功能区,其中,所述遥感影像特征包括所述遥感影像信息的均值、极大值以及极小值,所述地形信息包括高程信息、坡度信息以及曲率信息,所述地形特征分别包括所述高程信息、坡度信息以及曲率信息的均值、极大值以及极小值。
所述功能区分类模型为随机森林模型,所述功能区分类模型的训练方法,包括:获取初始森林模型;将历史时间段内的各功能区的样本遥感影像信息、样本地形信息输入至所述初始森林模型,确定各功能区的预测标签和预测标签的概率;基于所述预测标签、所述预测标签的概率和各功能区的样本标签,对所述初始森林模型的参数进行调整,直至所述预测标签与所述样本标签一致,且所述预测标签概率达到设定阈值,得到所述功能区分类模型。其中,所述概率指的是初始森林模型的输出的预测标签的可能性,通过迭代调整初始森林模型的参数,如果预测标签与所述样本标签一致,且预测标签概率达到设定阈值,则预测标签是样本标签的可能性较大,将该迭代次数下的初始森林模型作为功能区分类模型。所述设定阈值可以是一个较大值,例如设定阈值取值为0.9。可选地,在将所述样本遥感影像信息输入至初始森林模型之前,还可以对样本遥感影像信息进行预处理,例如,对样本遥感影像信息进行大气校正、辐射校正、边缘检测、假彩色合成处理等,以提高样本遥感信息的功能区分类模型的训练精度。
可以理解的是,功能区分类模型输出的各功能区包括目标区域的所有功能区,包括城镇建筑、城市绿地、水体、农田、裸土、山林等功能区,各功能区的滑坡发生概率不同。例如,城镇建筑和城镇绿地发生滑坡的概率极小,水体和山林发生滑坡的概率较大。需要说明的是,通过将各功能区对应的历史预测概率与历史概率阈值比较,筛选出目标功能区,对目标功能区进行针对性的滑坡预测,可以减小计算量,有利于提高目标区域的滑坡预测效率。
S320、根据第一栅格尺寸将第二栅格尺寸对应的初始栅格划分成多个栅格,并基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的特征信息。
若通过S310确定目标区域的各功能区,并根据各功能区对应的标签筛选目标功能区,则S320可以替换为:根据第一栅格尺寸将第二栅格尺寸对应的目标功能区的初始栅格划分成多个栅格,并基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的特征信息。
所述滑坡影响因子包括动态因子和静态因子,其中,所述动态因子包括降雨量和土壤湿度中的至少一项,所述降雨量可以通过将获取的目标区域的地理环境数据输入至降雨插值模型确定。所述地理环境数据包括地理位置数据、大气数据、地形数据和下垫面数据。所述降雨插值模型可以是一种反向传播模型(BP,back propagation neuralnetwork),基于误差反向传播算法的由非线性变换单元组成的多层前馈型网络,BP一般由输入层、隐含层和输出层组成,每层中又包括N个神经元,同层神经元之间互相独立,每层之间神经元的输出经过特定的激励函数后,只会影响下层神经元的输入。
具体地,降雨插值模型的训练方法为:获取初始模型,确定所述初始模型的初始权值矩阵和初始阈值,并计算所述初始模型的适应度;将历史时间段内的样本地理环境信息输入至所述初始模型,确定预测降雨量,根据所述预测降雨量和所述历史时间段的实测降雨量,并基于遗传算法迭代调整所述初始模型的初始权值矩阵和初始阈值;基于迭代调整后的权值矩阵阈值,调整所述初始模型并计算调整后的模型的适应度,直至重新计算的适应度达到期望值,将达到所述期望值的适应度对应的模型作为所述降雨插值模型。
具体地,所述初始权值矩阵可以包括输入层与隐含层之间的连接权值和隐含层与输出层之间的连接权值;初始阈值可以包括隐含层的阈值和输出层的阈值。初始模型的适应度的计算公式为:其中,n是样本个数,y(i)是初始模型的第i个样本的实际输出,t(i)是初始模型的第i个样本的期望输出。具体地,将历史时间段内的样本地理环境信息输入至初始模型后,得到预测降雨量,计算预测降雨量和实测降雨量之间的均方根误差,如果均方根误差大于预设误差阈值,迭代调整初始模型的权值编码和阈值编码,并基于调整后的权值编码和阈值编码确定调整后的权值矩阵和阈值,直至均方根误差不小于预设误差阈值,则说明初始模型调整至稳定状态,并采用上述适应度计算公式计算调整后的模型的适应度,如果重新计算的适应度达到期望值,将达到期望值的适应度对应的模型作为所述降雨插值模型。
可选地,采用二进制编码、实数编码或者格雷码编码方式对权值和阈值进行编码,并根据权值编码和阈值编码确定调整后的权值矩阵和阈值。可选地,还可以通过计算预测降雨量和实测降雨量之间的平均误差、平均绝对误差或线性相关系数,并根据所述平均误差、所述平均绝对误差或所述线性相关系数迭代调整初始模型的权值编码和阈值编码,直至平均误差、平均绝对误差或线性相关系数不小于预设误差阈值,说明初始模型调整至稳定状态,并进一步计算调整后模型的适应度以及确定降雨插值模型。
需要说明的是,历史时间段内的样本地理环境数据是多维度的,在训练降雨插值模型时,多维度的样本地理环境数据可以提高降雨插值模型的精度和鲁棒性,因而,在获取每天的地理环境数据后,基于降雨插值模型预测可以提高降雨概率的预测精度,有利于提高滑坡预测的预测精度。
S330、基于第一预设模型,根据特征信息和滑坡影响因子确定每天滑坡的发生概率。
与前述实施例不同的,第一预设模型可以是自步分类学习模型,所述第一预设模型根据历史滑坡数据、历史非滑坡数据迭代调整自步因子和下采样比例确定。所述自步因子根据所述历史非滑坡数据的分箱个数和迭代次数确定,所述下采样比例根据每个分箱的自步因子和所有分箱的自步因子确定。可选地,所述历史滑坡数据和历史非滑坡数据的确定方法,包括:确定历史区域的各栅格对应的标签,并确定各所述历史时间点下的滑坡栅格和非滑坡栅格分别对应的历史降雨量和特定相关因子;基于所述标签,生成包括滑坡栅格和非滑坡栅格分别对应的历史降雨量和特定相关因子的滑坡数据集,根据所述滑坡数据集中各个栅格对应的标签确定所述历史滑坡数据和历史非滑坡数据。
可选地,所述历史滑坡数据和历史非滑坡数据的确定方法,包括:确定历史区域的各栅格对应的标签,并确定各所述历史时间点下的滑坡栅格和非滑坡栅格分别对应的历史降雨量和特定相关因子;基于所述标签,生成包括滑坡栅格和非滑坡栅格分别对应的历史降雨量和特定相关因子的滑坡数据集,根据所述滑坡数据集中各个栅格对应的标签确定所述历史滑坡数据和历史非滑坡数据。
可选地,所述第一预设模型的训练方法,包括:获取初始预测模型,并获取各历史滑坡时间点下的滑坡栅格的降雨量和非滑坡栅格的降雨量,以及各所述标签下的特定相关因子;按天将所述历史滑坡时间点下的滑坡栅格的降雨量、所述非滑坡栅格的降雨量以及各所述标签下的特定相关因子输入至所述初始预测模型,基于初始预测模型输出的滑坡概率和所述历史滑坡时间点对应的滑坡概率调整所述初始预测模型的损失函数,直至所述损失函数达到设定阈值,得到所述第一预设模型。
其中,所述基于初始预测模型输出的滑坡概率和所述历史滑坡时间点对应的滑坡概率调整所述初始预测模型的损失函数,直至所述损失函数达到设定阈值,得到所述第一预设模型,包括:基于所述历史滑坡时间点下的滑坡栅格的降雨量、所述非滑坡栅格的降雨量、各所述标签以及各所述标签下的特定相关因子,确定所述初始预测模型的初始分类硬度;根据所述初始分类硬度确定所述非滑坡栅格的降雨量和非滑坡栅格的特定相关因子的分箱个数,基于所述分箱个数确定所述初始预测模型的各分箱的自步因子;基于所述自步因子确定每个分箱的下采样比例,并基于所述下采样比例确定各分箱内的下采样的非滑坡栅格的降雨量,以及确定下采样的非滑坡栅格的特定相关因子;将所述下采样的非滑坡栅格的降雨量、下采样的非滑坡栅格的特定相关因子以及滑坡栅格的降雨量输入至所述初始预测模型,基于所述初始预测模型输出的滑坡概率和所述历史滑坡时间点对应的滑坡概率确定所述损失函数,基于所述损失函数在单个样本上的训练分类硬度,迭代调整自步因子和每个分箱的下采样比例;根据迭代调整的所述自步因子和所述每个分箱的下采样比例调整所述初始预测模型,直至损失函数达到设定阈值,得到所述第一预设模型。
其中,每个分箱的下采样比例为每个分箱的自步因子与所有分箱的自步因子的比值,所述自步因子通过将分箱个数、分箱标签数的相反数以及迭代次数之和加1得到。
选择自步分类学习模型进行滑坡预测的好处是:根据历史非滑坡数据的分箱个数和迭代次数确定自步分类学习模型的自步因子,根据每个分箱的自步因子和所有分箱的自步因子确定下采样比例,因此,在基于历史滑坡数据、历史非滑坡数据迭代调整自步因子和下采样比例时,随着迭代次数的增加,每个分箱的采样比例由硬度低的下采样数量多到最后每个分箱内的下采样数量逐渐变化均匀,并且,在每个箱子的下采样比例均匀变化的同时,还可以保证硬度小的分箱的下采样数量始终高于硬度大的分箱的下采样数量,得到多样化、鲁棒性高且包容性强的第一预设模型。
S340、如果当前天滑坡的发生概率大于预设概率阈值,则基于第二预设模型,根据在特征信息和所述滑坡影响因子,确定各个栅格滑坡的发生概率。
可选地,所述第二预设模型也可以是自步集成学习模型。训练所述第二预设模型的历史滑坡数据的确定方式以及第二预设模型的训练方法可参考S330,该步骤不再详细阐述。
可选地,确定当前天的各个栅格滑坡的发生概率之后,还可以将所述当前天的各个栅格的发生概率与所述目标区域的预警阈值进行比较,基于得到的各个栅格的发生概率对应的预警等级确定所述目标区域的滑坡的发生等级,其中,所述预警阈值根据各区域栅格滑坡的第一系数和第二系数确定。
其中,所述第一系数和第二系数的确定方法包括:获取各区域栅格滑坡的设定时间段内的已发生滑坡的历史概率;基于所述历史概率的分布特征确定所述第一已确定概率区间和所述第二已确定概率区间;分别计算所述第一已确定概率区间和所述第二已确定概率区间内的同一风险等级的概率阈值与对应系数的乘积之和,得到中间已确定概率;如果所述中间已确定概率未达到设定评价指标,迭代调整各风险等级的概率阈值的对应系数,直至所述中间已确定概率达到所述设定评价指标,根据达到所述设定评价指标的对应系数确定第一系数和第二系数。
其中,确定所述预警阈值的方法,包括:基于所述各个栅格滑坡的发生概率的分布特征确定第一预测概率区间和第二预测概率区间;计算所述第一预测概率区间内各风险等级与所述第一系数的第一乘积,以及计算所述第二预测概率区间内各风险等级与所述第二系数的第二乘积,将所述第一乘积与所述第二乘积的和作为各风险等级的预警阈值。
其中,所述基于所述历史概率的分布特征确定所述第一已确定概率区间和所述第二已确定概率区间,包括:确定所述历史概率和所述历史概率对应的次数;根据所述历史概率和各所述历史概率对应的滑坡发生次数,确定所述历史概率的密集度分布特征和断点分布特征;基于所述密集度分布特征确定所述第一已确定概率区间,并基于所述断点分布特征确定所述第二已确定概率区间。
可以理解的是,上述方式基于各区域栅格滑坡的设定时间段内的已发生滑坡的历史概率,根据历史概率的分布特征确定第一已确定概率区间和第二已确定概率区间,根据第一已确定概率区间、第二已确定概率区间内的各风险等级的概率阈值以及设定评价指标确定各区域的第一系数和第二系数,当获取不同区域的发生概率后,可以根据不同区域对应的第一系数、第二系数和发生概率灵活确定不同区域的预警阈值,有利于后续根据各区域对应的预警阈值准确确定该区域的风险等级。
其中,所述基于得到的各个栅格的发生概率对应的预警等级确定所述目标区域的滑坡的发生等级,包括:筛选大于所述预警阈值中的第一等级的目标等级;计算所述目标等级对应的栅格的平均等级;基于所述平均等级、设定系数、目标等级对应的栅格数与目标区域内所有栅格数的比值,确定所述目标区域的滑坡的发生等级。
其中,目标区域的滑坡的发生等级的计算公式为:
其中,α为设定系数,通过对目标区域的滑坡的历史发生等级进行贝叶斯计算得到,为目标区域的栅格的平均等级,p为目标等级对应的栅格数与目标区域内所有栅格数的比值。
实施例四
图4为本发明实施例四提供的一种滑坡的预测装置的结果示意图,如图4所述,该滑坡的预测装置包括:栅格尺寸确定模块410、特征信息提取模块420和滑坡预测模块430。
其中,栅格尺寸确定模块410,用于获取目标区域的第一栅格尺寸,并确定由所述第一栅格尺寸重采样后的第二栅格尺寸,其中,所述第二栅格尺寸大于所述第一栅格尺寸;
特征信息提取模块420,用于根据所述第一栅格尺寸将所述第二栅格尺寸对应的初始栅格划分成多个栅格,并基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的特征信息,其中,所述滑坡影响因子为所述第二栅格尺寸对应的栅格下每个栅格的设定时间段内的每天的滑坡数据;
滑坡预测模块430,用于基于滑坡预测模型,根据所述特征信息和所述滑坡影响因子确定各栅格滑坡的发生概率,将所述发生概率作为所述目标区域的滑坡预测结果。
本发明实施例的技术方案,获取目标区域的第一栅格尺寸和由第一栅格尺寸重采样后的第二栅格尺寸,将第二栅格尺寸划分成多个栅格后并提取划分后的每个栅格的特征信息,可以同时兼顾栅格的平衡性、滑坡影响因子的丰富性以及计算量较小的优势,并便于在后续预测滑坡概率时,基于特征信息和滑坡影响因子,并以天为单位将滑坡影响因子输入至滑坡预测模型中,预测滑坡发生概率,将滑坡发生概率作为目标区域的滑坡预测结果,提高了滑坡预测的准确度和滑坡预测效率。
可选地,所述特征信息包括:所述滑坡影响因子的加权值、最大值、最小值、差分的均值和标准差特征中的至少一种。
可选地,所述特征信息提取模块420还用于,获取当前栅格的八邻域栅格内的滑坡影响因子的特征值;
根据所述八邻域栅格内的特征值,确定所述当前栅格的特征信息。
可选地,栅格尺寸确定模块410还用于,确定所述目标区域的每个栅格的初始尺寸和初始栅格数,以及确定滑坡栅格数和非滑坡栅格数,其中,所述滑坡栅格数是发生滑坡的栅格数量,所述非滑坡栅格数是未发生滑坡的栅格数量;
如果所述滑坡栅格数与非滑坡栅格数的比值处于设定平衡区间内,将所述初始尺寸作为所述第一栅格尺寸;
如果所述滑坡栅格数与所述非滑坡栅格数的比值未处于所述设定平衡区间内,调节所述目标区域的每个栅格的尺寸,重新确滑坡栅格数、非滑坡栅格数和当前栅格数;
如果重新确定的滑坡栅格数与非滑坡栅格数的比值处于所述设定平衡区间内,且所述当前栅格数与所述初始栅格数的比值小于设定阈值,将调节后的栅格的尺寸作为所述第一栅格尺寸。
可选地,栅格尺寸确定模块410还用于,获取所述目标区域的各雨量站的降雨数据、所述目标区域的每个栅格的初始尺寸和初始栅格数,以及确定滑坡栅格数和非滑坡栅格数,其中,所述滑坡栅格数是发生滑坡的栅格数量,所述非滑坡栅格数是未发生滑坡的栅格数量;
根据特定数量的雨量站的降雨数据计算任一当前栅格对应的雨量站的降雨预测值,基于所述降雨预测值与当前栅格对应的雨量站的实测雨量值计算降雨误差;
如果所述降雨误差小于设定误差阈值,且所述滑坡栅格数与所述非滑坡栅格数的比值处于设定平衡区间内,将所述初始尺寸作为所述第一栅格尺寸;
如果所述降雨误差大于设定误差阈值,调节所述目标区域的每个栅格的尺寸,根据调节后的所述任一栅格对应的雨量站的降雨预测值与实测雨量值重新计算降雨误差;
如果重新计算的降雨误差小于所述设定误差阈值,且重新确定的滑坡栅格数与重新确定的非滑坡栅格数的比值处于设定平衡区间内,且当前栅格数与所述初始栅格数的比值小于设定阈值,将调节后的栅格的尺寸作为所述第一栅格尺寸。
可选地,滑坡预测模块430还用于,基于第一预设模型,根据所述特征信息和所述滑坡影响因子确定每天滑坡的发生概率;
如果当前天滑坡的发生概率大于预设概率阈值,则基于第二预设模型,根据在所述特征信息和所述滑坡影响因子,确定各个栅格滑坡的发生概率。
可选地,所述滑坡影响因子包括动态因子和静态因子,其中,所述动态因子包括降雨量和土壤湿度中的至少一项,所述静态因子包括高程、坡度、坡向、平面曲率、剖面曲率、地形湿度指数、水流强度指数、沉积运输指数、地形粗糙度指数、与断层距离、与河流距离、与道路距离、岩性、土地利用和植被覆盖率中的至少一项。
本发明实施例所提供的滑坡的预测装置可执行本发明任意实施例所提供的滑坡的预测方法,具备执行方法相应的功能模块和有益效果。
实施例五
图5为本发明实施例五提供的一种滑坡的预测设备的结构示意图,如图5所示,该设备包括处理器510、存储器520、输入装置530和输出装置540;设备处理器510的数量可以是一个或多个,图5中以一个处理器510为例;设备中的处理器510、存储器520、输入装置530和输出装置540可以通过总线或其他方式连接,图5中以通过总线连接为例。
存储器520作为一种计算机可读存储介质,可用于存储软件程序、计算机可执行程序以及模块,如本发明实施例中的滑坡的预测方法对应的程序指令/模块(例如,滑坡的预测装置中的栅格尺寸确定模块410、特征信息提取模块420和滑坡预测模块430)。处理器510通过运行存储在存储器520中的软件程序、指令以及模块,从而执行设备的各种功能应用以及数据处理,即实现上述的滑坡的预测方法。
存储器520可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作***、至少一个功能所需的应用程序;存储数据区可存储根据终端的使用所创建的数据等。此外,存储器520可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中,存储器520可进一步包括相对于处理器510远程设置的存储器,这些远程存储器可以通过网络连接至设备/终端/服务器。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置530可用于接收输入的数字或字符信息,以及产生与设备的用户设置以及功能控制有关的键信号输入。输出装置540可包括显示屏等显示设备。
实施例六
发明实施例六还提供一种包含计算机可执行指令的存储介质,所述计算机可执行指令在由计算机处理器执行时用于执行一种滑坡的预测方法,该方法包括:
获取目标区域的第一栅格尺寸,并确定由所述第一栅格尺寸重采样后的第二栅格尺寸,其中,所述第二栅格尺寸大于所述第一栅格尺寸;
根据所述第一栅格尺寸将所述第二栅格尺寸对应的初始栅格划分成多个栅格,并基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的特征信息,其中,所述滑坡影响因子为所述第二栅格尺寸对应的栅格下每个栅格的设定时间段内的每天的滑坡数据;
基于滑坡预测模型,根据所述特征信息和所述滑坡影响因子确定各栅格滑坡的发生概率。
当然,本发明实施例所提供的一种包含计算机可执行指令的存储介质,其计算机可执行指令不限于如上所述的方法操作,还可以执行本发明任意实施例所提供的滑坡的预测方法中的相关操作。
通过以上关于实施方式的描述,所属领域的技术人员可以清楚地了解到,本发明可借助软件及必需的通用硬件来实现,当然也可以通过硬件实现,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如计算机的软盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、闪存(FLASH)、硬盘或光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
值得注意的是,上述用滑坡的预测装置的实施例中,所包括的各个单元和模块只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可;另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。
Claims (8)
1.一种滑坡的预测方法,其特征在于,包括:
获取目标区域的第一栅格尺寸,并确定由所述第一栅格尺寸重采样后的第二栅格尺寸,其中,所述第二栅格尺寸大于所述第一栅格尺寸;
根据所述第一栅格尺寸将所述第二栅格尺寸对应的初始栅格划分成多个栅格,并基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的特征信息,其中,所述滑坡影响因子为所述第二栅格尺寸对应的栅格下每个栅格的设定时间段内的每天的滑坡数据;
基于滑坡预测模型,根据所述特征信息和所述滑坡影响因子确定各栅格滑坡的发生概率;
第一栅格尺寸的确定方法,包括:
确定所述目标区域的每个栅格的初始尺寸和初始栅格数,以及确定滑坡栅格数和非滑坡栅格数,其中,所述滑坡栅格数是发生滑坡的栅格数量,所述非滑坡栅格数是未发生滑坡的栅格数量;
如果所述滑坡栅格数与非滑坡栅格数的比值处于设定平衡区间内,将所述初始尺寸作为所述第一栅格尺寸;
如果所述滑坡栅格数与所述非滑坡栅格数的比值未处于设定平衡区间内,调节所述目标区域的每个栅格的尺寸,重新确滑坡栅格数、非滑坡栅格数和当前栅格数;
如果重新确定的滑坡栅格数与非滑坡栅格数的比值处于所述设定平衡区间内,且所述当前栅格数与所述初始栅格数的比值小于设定阈值,将调节后的栅格的尺寸作为所述第一栅格尺寸;
所述基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的特征信息,包括:
获取当前栅格的八邻域栅格内的滑坡影响因子的特征值;
根据所述八邻域栅格内的特征值,确定所述当前栅格的特征信息。
2.根据权利要求1所述的预测方法,其特征在于,所述特征信息包括:所述滑坡影响因子的加权值、最大值、最小值、差分的均值和标准差特征中的至少一种。
3.根据权利要求1所述的预测方法,其特征在于,第一栅格尺寸的确定方法,包括:
获取所述目标区域的各雨量站的降雨数据、所述目标区域的每个栅格的初始尺寸和初始栅格数,以及确定滑坡栅格数和非滑坡栅格数,其中,所述滑坡栅格数是发生滑坡的栅格数量,所述非滑坡栅格数是未发生滑坡的栅格数量;
根据特定数量的雨量站的降雨数据计算任一当前栅格对应的雨量站的降雨预测值,基于所述降雨预测值与当前栅格对应的雨量站的实测雨量值计算降雨误差;
如果所述降雨误差小于设定误差阈值,且所述滑坡栅格数与所述非滑坡栅格数的比值处于设定平衡区间内,将所述初始尺寸作为所述第一栅格尺寸;
如果所述降雨误差大于设定误差阈值,调节所述目标区域的每个栅格的尺寸,根据调节后的任一栅格对应的雨量站的降雨预测值与实测雨量值重新计算降雨误差;
如果重新计算的降雨误差小于所述设定误差阈值,且重新确定的滑坡栅格数与重新确定的非滑坡栅格数的比值处于设定平衡区间内,且当前栅格数与所述初始栅格数的比值小于设定阈值,将调节后的栅格的尺寸作为所述第一栅格尺寸。
4.根据权利要求1所述的预测方法,其特征在于,所述基于滑坡预测模型,根据所述特征信息和所述滑坡影响因子确定各栅格滑坡的发生概率,包括:
基于第一预设模型,根据所述特征信息和所述滑坡影响因子确定每天滑坡的发生概率;
如果当前天滑坡的发生概率大于预设概率阈值,则基于第二预设模型,根据在所述特征信息和所述滑坡影响因子,确定各个栅格滑坡的发生概率。
5.根据权利要求1所述的预测方法,其特征在于,所述滑坡影响因子包括动态因子和静态因子,其中,所述动态因子包括降雨量和土壤湿度中的至少一项,所述静态因子包括高程、坡度、坡向、平面曲率、剖面曲率、地形湿度指数、水流强度指数、沉积运输指数、地形粗糙度指数、与断层距离、与河流距离、与道路距离、岩性、土地利用和植被覆盖率中的至少一项。
6.一种滑坡的预测装置,其特征在于,包括:
栅格尺寸确定模块,用于获取目标区域的第一栅格尺寸,并确定由所述第一栅格尺寸重采样后的第二栅格尺寸,其中,所述第二栅格尺寸大于所述第一栅格尺寸;
所述栅格尺寸确定模块,还用于确定所述目标区域的每个栅格的初始尺寸和初始栅格数,以及确定滑坡栅格数和非滑坡栅格数,其中,所述滑坡栅格数是发生滑坡的栅格数量,所述非滑坡栅格数是未发生滑坡的栅格数量;
如果所述滑坡栅格数与非滑坡栅格数的比值处于设定平衡区间内,将所述初始尺寸作为所述第一栅格尺寸;
如果所述滑坡栅格数与所述非滑坡栅格数的比值未处于所述设定平衡区间内,调节所述目标区域的每个栅格的尺寸,重新确滑坡栅格数、非滑坡栅格数和当前栅格数;
如果重新确定的滑坡栅格数与非滑坡栅格数的比值处于所述设定平衡区间内,且所述当前栅格数与所述初始栅格数的比值小于设定阈值,将调节后的栅格的尺寸作为所述第一栅格尺寸;
特征信息提取模块,用于根据所述第一栅格尺寸将所述第二栅格尺寸对应的初始栅格划分成多个栅格,并基于第一栅格尺寸提取第二栅格尺寸对应的栅格中滑坡影响因子的特征信息,其中,所述滑坡影响因子为所述第二栅格尺寸对应的栅格下每个栅格的设定时间段内的每天的滑坡数据;
所述特征信息提取模块,还用于获取当前栅格的八邻域栅格内的滑坡影响因子的特征值;
根据所述八邻域栅格内的特征值,确定所述当前栅格的特征信息;
滑坡预测模块,用于基于滑坡预测模型,根据所述特征信息和所述滑坡影响因子确定各栅格滑坡的发生概率,将所述发生概率作为所述目标区域的滑坡预测结果。
7.一种滑坡的预测设备,其特征在于,所述设备包括:
一个或多个处理器;
存储器,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如权利要求1-5中任一所述的滑坡的预测方法。
8.一种包含计算机可执行指令的存储介质,其特征在于,所述计算机可执行指令在由计算机处理器执行时用于执行如权利要求1-5中任一所述的滑坡的预测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011063199.4A CN112200355B (zh) | 2020-09-30 | 2020-09-30 | 一种滑坡的预测方法、装置、设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011063199.4A CN112200355B (zh) | 2020-09-30 | 2020-09-30 | 一种滑坡的预测方法、装置、设备和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112200355A CN112200355A (zh) | 2021-01-08 |
CN112200355B true CN112200355B (zh) | 2024-01-05 |
Family
ID=74013561
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011063199.4A Active CN112200355B (zh) | 2020-09-30 | 2020-09-30 | 一种滑坡的预测方法、装置、设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112200355B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112785078A (zh) * | 2021-02-02 | 2021-05-11 | 武汉中地云申科技有限公司 | 联合土壤水分信息和有效降雨的滑坡预测方法及终端 |
CN117974927B (zh) * | 2024-03-28 | 2024-06-25 | 煤炭科学研究总院有限公司 | 四维属性地质模型的煤矿灾害实时预警方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102915616A (zh) * | 2012-05-11 | 2013-02-06 | 新疆大学 | 一种区域栅格累加土地沙漠化预警方法 |
KR101721114B1 (ko) * | 2016-06-27 | 2017-03-30 | 서울대학교산학협력단 | 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법 |
CN110362949A (zh) * | 2019-07-23 | 2019-10-22 | 电子科技大学 | 一种基于神经网络的滑坡敏感性评估方法 |
CN111144651A (zh) * | 2019-12-26 | 2020-05-12 | 杭州鲁尔物联科技有限公司 | 一种地质灾害的预测方法、装置及设备 |
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/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102915616A (zh) * | 2012-05-11 | 2013-02-06 | 新疆大学 | 一种区域栅格累加土地沙漠化预警方法 |
KR101721114B1 (ko) * | 2016-06-27 | 2017-03-30 | 서울대학교산학협력단 | 위치정보가 포함된 포인트 데이터를 다축척의 웹 지도상에 클러스터링하기 위하여 격자의 크기를 결정하는 방법 |
CN110362949A (zh) * | 2019-07-23 | 2019-10-22 | 电子科技大学 | 一种基于神经网络的滑坡敏感性评估方法 |
CN111144651A (zh) * | 2019-12-26 | 2020-05-12 | 杭州鲁尔物联科技有限公司 | 一种地质灾害的预测方法、装置及设备 |
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 (zh) | 2021-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112200358B (zh) | 一种滑坡等级预测方法、装置、设备和存储介质 | |
CN112200362B (zh) | 一种滑坡的预测方法、装置、设备和存储介质 | |
Amini Parsa et al. | Spatio-temporal analysis of land use/land cover pattern changes in Arasbaran Biosphere Reserve: Iran | |
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 (zh) | 一种滑坡的预测方法、装置、设备和存储介质 | |
Zhao et al. | A Comparative Study of Landslide Susceptibility Mapping Using SVM and PSO‐SVM Models Based on Grid and Slope Units | |
CN111784044A (zh) | 一种滑坡的预测方法、装置、设备和存储介质 | |
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 (zh) | 一种基于rnn的滑坡灾害易发性预测方法 | |
CN115688404B (zh) | 基于svm-rf模型的降雨型滑坡预警方法 | |
CN116151483B (zh) | 一种区域级降雨滑坡概率性预测方法及预测终端 | |
CN112200356A (zh) | 一种滑坡的预测方法、装置、设备和存储介质 | |
Biard et al. | Automated detection of weather fronts using a deep learning neural network | |
CN112200355B (zh) | 一种滑坡的预测方法、装置、设备和存储介质 | |
Meena et al. | Assessing the importance of conditioning factor selection in landslide susceptibility for the province of Belluno (region of Veneto, northeastern Italy) | |
CN117556197B (zh) | 一种基于人工智能的台风涡旋初始化方法 | |
KR20200052398A (ko) | 머신 러닝 기법을 이용한 산사태 취약성 지도 작성 방법 및 장치 | |
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 | |
CN117708551B (zh) | 基于双精度gdp数据展布的洪涝灾害影响评估方法和*** | |
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 (zh) | 一种滑坡的预测方法、装置、设备和存储介质 | |
CN117233869B (zh) | 一种基于GRU-BiTCN的站点短期风速预测方法 | |
Momeni et al. | Pattern‐based calibration of cellular automata by genetic algorithm and Shannon relative entropy |
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 |