内蒙古大兴安岭林区，在内蒙古自治区的东北方向（119°36′26″～125°24′10″ E，47°03′26″～53°20′00″ N）^[8]，海拔425～1 760 m^[9]，是我国4大国有林区之一。林区地处高纬地带，南北跨7°，地域辽阔，是中国最大的集中连片的国有林区^[10]。内蒙古大兴安岭地区位于寒温带大陆性季风气候带，全年平均气温在−4 ℃至−2 ℃之间，平均气温小于10 ℃的时间长达9个月。全年降水量在300 mm以上，相对湿度为70%～75%，林区总面积约1.07 × 10⁷ hm²，森林覆盖率达79.82%。2009—2018年期间，研究区发生过多次森林火灾，火灾分布情况如图1。

图  1  研究区及森林火灾分布（2009—2017）

Figure 1.  Research area and forest fire distribution(2009—2017)

气温直接影响森林可燃物的含水量，高温且降水量少、强风且持续干旱的气候条件下森林火灾高发^[11-12]。地形条件对森林火灾的作用是间接的。地形条件影响局部小气候进而影响到森林可燃物含水率变化，并对森林火灾的发生产生影响^[13-14]。区域火灾历史数据可以反映区域的历史森林火灾状况，发生过森林火灾的地方有再发生森林火灾的潜力，因此需要关注。植被指数能够反映森林覆盖，是对地表植被状况的简单、有效和经验的度量，植被指数越高的地方，发生森林火灾的概率也相应较大。人为因素对森林火灾发生具有“增殖”效应^[15] ，通常考虑与道路的距离、与河流的距离等。

火点数据为MCD64A1月火点数据产品，获取自基于云的地理空间处理平台谷歌地球引擎（Google Earth Engine, https://earthengine.google.com/），其空间分辨率为500 m，时间分辨率为8 d，时间范围为2009—2018年，本研究同时选取低、中、高3个置信水平等级的火点数据。地形数据为ASTER GDEM数据，空间分辨率为30 m，下载自地理空间数据云（https://www.gscloud.cn/），在地理数据处理软件中进一步计算得到坡度、坡向、平面曲率数据。气象数据（平均气温、总降水量、平均风速）源自国家地球系统科学数据中心，空间分辨率为1 km，时间分辨率为月。NDVI来源于中国科学院资源环境科学数据中心（https://www.resdc.cn/DataSearch.aspx）。人为因素矢量数据（与建筑物的距离、与道路的距离、与水域的距离）获取自开源网站OpenStreetMap（https://www.openstreetmap.org/），时间分辨率为5 a。

在输入模型前，对每个森林火灾影响因子图层进行最大最小归一化处理，将图层像素点取值控制在[0,1]范围。在ArcGIS 10.8中，将论文中涉及的森林火点及火灾影响因素图层像素大小统一为1 km × 1 km，坐标系统一为WGS_1984_UTM_Zone_51N。最后，使用波段合成工具将1 a内所有的影响因子图层合并为1个栅格图层，最终建立11个多变量栅格数据集。为解决原有火点及非火点数据不平衡的问题，本研究通过创建随机点工具生成非火灾点，与已知火灾点合并生成5 560个训练样本（图2）。设定2009—2017年所有火点的森林火灾发生概率为1，重复出现的火点只记1次，非火点的森林火灾发生概率为0，作为后续模型的标签数据。

图  2  样本点分布情况

Figure 2.  The distribution of sample points

本研究通过计算容差值和VIF值来进行特征选择，同时计算Pearson相关系数来度量变量之间的线性关系强度，以及计算Spearman秩相关系数来度量变量之间的非线性关系(p ≤ 0.05)^[16]。

支持向量机（Support Vector Machines, SVM）有很强的数学基础和严谨的理论支持，是常用的监督学习模型，可运用于数据分类、回归分析和模式识别模型。支持向量机算法的最优分类面求解问题，本质是求样本分类间隔最大的二次函数的目标解。使用网格搜索法，最终确定γ为1，正则化参数设置为80。

随机森林（Random Forest, RF）是一种基于Bagging算法和决策树方法的集成学习方法。RF利用重抽样技术对训练样本进行随机采样，训练出多个不同的决策树，将多个决策树的分类结果进行组合，投票或取平均得出最终的预测结果。同样使用网格搜索法，最终确定生成树个数为450，树的最大深度为18。

卷积神经网络（Convolutional Neural Networks, CNN）是一种前馈神经网络，参数采用经典的随机梯度下降法（Stochastic Gradient Descent，SGD）不断调整。本研究所用的卷积神经网络（CNN）模型的框架是以Alexnet为参考，使用TensorFlow作为后端。为解决特征过多而产生的过拟合问题，加入Dropout算法和增加Batch Normalization层。

每个输入patch都是1个大小为15 × 15 × 11的三维数据（像元大小15，波段数为11）。本研究的主要架构包括3个卷积层、3个最大池化层和3个全连接层。卷积层的核数均为64，核大小均为3 × 3，步长为1，填充方式选择“same”。每个卷积层后都有一个ReLU激活函数、一个BN层和一个最大池化层。最大池化层的核大小为2 × 2，步长为2，3个全连接层均有32个神经元，在全连接层中调整权重衰减，并在全连接层后添加1个双向分类器的激活函数，即Softmax，可以计算出该点属于火点/非火点的概率值。参数总量为84706个，参数量设置较为合理。

本研究基于前人经验选择超参数取值并根据训练过程中的准确率和损失值图（图3），确定超参数的最优取值为：学习率为0.001，批处理大小为32，Dropout值为0.5，优化方法选择RMSprop，卷积核数量为64。

图  3  训练中的准确率和损失值变化曲线

Figure 3.  Accuracy and loss curve during training

本研究采用4个统计指标，包括总体准确率、特异性、召回率、精确率^[17]。这4个统计指标按以下方式计算:

式中TP表示实际为火点，预测结果也为火点，FP表示实际为非火点，预测结果为火点，TN表示实际为非火点，预测结果也为非火点，FN表示实际为火点，预测结果为非火点。在分类模型中，预测结果一般以概率值输出，通常设置阈值来将概率转化为具体类别。因此，阈值的设定与模型的准确率十分相关。受试者工作特征（Receiver Operating Characteristic, ROC）曲线和曲线下面积（Area Under ROC Curve, AUC）用于判断预测结果的准确性，并表达了模型的泛化能力。

基于林火影响因子相关性热力图（图4）结果表明，坡度与粗糙度呈现明显的正相关性，而平面曲率和剖面曲率则呈现较明显的负相关性。为了防止模型过拟合，将影响火灾发生的因子中的粗糙度和剖面曲率剔除，而保留坡度和平面曲率用于后续研究。

图  4  林火影响因子相关性热力图

Figure 4.  Correlation heat map of forest fire influencing factors

随机森林计算后，得到11个影响林火影响因子的重要性值。气温、海拔、与水域的距离、总降水量是影响林火发生最为重要的4个因素。此外，与建筑物的距离、月均风速和NDVI等因素也对森林火灾的发生产生了较为显著的影响，其重要性值均大于0.05。这些变量均与森林火灾有着明显的关系，因此可以用于后续的建模研究。

训练过程包括训练和验证2个阶段。训练样本和验证样本分别占训练集的80%、20%。在每个epoch结束后，使用验证样本数据集对超参数进行微调。结果表明（图5），每次训练后，训练样本数据集和验证样本数据集的准确性均达到95%以上，验证集的损失值低于0.05，而验证集的损失值在0.15左右。最终，选择验证集的损失值最小时（epoch=300）对应的模型作为最终的分类模型。

图  5  训练过程中准确率和损失值变化

Figure 5.  The change of accuracy and loss value during graph training

本研究将CNN与RF及SVM这两种常见的机器学习方法进行比较，并通过绘制评价指标雷达图（图6）对这3种方法进行比较。在预测样本数据集中，CNN精确率和召回率远高于RF及SVM，达到90%以上，CNN特异性和总体准确率略低于RF及SVM，但可明显观察到CNN总体性能优于RF及SVM，故CNN预测结果最为可靠。

图  6  3种模型评价指标雷达

Figure 6.  Three kinds of model evaluation index radar map

AUC值是评价模型分类性能的一个重要指标，用来评估模型的平均性能。其中，CNN的AUC值为0.838，表示CNN对2018年大兴安岭林区森林火灾预测的整体拟合度达到了83.8%，高于RF（0.788）和SVM（0.794）（图7）。

图  7  3种模型的ROC曲线

Figure 7.  ROC curves of the three models

经过训练，最终得到分类模型，使用测试数据集（2018年的多变量栅格数据集）来评估最终模型性能。因为滑动窗口在预处理阶段密集重叠并覆盖了整个栅格数据集，故模型可以对所需预测图中的每个像素得出火灾和非火灾类别的概率，选择火灾类别的概率作为最终的预测概率值，可视化生成森林火灾敏感性图（图8）。3种模型均识别出在大兴安岭东南方向森林火灾易感性的值较大，中部和西部位置火灾发生概率较小，这表明大兴安岭东南地区更易发生森林火灾，中部和西部地区较不易发生森林火灾。

图  8  RF、SVM、CNN模型预测2018年森林火灾敏感性

Figure 8.  RF, SVM and CNN model prediction of forest fire sensitivity in 2018

通过分析图9可以发现，CNN模型中极高和极低的森林火灾类别占总研究区森林面积的90%以上，在3种模型中所占比例最高，高、中、低3类分别占总面积的2.12%、1.89%和2.74%。因此， CNN模型能够有效划分出森林火灾易感性极高及极低的区域，有利于在实践工作中加强对森林火灾易感区的火灾预防。相比之下，RF和SVM的预测图中存在较大范围的中高易感区，对林火高易感性区域没有明确的判定。因此，CNN模型的结果更适合该地区森林火灾的预测。

图  9  3种模型下林火敏感性等级占比情况

Figure 9.  The proportion of forest fire sensitivity grades under the three models

针对变量选择，本研究将粗糙度和剖面曲率这两个因子剔除，保留剩余的11个特征用于后续建模。气温、海拔、与水域的距离、总降水量是影响林火发生较为重要的4个因素，这与Najafabadi、Pourtaghi、马文苑、Nami、Abdi等的研究结果基本一致^[18-22]。随着海拔的增加，林火发生概率明显下降，这是因为海拔的增加使得植被和土壤的湿度增加^[22]，且受人为干扰较少^[23]，不利于森林火灾的发生。气象条件大多是通过影响可燃物的含水量和温度对林火的发生时间和区域产生影响^[24]，气温升高会使林区内可燃物的含水量减少，也会使可燃物本身的温度升高，减少了外部因素达到其燃点所需的热量，可燃物容易被点燃^[25-27]；降水可以反映出可燃物湿度的变化情况，降水量较多时林区内可燃物含水量处于饱和状态，不易被点燃，发生森林火灾的可能性和严重程度都会较低^[26-27]。人为因素也会对森林火灾的发生有一定影响，距离人为活动的距离越远，森林火灾发生的概率也就越小^[28-29]。而坡向对森林火灾发生的解释度较小，这与li等^[30]的研究结果一致。

本研究证实了CNN在森林火灾预测方面具有比RF和SVM更为突出的优势，能够有效地划分出森林火灾易感性极高和极低的区域，因此CNN模型的结果更加适合森林火灾预测。CNN模型AUC值达到0.8，其精度介于Dimuccio等^[1]、Zhang等^[3]研究的AUC值之间，基本满足预测精度。在Bisquert^[31]关于西班牙西北部加利西亚森林火灾概率的预测研究中也证实，人工神经网络的准确性和精度都高于逻辑回归。Bergado^[32]运用29个量化特征，采用深度学习方法生成2006—2017年澳大利亚维多利亚州未来7天内野火燃烧概率的每日地图，该地图更平滑、更正则化，在各项定量和定性指标中都显示出很大的优越性。Muhammad^[33]探索了CNN并设计了一种微调架构，用于在有效灾害管理系统的监视期间进行早期火灾检测，研究证实该框架有高准确性和低火灾误报率。总的来看，与机器学习SVM、RF相比，CNN具有以下优势：首先，CNN可以考虑相邻空间区域的相关性，因此适用于空间和地理相关的研究。其次，深度学习可以揭示更深层的特征。通过多次卷积和池化操作，CNN可以提取更高级和抽象的特征，这些特征对于森林火灾的发生起着决定性的作用。最后，CNN通过权值共享减少了训练过程中的超参数数量，从而降低了CNN的结构复杂度。

本研究选取地形、气候、植被和人为方面共11个因素形成多变量栅格数据集，构建出3个林火敏感性模型，并将大兴安岭森林火灾的风险进行可视化展示。对大兴安岭森林火灾进行预测，有助于相关部门重点盯防，做到早预防、早发现，对大兴安岭森林火灾预防工作有重要意义。

（1）通过特征选择方法，最终确定11个森林火灾潜在影响因子，分别是地形、平面曲率、坡度、坡向、平均气温、总降水量、平均风速、与道路的距离、与水域的距离、与建筑的距离、NDVI等。其中，气温、海拔、与水域的距离、总降水量是影响林火发生最为重要的4个因素。

（2）CNN模型的AUC值高于RF、SVM，说明CNN对林火的预测更准确，模型可靠性更高。

（3）空间上，大兴安岭森林地区东南地区森林火险等级较高，较易发生火灾。