一种结合树木阴影特征的遥感影像毛白杨识别方法
【技术领域】
[0001]本发明是一项对地观测领域的技术，针对高分辨率多光谱遥感影像光谱及纹理特征识别树种信息存在精度低、难度大的问题，本发明能够有效修正树种识别结果，大幅提高识别精度，是一项具有研究及应用价值的空间分析技术。
【背景技术】
[0002]树木在生态系统中扮演着重要的角色，具有涵养水源、保持水土、净化空气、卫生保健、美化环境等功能，对于保护自然界生态平衡具有非常重要的作用。
[0003]传统的森林监测及城市树种精细调查主要通过基层单位实地抽样调查完成，存在投入的资金和人力大、完成的时间周期长、条件恶劣区作业困难、数据精度低、缺乏空间统计分析功能等问题。
[0004]航天遥感是一项综合性探测技术，随着其分辨率不断提高，高分辨率多光谱遥感技术在森林类型解译、蓄积测定、健康评价等方面得到越来越广阔应用，通过遥感技术快速提取树木信息、识别树种对城市树种调查、配置、管理及森林生态系统监测具有重要意义。
[0005]目前，用于树种精细识别的商用卫星遥感空间分辨率可达0.4m(W0rldvieW-3)，完全能分辨树丛甚至单颗树木，但是要准确分辨具体的树种还存在较大的难度。在树种识别方法方面，相对于传统的基于像元的处理方法而言，随着高空间分辨率遥感影像的出现，面向对象方法受到越来越广的应用。研究表明，面向对象方法能够通过具体的实体对象来分析地物的特征，更接近提取信息的思维逻辑，能够极大地提高高分辨率影像的信息提取精度。Le Wang等用面向对象分类技术识别不同的红树林树种，并和传统的GLCM纹理特征比较，得到了更高的精度；Su等把对象的光谱和GLCM纹理信息结合起来，对复杂城市区域进行实验，结果表明纹理和光谱信息结合能够有效改善面向对象分类的效果。
[0006]以上方法都是针对对象的光谱及纹理特征进行树种的识别和提取，然而不同种类的树木可能存在光谱及纹理特征相似的情况，给树种的识别和提取带来很大困难。阴影常被视为遥感影像地物识别的干扰信息，然而，从另一个角度，阴影包含着物体的三维信息，通过阴影特征分析可以反映不同地物的高度差异。因此，根据生态系统群落特征，结合各树种阴影差异来辅助遥感影像光谱、纹理特征进行树种高精度识别具有重大的研究价值和推广空间。
[0007]毛白杨是城乡绿化、速生用材以及防护林的优质树种，其品种适应性强、枝叶茂密，广泛分布于我国以黄河流域中下游为中心的十多个省。毛白杨高可达30米，在植被生态系统群落中往往处于顶层，在高分辨率多光谱遥感影像上可体现出比其它树种更明显的阴影特征。基于此特征，本发明试图建立一种结合树木阴影特征的遥感影像毛白杨识别方法。该方法通过引入树木的阴影特征，可有效弥补遥感影像光谱及纹理特征相似的树种难以区别的技术难题，大幅度修正毛白杨识别结果，对处于生态系统上游的高大树种识别具有重大的借鉴意义。
【发明内容】
[0008]本发明提出一种结合树木阴影特征的遥感影像毛白杨识别方法，通过将树木阴影特征作为树种识别的辅助特征，有效解决了遥感影像光谱及纹理特征识别树种的技术难题，有益于研究、应用和进一步推广。该度量方法包括如下步骤:
步骤1)对研究区遥感影像进行多尺度分割，形成分割结果图；步骤2)基于所述分割结果图中每个对象的光谱及纹理特征进行毛白杨初步识别，形成毛白杨初步识别结果图；步骤3)基于所述分割结果图中每个对象的光谱特征进行阴影提取，形成阴影提取结果图；步骤4)从所述阴影提取结果图中选取所述毛白杨初步识别结果图中各毛白杨斑块所对应的阴影，并估算其阴影长度；步骤5)基于所述阴影长度进行树木高度的估算，作为度量各毛白杨斑块平均树高的参考值；步骤6)基于所述毛白杨斑块平均树高信息，对所述毛白杨初步识别结果进行修正，形成毛白杨精确识别结果图。
[0009]进一步地，所述步骤1)中研究区遥感影像分割的方法如下:采用多尺度分割算法，确定最大分割尺度130，尺度数3，尺度间隔0.6，将遥感影像分割成高度同质性的分割对象，形成分割结果图。
[0010]进一步地，所述步骤2)中毛白杨初步识别方法如下:利用步骤1)中形成的分割结果图中每个对象的光谱特征亮度Brightness及纹理特征灰度共生矩阵同质性GLCM-H0M0，采用多参数阈值法，即使用如下计算公式(a ( Brightness ( β )and( γ彡GLCM-HOMO彡δ )，形成毛白杨初步识别结果图，其中，α、β为Brightness的阈值，γ、δ为GLCM-HOMO的阈值。
[0011]进一步地，所述步骤3)中阴影提取方法如下:利用步骤1)中形成的分割结果图中每个对象的光谱特征亮度Brightness，采用参数阈值法，即使用如下计算公式Brightness ^ ε，形成阴影提取结果图，其中，ε为Brightness进行阴影提取的阈值。
[0012]进一步地，所述步骤4)中毛白杨斑块阴影长度的估算方法如下:以步骤2)中毛白杨初步识别结果图各毛白杨斑块为基础，从步骤3)阴影提取结果图中选取各毛白杨斑块所对应的阴影，形成毛白杨及其阴影提取结果图，沿着树木阴影方向在遥感影像上绘制一系列平行直线，在阴影边界与平行直线的交点处对平行直线进行打断，求取每个阴影内所有平行线段长度均值作为相应毛白杨斑块的阴影长度。
[0013]进一步地，所述步骤5)中毛白杨斑块平均树高的估算方法如下:基于各毛白杨斑块的阴影长度，根据遥感影像成像时刻的太阳高度角进行树高Height估算，即使用如下计算公式Height=(Length+r).tan Θ，得到各毛白杨斑块平均树高信息，其中，Length为阴影长度，Θ为太阳高度角，r为毛白杨冠幅平均半径。
[0014]进一步地，所述步骤6)中毛白杨精确识别方法如下:根据步骤5)中所述的毛白杨斑块的高度特征，对所述的步骤2)中毛白杨初步识别结果进行修正，采用阈值法，即使用如下计算公式Height多σ，剔除初步识别结果图中树高小于σ的对象，形成毛白杨精确识别结果图，其中，σ为毛白杨精确识别的高度阈值。
【附图说明】
[0015]图1为分割结果图；
图2为毛白杨初步识别结果图； 图3为阴影提取结果图；
图4为毛白杨及其阴影提取结果图；
图5为辅助阴影长度计算的平行线段分布图；
图6为阴影长度计算结果图；
图7为树高、阴影长度、树冠幅半径及太阳高度角关系图；
图8为毛白杨斑块平均树高计算结果图；
图9为毛白杨精确识别结果图；
图10为误剔除毛白杨斑块位置图。
【具体实施方式】
[0016]下面结合附图对本发明“一种结合树木阴影特征的遥感影像毛白杨识别方法”作进一步说明。
[0017]本发明试图引入树木的阴影特征来辅助高分辨率多光谱遥感影像光谱和纹理特征进行毛白杨精细识别，大幅度提高了毛白杨的识别精度，为高分辨率多光谱遥感影像树种识别和提取提供参考。
[0018]1)研究区遥感影像分割
研究区位于奥林匹克森林公园，覆盖0.72kmX0.48km区域范围，研究区分布着多个毛白杨林地斑块，各斑块内毛白杨分布密集，数量从数颗到上百颗不等，平均树高15m以上。除毛白杨外，研究区分布着柏树、松树、白蜡、柳树、天目琼花、鞑靼忍冬、郁香忍冬等多种其它植被景观，高度由lm至10m不等。各毛白杨斑块与其它树木交相辉映，形成高低错落有致的生态系统群落。在0.4m高空间分辨率Worldview-3标准假彩色遥感影像上，毛白杨斑块色调暗红、纹理略粗。
[0019]根据研究区毛白杨影像特征，采用多尺度分割算法进行遥感影像面向对象分割。多尺度分割算法是较为常用的一种分割方法，能够生成高度同质性的影像分割区域，从而以最佳的尺度分离和表示地物目标。在某指定的尺度下分割时，通过合并相邻的像素或小的分割对象，在保证对象内部像元之间同质性最大的前提下，基于区域合并技术实现影像分割。
[0020]将最大分割尺度设为130，尺度数设为3，尺度间隔设为0.6，即最大尺度、中间尺度和最小尺度分别为130、130X0.6=78,78X0.6=47，分割中颜色因子占的比例设为0.7，分割中形状因子中平滑度因子所占的比例设为0.5，将研究区遥感影像分割成高度同质性的影像分割单元，将每个分割单元作为一个对象。
[0021]图1为分割结果图，分割结果表明，分割尺度适中，未造成地物景观过度破碎，且各对象内部树木种类均一，为毛白杨精细识别奠定了良好基础。
[0022]2)毛白杨初步识别
本发明充分综合毛白杨的光谱特征亮度Brightness及纹理特征灰度共生矩阵同质性GLCM-H0M0,采用多参数阈值法进行面向对象毛白杨初步识别。
[0023]亮度是最基本的光谱特征，反映了各波段像元亮度的平均值。
[0024]纹理特征是细小物体在影像上大量重复出现所形成的规律和特征，它是大量个体的大小、形状、阴影和色彩的综合反映，它描述像元亮度的空间变化特性。灰度共生矩阵(GLCM)是广泛使用的纹理特征提取方法，其基本原理是计算局域范围内像元灰度级共同出现的频率，不同的空间关系和纹理会产生不同的共生矩阵，以此来区分不同的纹理和结构特性。同质性是灰度共生矩阵常用的统计测度方法，是影像纹理相似性的度量，其值越高代表局部区域缺乏变化，具有较小的灰度差异。遥感影像上道路、草地、屋顶的特征值较大，而树木、小路、阴影的特征值较小。
[0025]对每个分割对象进行光谱特征亮度Brightness及纹理特征灰度共生矩阵同质性GLCM-H0M0 的计算。
[0026]采用多参数阈值法，即Brightness、GLCM-HOMO分别取阈值，取二者的交集(见公式1 )，完成毛白杨的初步识别。
[0027]( α Brightness β ) and ( γ GLCM-HOMO δ ) (1)
其中，α、β为Brightness的阈值，γ、δ为GLCM-H0M0的阈值。
[0028]在确保识别结果中最大程度涵盖毛白杨信息的前提下，主要根据毛白杨在影