CN111191628A - 基于决策树与特征优化的遥感影像震害建筑物识别方法 - Google Patents

Abstract

本发明针对在震前参考信息缺乏情况下，仅依赖震后遥感影像对震害建筑物进行特征建模和随机森林分类中所存在的困难和局限，提出了一种基于决策树与特征优化的遥感影像震害建筑物识别方法。首先，结合影像分割及非建筑物筛选规则提取潜在建筑物对象集合；在此基础上，提出了一种基于分类正确率曲线波动判别的自适应决策树数量提取策略；同时，在特征重要度指标指引下对光谱、纹理、几何形态学这三大类特征进行筛选，获得具有代表性的震害特征集合；最后，基于所构建的优化随机森林模型对震害建筑物进行识别。通过对四组不同的遥感图像进行实验表明，本发明方法在震后复杂场景下的震害建筑物识别中展现出了出色的性能，其总体精度可达到85％以上。

Description

基于决策树与特征优化的遥感影像震害建筑物识别方法

技术领域

本发明公开了一种基于决策树与特征优化的遥感影像震害建筑物识别方法，属于图像识别技术领域。

背景技术

作为一种严重的自然灾害，地震的发生往往伴随着巨大的人员伤亡以及经济财产损失。震后及时、准确地识别震害建筑物，对快速评估灾情、开展应急救援响应及灾后重建等具有重要意义。相较于传统震后人工实地探查的方式，基于遥感影像的震害建筑物识别具有数据获取迅速、覆盖范围广等优势，已经成为震后应急响应中一种重要的技术手段。

伴随着卫星和传感器技术的不断发展，高分辨率遥感影像的广泛应用带来了更加丰富的空间细节信息，从而有利于对震害建筑物的精细刻画。在已有震害建筑物识别方法中，仅依赖于震后影像的震害建筑物识别方法突破了对震前影像的依赖，因而更具应用推广性。尽管如此，此类方法不仅无法提取震前、震后变化信息作为识别的震害建筑物识别的依据，同时还需要面对遥感影像空间分辨率的提高所产生的更加突出的“同物异谱”和“同谱异物”现象，即完好建筑物、震害建筑物及其他地物的类内方差增大而类间方差减小所带来的严峻挑战。为此，能否构建更加高效的特征空间以准确描述建筑物的震害信息，是实现高精度震害建筑物识别的前提和基础。目前震害建筑物识别方法中采用的特征主要包括光谱、纹理以及几何形态学三类。例如，Liu Yu等人利用形态学属性剖面(MAP)和与局部二值模式(LBP)算子提取影像的几何特征与纹理特征，进而利用随机森林分类器提取震害建筑物[刘宇，曹国.基于多特征结合的损毁建筑物检测[J].计算机应用,2015,35(9):2652-2655.]；Asli Sabuncu等人综合利用光谱、compactness以及smoothness特征，提出了一种针对震后场景坍塌建筑物的识别方法[Sabuncu A.A study of earthquake-inducedbuilding detection by object oriented classification approach[C]//Egu GeneralAssembly Conference.2017.]。尽管结合不同种类特征有利于对震害建筑物的多维度描述，但由于不同特征之间存在冗余信息，不仅会增加计算复杂度，同时还可能由于不同特征之间作为震害证据时的相互冲突而降低检测精度。为此，需要多特征的筛选和优化组合策略，从而构建精炼的特征集合。在此基础上，特征集合还必须结合合适的分类方法以获得可靠的震害建筑物识别结果。随机森林(Random Forest，RF)是目前比较流行的集成分类器，具有模型参数少、利用样本和特征选择的双重随机性可以避免过拟合的优点。但合理选择决策树的数量是提高随机森林分类器性能的一个关键因素，决策树数量过少时难以获得可靠的检测结果；数量过多则会降低算法运行效率。特别是当决策树超过一定数量时，分类的精度就在一定范围内上下浮动甚至有下降的趋势。因此要选择合适的决策树棵树。尽管如此，在RF理论中并未给出明确的决策树数量标准，而通常采用的人工赋值方式，不仅容易受到主观因素影响而陷入局部最优，同时又降低了分类过程的自动化程度。

发明内容

发明目的：针对现有技术的不足，本发明的目的在于提出了一种基于决策树与特征优化的遥感影像震害建筑物识别方法，以实现震后复杂场景下的震害建筑物快速、准确地识别，提升总体识别精度。

技术方案：为实现上述发明目的，本发明采用如下技术方案：

一种基于决策树与特征优化的遥感影像震害建筑物识别方法，包括如下步骤：

(1)对卫星遥感影像进行分割提取初始对象集合，并结合几何形态学特征提取潜在的建筑物对象集合；

(2)基于分类正确率曲线波动判别规则自适应确定决策树数量；

(3)在特征重要度指标指引下对光谱、纹理、几何形态学三大类特征进行筛选，获得具有代表性的震害特征集合，进而构建优化的随机森林模型；

(4)根据随机森林判别结果将潜在建筑物对象集合进一步进行分类，识别出震害建筑物。

作为优选，所述步骤(1)中采用高分辨率遥感影像分割算法WJSEG提取遥感影像的初始对象集合。

作为优选，所述步骤(1)中结合几何形态学特征提取潜在的建筑物对象集合的具体规则包括面积规则、矩形度规则和MBI规则；所述面积规则是根据对象所包含的像素个数确定是否剔除，若像素个数不超过设定的阈值，则剔除；所述矩形度规则是根据对象所包含的像素个数与最小外接矩形所包含的像素个数的比值和最小外接矩形的长宽比确定是否剔除，若比值小于设定的阈值且长宽比大于设定的阈值，则剔除；所述MBI规则是基于影像中所有像素的MBI值利用最大类间方差法自适应确定分离阈值，获得对象内属于非建筑物像素的比例，若比例大于设定的阈值则剔除对象。

作为优选，所述步骤(2)中基于分类正确率曲线波动判别规则自适应确定决策树数量，包括：

(2.1)在潜在的建筑物对象集合中，采用人工标记的方式提取样本，构成训练样本集；

(2.2)确定初始建筑物震害特征集合，包括光谱特征、纹理特征和几何形态学特征；

(2.3)构建初始随机森林模型，输入训练样本集，定义分类正确率为每棵树的分类正确的样本数与输入的样本数的比值，比值越大说明分类结果与真实情况越接近，计算每棵树的分类正确率并求所有决策树的正确率均值；

(2.4)以设定的步长，迭代增加决策树后构建新的随机森林模型，获得迭代次数为n时的正确率均值Acc_n；

(2.5)若某Acc_n满足Acc_n-1＜Acc_n＜Acc_n+1，则认为Acc_n是一个子区间内的峰值点，记Acc_opt0＝Acc_n；

(2.6)继续计算Acc_n后连续出现的3个峰值点，分别计为Acc_opt1、Acc_opt2、Acc_opt3，若满足Acc_opt0≥argmax{Acc_opt1,Acc_opt2,Acc_opt3}，则记Acc_opt0对应的决策树颗数为I_opt0；否则，进入下一步；

(2.7)重复步骤(2.5)和(2.5)，直到取得满足

Acc_opt0≥argmax{Acc_opt1,Acc_opt2,Acc_opt3}的Acc_opt0，停止迭代并提取I_opt0；若无法确定Acc_opt0，则记已经获得的Acc_n最大值对应的决策树颗数为I_opt0；

(2.8)继续比较I_opt0与其前、后相邻的N-1棵数的分类正确率，根据分类正确率最大值确定最终提取的优化决策树棵数I_opt，N为步骤(2.4)中设定的步长。

作为优选，所述步骤(3)中在特征筛选的方法是：在包含I_opt颗决策树的随机森林模型中，首先计算所有特征的重要度，再分别求取光谱、纹理和几何形态学三类特征的类间重要度，再分别对三类特征中包含的子特征进行类内重要度排序，最后在重要度指引下获得优化后的震害建筑物分类特征集合；其中I_opt为步骤(2)中优化决策树棵数。

作为优选，任意特征f_t对于第i单棵决策树的重要度由下式计算得到：

式中Φ_B表示OOB样本集合，x_j和l_j分别表示袋外数据中任意样本与其被赋予的类别标签，c_j ⁽ⁱ⁾表示样本x_j所获得的类别标签，

表示将特征f_t的值替换为其他随机值后样本x_j所得的类别标签，N(·)为计数函数；特征f_t对于随机森林分类器的重要度

作为优选，将光谱、纹理和几何形态学三类特征的重要度分别求和得到三类特征的类间重要度，分别记为W_L，W_T和W_G，并分别计算类间冗余归一化惩罚因子；对各个特征子集内的特征按照类内重要度从高到低排列；对各个特征子集，依据冗余归一化惩罚因子的比例取整后，分别仅保留类内重要度相对较高的v_L，v_T和v_G个特征，从而获得共包含V＝v_L+v_T+v_G维特征的优化特征集合F_opt。

作为优选，所述步骤(3)中基于提取的优化决策树棵数I_opt及优化特征集合F_opt，构建的优化随机森林模型如下：

式中，P(x)表示分类结果；p_i(x)表示单棵决策树的分类结果；c表示分类标签，c∈{Undamaged Building,Damage Building,Ruin,The Others}。

有益效果：本发明针对在震前参考信息缺乏情况下，仅依赖震后遥感影像对震害建筑物进行特征建模和随机森林分类中所存在的困难和局限，提出了一种基于决策树与特征优化的震害建筑物识别方法。首先，结合影像分割及非建筑物筛选规则集合提取潜在建筑物对象集合；在此基础上，提出了一种基于分类正确率曲线波动判别的自适应决策树数量提取策略；同时，在特征重要度指标指引下对光谱、纹理、几何形态学这三大类特征进行筛选，获得具有代表性的震害特征集合；最后，基于所构建的优化随机森林模型将建筑物划分为完好建筑物、部分震损建筑物、废墟及其他。通过对四组不同传感器的遥感图像进行实验表明，本发明方法在震后复杂场景下的震害建筑物识别中展现出了出色的性能，其总体精度可达到85％以上。

附图说明

图1为本发明实施例的方法流程图。

图2为本发明实施例中使用的实验数据集的卫星遥感影像图；包括(a)Dataset 1、(b)Dataset 2、(c)Dataset 3和(d)Dataset 4。

图3为对图2中的遥感影像进行分割及初始对象提取结果图。

图4为对图3中提取的对象进行非建筑物对象筛选后的结果图。

图5为本发明实施例中Dataset 1决策树棵数与分类正确率关系图。

图6为本发明实施例与现有方法对Dataset 1的震害建筑物分类结果对比图。

图7为本发明实施例与现有方法对Dataset 2的震害建筑物分类结果对比图。

图8为本发明实施例与现有方法对Dataset 3的震害建筑物分类结果对比图。

图9为本发明实施例与现有方法对Dataset 4的震害建筑物分类结果对比图。

图10为Dataset 1～4震害建筑物分类的参考图。

图11为决策树数量与总体精度的关系曲线图。

具体实施方式

下面结合具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。

本发明实施例公开的一种基于决策树与特征优化的遥感影像震害建筑物识别方法，主要包括四个步骤：潜在建筑物集合提取、随机森林决策树颗数自适应选择、特征重要度指引下的特征集合优化及基于优化随机森林模型的影像分类。具体实现流程如图1所示，下面对各步骤进行详细说明：

(1)潜在建筑物集合提取

(1.1)初始对象集合提取

首先通过影像分割将离散像素划分为具有语义信息的地理对象集合，从而为后续震害建筑物识别提供有效的基本分析单元。为此，本发明采用高分辨率遥感影像分割算法WJSEG，其优势在于：与广泛采用的著名商业软件eCognition相比，WJSEG保持包括建筑物在内的地理对象的轮廓更加完整，且不存在狭长的“虚假单元”，同时有助于提高所提出算法的透明度。

WJSEG的具体实现流程可参见文献[Wang C et al.,“A novel multi-scalesegmentation algorithm for high resolution remote sensing images based onwavelet transform and improved JSEG algorithm,”Optik-International Journalfor Light and Electron Optics.125(19),5588-5595(2014).]，分割参数的设定均采用原文的建议值。经过影像分割后所提取的初始对象集合本发明采用R_origin表示。

(1.2)非建筑物对象筛选

在R_origin基础上，本发明设计了一组基于几何形态学特征的非建筑物对象筛选规则集合，以剔除与建筑物有明显差异的对象，从而在减少计算量的同时避免此类地物在后续处理中产生误检。对于R_origin中的任意对象R_b，具体筛选规则如下：

(1)面积规则。统计R_b所包含的像素个数N_pixels，若N_pixels≤80，则认为R_b为车辆、噪声等弱小目标，予以剔除；

(2)矩形度规则。矩形度是衡量一个对象与其最小外接矩形充满程度的参数，通过计算R_b的最小外接矩形所包含的像素个数N_rectangle，矩形度可表示为Rd＝N_pixels/N_rectangle。令R_b的最小外接矩形的长宽比为Ar。若R_b满足Rd＜0.8且Ar＞5，则认为该对象为道路、河道等狭长目标，予以剔除。

(3)形态学建筑物指数(MBI，Morphological Building Index)规则。MBI利用了属于建筑物的像素在高帽变换后的灰度影像中多具有高亮特征的特点，通过计算多尺度差分序列获得某像素对应的指标值[Huang X,Zhang L.A Multidirectional and MultiscaleMorphological Index for Automatic Building Extraction from MultispectralGeoEye-1 Imagery[J].Photogrammetric Engineering&Remote Sensing,2011,77(7):721-732.]。当MBI值越大时，该像素属于建筑物的可能性越高，计算公式如下：

其中，D和S分别代表线性结构元素的方向和尺度；DMP(d,s)为多尺度差分形态学特征值，依据文献[付乾坤,吴波,汪小钦,et al.基于形态学建筑物指数的城市建筑物提取及其高度估算[J].遥感技术与应用,2015,Issue(1).]的建议，本文设定D＝8，S＝50，d的取值范围为1到8且间隔为1，s的取值范围为1到50且间隔为1。基于影像中所有像素的MBI值，利用最大类间方差法(Otsu)自适应确定分离阈值T_MBI，从而获得R_b内属于非建筑物像素的比例Ratio_false，若满足Ratio_false＞0.8，则剔除该对象。

根据规则(1)～(3)遍历R_origin中所有对象，被剔除的非建筑物对象将不再参与后续的分析和判别。剩余对象将构成潜在建筑物集合R_BLD，用于进一步的震害建筑物识别。

(2)随机森林决策树颗数自适应选择

基于机器学习的高分辨率遥感影像震害建筑物识别本质上是通过对震后影像的特征提取，将目标识别问题转化为影像分类为题。本发明采用的随机森林分类器是以无剪枝决策树为基分类器的一种集成算法，相较于其他机器学习算法具有精度高、泛化能力强、参数少等优势。其中，随机森林中的决策树颗数I是一个显著影响分类精度的关键参数。通常情况下，随着决策树颗数的不断增加，分类精度会随之逐渐上升直至稳定于某一常数。此时，若继续增加决策树颗数，不仅会无谓增加计算量，甚至可能导致分类精度的降低。为此，本发明提出了一种基于分类正确率曲线波动趋势判别的决策树数量自适应确定策略，主要包括以下五个步骤：

Step1:在潜在建筑物对象集合R_BLD中，采用人工标记的方式提取20个完好建筑物、20个震害建筑物、20个废墟和20个非建筑物样本，构成训练样本集H。

Step2：确定初始建筑物震害特征集合F。本发明选择包括光谱、纹理和几何形态学在内的30种震害建筑物识别领域的常用特征构建初始建筑物震害特征集合F[Ximeng C,Zhanfeng S,Tingyan X,et al.Damaged building extraction and rapid assessmentfor earthquake disasters based on high-resolution remote sensing images[J].Journal of Natural Disasters,2016.//李强[1,2],焦其松[2],et al.基于地面激光雷达技术点云的北川县城震害建筑物提取研究[J].科学技术与工程,2016,16(19):244-249.//Qihao C,Yuliang N,Linlin L I,et al.Buildings damage assessment usingtexture features of polarization decomposition components[J].Journal ofRemote Sensing,2017.]。其中，光谱特征包括R波段均值、G波段均值、B波段均值、R波段标准差、G波段标准差、B波段标准差、R波段贡献率、G波段贡献率、B波段贡献率、亮度；几何形态学特征包括Area(区域面积)、MajorAxisLength(与区域具有相同标准二阶中心矩的椭圆的长轴长度)、Eccentricity(与区域具有相同标准二阶中心矩的椭圆的离心率)、Orientation(与区域具有相同标准二阶中心矩的椭圆的长轴与x轴的交角)、MinorAxisLength(与区域具有相同标准二阶中心矩的椭圆的短轴长度)、Extent(同时在区域和其最小边界矩形中的像素比例)、Perimeter(区域的周长)、ConvexArea(填充区域凸多边形图像中用于填充空洞的所有像素的个数)、EquivDiameter(与区域具有相同面积的圆的直径)、Solidity(同时在区域和其最小凸多边形中的像素比例)；纹理特征包括灰度共生矩阵对比度、同质性、相关性、熵、Jvalue、Roberts算子、Sobel算子、Prewitt算子、Laplacian算子、Canny算子。

Step3：以5棵决策树构建初始随机森林模型，输入包含30种特征的训练样本集H。定义分类正确率为每棵树的分类正确的样本数与输入的样本数的比值，比值越大说明分类结果与真实情况越接近。计算每棵树的分类正确率T_tree1～T_tree5。，并求所有决策树的正确率均值

其中i表示决策树的序号；

Step4:以5棵树为步长，迭代增加决策树后构建新的随机森林模型。采用与Step3相同的步骤，获得迭代次数为n时的正确率均值

(n∈[1,99])；

Step5：若某Acc_n满足Acc_n-1＜Acc_n＜Acc_n+1，则认为Acc_n是一个包含10颗决策树的子区间内的峰值点，记Acc_opt0＝Acc_n；

Step6：继续计算Acc_n后连续出现的3个峰值点，分别计为Acc_opt1、Acc_opt2、Acc_opt3。若满足：

Acc_opt0≥argmax{Acc_opt1,Acc_opt2,Acc_opt3} (2)

则记Acc_opt0对应的决策树颗数为I_opt0。否则，进入下一步；

Step7：在n∈[1,99]区间内，重复Step5和Step6，n为迭代次数，直到取得满足公式(2)的Acc_opt0，停止迭代并提取I_opt0。若在n∈[1,99]内无法确定Acc_opt0，则记已经获得的Acc_n最大值对应的决策树颗数为I_opt0；

Step8：继续比较I_opt0与其前、后相邻的4棵数(共9种模型)的分类正确率，根据分类正确率最大值确定最终提取的优化决策树棵数I_opt。

(3)特征重要度指引下的特征集合优化

在对随机森林的决策树颗数进行的优化选择基础上，为进一步减少候选特征集合F中特征间的冗余和证据冲突，本发明定义了一种特征的重要度指标，进而提出了特征集合优化策略：在包含I_opt颗决策树的随机森林模型中，首先计算所有特征的重要度，再分别求取光谱、纹理和几何形态学三类特征的类间重要度，再分别对三类特征中包含的子特征进行类内重要度排序，最后在重要度指引下获得优化后的震害建筑物分类特征集合。具体步骤如下：

Step1:计算所有特征的重要度。在训练样本集H中，经过有放回的随机采样(bootsrap)后未被选中的样本组成袋外数据集(out-of-bag，OOB)。任意特征f_t(f_t∈F,t∈[1,30])对于第i单棵决策树的重要度由公式(3)计算得到：

式中Φ_B表示OOB样本集合，x_j和l_j分别表示袋外数据中任意样本与其被赋予的类别标签，c_j ⁽ⁱ⁾表示样本x_j所获得的类别标签，

表示将特征f_t的值替换为其他随机值后样本x_j所得的类别标签，N(·)为计数函数。遍历所有决策树，可求得特征f_t对于随机森林分类器的重要度

Step2:计算特征的类间重要度。由于F中的特征分别归属于光谱、纹理和几何形态学(下文中分别用下标L，T，G表示)三类，因此根据每个f_t的重要度W(f_t)及其所属的类别分别求和，可得到光谱、纹理和几何形态学三类特征的类间重要度，分别记为W_L，W_T和W_G。在此基础上，定义光谱特征的类间冗余归一化惩罚因子为

以此类推，可得到纹理、几何形态学特征的类间冗余归一化惩罚因子分别为

和

Step3:特征的类内重要度排序。将30种特征分别按照其所属的光谱、纹理或几何形态学分为三类特征子集，分别仅采用这三类特征子集单独进行模型训练，采用Step1的步骤计算特征重要度。在此基础上，对各个特征子集内的特征按照重要度从高到低排列从而得到类内重要度排序。

Step4：对各个特征子集，依据冗余归一化惩罚因子的比例取整后，分别仅保留类内重要度相对较高的v_L，v_T和v_G个特征，从而获得共包含V＝v_L+v_T+v_G维特征的优化特征集合F_opt：

其中F_opt是提取的特征集合；δ，β，γ分别代表按各自类内重要度排序并抽取后的三大类特征。若加权后的抽取个数大于某类特征总数，则取该类全部特征作为分类特征。

(4)基于优化随机森林模型的影像分类

基于所提取的F_opt及I_opt，构建优化随机森林模型如下：

式中，P(x)表示分类结果；p_i(x)表示单棵决策树的分类结果；x为样本；c表示分类标签，c∈{Undamaged Building,Damage Building,Ruin,The Others}，分别对应完好建筑物、震害建筑物、废墟和其他；E(·)为示性函数，即当结果为Undamaged Building时，取值为1，Damage Building、Ruin、The Others依次为2,3,4。此基础上，根据每棵决策树给出的分类标签进行投票表决，以票数多作为该判别样本最终分类标签的标准。若属于某几类的投票结果出现平局情况，则根据公式(6)判断该样本与此几类训练样本之间的距离Dist，取Dist较小的类别作为该样本的最终分类结果。

式中,x_test(v)和x_train(v)分别为测试样本和训练样本中第v个特征的值；

是第v个特征的方差。

为了说明本发明的有益效果，实验采用四组不同地区、不同传感器的震后高分辨率遥感影像，通过目视分析和定量精度评价，并与多种先进的震害建筑物识别方法进行了对比，验证了所提出方法的性能。

实验数据：Dataset 1和Dataset 2为中国青海省玉树地区的GE01卫星遥感影像，采集时间为2010年5月6日。地震发生时间为2010年4月4日，地震最高震级为7.1级。影像包括全色波段和多光谱(蓝、绿、红、近红外)波段，空间分辨率分别为0.41m和1.65m，尺寸为1024×1024像素。实验采用ENVI软件融合后的空间分辨率为0.41m的pan-sharpened RGB影像。如图2(a)，(b)所示。Dataset 3和Dataset 4为中国四川省汶川地区的QuickBird卫星遥感影像，采集时间为2008年6月3日。地震发生时间2018年5月12日，地震最高震级为8.0级。影像包括全色波段和多光谱(蓝、绿、红、近红外4个波段)波段，空间分辨率分别为0.6m和2.4m。实验采用ENVI软件融合后的空间分辨率为0.6m的pan-sharpened RGB影像，如图2(c)，(d)所示

采用WJSEG方法对影像进行分割，所获得的联通区域构成初始对象集合R_origin。为便于观察，首先在原始影像中叠加一个半透明白色图层，再将分割结果用黑色像素表示并投影到原始影像中，如图3所示。

由图可以看出，分割结果较为完整地提取了完好建筑物与震害建筑物，且几乎不存在欠分割现象；过分割现象则主要存在于大尺寸的河道，草地区域中，并非本研究的感兴趣对象Region of Interesting(ROI)。因此，所提取的初始对象集合可为后续建筑物震害识别提供有效的分析基元。

根据前述规则集合对影像非建筑物对象进行筛选，结果如图4所示。其中，被剔除的非建筑物对象用黑色像素表示，剩余白色对象构成潜在建筑物集合R_BLD。同时，为便于分析，本研究采用字母符号(alphabetic symbol)对影像中的一些位置和对象进行了标注，后续内容中也采用了同样的方式。

如图4所示，Dataset 1和Dataset 2中狭长的道路(如位置A、C)，Dataset1中面积较小的帐篷(如位置B)，Dataset3和Dataset4中的河道(如位置G、E)，以及Dataset4中大面积的植被区域(如位置F)均得到了有效剔除。尽管如此，仍有一些非建筑物对象没有被剔除，如Dataset2中的裸地(如位置D)，Dataset3中的植被、裸地等(如位置H)。因此，所提出的非建筑物筛选策略是可行且有效的，但同时在后续的分类结果中也有必要保留其他地物这一类别。

Dataset 1～Dataset 4所确定的优化决策树棵数分别为131,64,91和63。下文以Dataset 1为例进行详细说明。如图5所示，依据分类正确率曲线波动判别规则，所选定满足条件的峰值点为I_opt0＝130。在此基础上，以I_opt0为中心构造一个区间，分别计算I_opt0前后4棵树的分类正确率，并比较区间内9种模型的分类正确率。其中，当决策树颗数为131时分类正确率取得最大值，因此取I_opt＝131。

对候选特征集合F进行筛选后，四组实验数据集合对应的优化特征集合如表1所示：

表1特征优化集合

为客观分析和验证本发明方法的性能，本发明采用传统随机森林以及两种先进的震害建筑物识别方法进行了对比实验。方法1[Breiman L.Random Forest.MachineLearning,2001,45(1):5-32.]采用传统随机森林方法，其输入特征采用初始建筑物震害特征集合F，决策树棵树为500棵。通过与方法1对比，以分析本发明所提出的特征及决策树颗数优化策略的有效性；方法2采用文献[Dou J,Yunus A P,Tien Bui D,et al.Assessmentof advanced random forest and decision tree algorithms for modeling rainfall-induced landslide susceptibility in the Izu-Oshima Volcanic Island,Japan[J].Science of The Total Environment,2019,662:332-346.]中提出的震害建筑物识别方法。此方法采用十折交叉验证法对训练样本集进行优化选择，再利用传统随机森林进行分类。其中，输入特征为本发明所提取的优化特征集合F_opt，决策树棵树则按照该文献中所使用的500棵。通过与方法2对比，有助于分析本发明方法自适应提取的决策树颗数对提高识别精度的作用。方法3采用文献[王之,刘超,刘秀菊,et al.基于SEaTH算法的芦山地震无人机低空遥感影像信息对象级分类[J].地震研究,2018,v.41；No.186(02):19-25+189.]中提出的震害建筑物识别方法。该方法利用改进的SEaTH算法进行了特征优化，再采用隶属度分类法对影像进行分类。由于方法3同样从特征优化角度出发，同时采用了不同的分类器，因此通过与该方法对比有助于客观评价所提出方法的总体精度。

实验中，由于方法2属于像素级分类方法，难以直接与本发明对象级分类方法结果进行比较，因此直接利用本发明提取的初始对象集合R_origin中的对象代替像素作为后续分类的基本单元，其他实现步骤与原文一致；由于方法3为对象级分类方法，为保证分类对象的一致性，同样采用R_origin中的对象作为震害建筑物识别的基本单元，其他实现步骤均与原文一致。此外，为了避免由于是否进行了非建筑物初筛所导致的震害建筑物识别结果差异，四种对比方法均在本发明所提取的潜在建筑物集合R_BLD基础上进行特征提取和分类。此外，实验结果将建筑物划分为完好建筑物、震害建筑物、废墟和其他四类，并分别用不同颜色进行表示。本发明方法和三种对比方法的震害建筑物识别结果如图6～图9所示。

为了开展对实验结果的精度评价，本发明依据目视解译和实地考察绘制了四组数据的参考map，并以此作为后续精度评价的依据，如图10所示。

目视分析：通过将实验结果图和参考图对比，本发明方法对震害建筑物的识别效果优于其他三种对比方法。主要体现在：(1)在四组实验结果中，位于位置I、Q的完好建筑物、位置O的非建筑物、位置R的废墟只有本发明方法正确识别出来，其他方法均出现了漏检或者错检；(2)对于形状、纹理规则的完好建筑物如位置S、W、Z、C*，四种方法都做出了正确的判别结果；但对于Partly damaged building，如位置K、P的只有本发明方法和方法3做出了正确判别，位置L、J只有本发明方法分别和方法2、3做出了正确判别；(3)对于废墟，如位置U、A只有本发明和方法2做出了正确判别；(4)在构建潜在建筑物集合时未被筛掉的，与建筑物有相似光谱、形状特征的非建筑物，如位置V、X只有本发明方法和方法2做出正确判别；对于一些面积较小的对象，如B*、Z只有本发明方法和方法做出正确判别。

定量分析：在目视分析的基础上，本发明采用总体精度(Overall accuracy)、误检率(False positive)、漏检率(False negative)、Kappa系数四个指标进行定量精度评价。本发明方法和三种对比方法的精度评价表分别见表2～表5。

表2 Dataset1震害建筑物识别精度评价表

表3 Dataset2震害建筑物识别精度评价表

表4 Dataset3震害建筑物识别精度评价表

表5 Dataset4震害建筑物识别精度评价表

由表2～表5可以看出，本发明方法总体精度能够达到85％以上，且四种精度指标均优于另外三种对比方法，与目视分析结果一致。在数据集1和数据集4实验中，4种对比方法的精度均显著高于其他两组数据集，这主要是由于非建筑物初筛结果不同所导致的，具体体现在：数据集1和数据集4中被筛选掉的非建筑物对象占初始对象集合R_origin中的比例分别为67.84％和66.6％,显著高于数据集2和数据集3的34.66％和38.35％。这也就意味着，数据集2和数据集3中的精度指标能够更加真实的反映不同方法之间的性能差异。因此，我们重点对这两个数据集的精度指标展开分析和讨论：与方法1对比，所提出方法的总体精度提高了5％以上，从而证明了所提出的特征及决策树颗数优化策略是必要、可行且有效的。方法2采用不仅采用了优化训练样本集合，同时特征集合也与所提出方法一致，但由于决策树颗数均设定为500颗，其总体精度依然低于所提出方法5％以上。这也再次证明了确定合理的决策树颗数对于震害建筑物识别精度的重要作用。本发明方法和方法3都是在分类前对初始特征集合进行了优化，但总体精度低于所提出方法超过8％，这也在一定程度上证明了本发明所构建的优化随机森林模型的合理性和有效性。

决策树数量对总体精度影响分析：为进一步分析决策树数量对总体精度的影响，以及评价本发明所自适应提取的决策树颗数的合理性。本发明依据参考图及所提出方法，在[50,200]区间内以10为间隔，统计四组实验中总体精度与决策树数量的变化曲线，并对对统计获得的最大、最小及平均总体精度及本发明获得的总体精度分别用不同样式的直线表示，如图11所示。

如图所示，决策树数量的变化对总体精度具有显著影响，因此选择合理的决策树颗数是必要的。另一方面，尽管所提取的决策树颗数对应的总体精度均非[50,200]区间内对应的最大值，但差异均小于1％，且远高于平均总体精度。因此，所提出的决策树颗数自适应选择策略是可行且有效的，不仅提高了算法自动化程度，同时取得了理想的分类精度。

特征组合对总体精度影响分析

除决策树棵树外，本发明进一步分析了不同类别特征及特征组合对震害建筑物检测精度的影响。在初始建筑物震害特征集合F中，依据特征所属类别不同，统计获得的总体精度如表6所示。

表6不同特征及特征组合所获得的总体精度

如上表所示，本发明所提取的特征集合效果最好。在单独使用光谱、几何形态、纹理特征进行分类时，光谱特征的总体精度明显高于另外两种特征，其原因在于震害建筑物的几何形态和纹理特征均遭到破坏，增加了分类的不确定性，而光谱特征不易发生大的改变因此在震害建筑物识别这一特定应用领域中更加可靠。而采用三大类特征所取得的总体精度甚至低于单独采用光谱特征对应的总体精度，其原因就在于存在特征冗余和证据冲突，而本发明所给出的特征集合优策略给出了一种有效的解决方案。

综上，在震前参考信息缺乏前提下，本发明提出了一种基于决策树与特征优化的震害建筑物识别方法。在对不同地区、不同传感器的多组震后高分辨率遥感影像的开展的实验中，其总体精度均可达到85％以上，且误检率小于6％，可为震后应急响应及灾后重建提供关键而可靠的决策支持信息。其理论主要贡献在于，(1)提出的基于特征重要性指标的特征集优化策略，为地震破坏建筑特征空间的自动构建提供了可行的解决方案；(2)提出的特征集筛选策略，结合决策树数的自适应提取策略，构建了一种新的高效的地震损伤建筑识别优化RF模型。

Claims (8)

1.一种基于决策树与特征优化的遥感影像震害建筑物识别方法，其特征在于，包括如下步骤：

(1)对卫星遥感影像进行分割提取初始对象集合，并结合几何形态学特征提取潜在的建筑物对象集合；

(2)基于分类正确率曲线波动判别规则自适应确定决策树数量；

(3)在特征重要度指标指引下对光谱、纹理、几何形态学三大类特征进行筛选，获得具有代表性的震害特征集合，进而构建优化的随机森林模型；

(4)根据随机森林判别结果将潜在建筑物对象集合进一步进行分类，识别出震害建筑物。

2.根据权利要求1所述的基于决策树与特征优化的遥感影像震害建筑物识别方法，其特征在于，所述步骤(1)中采用高分辨率遥感影像分割算法WJSEG提取遥感影像的初始对象集合。

3.根据权利要求1所述的基于决策树与特征优化的遥感影像震害建筑物识别方法，其特征在于，所述步骤(1)中结合几何形态学特征提取潜在的建筑物对象集合的具体规则包括面积规则、矩形度规则和MBI规则；所述面积规则是根据对象所包含的像素个数确定是否剔除，若像素个数不超过设定的阈值，则剔除；所述矩形度规则是根据对象所包含的像素个数与最小外接矩形所包含的像素个数的比值和最小外接矩形的长宽比确定是否剔除，若比值小于设定的阈值且长宽比大于设定的阈值，则剔除；所述MBI规则是基于影像中所有像素的MBI值利用最大类间方差法自适应确定分离阈值，获得对象内属于非建筑物像素的比例，若比例大于设定的阈值则剔除对象。

4.根据权利要求1所述的基于决策树与特征优化的遥感影像震害建筑物识别方法，其特征在于，所述步骤(2)中基于分类正确率曲线波动判别规则自适应确定决策树数量，包括：

(2.1)在潜在的建筑物对象集合中，采用人工标记的方式提取样本，构成训练样本集；

(2.2)确定初始建筑物震害特征集合，包括光谱特征、纹理特征和几何形态学特征；

(2.3)构建初始随机森林模型，输入训练样本集，定义分类正确率为每棵树的分类正确的样本数与输入的样本数的比值，比值越大说明分类结果与真实情况越接近，计算每棵树的分类正确率并求所有决策树的正确率均值；

(2.4)以设定的步长，迭代增加决策树后构建新的随机森林模型，获得迭代次数为n时的正确率均值Acc_n；

(2.5)若某Acc_n满足Acc_n-1＜Acc_n＜Acc_n+1，则认为Acc_n是一个子区间内的峰值点，记Acc_opt0＝Acc_n；

(2.6)继续计算Acc_n后连续出现的3个峰值点，分别计为Acc_opt1、Acc_opt2、Acc_opt3，若满足Acc_opt0≥argmax{Acc_opt1,Acc_opt2,Acc_opt3}，则记Acc_opt0对应的决策树颗数为I_opt0；否则，进入下一步；

(2.7)重复步骤(2.5)和(2.5)，直到取得满足

Acc_opt0≥argmax{Acc_opt1,Acc_opt2,Acc_opt3}的Acc_opt0，停止迭代并提取I_opt0；若无法确定Acc_opt0，则记已经获得的Acc_n最大值对应的决策树颗数为I_opt0；

(2.8)继续比较I_opt0与其前、后相邻的N-1棵数的分类正确率，根据分类正确率最大值确定最终提取的优化决策树棵数I_opt，N为步骤(2.4)中设定的步长。

5.根据权利要求1所述的基于决策树与特征优化的遥感影像震害建筑物识别方法，其特征在于，所述步骤(3)中在特征筛选的方法是：在包含I_opt颗决策树的随机森林模型中，首先计算所有特征的重要度，再分别求取光谱、纹理和几何形态学三类特征的类间重要度，再分别对三类特征中包含的子特征进行类内重要度排序，最后在重要度指引下获得优化后的震害建筑物分类特征集合；其中I_opt为步骤(2)中优化决策树棵数。

6.根据权利要求5所述的基于决策树与特征优化的遥感影像震害建筑物识别方法，其特征在于，任意特征f_t对于第i单棵决策树的重要度由下式计算得到：

式中Φ_B表示OOB样本集合，x_j和l_j分别表示袋外数据中任意样本与其被赋予的类别标签，c_j ⁽ⁱ⁾表示样本x_j所获得的类别标签，

表示将特征f_t的值替换为其他随机值后样本x_j所得的类别标签，N(·)为计数函数；特征f_t对于随机森林分类器的重要度

7.根据权利要求6所述的基于决策树与特征优化的遥感影像震害建筑物识别方法，其特征在于，将光谱、纹理和几何形态学三类特征的重要度分别求和得到三类特征的类间重要度，分别记为W_L，W_T和W_G，并分别计算类间冗余归一化惩罚因子；对各个特征子集内的特征按照类内重要度从高到低排列；对各个特征子集，依据冗余归一化惩罚因子的比例取整后，分别仅保留类内重要度相对较高的v_L，v_T和v_G个特征，从而获得共包含V＝v_L+v_T+v_G维特征的优化特征集合F_opt。

8.根据权利要求1所述的基于决策树与特征优化的遥感影像震害建筑物识别方法，其特征在于，所述步骤(3)中基于提取的优化决策树棵数I_opt及优化特征集合F_opt，构建的优化随机森林模型如下：

式中，P(x)表示分类结果；p_i(x)表示单棵决策树的分类结果；c表示分类标签，c∈{Undamaged Building,Damage Building,Ruin,The Others}。

Images