[发明专利]基于概率编程技术的山火风险预测方法

摘要

本发明公开了一种基于概率编程技术的山火风险预测方法。该基于概率编程技术的山火风险预测方法以地理遥感、气象数据作为输入，实现对山火发生风险的预测，而地理遥感、气象数据存在多源、异构、时空分辨率多样等的特征，利用深度学习、概率编程技术对山火发生风险进行预测，解决了山火标记数据少、山火发生机理复杂、信息不完备情况下的预测建模，该预测方法能够对预测的不确定性进行分析，更好的满足山火发生风险预测这种机理复杂、数据不完备的应用需要且能够快速完成山火发生风险的预测预警，并且该预测方法预测的整体效果较好；通过实验验证AUC(Area Under Curve，ROC曲线下的面积)指标达到了0.8。适合在山火短期发生风险预测技术领域推广应用。

主权项

1.基于概率编程技术的山火风险预测方法，其特征在于包括以下步骤：1)、建立山火风险预测模型，所述山火风险预测模型的建立过程如下所述：11)、采集一段历史时间内的地理遥感数据、气象数据、历史山火点数据；所述地理遥感数据包括植被指数NVDI、地表覆盖LC、地形地貌DEM、水汽数据；所述气象数据包括降雨、温度、湿度、风速、风向数据；12)、采用时空数据融合技术分别对11)步骤采集的地理遥感数据、气象数据进行融合处理得到在时间、空间上分辨率均一样的融合数据；13)、将经过12)步骤处理得到的地理遥感数据、气象数据按照数据类型分为数值型数据和字符型数据，所述数值型数据包括植被指数NVDI、降雨、温度、湿度、风速、水汽数据；所述字符型数据地表覆盖LC、地形地貌DEM、风向数据；14)、对数值型数据采用z‑score标准化方式进行数据处理；对字符型数据采用one‑hot编码方式进行数据处理；15)、将历史山火点数据进行山火点风险标记，所述标记规则为：未来n天中发生山火则该点的目标变量标记为1，否则为0；目标变量标记为1的历史山火点数据集为山火样本集，目标变量标记为0的历史山火点数据集为非山火样本集；16)、将15)步骤得到的山火样本集和非山火样本集采用如下算法进行处理：A、输入：空间半径radius、训练集d、目标变量y、多数类抽样数量n、采样结果Sample_result；B、对目标y中每个取值进行统计计数，记为为y_i的样本数；C、确定多数类少数类对应的多数类和少数类的数据集分别记为和d_min，将d_min追加到Sample_result中；D、依次对多数类中每一个类型进行欠抽样，算法流程为：D‑1)、获取和y_min所对应样本的全部空间点集合分别记为p_ma、p_min，并将样本点的WGS84坐标转换为UTM平面投影坐标系；D‑2)、计算p_ma与p_min点集合中两两之间的欧式距离，距离计算公式为：D‑3)、若d_{p1_p2}＞radius则p₁保留到p_ma中，否则将p₁从p_ma剔除，最终得到sp_ma；D‑4)、从中选择空间点sp_ma对应的样本数据D‑5)、在中采用有放回随机抽样方式，抽取n个样本追加到Sample_result中；17)、构建神经网络模型，所述神经网络模型为贝叶斯神经网络结构，所述贝叶斯神经网络结构包括输入层、隐藏层、输出层，所述输入层的节点数与转换处理后的特征数一致；所述隐藏层对数据进行深度的处理和学习，隐藏层的参数包括隐藏层数n_levels∈[3,…,10]，第i层的神经元数量为n_cells_i∈random(n_min,n_max)，激活函数为act_fun∈[tanh,ReLU,…]；所述输出层的输出节点数为2，输出层的激活函数采用sigmoid函数；18)、采用概率建模对贝叶斯神经网络的神经元权重参数进行设置，所有神经元的初始权重参数采用高斯分布，即：w_ij～N(μ,σ²)，其中μ＝0,σ²＝1；将14)步骤和16)步骤得到数据通过贝叶斯神经网络结构的输入层输入并以KL散度最小化为目标，采用变分推理(ADVI)算法实现神经元权重参数的推理及学习，得到神经元的最终权重参数分布；2)、利用山火风险预测模型进行未来山火风险预测，具体预测过程如下所述：21)、针对山火风险预测任务的需要，采集从当前时间开始往前一段历史时间内的地理遥感数据、气象数据、历史山火点数据；所述地理遥感数据包括植被指数NVDI、地表覆盖LC、地形地貌DEM、水汽数据；所述气象数据包括降雨、温度、湿度、风速、风向数据；22)、采用时空数据融合技术分别对21)步骤采集的地理遥感数据、气象数据进行融合处理得到在时间、空间上分辨率均一样的融合数据；23)、将经过22)步骤处理得到的地理遥感数据、气象数据按照数据类型分为数值型数据和字符型数据，所述数值型数据包括植被指数NVDI、降雨、温度、湿度、风速、水汽数据；所述字符型数据地表覆盖LC、地形地貌DEM、风向数据；24)、对数值型数据采用z‑score标准化方式进行数据处理；对字符型数据采用one‑hot编码方式进行数据处理；25)、对18)步骤得到的神经元的最终权重参数分布进行抽样，抽样得到一个神经元的权重参数，将该神经元权重参数代入山火风险预测模型中，并将24)步骤得到的数据输入确定了神经元权重参数的山火风险预测模型中得到山火发生风险值；26)、重复步骤25)N次，得到N个山火发生风险值，进而获得山火发生概率分布的风险期望和方差。