王雅慧，陈尔学，郭 颖＊，李增元，金玉栋，赵俊鹏，周 瑶

(1. 中国林业科学研究院资源信息研究所，北京 100091；2. 喀喇沁旗旺业甸实验林场，内蒙古 赤峰 024423；3. 赤峰市红山区棚户区改造办公室，内蒙古 赤峰 024000)

深度学习模型是近年来新兴的一类机器学习算法，已被有效的应用于许多领域当中。其中，在遥感领域，该类方法的应用效果也得到了有效的证明。对于遥感地物分类的应用问题，研究表明全卷积类深度神经网络方法（Fully Convolutional Networks，FCN）可有效的解决该问题。方旭等[1]证明通过将均值漂移分割算法与FCN-8s方法相结合能够有效提高高分辨率遥感影像的分类精度，与传统机器学习分类方法SVM相比提高了10%以上的精度。邓国徽等[2]将FCN-8s应用于高分遥感数据的语义分割，以分割图像中的施工场地为例，结果表明：提出的改进全卷积网络模型能够较为显著的提高了施工场地的识别率。Piramanayagam[3]利用FCN将高分辨率遥感影像的光谱信息与DSM的几何信息有效融合，从而提高了高分辨遥感图像的语义分割能力。条件随机场方法（Conditional Random Field，CRF）作为图像处理中常用的能量优化算法，也已被学者证明可有效的改善FCN网络的分类效果[4-5]。

U-net网络是全卷积神经网络模型中常用的一种，该模型通过一个收缩网络以及一个扩张网络，构成了一个U型结构，与经典的FCN网络相比，该网络在上采样过程中可以结合更多层的不同的特征，能够最大程度地保留下采样过程中的一些重要的特征信息[6]。研究表明，该类网络可有效的进行遥感图像信息提取。陈睿敏等[7]使用UNet深度学习模型对红外遥感影像进行5类地物信息提取，其分类效果与传统机器学习方法相比提高了6%。Li等[8]提出了一种卷积块代替卷积层的深度卷积神经网络DeepUNet，提高了海岸线的分割精度。

有研究表明，通过少量的调整相关遥感特征，可以优化FCN深度学习模型的效果[2，9]。邓国徽等[2]提出了融合蓝板房光谱识别指数和裸土光谱识别指数的改进全卷积网络模型，结果表明：添加的指数可以较为显著的提高了施工场地的识别率；Cao Q等[9]将NDVI与CNN的高级特征进行了融合，该算法对土地类型分类效果良好。

本研究借鉴前人对FCN-8s模型的优化思路，基于U-net网络模型，在模型训练过程中通过在原始波段的基础上加入NDVI特征，并增加CRF后处理过程，以期改善U-net深度学习模型用于森林类型分类的效果。

1 实验区概况和数据

1.1 实验区概况

旺业甸林场位于赤峰市喀喇沁旗西南部，地理坐标为 118°09′～118°30′ E，41°35′～41°50′ N，属于生态主导型多功能林业区。山体相对高度一般为200～ 400 m， 坡 度 15°～ 35°。 林 场 林 地 面 积2.57万hm2，其中有林面积是2.33万 hm2。人工林以油松（Pinus tabulaeformis Carr.）、落叶松（Larix principis-rupprechtii May）为主，天然林的树种以白桦（Betula platyphylla Suk），山杨（Populus davidiana Dode）等为主。

1.2 遥感数据

采用的遥感数据为1景GF-2 PMS多光谱遥感影像（空间分辨率为3.2 m），完全覆盖实验区。成像时间为2015年4月26日，这时地表土壤尚未解冻，植被尚未开始萌芽，油松、落叶松和白桦等森林类型在遥感影像上的区分度较高。

对GF-2多光谱影像进行辐射定标、大气校正、正射校正、影像配准、地形辐射校正等预处理。首先，对GF-2多光谱影像进行辐射定标和FLAASH大气校正，得到地表反射率影像[10]。其次，采用GF-2多光谱影像自带的RPB文件对地表反射率影像进行基于RPC模型的正射校正，所采用的数字高程模型（DEM）为5 m空间分辨率的ZY-3 DEM；然后，以2 m空间分辨率的ZY-3数字正射影像产品（DOM）为参考，采用自动匹配算法对上一步产生的“正射校正”结果影像进行影像到影像的配准处理，并手工采集30个控制点对校正精度进行检验，东西向误差为4.05 m，南北向误差为1.35 m；进一步采用VECA地形辐射校正模型[11]和ZY-3 DEM数据进行了地形辐射校正处理。

1.3 小班数据

收集到研究区2007年的矢量小班数据，其数据库包括林种、优势树种、小班面积、地类等属性字段，其中优势树种类型包括油松、落叶松、白桦等，基于该属性字段制作了小班优势树种类型分布图（图1），这里的小班优势树种类型和文本的森林类型具有相同含义，可以作为获取森林类型遥感分类模型训练和精度检验样本的参考数据之一。

图 1 旺业甸林场优势树种类型分布图Fig. 1 Distribution map of dominant tree species in Wangyedian Forest Farm

1.4 外业调查数据

外业调查于2017年9月开展，通过实地调查，共采样旺业甸林场内311个样点，包含的主要森林类型有落叶松、油松、白桦、杨树、榆树等，其中分布最为广泛的为落叶松和油松。图2给出了所调查样点的空间分布，样点基本覆盖整个林场，其中油松样点74个，落叶松69个，白桦46个，其他林地46个，耕地36个以及建设用地40个。外业调查过程中，对所调查样点在2015年的实际地类进行了分析判断，确保不会因调查时间和遥感影像成像时间不一致导致地物类型判断出现偏差。

图 2 主要地类的外业调查样点分布图Fig. 2 Distribution map of field sample points of the major land types

1.5 样本数据

深度学习网络模型的训练需要输入块（Patch）状的样本，本研究参考小班优势树种分布图（图1）、外业调查数据，通过目视解译制作样本块。共制作了样本块142块，其中114块作为训练样本，28块作为验证样本。从图3可以看出，训练样本块在空间上的分布较均匀，而且作者还特地调整了各地类样本像元的数量，尽量使各类别的样本像元数相差不大，但由于研究区建设用地太少，其样本数量和其他类别差别仍然较大。图4为各地类样本像元数对比图。所选取的每个样本块的大小并不一致，为了训练深度学习网络模型的需要，采用最临近值法将不同大小的样本块统一重采样为256像元×256像元大小。

图 3 样本块的空间分布图Fig. 3 Spatial distribution of the sample patchs

图 4 各地类像元数量Fig. 4 Total pixel number of each land types

为了能够采用统一的标准比较评价传统机器学习模型和深度学习模型，SVM、RF模型的训练和精度检验也采用上述的同一套样本块。

2 研究方法

2.1 森林类型分类

依据研究区出现的地类及其空间分布特点，在外业调查数据和已有小班数据的支持下，通过对GF-2遥感影像进行预分类，对能够有效区分的潜在地类进行了分析，最终确定了本研究的待分类类别包括：落叶松、油松、白桦、耕地、建设用地和其他。其中前3个地类属于有林地的细分类别，是以林地上林分的优势树种类型名称作为类别名称。其他地类主要包括在该研究区面积占比很小的樟子松、云杉、榆树等。

2.2 分类总技术方法

本文研究从3个部分展开分类方法研究：首先，提取NDVI特征，其次，基于NDVI特征和原始影像4个波段，共计5个特征的数据构建深度U-net模型并进行分类；最后，对分类结果进行CRF后处理。在本研究中，将上述流程简称为U-net-NDVI-CRF方法。2.2.1 U-net网络 U-net网络由一个收缩网络以及一个扩张网络构成，一共包含23个卷积层。其网络的结构具体为：收缩网络主要用于重采样的过程，随着网络深度的增加，提取高维抽象特征，每一次重采样都包含两个3×3的卷积操作，使用修正线性单元（rectified linear unit，ReLU）作为激活函数。每一次重采样的卷积操作将图片的尺寸变为原来图片尺寸大小的1/2，特征通道变为原来的两倍；扩张网络主要对收缩网络中的特征进行过采样的工作，每一次过采样包含两个3×3的卷积操作，同样使用RuLU作为激活函数。在每一次上采样后，图片大小变为原来的2倍，特征数量变为原来的1/2[6]。在上采样操作中，将每一次输出的特征与对应的收缩网络特征合并在一起，补全中间丢失的特征信息。最后，加入1×1的卷积操作将之前所获得的特征映射到所属类别上[6]。

2.2.2 网络训练 U-net模型基于Tensorflow框架训练，使用Adam方法作为优化器来训练网络，学习率为0.001。网络的权重共更新10 000次。网络在一块NVIDIA Quadro K5200显卡上训练，GPU内存为8G。

本研究使用SVM和RF分类器作对比实验[12-13]。针对SVM模型，设置核函数为GRBF（Gaussian radial basis function），使用的目标函数类型是C-支持向量分类（C-SVC），最优惩罚因子C和核函数参数γ进行自动寻优。针对RF模型，实验中设置的决策树的数量为100，每个决策树随机选择的5个特征构建模型，并通过投票法得到分类类别。

2.3 分类精度检验

以外业调查样本地类为各地类真值，对分类的正确性进行检验，采用的评价指标包括分类混淆矩阵，得到制图精度、用户精度、总体分类精度及Kappa系数。

3 结果与分析

3.1 基于U-net-NDVI-CRF方法的分类结果

基于U-net-NDVI-CRF方法的分类结果混淆矩阵如表1所示，总体分类精度达84.89%，Kappa系数为0.82。除落叶松和其他地类外，各地类制图精度均在83%以上。其中，其他地类精度较低，易错分为白桦，主要原因在于GF-2多光谱数据成像时间为冬季，其他地类（主要包括榆树、杨树、柞树等植被类型）均属于落叶性植被，他们的光谱曲线和NDVI特征与白桦差异较小，容易导致误分。落叶松存在部分错分为油松的情况，主要原因在于此季节的落叶松已经开始萌芽，且与油松同属于针叶林、人工林，在光谱和纹理上具有一定的相似性，造成部分落叶松错分为油松。

表 1 U-net-NDVI-CRF方法的分类结果混淆矩阵Table 1 Confusion matrix of U-net-NDVI-CRF classification result

3.2 NDVI对改善分类效果的贡献分析

加入NDVI特征前后U-net模型分类精度对比如表2所示。从表中可以看出，加入NDVI特征之后，分类精度得到有效的提高，总体精度提高了2.25%，Kappa系数提高了0.02，其他地类、耕地、建设用地的分类精度有显著提高，其中，其他地类提高了15.22%，耕地提高了13.89%。从分类细节图（图5）中可以看出，加入NDVI特征可有效改善林地分布区域的分类效果。

3.3 CRF对U-net分类结果的优化效果分析

加入CRF前后U-net分类结果的精度对比如表2所示。在加入CRF后，总体精度较之前提高了0.65%，Kappa系数提高了0.01，其中，落叶松和其他地类的精度有所提高。从定量上的分析可见，分类精度提升幅度并不大，但从目视角度分析，各类地物之间的边缘得到细化。从分类细节图（图5）可以看出，在其他和白桦混交的区域，采用CRF后处理分类效果提升的贡献明显。

表 2 加入NDVI特征以及CRF前后U-net分类结果的精度对比Table 2 Accuracy comparison of U-net classification results before and after adding NDVI features and CRF

图 5 部分油松、落叶松分类结果细节图Fig. 5 The detailed classification results

3.4 与其他3种相关分类方法的对比分析

为了验证本研究所设计的U-net-NDVI-CRF分类方法的有效性，将其与加入NDVI特征的传统机器学习SVM、RF的分类结果，加入NDVI特征和CRF的FCN-8s的分类结果进行了比较分析。这4种分类方法的分类结果如图6所示，分类精度对比如表3所示。深度学习U-net、FCN-8s模型的总体精度和Kappa系数均高于SVM、RF。U-net的总体精度和Kappa系数明显高于其他模型，对油松、白桦、耕地、建设用地和其他的分类精度均较高，尤其是其他地类的分类表现更为突出。对于落叶松，U-net的分类精度较低，主要原因为落叶松与油松同属于人工林，在空间特征上存在一定的相似性，使得油松和落叶松存在一定的混淆。对建设用地，U-net和FCN-8s比SVM或RF改进很多，主要是因为此季节的建设用地和耕地在光谱上存在一定的相似性，且耕地中的大棚与建设用地十分相似，U-net和FCN-8s可以自动提取特征以区分两个地类，而SVM和RF仅靠光谱特征不能完全将建设用地与耕地等地类分出。从分类细节图（图7）可以看出，U-net-NDVI-CRF对油松和落叶松的分类效果要好明显于其他分类方法。

4 讨论

（1）本研究采用U-net网络模型结合CRF对5个特征（原始四波段以及NDVI）的GF-2 PMS 3.2 m的遥感数据进行森林类型分类，总体分类精度为84.89%，与传统的机器学习分类方法RF、SVM在同等五波段时的分类结果相比，提高了8.04%至12.54%，提升效果显著，该结论与已有的将深度学习网络模型用于土地覆盖/利用遥感分类的文献所得到的结论一致[2,5,9]。

（2）虽然深度学习方法，可以自动提取遥感特征，但是现有文献[2]、[9]表明，可以通过加入NDVI等少量特征可以提高分类精度。本研究为与机器学习分类方法在数据上形成对等性，不进行过多的手动提取和选择特征的步骤，选取NDVI与影像原始4个波段组成5个特征的数据进行分类，机器学习SVM、RF分类精度低于U-net的分类精度。

图 6 分类结果 (a) RF；(b) SVM；(c) FCN-8s；(d) U-netFig. 6 Classification result (a) RF; (b) SVM; (c) FCN-8s; (d) U-net

表 3 4种相关分类方法的分类精度对比Table 3 Comparisons of classification accuracy of four relevant methods

（3）由于采用的遥感影像成像时间和光谱限制，其他地类（主要包括榆树、杨树、柞树等植被类型）均属于落叶性植被，他们的光谱曲线和NDVI特征与白桦差异较小，造成目视解译困难，对分类结果造成了一定的影响。因此，考虑使用生长季影像对其他地类进行细化之后进行深度学习精细分类。

（4）本研究实验采用的数据是GF-2 PMS 3.2 m分辨率的遥感数据，而全色数据的分辨率可达0.8 m，结合纹理信息丰富的全色数据进行实验，对比在更高分辨率下地物的边缘信息是否更好的保留为下一步的研究方向。且本研究的GF-2 PMS遥感数据的时相为4月，使得其他地类的光谱差异不明显，使用生长季的数据进行森林类型的精细分类是下一步的研究方向。

图 7 局部区域油松、落叶松4种分类方法的分类结果对比：(a) GF-2 PMS；(b) 优势树种分布图；(c) RF；(d) SVM；(e) FCN-8s；(f) U-net；(g) 图例Fig. 7 Comparison of classification result of the four classification method for some region: (a) GF-2 PMS;(b) Distribution map of dominant tree species; (c) RF; (d) SVM; (e) FCN-8s; (f) U-net; (g)Legend

5 结论及展望

针对国家森林资源规划设计调查业务对高分辨率多光谱遥感分类技术的应用需求，发展了一种深度学习森林类型遥感分类方法。研究结果表明，通过加入NDVI特征可有效提高U-net网络森林类型分类结果的精度，采用CRF后处理则可以有效改善分类效果，且与FCN-8s、SVM、RF相比，优化后的U-net分类效果有较大的提升。

与传统机器学习分类方法相比，本研究所发展的方法可以自动提取遥感特征，从而可有效避免传统方法由于手工提取和选择特征等而导致的分类结果容易受人为因素影响的问题，具有更高的可操作性。当然，目前深度学习方法还存在对高性能GPU计算机依赖性强，模型训练耗时，技术门槛较高等问题。随着我国高性能GPU计算机的快速普及和计算能力的快速提高，这些制约深度学习分类方法林业业务化应用的障碍将很快得到克服，在本研究研究的基础上，针对林业遥感业务特点，通过深度学习网络结构的创新设计，进一步提高分类精度和效率，具有重要的研究和应用价值。