一种滑坡的预测方法、装置、设备和存储介质

技术领域

本发明实施例涉及滑坡监测技术领域，尤其涉及一种滑坡的预测方法、装置、设备和存储介质。

背景技术

山体滑坡是最常见的灾难性自然灾害之一，其具有分布范围广、发生频次高、多发性、区域性和严重性等特点，因滑坡每年会造成大量的人员伤亡和重大的环境和基础设施损失。对滑坡的易发性进行评估具有重要意义。

现有的滑坡易发性预测根据其所依据的理论基础的差异，可以分为确定性方法和非确定性方法。确定性方法主要是基于专家经验和知识的定向分析和基于滑坡过程或者物理学模型的进行分析的方法，预测准确度较差。近年来随着计算机技术和3S技术的高速发展，非确定性方法得到了广泛的应用，主要包括模糊逻辑法、层次分析法、决策树等。但是上述方法所处理的滑坡因子的时间精度较差，尤其是降雨量因子，通常仅考虑一月或一年的降雨量，因此预测精度也不太理想，同时，由于没有进行栅格尺寸确定，不能保证基于最优栅格下的滑坡数据进行滑坡预测，进而不能保证预测精度。

发明内容

本发明提供了一种滑坡的预测方法、装置、设备和存储介质，实现了基于目标区域的实际情况确定最优栅格，按每天的各个栅格的滑坡数据进行滑坡预测，提高了预测精度，同时基于滑坡预测模型进行滑坡预测，提高了预测效率。

第一方面，本发明实施例提供了一种滑坡的预测方法，该滑坡的预测方法包括：

获取目标区域的第一栅格尺寸和第二栅格尺寸，分别在所述第一栅格尺寸下和所述第二栅格尺寸下获取所述目标区域的滑坡影响因子，其中，所述第二栅格尺寸大于所述第一栅格尺寸；

分别基于所述第一栅格尺寸下的滑坡影响因子和第二栅格尺寸下的滑坡影响因子的分布特征，确定所述第二栅格尺寸是否为期望栅格尺寸，其中，所述期望栅格尺寸小于所述第二栅格尺寸且大于或等于所述第一栅格尺寸；

如果否，基于特定尺度因子，迭代调整所述第二栅格尺寸；

如果当前迭代次数对应的当前栅格尺寸为所述期望栅格尺寸，将所述期望栅格尺寸对应的各栅格的每天的滑坡影响因子输入至滑坡预测模型，基于所述滑坡预测模型确定各栅格滑坡的发生概率。

第二方面，本发明实施例还提供了一种滑坡的预测装置，该滑坡的预测装置包括：

数据获取模块，用于获取目标区域的第一栅格尺寸和第二栅格尺寸，分别在所述第一栅格尺寸下和所述第二栅格尺寸下获取所述目标区域的滑坡影响因子，其中，所述第二栅格尺寸大于所述第一栅格尺寸；

第二栅格尺寸确定模块，用于分别基于所述第一栅格尺寸下的滑坡影响因子和第二栅格尺寸下的滑坡影响因子的分布特征，确定所述第二栅格尺寸是否为期望栅格尺寸，其中，所述期望栅格尺寸小于所述第二栅格尺寸且大于或等于所述第一栅格尺寸；

第二栅格尺寸调整模块，用于如果所述第二栅格尺寸不是期望栅格尺寸，基于特定尺度因子，迭代调整所述第二栅格尺寸；

滑坡预测模块，用于如果当前迭代次数对应的当前栅格尺寸为所述期望栅格尺寸，将所述期望栅格尺寸对应的各栅格的每天的滑坡影响因子输入至滑坡预测模型，基于所述滑坡预测模型确定各栅格滑坡的发生概率。

第三方面，本发明实施例还提供了一种滑坡的预测设备，该设备包括：

一个或多个处理器；

存储器，用于存储一个或多个程序；

当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现本发明任意实施例提供的滑坡的预测方法。

第四方面，本发明实施例还提供了一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行本发明任意实施例提供的滑坡的预测方法。

本发明实施例的技术方案，获取目标区域的第一栅格尺寸和第二栅格尺寸，分别在所述第一栅格尺寸下和所述第二栅格尺寸下获取所述目标区域的滑坡影响因子，其中，所述第二栅格尺寸大于所述第一栅格尺寸，分别基于所述第一栅格尺寸下的滑坡影响因子和第二栅格尺寸下的滑坡影响因子的分布特征，确定所述第二栅格尺寸是否为期望栅格尺寸；如果否，基于特定尺度因子，迭代调整所述第二栅格尺寸；如果当前迭代次数对应的当前栅格尺寸为所述期望栅格尺寸，将所述期望栅格尺寸对应的各栅格的每天的滑坡影响因子输入至滑坡预测模型，基于所述滑坡预测模型确定各栅格滑坡的发生概率。可以根据不同栅格尺寸下的滑坡影响因子的分布特征进行栅格尺寸划分，使不同栅格的栅格尺寸不同，使滑坡影响因子在划分后的栅格尺寸下同时兼顾保留丰富的特征信息并精简数据的冗余度的优势；在后续预测滑坡概率时，根据不同栅格尺寸对应的栅格内的滑坡影响因子，以天为单位将滑坡影响因子输入至滑坡预测模型中，预测滑坡发生概率，提高了滑坡预测的准确度。

附图说明

图1是本发明实施例一中的一种滑坡的预测方法的流程图；

图2是本发明实施例二中的一种滑坡的预测方法的流程图；

图3是本发明实施例三中的一种滑坡的预测方法的流程图；

图4是本发明实施例四中的一种滑坡的预测装置的结构示意图；

图5是本发明实施例五中的一种滑坡的预测设备的结构示意图。

具体实施方式

下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。

实施例一

图1是本发明实施例一提供的一种滑坡的预测方法的流程图，本实施例可适用于对滑坡易发性进行评估的情况，该方法可以由滑坡的预测装置来执行，如图1所示，该方法包括如下步骤：

S110、获取目标区域的第一栅格尺寸和第二栅格尺寸，分别在第一栅格尺寸下和第二栅格尺寸下获取目标区域的滑坡影响因子。

其中，所述第二栅格尺寸大于所述第一栅格尺寸。目标区域通常是发生过滑坡的区域，也可以是指定的任意一个区域。栅格就是将空间分割成有规律的网格，每一个网格成为一个栅格或单元，并在各个单元上赋予相应的属性值来表示实体的一种数据形式。所述滑坡影响因子即为栅格数据。所述滑坡影响因子可以是目标区域的各个时间段内每天的数据。设定时间段可以是一天、三天、一周、一个月或者其他时间段。具体的，可以在目标区域设置多个监测点，以实时获取各个监测点的滑坡影响因子。进一步结合预设部门收集的数据，形成该目标区域的各个栅格的设定时间段内每天的滑坡影响因子。

所述滑坡影响因子包括动态因子和静态因子，其中，所述动态因子包括降雨量和土壤湿度中的至少一项，所述静态因子包括高程、坡度、坡向、平面曲率、剖面曲率、地形湿度指数、水流强度指数、沉积运输指数、地形粗糙度指数、与断层距离、与河流距离、与道路距离、岩性、土地利用、植被覆盖率中的至少一项。其中，所述第一栅格尺寸可以是目标区域的最小尺寸，在第一栅格尺寸下的滑坡栅格可以很好的保留滑坡影响因子的丰富性。所述第二栅格尺寸可以是目标区域的最大尺寸，将目标区域的栅格划分成第二栅格尺寸可以减小数据冗余度并降低计算量。例如，所述第一栅格尺寸为30m*30m，所述第二栅格尺寸为1km*1km。

S120、分别基于第一栅格尺寸下的滑坡影响因子和第二栅格尺寸下的滑坡影响因子的分布特征，确定第二栅格尺寸是否为期望栅格尺寸。

可以理解的是，滑坡影响因子在不同尺寸下的栅格的表征程度不同，尺寸较小的栅格能够保留滑坡影响因子的丰富性，尺寸较大的栅格能够精简数据的冗余度。目标区域包括的各细分区域的地形、经纬度信息以及地势信息均不同(例如高程与坡度均呈块状分布，与湖泊距离呈圈状分布，与铁路距离呈线条分布)，将目标区域的栅格划分为统一尺寸，不能同时兼顾保留丰富的特征信息并精简数据的冗余度的优势。为了降低滑坡预测过程的计算量，同时兼顾滑坡影响因子的丰富性，本实施例可以根据目标区域的第一栅格尺寸和第二栅格尺寸下的滑坡影响因子的分布特征，以第一栅格尺寸为基础对第二栅格尺寸进行调节，以确定目标区域的最优栅格尺寸。

其中，所述期望栅格尺寸小于第二栅格尺寸且大于或等于第一栅格尺寸。可选地，确定所述第二栅格尺寸是否为期望栅格尺寸的第一种方法包括如下步骤：基于所述第一栅格尺寸下的滑坡影响因子的分布特征，对所述第二栅格尺寸下的滑坡影响因子的分布特征进行显著性检验；如果所述第二栅格尺寸下的滑坡影像因子通过显著性检验，确定所述第二栅格尺寸为所述期望栅格尺寸，否则，所述第二栅格尺寸不是所述期望栅格尺寸。

具体地，可以基于单因素非参数方差分析方法和双因素秩方差分析方法，第一栅格尺寸下的滑坡影响因子的分布特征为基础，对第二栅格尺寸下的滑坡影响因子的分布特征进行非参数显著性检验，以判断第二栅格尺寸下的滑坡影响因子的分布特征与第一栅格尺寸下的滑坡影响因子的分布特征是否具有显著性差异，如果第二栅格尺寸下的滑坡影响因子的分布特征与第一栅格尺寸下的滑坡影响因子的分布特征具有明显差异，则确定第二栅格尺寸不是期望栅格尺寸，如果第二栅格尺寸下的滑坡影响因子的分布特征与第一栅格尺寸下的滑坡影响因子的分布特征不具有明显差异，则确定第二栅格尺寸是期望栅格尺寸。

可选地，确定所述第二栅格尺寸是否为期望栅格尺寸的第二种方法包括如下步骤：基于所述第一栅格尺寸下的滑坡影响因子的分布特征计算第一特征数据，基于所述第二栅格尺寸下的滑坡影响因子的分布特征计算第二特征数据；其中，所述第一特征数据和所述第二特征数据包括方差、均值中的至少一个；将所述第一特征数据与所述第二特征数据之间的差值与设定阈值比较；如果所述差值小于所述设定阈值，确定所述第二栅格尺寸为所述期望栅格尺寸，否则，所述第二栅格尺寸不是所述期望栅格尺寸。

其中，设定阈值可以是一个极小的百分数，例如1％、2％。如果所述差值小于所述设定阈值，则第二栅格尺寸下的滑坡影响因子的分布特征与第一栅格尺寸下的滑坡影响因子的分布特征接近，无需调整第二栅格尺寸，直接将第二栅格尺寸作为期望栅格尺寸；如果所述差值不小于所述设定阈值，则第二栅格尺寸下的滑坡影响因子的分布特征与第一栅格尺寸下的滑坡影响因子的分布特征差异较大，需要将第二栅格尺寸调整至期望栅格尺寸。

S130、如果第二栅格尺寸不是期望栅格尺寸，基于特定尺度因子，迭代调整第二栅格尺寸。

S140、如果当前迭代次数对应的当前栅格尺寸为期望栅格尺寸，将期望栅格尺寸对应的各栅格的每天的滑坡影响因子输入至滑坡预测模型，基于滑坡预测模型确定各栅格滑坡的发生概率。

其中，所述特定尺度因子可以是栅格尺寸的缩小梯度，例如，特定尺度因子是2；可选地，所述基于特定尺度因子，迭代调整所述第二栅格尺寸，得到当前迭代次数下的当前栅格尺寸，包括：按照所述尺度因子迭代缩小所述第二栅格尺寸，得到所述当前迭代次数下的当前栅格尺寸。

具体地，按照所述尺度因子迭代缩小所述第二栅格尺寸，并确定当前迭代次数下的第一栅格尺寸下的滑坡影响因子和第二栅格尺寸下的滑坡影响因子的分布特征，根据两个尺寸下的分布特征确定当前迭代次数下的第二栅格尺寸是否为期望栅格尺，直至当前迭代次数对应的当前栅格尺寸为期望栅格尺寸，将期望栅格尺寸对应的各栅格的每天的滑坡影响因子输入至滑坡预测模型，基于所述滑坡预测模型确定各栅格滑坡的发生概率。

示例性地，第一栅格尺寸是1km*1km，第二栅格尺寸是30m*30m，特定尺度因子是2。以30m*30m下的滑坡影响因子的分布特征为基础，对1km*1km下的滑坡影响因子的分布特征进行非参数显著性检验，如果1km*1km下的滑坡影响因子的分布特征与30m*30m下的滑坡影响因子的分布特征具有显著性差异，基于特定尺度因子缩小1km*1km的栅格，缩小第一次得到的当前栅格尺寸为500m*500m；再以30m*30m下的滑坡影响因子的分布特征为基础，对500m*500m下的滑坡影响因子的分布特征进行非参数显著性检验，如果500m*500m下的滑坡影响因子的分布特征与30m*30m下的滑坡影响因子的分布特征具有显著性差异，继续基于特定尺度因子缩小500m*500m的栅格，直至得到的当前迭代次数对应的当前栅格尺寸为所述期望栅格尺寸，结束调整第二栅格尺寸的过程。其中，所述期望栅格尺寸可以是30m*30m或者是小于500m*500m且大于30m*30m的尺寸。

通过上述方式，根据不同栅格尺寸下的滑坡影响因子的分布特征进行栅格尺寸划分，使不同栅格的栅格尺寸不同，使滑坡影响因子在划分后的栅格尺寸下表征效果最好，同时兼顾保留丰富的特征信息并精简数据的冗余度的优势，有利于提高滑坡预测的准确度和效率。

其中，所述滑坡预测模型可以是单一神经网络模型，也可以是输入端的和输出端连接形成的复杂神经网络模型。具体地，滑坡预测模型的训练过程为：从滑坡野外勘察报告、典型滑坡监测报告等文件提取目标区域或全部区域的滑坡的基本信息包括滑坡发生的时间、经纬度、受灾规模等滑坡信息，并根据该滑坡信息确定各个滑坡发生地的历史时间段的滑坡影响因子，以设定比例组成训练集和验证集，如8:2或7:3；将各个滑坡影响因子进行数据校正与配准，并进行坐标统一处理和栅格统一处理，从而得到栅格大小一致的滑坡影响因子，并将通过上述步骤处理后的每天的滑坡影响因子训练滑坡预测模型。其中，该滑坡影响因子通过特征工程提取各个滑坡影响因子的特征，组成滑坡预测模型的输入特征矩阵；初始化滑坡预测模型的参数，并将该输入特征矩阵输入滑坡预测模型中，进行模型训练，得到历史时间段的历史滑坡信息，并基于F1-值(F1-Score)与ROC(Receiver OperatingCharacteristic，接收者操作特征)作为评价指标，根据评价结果进行滑坡预测模型参数的调整，当满足时，训练初步完成，通过验证集进行模型验证，验证通过则得到训练好的滑坡预测模型。其中，所述历史滑坡信息包括各栅格发生滑坡的概率。

本发明实施例的技术方案，获取目标区域的第一栅格尺寸和第二栅格尺寸，分别在所述第一栅格尺寸下和所述第二栅格尺寸下获取所述目标区域的滑坡影响因子，其中，所述第二栅格尺寸大于所述第一栅格尺寸，分别基于所述第一栅格尺寸下的滑坡影响因子和第二栅格尺寸下的滑坡影响因子的分布特征，确定所述第二栅格尺寸是否为期望栅格尺寸；如果否，基于特定尺度因子，迭代调整所述第二栅格尺寸；如果当前迭代次数对应的当前栅格尺寸为所述期望栅格尺寸，将所述期望栅格尺寸对应的各栅格的每天的滑坡影响因子输入至滑坡预测模型，基于所述滑坡预测模型确定各栅格滑坡的发生概率。可以根据不同栅格尺寸下的滑坡影响因子的分布特征进行栅格尺寸划分，使滑坡影响因子在划分后的栅格尺寸下同时兼顾保留丰富的特征信息并精简数据的冗余度的优势；在后续预测滑坡概率时，根据不同栅格尺寸对应的栅格内的滑坡影响因子，以天为单位将滑坡影响因子输入至滑坡预测模型中，预测滑坡发生概率，，提高了滑坡预测的准确度。

实施例二

图2为本发明实施例二提供的一种滑坡的预测方法的流程图，本实施例是对上一实施例的进一步细化，所述将所述期望栅格尺寸对应的各栅格的每天的滑坡影响因子输入至滑坡预测模型，基于所述滑坡预测模型确定各栅格滑坡的发生概率，包括：基于第一预设模型，根据当前栅格尺寸对应的各栅格的每天的滑坡影响因子，确定每天滑坡的发生概率；如果当前天滑坡的发生概率大于预设概率阈值，则基于第二预设模型，根据当前栅格尺寸对应的各栅格的每天的滑坡影响因子，确定各个栅格滑坡的发生概率。在该方法实施例中未详尽描述的部分请参考上述实施例。具体参见图2所示，该方法可以包括如下步骤：

S210、获取目标区域的第一栅格尺寸和第二栅格尺寸，分别在第一栅格尺寸下和第二栅格尺寸下获取目标区域的滑坡影响因子。

S220、分别基于第一栅格尺寸下的滑坡影响因子和第二栅格尺寸下的滑坡影响因子的分布特征，确定第二栅格尺寸是否为期望栅格尺寸。

S230、如果第二栅格尺寸不是期望栅格尺寸，基于特定尺度因子，迭代调整第二栅格尺寸。

S240、如果当前迭代次数对应的当前栅格尺寸为期望栅格尺寸，基于第一预设模型，根据当前栅格尺寸对应的各栅格的每天的滑坡影响因子，确定每天滑坡的发生概率。

其中，第一预设模型可以是神经网络模型，或者其他学习算法。示例性的，第一预设模型可以是支持向量机算法(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)等。在对第一预设模型进行训练时，获取各个滑坡发生地的历史时间段的历史滑坡数据，历史滑坡数据中包括各个滑坡预测因子并提取各个滑坡预测因子的特征，组成第一预设模型的输入特征矩阵；初始化第一预设模型的参数，并将该输入特征矩阵输入第一预设模型中，进行模型训练，得到历史时间段对应的每天的滑坡预测概率，并基于F1-值(F1-score)与ROC(receiver operating characteristic，接收者操作特征)作为评价指标，根据评价结果进行第一预设模型参数的调整，直至第一预设模型训练完成。所述第一预设模型的训练方法与前述实施例中的滑坡预测模型的训练方法一样，此处不再详细描述。

可选地，所述基于第一预设模型，根据当前栅格尺寸对应的各栅格的每天的滑坡影响因子，确定每天滑坡的发生概率，包括：计算所述当前栅格尺寸对应的栅格的动态因子的构造特征，其中，所述构造特征包括各个栅格的每个动态因子的总和、平均值、最大值、最小值、极差、四分位数和雨量时间，并根据每天的所述构造特征和静态因子得到所述第一预设模型的第一特征矩阵；将所述第一特征矩阵输入所述第一预设模型中，以根据所述第一预设模型确定每天滑坡的发生概率。

具体地，所述滑坡影响因子包括：土壤湿度、高程、坡度、坡向、平面曲率、剖面曲率、地形湿度指数、水流强度指数、沉积运输指数、地形粗糙度指数、与断层距离、与河流距离、与道路距离、岩性、土地利用和植被覆盖率中的至少一项。所述降雨量和所述土壤湿度可以作为预测滑坡的发生概率的动态因子，所述高程、坡度、坡向、平面曲率、剖面曲率、地形湿度指数、水流强度指数、沉积运输指数、地形粗糙度指数、与断层距离、与河流距离、与道路距离、岩性、土地利用和植被覆盖率可以作为预测滑坡的发生概率的静态因子，对上述动态因子构造特征并结合静态因子确定第一特征矩阵，基于所述第一特征矩阵和第一预设模型确定每天滑坡的发生概率。

示例性地，所获取的每天的各个栅格的滑坡影响因子X为：X＝{X

S250、如果当前天滑坡的发生概率大于预设概率阈值，则基于第二预设模型，根据当前栅格尺寸对应的各栅格的每天的滑坡影响因子，确定各个栅格滑坡的发生概率。

其中，预设概率阈值可以是0.5、0.6或者其他值，当然预设概率阈值也可以采用分数或者百分数进行表示。当滑坡的发生概率大于该预设概率阈值时，表示当前天滑坡的易发性较高。通过设置预设概率阈值对滑坡影响因子进行筛选，只有当前天滑坡的发生概率大于预设概率阈值时，才将当前天的数据传输至第二预设模型进行进一步预测，大大减少了模型输入的数据量，提高了处理效率，同时，提高了预测的精度。当前天的个数可以是1个或多个。即通过第一预设模型预测设定时间段内每天的滑坡发生概率，当某一天的滑坡发生概率大于预设概率阈值时，则将该天对应的滑坡预设数据或第一特征矩阵发送至第二预设模型，以进行该天的各个栅格的滑坡发送概率的预测。

可选地，所述基于第二预设模型，根据当前栅格尺寸对应的各栅格的每天的滑坡影响因子，确定各个栅格滑坡的发生概率，包括：对所述当前栅格尺寸对应的各栅格的每天的滑坡影响因子进行特征提取，以生成所述第二预设模型的第二特征矩阵；对所述第二特征矩阵进行归一化处理；将归一化处理后的所述第二特征矩阵输入所述第二预设模型中，以根据所述第二预设模型确定当前天的各个栅格滑坡的发生概率。

具体地，假设第c天的滑坡发生概率大于预测概率阈值，则将第c天的滑坡影响因子X

进一步地，可以基于max-min归一化算法对第二特征矩阵的各个特征进行归一化处理。当然也可以选取其他归一化算法进行归一化处理。本发明实施例对第一特征矩阵和第二特征矩阵归一化算法不进行限定。在归一化处理之后，可以基于多重共线性方法对第二特征矩阵进行特征选择，以根据共线性程度进行特征筛选，得到筛选的后的第二特征矩阵。

具体的，第二预设模型可以是一维卷积神经网络模型，包括卷积层、批标准化层、激活函数和优化层，其中，卷积层的卷积核个数的取值范围为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。可选地，在将所述样本遥感影像信息输入至初始森林模型之前，还可以对样本遥感影像信息进行预处理，例如，对样本遥感影像信息进行大气校正、辐射校正、边缘检测、假彩色合成处理等，以提高样本遥感信息的功能区分类模型的训练精度。

可以理解的是，功能区分类模型输出的各功能区包括目标区域的所有功能区，包括城镇建筑、城市绿地、水体、农田、裸土、山林等功能区，各功能区的滑坡发生概率不同。例如，城镇建筑和城镇绿地发生滑坡的概率极小，水体和山林发生滑坡的概率较大。需要说明的是，通过将各功能区对应的历史预测概率与历史概率阈值比较，筛选出目标功能区，对目标功能区进行针对性的滑坡预测，可以减小计算量，有利于提高目标区域的滑坡预测效率。

若通过S310确定目标区域的各功能区，并根据各功能区对应的标签筛选目标功能区，则S320可以替换为：获取目标区域的目标功能区的第一栅格尺寸和第二栅格尺寸，分别在第一栅格尺寸下和第二栅格尺寸下获取目标区域的滑坡影响因子。

S320、分别基于第一栅格尺寸下的滑坡影响因子和第二栅格尺寸下的滑坡影响因子的分布特征，确定第二栅格尺寸是否为期望栅格尺寸。

S330、如果第二栅格尺寸不是期望栅格尺寸，基于特定尺度因子，迭代调整第二栅格尺寸。

S340、如果当前迭代次数对应的当前栅格尺寸为期望栅格尺寸，基于第一预设模型，根据当前栅格尺寸对应的各栅格的每天的滑坡影响因子，确定每天滑坡的发生概率。

S350、如果当前天滑坡的发生概率大于预设概率阈值，则基于第二预设模型，根据当前栅格尺寸对应的各栅格的每天的滑坡影响因子，确定各个栅格滑坡的发生概率。

如前述实施例所述，滑坡影响因子包括动态因子和静态因子。所述动态因子包括降雨量和土壤湿度中的至少一项。所述降雨量可以通过将获取的目标区域的地理环境数据输入至降雨插值模型确定。所述地理环境数据包括地理位置数据、大气数据、地形数据和下垫面数据。所述降雨插值模型可以是一种反向传播模型(BP，back propagation neuralnetwork)，基于误差反向传播算法的由非线性变换单元组成的多层前馈型网络，BP一般由输入层、隐含层和输出层组成，每层中又包括N个神经元，同层神经元之间互相独立，每层之间神经元的输出经过特定的激励函数后，只会影响下层神经元的输入。

具体地，降雨插值模型的训练方法为：获取初始模型，确定所述初始模型的初始权值矩阵和初始阈值，并计算所述初始模型的适应度；将历史时间段内的样本地理环境信息输入至所述初始模型，确定预测降雨量，根据所述预测降雨量和所述历史时间段的实测降雨量，并基于遗传算法迭代调整所述初始模型的初始权值矩阵和初始阈值；基于迭代调整后的权值矩阵阈值，调整所述初始模型并计算调整后的模型的适应度，直至重新计算的适应度达到期望值，将达到所述期望值的适应度对应的模型作为所述降雨插值模型。

具体地，所述初始权值矩阵可以包括输入层与隐含层之间的连接权值和隐含层与输出层之间的连接权值；初始阈值可以包括隐含层的阈值和输出层的阈值。初始模型的适应度的计算公式为：

可选地，采用二进制编码、实数编码或者格雷码编码方式对权值和阈值进行编码，并根据权值编码和阈值编码确定调整后的权值矩阵和阈值。可选地，还可以通过计算预测降雨量和实测降雨量之间的平均误差、平均绝对误差或线性相关系数，并根据所述平均误差、所述平均绝对误差或所述线性相关系数迭代调整初始模型的权值编码和阈值编码，直至平均误差、平均绝对误差或线性相关系数不小于预设误差阈值，说明初始模型调整至稳定状态，并进一步计算调整后模型的适应度以及确定降雨插值模型。

需要说明的是，历史时间段内的样本地理环境数据是多维度的，在训练降雨插值模型时，多维度的样本地理环境数据可以提高降雨插值模型的精度和鲁棒性，因而，在获取每天的地理环境数据后，基于降雨插值模型预测可以提高降雨概率的预测精度，有利于提高滑坡预测的预测精度。

为了提高滑坡预测的准确性，还可以对当前栅格尺寸进行重采样，结合重采样的栅格尺寸的特征信息进行滑坡预测。具体方法包括：将当前栅格尺寸作为第一尺寸，并确定第一尺寸重采样后的第二尺寸，所述第二尺寸大于所述第一尺寸；根据第一尺寸将第二尺寸对应的栅格划分成多个栅格，并基于第一尺寸提取第二尺寸的栅格下的滑坡影响因子的特征信息，将特征信息和滑坡影响因子输入至第一预设模型以及第二预设模型进行滑坡预测。其中，特征信息的提取方法为：获取当前栅格的八邻域栅格内的滑坡影响因子的特征值；根据所述八邻域栅格内的特征值，确定所述当前栅格的特征信息。

本实施例通过将特征信息和滑坡影响因子结合确定各栅格滑坡的发生概率，可以增加滑坡影响因子的数据量，根据特征信息可以更准确的确定各栅格滑坡的发生概率，提高了滑坡预测的准确性。

与前述实施例不同的，第一预设模型和第二预设模型均可以是自步分类学习模型，所述第一预设模型和第二预设模型均可以根据历史滑坡数据、历史非滑坡数据迭代调整自步因子和下采样比例确定。所述自步因子根据所述历史非滑坡数据的分箱个数和迭代次数确定，所述下采样比例根据每个分箱的自步因子和所有分箱的自步因子确定。可选地，所述历史滑坡数据和历史非滑坡数据的确定方法，包括：确定历史区域的各栅格对应的标签，并确定各所述历史时间点下的滑坡栅格和非滑坡栅格分别对应的历史降雨量和特定相关因子；基于所述标签，生成包括滑坡栅格和非滑坡栅格分别对应的历史降雨量和特定相关因子的滑坡数据集，根据所述滑坡数据集中各个栅格对应的标签确定所述历史滑坡数据和历史非滑坡数据。

可选地，所述历史滑坡数据和历史非滑坡数据的确定方法，包括：确定历史区域的各栅格对应的标签，并确定各所述历史时间点下的滑坡栅格和非滑坡栅格分别对应的历史降雨量和特定相关因子；基于所述标签，生成包括滑坡栅格和非滑坡栅格分别对应的历史降雨量和特定相关因子的滑坡数据集，根据所述滑坡数据集中各个栅格对应的标签确定所述历史滑坡数据和历史非滑坡数据。

可选地，所述第一预设模型的训练方法，包括：获取初始预测模型，并获取各历史滑坡时间点下的滑坡栅格的降雨量和非滑坡栅格的降雨量，以及各所述标签下的特定相关因子；按天将所述历史滑坡时间点下的滑坡栅格的降雨量、所述非滑坡栅格的降雨量以及各所述标签下的特定相关因子输入至所述初始预测模型，基于初始预测模型输出的滑坡概率和所述历史滑坡时间点对应的滑坡概率调整所述初始预测模型的损失函数，直至所述损失函数达到设定阈值，得到所述第一预设模型。

其中，所述基于初始预测模型输出的滑坡概率和所述历史滑坡时间点对应的滑坡概率调整所述初始预测模型的损失函数，直至所述损失函数达到设定阈值，得到所述第一预设模型，包括：基于所述历史滑坡时间点下的滑坡栅格的降雨量、所述非滑坡栅格的降雨量、各所述标签以及各所述标签下的特定相关因子，确定所述初始预测模型的初始分类硬度；根据所述初始分类硬度确定所述非滑坡栅格的降雨量和非滑坡栅格的特定相关因子的分箱个数，基于所述分箱个数确定所述初始预测模型的各分箱的自步因子；基于所述自步因子确定每个分箱的下采样比例，并基于所述下采样比例确定各分箱内的下采样的非滑坡栅格的降雨量，以及确定下采样的非滑坡栅格的特定相关因子；将所述下采样的非滑坡栅格的降雨量、下采样的非滑坡栅格的特定相关因子以及滑坡栅格的降雨量输入至所述初始预测模型，基于所述初始预测模型输出的滑坡概率和所述历史滑坡时间点对应的滑坡概率确定所述损失函数，基于所述损失函数在单个样本上的训练分类硬度，迭代调整自步因子和每个分箱的下采样比例；根据迭代调整的所述自步因子和所述每个分箱的下采样比例调整所述初始预测模型，直至损失函数达到设定阈值，得到所述第一预设模型。

其中，每个分箱的下采样比例为每个分箱的自步因子与所有分箱的自步因子的比值，所述自步因子通过将分箱个数、分箱标签数的相反数以及迭代次数之和加1得到。

选择自步分类学习模型进行滑坡预测的好处是：根据历史非滑坡数据的分箱个数和迭代次数确定自步分类学习模型的自步因子，根据每个分箱的自步因子和所有分箱的自步因子确定下采样比例，因此，在基于历史滑坡数据、历史非滑坡数据迭代调整自步因子和下采样比例时，随着迭代次数的增加，每个分箱的采样比例由硬度低的下采样数量多到最后每个分箱内的下采样数量逐渐变化均匀，并且，在每个箱子的下采样比例均匀变化的同时，还可以保证硬度小的分箱的下采样数量始终高于硬度大的分箱的下采样数量，得到多样化、鲁棒性高且包容性强的第一预设模型。

可选地，确定当前天的各个栅格滑坡的发生概率之后，还可以将所述当前天的各个栅格的发生概率与所述目标区域的预警阈值进行比较，基于得到的各个栅格的发生概率对应的预警等级确定所述目标区域的滑坡的发生等级，其中，所述预警阈值根据各区域栅格滑坡的第一系数和第二系数确定。

其中，所述第一系数和第二系数的确定方法包括：获取各区域栅格滑坡的设定时间段内的已发生滑坡的历史概率；基于所述历史概率的分布特征确定第一已确定概率区间和第二已确定概率区间；分别计算所述第一已确定概率区间和所述第二已确定概率区间内的同一风险等级的概率阈值与对应系数的乘积之和，得到中间已确定概率；如果所述中间已确定概率未达到设定评价指标，迭代调整各风险等级的概率阈值的对应系数，直至所述中间已确定概率达到所述设定评价指标，根据达到所述设定评价指标的对应系数确定第一系数和第二系数。

其中，确定所述预警阈值的方法，包括：基于所述各个栅格滑坡的发生概率的分布特征确定第一预测概率区间和第二预测概率区间；计算所述第一预测概率区间内各风险等级与所述第一系数的第一乘积，以及计算所述第二预测概率区间内各风险等级与所述第二系数的第二乘积，将所述第一乘积与所述第二乘积的和作为各风险等级的预警阈值。

其中，所述基于所述历史概率的分布特征确定所述第一已确定概率区间和所述第二已确定概率区间，包括：确定所述历史概率和所述历史概率对应的次数；根据所述历史概率和各所述历史概率对应的滑坡发生次数，确定所述历史概率的密集度分布特征和断点分布特征；基于所述密集度分布特征确定所述第一已确定概率区间，并基于所述断点分布特征确定所述第二已确定概率区间。

可以理解的是，上述方式基于各区域栅格滑坡的设定时间段内的已发生滑坡的历史概率，根据历史概率的分布特征确定第一已确定概率区间和第二已确定概率区间，根据第一已确定概率区间、第二已确定概率区间内的各风险等级的概率阈值以及设定评价指标确定各区域的第一系数和第二系数，当获取不同区域的发生概率后，可以根据不同区域对应的第一系数、第二系数和发生概率灵活确定不同区域的预警阈值，有利于后续根据各区域对应的预警阈值准确确定该区域的风险等级。

其中，所述基于得到的各个栅格的发生概率对应的预警等级确定所述目标区域的滑坡的发生等级，包括：筛选大于所述预警阈值中的第一等级的目标等级；计算所述目标等级对应的栅格的平均等级；基于所述平均等级、设定系数、目标等级对应的栅格数与目标区域内所有栅格数的比值，确定所述目标区域的滑坡的发生等级。

其中，目标区域的滑坡的发生等级的计算公式为：

其中，α为设定系数，通过对目标区域的滑坡的历史发生等级进行贝叶斯计算得到，

实施例四

图4为本发明实施例四提供的一种滑坡的预测装置的结果示意图，如图4所述，该滑坡的预测装置包括：数据获取模块410、第二栅格尺寸确定模块420、第二栅格尺寸调整模块430和滑坡预测模块440。

其中，数据获取模块410，用于获取目标区域的第一栅格尺寸和第二栅格尺寸，分别在所述第一栅格尺寸下和所述第二栅格尺寸下获取所述目标区域的滑坡影响因子，其中，所述第二栅格尺寸大于所述第一栅格尺寸；

栅格尺寸确定模块420，用于分别基于所述第一栅格尺寸下的滑坡影响因子和第二栅格尺寸下的滑坡影响因子的分布特征，确定所述第二栅格尺寸是否为期望栅格尺寸，其中，所述期望栅格尺寸小于所述第二栅格尺寸且大于或等于所述第一栅格尺寸；

第二栅格尺寸调整模块430，用于如果所述第二栅格尺寸不是期望栅格尺寸，基于特定尺度因子，迭代调整所述第二栅格尺寸；

滑坡预测模块440，用于如果当前迭代次数对应的当前栅格尺寸为所述期望栅格尺寸，将所述期望栅格尺寸对应的各栅格的每天的滑坡影响因子输入至滑坡预测模型，基于所述滑坡预测模型确定各栅格滑坡的发生概率。

本发明实施例的技术方案，获取目标区域的第一栅格尺寸和第二栅格尺寸，分别在所述第一栅格尺寸下和所述第二栅格尺寸下获取所述目标区域的滑坡影响因子，其中，所述第二栅格尺寸大于所述第一栅格尺寸，分别基于所述第一栅格尺寸下的滑坡影响因子和第二栅格尺寸下的滑坡影响因子的分布特征，确定所述第二栅格尺寸是否为期望栅格尺寸；如果否，基于特定尺度因子，迭代调整所述第二栅格尺寸；如果当前迭代次数对应的当前栅格尺寸为所述期望栅格尺寸，将所述期望栅格尺寸对应的各栅格的每天的滑坡影响因子输入至滑坡预测模型，基于所述滑坡预测模型确定各栅格滑坡的发生概率。可以根据不同栅格尺寸下的滑坡影响因子的分布特征进行栅格尺寸划分，使不同栅格的栅格尺寸不同，使滑坡影响因子在划分后的栅格尺寸下同时兼顾保留丰富的特征信息并精简数据的冗余度的优势；在后续预测滑坡概率时，根据不同栅格尺寸对应的栅格内的滑坡影响因子，以天为单位将滑坡影响因子输入至滑坡预测模型中，预测滑坡发生概率，提高了滑坡预测的准确度。

可选地，栅格尺寸确定模块420还用于，基于所述第一栅格尺寸下的滑坡影响因子的分布特征，对所述第二栅格尺寸下的滑坡影响因子的分布特征进行显著性检验；

如果所述第二栅格尺寸下的滑坡影像因子通过显著性检验，确定所述第二栅格尺寸为所述期望栅格尺寸，否则，所述第二栅格尺寸不是所述期望栅格尺寸。

可选地，栅格尺寸确定模块420还用于，基于所述第一栅格尺寸下的滑坡影响因子的分布特征计算第一特征数据，基于所述第二栅格尺寸下的滑坡影响因子的分布特征计算第二特征数据；其中，所述第一特征数据和所述第二特征数据包括方差、均值中的至少一个；

将所述第一特征数据与所述第二特征数据之间的差值与设定阈值比较；

如果所述差值小于所述设定阈值，确定所述第二栅格尺寸为所述期望栅格尺寸，否则，所述第二栅格尺寸不是所述期望栅格尺寸。

可选地，第二栅格尺寸调整模块430还用于，按照所述尺度因子迭代缩小所述第二栅格尺寸，得到所述当前迭代次数下的当前栅格尺寸。

可选地，滑坡预测模块440还用于，基于第一预设模型，根据当前栅格尺寸对应的各栅格的每天的滑坡影响因子，确定每天滑坡的发生概率；

如果当前天滑坡的发生概率大于预设概率阈值，则基于第二预设模型，根据当前栅格尺寸对应的各栅格的每天的滑坡影响因子，确定各个栅格滑坡的发生概率。

可选地，第二栅格尺寸调整模块430还用于，计算所述当前栅格尺寸对应的栅格的动态因子的构造特征，其中，所述构造特征包括各个栅格的每个动态因子的总和、平均值、最大值、最小值、极差、四分位数和雨量时间，并根据每天的所述构造特征和静态因子得到所述第一预设模型的第一特征矩阵；

将所述第一特征矩阵输入所述第一预设模型中，以根据所述第一预设模型确定每天滑坡的发生概率。

可选地，第二栅格尺寸调整模块430还用于，对所述当前栅格尺寸对应的各栅格的每天的滑坡影响因子进行特征提取，以生成所述第二预设模型的第二特征矩阵；

对所述第二特征矩阵进行归一化处理；

将归一化处理后的所述第二特征矩阵输入所述第二预设模型中，以根据所述第二预设模型确定当前天的各个栅格滑坡的发生概率。

本发明实施例所提供的滑坡的预测装置可执行本发明任意实施例所提供的滑坡的预测方法，具备执行方法相应的功能模块和有益效果。

实施例五

图5为本发明实施例五提供的一种滑坡的预测设备的结构示意图，如图5所示，该设备包括处理器510、存储器520、输入装置530和输出装置540；设备处理器510的数量可以是一个或多个，图5中以一个处理器510为例；设备中的处理器510、存储器520、输入装置530和输出装置540可以通过总线或其他方式连接，图5中以通过总线连接为例。

存储器520作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块，如本发明实施例中的滑坡的预测方法对应的程序指令/模块(例如，滑坡的预测装置中的数据获取模块410、第二栅格尺寸确定模块420、第二栅格尺寸调整模块430和滑坡预测模块440)。处理器510通过运行存储在存储器520中的软件程序、指令以及模块，从而执行设备的各种功能应用以及数据处理，即实现上述的滑坡的预测方法。

存储器520可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据终端的使用所创建的数据等。此外，存储器520可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储器520可进一步包括相对于处理器510远程设置的存储器，这些远程存储器可以通过网络连接至设备/终端/服务器。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。

输入装置530可用于接收输入的数字或字符信息，以及产生与设备的用户设置以及功能控制有关的键信号输入。输出装置540可包括显示屏等显示设备。

实施例六

发明实施例六还提供一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行一种滑坡的预测方法，该方法包括：

获取目标区域的第一栅格尺寸和第二栅格尺寸，分别在所述第一栅格尺寸下和所述第二栅格尺寸下获取所述目标区域的滑坡影响因子，其中，所述第二栅格尺寸大于所述第一栅格尺寸；

分别基于所述第一栅格尺寸下的滑坡影响因子和第二栅格尺寸下的滑坡影响因子的分布特征，确定所述第二栅格尺寸是否为期望栅格尺寸，其中，所述期望栅格尺寸小于所述第二栅格尺寸且大于或等于所述第一栅格尺寸；

如果否，基于特定尺度因子，迭代调整所述第二栅格尺寸；

如果当前迭代次数对应的当前栅格尺寸为所述期望栅格尺寸，将所述期望栅格尺寸对应的各栅格的每天的滑坡影响因子输入至滑坡预测模型，基于所述滑坡预测模型确定各栅格滑坡的发生概率。

当然，本发明实施例所提供的一种包含计算机可执行指令的存储介质，其计算机可执行指令不限于如上所述的方法操作，还可以执行本发明任意实施例所提供的滑坡的预测方法中的相关操作。

通过以上关于实施方式的描述，所属领域的技术人员可以清楚地了解到，本发明可借助软件及必需的通用硬件来实现，当然也可以通过硬件实现，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如计算机的软盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、闪存(FLASH)、硬盘或光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。

值得注意的是，上述用滑坡的预测装置的实施例中，所包括的各个单元和模块只是按照功能逻辑进行划分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能单元的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。

注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。