基于多特征融合的遥感影像变化检测算法.doc

基于多特征融合的遥感影像变化检测算法 第30卷，第2期国土资源遥感V30，N.22018年6月REMOTESENSINGFORLAND&RESOURCES Jun.，2018doi10.6046/gtzyyg.2018.02.13引用格式王光辉，李建磊，王华斌，等.基于多特征融合的遥感影像变化检测算法J.国土资源遥感，2018,3 (2):93-99.Wang GH，Li JL，Wang HB，et al.Change detection based on adaptive fusionof multijDlefeaturesJ.Remote Sensingfor Land andResources，2018，30 (2)93-99.基于多特征融合的遥感影像变化检测算法王光辉李建磊2，王华斌杨化超2(1.国家测绘地理信息局卫星测绘应用中心，北京100830;2.中国矿业大学（徐州），徐州221116)摘要针对传统的变化检测算法主要依赖像斑的光谱信息，未能有效地利用影像多特征检测优势的问题，基于面向对象的分析思想，提出一种多特征融合的遥感影像变化检测算法。 首先，以多尺度分割的影像对象为基础，统计各对象的颜色直方图和边缘直线梯度直方图；然后，利用推土机距离计算不同时期对象之间的颜色距离和边缘直线特征距离，采用自适应加权方法将颜色距离和边缘直线特征距离组合构建对象的异质性；最后，采用直方图曲率分析获得像斑的变化检测结果。 实验结果表明，该方法能够充分融合颜色和边缘直线特征，提高变化检测的精度。 关键词变化检测；颜色直方图；直线梯度直方图；直方图距离；直方图曲率中图法分类号#TP79文献标志码A#1001-070X (2018)02-0093-070引言随着空间信息技术的飞速发展，遥感影像的空间分辨率、光谱分辨率和时间分辨率不断提高，其获取的地表信息更加丰富，应用领域越来越广。 如何利用高空间分辨率遥感图像丰富的地物细节信息，并抑制成像条件与自然条件带来的信息干扰，是目前遥感影像变化检测研究中倍受关注的问题。 传统的单像素、单一特征的变化检测算法已经不能满足变化目标检测的准确性和完整性1，面向对象多特征融合的变化检测算法2-3&已成为高空间分辨率遥感变化检测研究领域的热点。 由影像派生的多特征是变化检测信息的载体，主要包括光谱、边缘和纹理等特征。 利用影像对象特征进行多时相变化检测时传统的变化检测算法选用单一特征，如Bovlo等4&采用像斑均值信息进行变化检测；庄会富等5提出基于纹理特征的影像变化检测算法，通过多方向、多个纹理特征构建差异影像提高了变化检测的精度；李松等6提出结合纹理分析与比率变换的变化检测算法，通过分析纹理特征，统计纹理的指标，有效地识别出滑坡灾害信息;钟家强等7提出线特征变化检测算法，通过提取影像的边缘梯度信息、压缩、拟合等处理，获得变化检测结果，避免了线匹配的复杂过程，具有较强的实用性。 虽然以上的算法都从不同的角度验证特征的有效性和可用性，但单一特征模型的适用性和精度并不理想。 刘臻等8&利用梯度和纹理相似度进行变化检测；李亮等9&将光谱信息与纹理信息结合提出多特征融合的变化检测算法；魏立飞等10&提出应用多特征流行嵌入模型，充分综合多特征的信息，获得了完整的变化信息。 但上述多特征融合方法多采用光谱与纹理的特征，未能考虑颜色的空间分布信息以及影像的边缘特征。 鉴于上述情况，本文通过研究不同特征的特点，提出多特征自适应融合的变化检测方法，以影像对象的颜色直方图为依据，通过统计颜色直方图、边缘直线梯度直方图的最大概率差别自适应确定权值，实现颜色特征与边缘直线特征的自适应加权组合，为下一步变化强度阈值的选取做准备。 1基于多特征融合的变化检测为了充分利用遥感图像的颜色空间分布信息以及边缘直线梯度特征，本文提出一种基于面向对象思想的多特征融合变化检测方法。 首先，运用多尺度分割算法（multi-solution segmentation)获得影像对象，并统计影像对象的颜色直方图和边缘直线梯#xx-10-08;修订日期#xx-12-03基金项目高分辨率对地观测系统重大专项项目“高分遥感测绘应用示范系统（一期）”（编号AH1601)资助。 第一作者王光辉（1982-),男，工程师，主要从事遥感影像数据处理、信息提取与变化检测等技术研究和软件开发工作。 Emal:iangghsasmac.。 +94+国土资源遥感2018年Fig.1.1多尺度影像分割影像分割是面向对象影像分析的基础，其结果将对后续的分析产生重要的影响11。 本文选用多尺度分割算法12，通过计算各波段的形状与光谱异质性的特征值，根据各个波段信息自适应确定权重，进行迭代运算，直到影像对象的形状和光谱的综合加权值大于给定的阈值时停止迭代，完成影像的多尺度分割。 对于任意波段的影像&其异质性/&为的权重）为形状异质性，通过形状平滑度的权重2*形状平滑度；s m和形状紧凑度；。 m确定。 1.2颜色和边缘直线特征提取直方图是一种空间统计特征，颜色和边缘直线梯度直方图能较好地表达对象的空间特征分布情况，为下一步计算前后2期差异提供依据。 本文采用影像对象的颜色直方图和边缘直线梯度直方图来描述不同时相影像对象的光谱特征和空间特征。 1.2.1颜色特征提取颜色作为影像上一种直观显著的重要视觉特征，对影像的形态、大小、方向和视角依赖较小。 窦建军等13&曾将全局颜色直方图应用于图像检索领域。 常用的RGB模型都存在以下不足不符合人对颜色的感知心理）在颜色感知上存在较小差异的2个点，难以在非均勻分布的RGB颜色空间中以2个点来表示。 HSV颜色空间是面向视觉感知的模型，符合人眼感知的色调、亮度和饱合度3要素，能更好地模拟人工变化检测的过程，因此能更好地用于提取颜色特征。 颜色直方图的提取主要包括2步RGB模型与HSV模型转换；HSV的量化。 N!丄（#+K+H)槡Vmin(#，H，K)Oe2$H(K H (3) (5)， (8)式中#，H和K分别为红光、绿光和蓝光波段影像。 经过HSV颜色变换，影像上的每一个像素点都可以用O,9和V来表示；然后，依据人眼对颜色感知的特性，对HSV3分量均勻量化，将O,9和V分别空间量化为8,5和5份；最后，统计获得影像对象的颜色直方图。 1.2.2边缘直线特征提取高空间分辨率遥感图像纹理比较丰富，采用canny，soble和lapace等算子可以直接提取低矮建筑物区的边缘，但边缘比较琐碎，难以表达影像对象的特征。 本文采用linesegment detector算子14检测遥感图像中几何形状明显的目标。 由于直线区域的灰度通常是不连续的，直线上每一像素对应一个梯度方向，基于直线像素点的梯度方向建立的直方图能够表示影像对象的形状特征，对于发生变化的对象能很好地反映前后2期的差异。 对于影像上点(/G处/(/G的梯度定义为=Gx,GyT=,(V)梯度是一个矢量，其大小和方向为G(f(x,y)=V f(x,y)|=G2x+G2y, (8)$(/，j)=arctan(H/H x)， (9)度直方图信息；然后，利用推土机距离（earth movers distance，E MD)计算不同时相影像对象之间的颜色距离和边缘直线距离，并自适应加权颜色距离和边缘直线距离，获得对象的异质性，对越接近于0变化可能性越低，远离0则变化可能性越高；最后，基于以上理论采用直方图曲率分析获得像斑的变化检测结果。 的异质性具体变化检测流程见图I图1变化检测流程,1Workflow ofchange detection首先，R GB模型转换成H SV模型，具体公式为一前后期矢量文件叠加一+-_I I_多尺度分割一|多尺度分割第！期王光辉，等基于多特征融合的遥感影像变化检测算法+95+式中为/(/G在该点的灰度变化最大的方向。 根据直线像素点方向统计直线梯度方向直方图(越远离0值地区，则变化的可能性越大。 影像上大部分是没有发生变化的，则异质性直方图是逐渐递减的一种形态。 本文采用异质性直方图曲率分析histogramcurvature analysis，HCA)方法确定变化检测阈值，该方法自适应能力较强。 影像对象的异质性值呈现出单侧正态分布趋势，且对象中变化区域和非变化区域的异质性往往会有很明显的差异，通过分析直方图的曲率可以确定变化检测的阈值。 依据微积分基本原理，曲率值越小表示曲线的弯曲程度比较缓，曲线发生突变的可能性就越小；而对于离散数据其曲率值7可近似表达为式中/&为第i个位置对应的频率值。 取7最大值对应的&乍为变化检测的最优阈值，实现最终变化区域的提取。 2结果与讨论2.1实验数据实验数据采用苏州市某城区的xx年和xx年资源三号卫星遥感影像，尺寸大小为1800像素$1420像素，空间分辨率为2.1m，由蓝光、绿光、红光和近红外4个波段组成；影像已经经过影像配准和直方图匹配等预处理。 实验数据如图4所示。 a)xx年b)xx年图4苏州市某城区影像Fig.4Su zhou Cityimages2.2实验分析实验采用多尺度分割算法对配准和直方图匹配后的多时相遥感影像进行分割，获得具备光谱和形状同质性的对象。 最优分割尺度根据经验和知识确定；形状因子参与形状一致度的权重，最优的形状因子能够得到合理边界的影像对象；紧致度表示影像对象的规则程度，紧致度过大则分割后影像地物不唯一，出现不纯净的对象。 最终确定分割参数分别为分割尺度为 100、形状因子为0. 45、紧致度为0.5；分割尺度为 80、形状因子为0. 5、紧致度为0.64。 然后将前后2个时期的分割结果进行叠加生成最终的矢量文件如图5a)所示。 图5b)为标准变化检测结果，其中黑色表示未变化区域，白色表示变化区域。 为了验证多特征融合算法的有效性，将单一特征、固定权重算法同本文算法进行了比较，检测结果如图6所示。 图6中椭圆为虚检区域，矩形为漏检区域。 6种变化检测结果的精度评定如表1所示。 第！期王光辉，等基于多特征融合的遥感影像变化检测算法+97+分割结果（b)标准变化检测结果图5分割结果及标准检测结果Fig.5Segmentation andstandard detectionresults(a)w?!1，lme!6b)W?!67，lme!63C(h SV!67，lme!67(d)!3，lme!?7(e)h?!6，lme!1(本文检测结果图J不同权重变化检测结果Fig.6Different weightresults ofchange detection表1不同权重变化检测结果精度比较Tab.1Auracy parisonby methodswith differentweights参数h?E!1lme!6hsv!6U7lme!63hsv!6?5lme!6?5变化未变化变化未变化变化未变化实际变化/m2684947.40781227.66926636.56575792.46814214.35739436.11实际未变化/m2319829.909796962.46391942.469682541.86393486.869629782.65合计/m21664777.3616572136.661318572.9616258334.261267695.1516369212.16虚检率6.326.366.33漏检率6.536.386.46正确率6.966.926.96参数h?!63lme!67hsv!16lme!1本文算法变化未变化变化未变化变化未变化实际变化/m21696269.13348681.261288899.66177276.661174865.98327581.62实际未变化/m2583248.279554768.76867573.669243159.26374527.629699933.8合计/m21673457.469963449.962156472.669426435.261549393.6616627514.16虚检率6.356.466.24漏检率2正确率6.926.916.94分析图8可以得出如下结论图5(a)仅利用漏检现象比较严重。 利用颜色特征与边缘直线特征颜色特征进行变化检测，存在较多的虚检和一定的加权进行变化检测，当颜色权重为.7时，虽然部分漏检。 图5(+仅利用边缘直线特征进行变化检测，图5(b)中虚检部分被标记为未变化，但仍存在一定?98?国土资源遥感2018年的漏检现象；当颜色权重与边缘直线权重都为0.5时，图5(c)也存在较多的漏检和虚检；当颜色权重为0.3时，与图5(A)和（=相比虚检有所下降，但是漏检增加；利用本文算法进行变化检测，图5(变化区域的识别仍然存在虚检和漏检，但虚检率和漏检率比固定权重有所降低。 文献9与文献16都采用了光谱特征与纹理特征融合，依据像斑灰度的标准差自适应确定阈值，其精度与本文精度的对比如表2所示。 由于不同时期遥感影像存在太阳高度角、光照条件和季节等因素不同程度的影响，使得不同地物在影像上具有不同的灰度值，而纹理特征仅反映地物结构信息，无法作为区分不同地物的主要依据，故实验结果表明，本文采用光谱特征与几何特征的变化检测精度最高。 表2不同变化检测算法精度比较Tab.2Auracy parisonof differentmethods ofchange detection参数文献9方法文献16方法本文算法虚检率0&320.320&24漏检率0&210&250&22正确率0&930.920&943结论基于面向对象的分析思想，本文提出了一种多特征融合的遥感影像变化检测算法。 实现了颜色特征与边缘直线特征的自适应加权组合，为融合HSV颜色直方图与边缘直线梯度直方图进行变化检测提供了一种思路。 采用直方图曲率分析获得像斑的变化检测结果，实现了变化检测最优阈值的自动确定。 实验结果表明，该方法依据像斑的光谱信息，充分、有效地挖掘和融合了影像多特征的检测优势，降低了虚检率和漏检率，提高了变化检测的正确率。 为了进一步提高变化检测的精度，后继研究将在本文基础上，引入纹理和形状等特征信息。 参考文献（References):1魏立飞，钟燕飞,张良培，等.遥感影像融合的自适应变化检测(英文）J.遥感学报,xx,14 (6):1196-1211.Wei LF,Zhong YF,Z hangLP,et al.Adaptive changemethod ofremote sensing image fusionJ.Journal of Remote Sensing,xx,14 (6)1196-1211.2Gong JY，Sui HG，Sun KM，et al.Object-level change detection based on full-scale imagesegmentation andits applicationto WenchuanEartliquakeJ.Science inChina SeriesETechnologi-cal Sciences，xx,51(S2)110-122.3Zhang G，Li Y，Li ZJ.A newapproach towardobject-basedchange detectionJ.Science ChinaTechnological Sciences，xx，53(S1):105-110.4Bovolo F，Bruzzone L，Marconcini M.A novelapproach tounsupervised change detection based onasemisiaperised SVMandasimilarity measiareJ.IEEE Transactionson Geoscienceand Remote Sensing，xx,46 (7)2070-2082.5庄会富，邓喀中，范洪冬.纹理特征向量与最大化熵法相结合的SAR影像非监督变化检测J.测绘学报，xx,45 (3):339-346.Zhuang HF，Deng KZ，Fan HD.SAR imagesunsupervised change detectionbased on binationof texturefeature vectorwith maximumentropy principleJ.ActaGeodaeticaetCarto-graphica Sinica，xx，45 (3)339-346.6李松，李亦秋，安裕伦.基于变化检测的滑坡灾害自动识别J?遥感信息，xx (1)27-31.Li S，Li YQ，An YL.Automatic recognitionof landslidesbasedonchange detectionJ.Remote sensingInformation，xx (1)27-31.7钟家强，王润生.基于自适应参数估计的多时相遥感图像变化检测J?测绘学报，xx,34 (4)331-336.Zhong JQ，Wang RS.Multitemporal remote sensing image change detectionbasedonadaptive parameterestimationJ.Acta Geo-daetica et Cartographica Sinica，xx，34 (4)331-336.8刘臻，宫鹏，史培军，等.基于相似度验证的自动变化探测研究J?遥感学报，xx,9 (5)537-543.Liia Z，Gong P，Shi PJ，et al.Study onchange detectionautomatically basedon similaritycalilDrationJ.Journal ofRemoteSensing，xx，9 (5)537-543.9李亮，舒宁，王凯，等.融合多特征的遥感影像变化检测方法J.测绘学报，xx,43 (9)945-953.Li L，Shu N，Wang K，et al.Change detectionmethod for remote sensing images basedon multi-features fusionJ?Acta Geodaeticaet CartographicaSinica，xx，43 (9)945-953.10魏立飞，王海波.基于MSE模型的高分辨率遥感图像变化检测J.光谱学与光谱分析，xx，33 (3)728-732.Wei LF，Wang HB.Change detectionfrom high-resolution remotesensingimage basedonMSE modelJ.Spectroscopy andSpectral Analysis，xx，33 (3)728-732.11巫兆聪，胡忠文，张谦，等.结合光谱、纹理与形状结构信息的遥感影像分割方法J?测绘学报，xx，42 (1)44-50.Wu ZC，Hu ZW，Zhang Q，et al.On biningspectral，textural andshape featuresforremotesensingimagesegmentationJ?Acta GeodaeticaetCartographicaSinica，xx，42 (1)44-50.12陈云浩，冯通，史培军，等.基于面向对象和规则的遥感影像分类研究J武汉大学学报(信息科学版），xx,31 (4)316-320.Chen YH，Feng T，Shi PJ，et al.Classification ofremot sensingimagebasedon objectoriented andclass rulesJ.Geomatics andInformation Scienceof WuhanUniversity，xx,31 (4)316-320.13窦建军，文俊，刘重庆?基于颜色直方图的图像检索技术J.红外与激光工程，xx，34 (1)84-88.Dou JJ，Wen J，Liia CQ.Histogram-based colorimage retrievalJ.Infrared andLaser Engineering，xx，34 (1)84-88.14Von GioiR G，Jakubowicz J，Morel JM，et al.LSDA linesegment第2期王光辉，等基于多特征融合的遥感影像变化检测算法?99?detectorJ.Image Processingon Line，xx，235-55.15韩华，曹伟，龚涛.目标再确认中基于推土机距离的关联度建立J?华中科技大学学报（自然科学版），xx,43(S1):435-439.Han H，Cao W，Gong T.The establishmentof correlativedegree intarget rr-identification basedon earth movers distanceJ.H ua-zhong Univerityof Scienceand TechnologyNatiaral ScienceEdition)，xx，43(si)435-439.16Wang AP，Wang SG，Lucieer A.Segmentation ofmultis-Dectral high-resoliation satelliteimagety basedon integratedfeatiare dis-trilDutionsJ.International JournalofRemoteSensing，xx，31 (6)1471-1483.Changedetectionbasedonadaptivefusionof multiplefeaturesWANG Guanghui1，LI Jianlei，WANG Huabin1，YANG Huachao2(1.Satelite Surveyingand MappingApplication，NASG，Beijing100830，China)2.China Universityof Miningand Technology(Xuz