专利名称：：一种基于信息量维数序列的高光谱数据监督分类方法
技术领域：
：本发明涉及一种基于信息量维数序列的高光谱数据监督分类方法，属于高光谱数据处理方法与应用
技术领域：
，适用于高光谱数据目标识别及监督分类的理论方法和应用技术研究。
背景技术：
：目标识别与地物分类是高光谱遥感数据应用的主要方向之一。高光谱图像提供了近乎连续光谱的方式实现对目标的精确描述，和其它遥感方式相比，高光谱遥感科学技术最大的优点体现在对地物的分类和识别能力上，该类技术的发展扩展了高光语数据应用的深度和广度。基于光i普匹配的分类方法，利用光谱库中已知光谱数据，采用匹配算法来识别图像中的地面覆盖类型，该类方法主要通过全波段特征和局部量化特征两种方式实现。全波段特征匹配方式利用了整个光i普的形状特征，受照度、光谱定标和光i普重建精度等的影响较小，但对光镨的微小差异不够敏感。局部量化特征匹配方式虽然对光谱微小差异敏感，但仅仅利用了特定的一些特征，受图像信噪比、光谱定标和光谱重建精度等因素的影响较大。分形体具有精细的结构，形状不规则并具有自相似性，分形理论在遥感图像解译中，常常被用来提取遥感影像的空间分形特征。研究表明，高光谱图像光谱维数据的具有分形特征，每一个像元可以得到一个分形维数，分形维数可以构成一幅新的图像，可以充分的把高光镨图像的光谱特征与图像特征结合起来进行分析，用于分类和识别。但是实际应用中的光谱曲线，只有在某种被限制的观测尺度的范围内，自相似性才成立，仅利用分形维数一个数字去描述光谱曲线所有波段包含的信息无法满足应用需要。信息量维数用来表示仅用单一分形维数不能完全描述的光i普局部分形特征，从光谱的局部出发来研究其最终的整体光谱特征。目前，基于光谱维的分形分析已经受到光谱识别和地物分类等研究领域的关注。
发明内容本发明的一种基于信息量维数序列的高光谱数据监督分类方法的目的及其所要解决的问题是将信息量维数引入到光谱域分析，应用于高光谱数据的目标识别及监督分类，综合全波段匹配和局部量化特征匹配的优点，克服采用单一匹配方式的或对微小差异不敏感，或受信噪比和光语重建精度等影响较大的不足，以获得较高的分类效率和分类精度。本发明一种基于信息量维数序列的高光谱数据监督分类方法，包括如下步骤(一)读入某地区高光谱图像数据；(二)根据先验信息从光谱库中选择该区域可能存在的地物光语作为参考光谱，或从图像中选取训练样本进行波段平均获取参考光语；(三)获得参考光谱后，逐个计算参考光谱及所有测试光谱的信息量维数序列；(四)参考光谱和测试光谱的信息量维数序列计算完毕，测试光语的信息量维数序列与所有参考光语的信息量维数序列逐个进行向量角匹配，采用最小距离分类器确定与参考光i普最接近的一个或多个像元；(五)所有测试像元匹配完成，二值化匹配结果，将每类地物的匹配结果用二值图像表示，每幅图像中只包含一类地物。其中，在步骤(一)中所述的"读入某地区高光谦图像数据"是指在遥感影像处理平台ENVI4.2上打开高光谱数据文件来读入某地区高光谱图像数据；其中，在步骤(二)中所述的"波段平均获取参考光谱"，其具体实现方式》口下<formula>formulaseeoriginaldocumentpage6</formula>(8)式(8)中，x,;^为训练样本，M为训练样本数目，7V5为波段数目，i为A^维参考光谱。i第Z个波段的值由x,~&第/个波段值平均得到。其中，在步骤(三)中所述的"逐个计算参考光镨及所有测试光语的信息量维数序列"，其具体实现方式如下先给定光i普曲线的分段数目w以决定信息量维数序列的维数，光谱曲线信息量维数序列采用如下方法计算将光谱曲线等分为W个小区域，每个区域包含的波段数目为層S-y離(，)(9)式(9)中W。wOO表示取小于x的最大整数。当W5不能被iV整除时，前yv-i个区域包含的波段数目为,最后一个区间包含波段数目为AB-(iV-l).MS。每个小区域的线度(小区域长度)大小为A,光傳曲线在该小区域的生长几率为S,不同小区域生长几率不同，用不同标度指数",来表征《=Z,'，/=1,2,3,-..,TV(10)若线度A的大小趋于零，用下式表示局部分维lnPa,=lim"^，/=l，2,3,"，iV(11)计算每个小区域上的局部分维，得到一个由不同",.所组成的序列构成的谱，用/(")表示。局部分维计算方法如下用线度^将线小区域丄,划分成一系列更小的区间，改变&的大小，得到一系列的生长几率A，对風),ln(《2)，…風))和風),ln(。，…,ln(^"U为不同的&的个数)用最小二乘拟合直线，斜率就为该小区域上的局部分维，该局部分维就是对应小区间上的信息量维数，所有小区域的信息量维数构成光谱曲线的信息量维数序列光谱曲线在小区域生长的几率《采用如下方法计算设波段数目为7\^的光谱曲线义=(^12，.."^),每个分量x,为每个波段分量，定义x的光谱信息度量"化rw，其中<formula>formulaseeoriginaldocumentpage8</formula>用线度为A的区域包含了光谱曲线从x,到;c,.Oy)的_/一+1个波段，则该区域中光谱曲线生长的几率为i^(tA).ln(^>》，,=1,2,…,iV(13)A为光谱曲线X的概率分布，TV为线度A下小区域的个数。其中，在步骤(四)中所述的"测试光谱的信息量维数序列与所有参考光谱的信息量维数序列逐个进行向量角匹配，采用最小距离分类器确定与参考光镨最接近的一个或多个像元"，其具体实现方式如下将信息量维数序列当作向量，测试光谱的信息量维数序列与所有参考光i普得信息量维数序列计算向量角，向量角计算公式如下)=arccos-"(14)A和A分别为参考光语和测试光谱的信息量维数序列，w为其维数，6>(Z)r,D,)为Dr和D,所成的向量角度。每一测试光谱的信息量维数序列与所有参考光谱的信息量维数序列计算得到一个向量角，与测试光谱的信息量维数序列所成角度最小的参考光语代表的类别为测试光i普所属类别，即如果0(Z^,A)-min讽Dn,D,)刀(Z^,A)，…e(A，Z),),…0(Z^，A)〉，则测试光语属于类别表示参考光语类别数目。其中，在步骤(五)中所述的"二值化匹配结果，将每类地物的匹配结果用二值图像表示"，其具体实现方式如下在匹配结果图像中，将每个感兴趣类别设为"1"，其他类别或背景设为"0"，得到最终分类结果图，每一幅分类结果中只包含一类地物，其它均为背景。本发明一种基于信息量维数序列的高光谱数据监督分类方法，其优点是综合了全波段匹配和局部量化特征匹配的优点，一定程度上克服釆用单一匹配方式的对微小差异不敏感、受信噪比和光谱重建精度等影响较大的不足，能获得较高的分类效率和分类精度；将信息量维数引入到光谱分析，为高光谱数据监督分类和目标识别提供了一种全新的途径。图1为华盛顿地区原始影像；图2为本发明的流程图3为参考光谱辐亮度曲线；图4为参考光谱的信息量维数序列，维数为5;图5为草坪的分类结果；图6为人工建筑的分类结果；图7为特殊类别屋顶的分类结果；图8为阴影的分类结果；图9为街道的分类结果；图10为树木的分类结果。具体实施例方式采用的华盛顿地区都市影像如图l所示，为飞机搭载高光谱仪低空拍摄而成。此影像具有220个波段，去除水气吸收波段后剩余191个波段，图像大小150xl95像素。图像中包含草碎，建筑，屋顶，阴影，街道和树木6个类别。本发明一种基于信息量维数序列的高光谱数据监督分类方法，其流程图如图2所示，以对华盛顿地区的都市影像进行监督分类为例，本发明所述的方法，其具体实现步骤如下(一)读入某地区高光谱图像数据；本发明的一种基于信息量维数序列的高光谱数据监督分类方法利用的是华盛顿地区高光谱数据，在遥感影像处理平台ENVI4.2打开，原始影像如图1所示，可以得到地物分布化息5(二)根据先验信息从光谘库中选择该区域可能存在的地物光谱作为参考光谱，或从图像中选取训练样本进行波段平均获取参考光谱；从图像中随机提取六类地物的训练光谱和测试光语，训练样本不包含在测试样本中；下列表1给出了用于分类的地物类别及样本个数表1实验样本个数<table>tableseeoriginaldocumentpage10</column></row><table>所有样本从图1所示图像中随机选取；对训练光谱按波段平均得到六个类别的参考光谱如图3所示，其中，波段平均获取参考光谱实现方式如下<formula>formulaseeoriginaldocumentpage10</formula>式(15)中，jc,j^为训练样本，M为训练样本数目，A^为波段数目，及为维参考光谱。/第Z个波段的值由x,~xM第/个波段值平均得到。(三)获得参考光谱后，逐个计算参考光谱及所有测试光谱的信息量维数序列；给定光谱分段数目iV=5，以决定信息量维数序列的维数；将光谱波段等分为5个区间，光谱曲线信息量维数序列采用如下方法计算将光谱曲线等分为5个小区域，每个区域包含的波段数目为層s=y,(逻)=y縱(^)(16)5式(16)中，y/o。K"表示取小于x的最大整数；当A^不能被w整除时，前w-1个区域包含的波段数目为iVAS,最后一个区间包含波段数目为A-(iV-l).ASS,因为该数据为191个波段，划分结果为前4个小区域波段数目为38，第5个区域波段数目为39;每个小区域的线度(小区域长度)大小为A，光谱曲线生长界面在该小区域的生长几率为f，不同小区域生长几率不同，用不同标度指数",来表征《=丄,、/=l,2，3,.'-,iV(17)若线度z,的大小趋于零，用下式表示局部分维",=lim^,/=l，2，3，-.-，iV(18)计算每个小区域上的局部分维，到一个由不同",所组成的序列构成的谱，用/(")表示。局部分维计算方法如下用线度&将线小区域Z,划分成一系列更小的区间，改变&的大小，得到一系列的生长几率几，对(lnOUln(《2)，…,ln(&)〉和{lri(。,ln(《2),..，ln(。}U为不同的&的个数)用最小二乘拟合直线，斜率就为该小区域上的局部分维",，该局部分维就是对应小区间上的信息量维数，所有小区域的信息量维数构成光谱曲线的信息量维数序列/(");光谱曲线在小区域生长的几率^采用如下方法计算设波段数目A^为191的光谱曲线X-(x,，;c2,...,^),每个分量x,为每个波段分量，定义x的光谱信息度量分布/^化r^，其中191A=V2>'，X2，…,iV5(19)用线度为丄,的区域包含了光谱曲线从;c,到、(y^')的y-Z+l个波段，则该区域中光谱曲线生长的几率为A为光谱曲线义的概率分布，iv为线度丄,下小区域的个数；利用上述方法，得到的参考光谱的信息量维数序列如图4所示；然后计算华盛顿地区高光谱图像中每个像元的信息量维数序列；(四)参考光谱和测试光谱的信息量维数序列计算完毕，测试光谱的信息量维数序列与所有参考光谱的信息量维数序列逐个进行向量角匹配，采用最小距离分类器确定与参考光谱最接近的一个或多个像元；将信息量维数序列当作向量，测试光谱的信息量维数序列与所有参考光镨得信息量维数序列计算向量角，向量角计算公式如下A和D,分别为参考光i普和测试光谱的信息量维数序列，为其维数，A,)为Dr和所成的向量角度。每一测试光谱的信息量维数序列与所有参考光谱的信息量维数序列计算得到一个向量角，与测试光谱的信息量维数序列所成角度最小的参考光谱代表的类别为测试光谱所属类别，即如果^^,D,^min讽Z^，D,)^(A2,D,)，…^(Z)。，Z),),…e(Z^，A)〉，则测试光谱属于类别y。^表示参考光谱类别数目。(五)所有测试像元匹配完成，二值化匹配结果，将每类地物的匹配结果用二值图像表示，每幅图像中只包含一类地物；在匹配结果图像中，将每个感兴趣类别设为"1",其他类别或背景设为"0",得到最终分类结果图，每一幅分类结果中只包含一类地物，其它均为背景。釆用本发明的一种基于信息量维数序列的高光谱数据监督分类方法得到的分类结果如图5~图10所示；不同地物以不同灰度表示，图6图10中各包含一类地物，白色表示地物，黑色表示背景；图5表示的草坪和图10的树木在光谱波形上比较相似，两类地物依然较好的被分开，草坪在图像中占大部分；图6的建筑实际上包含了图像中大部分的人工建筑，图7中屋顶单独作为一类也成功的和其它人工建筑分开；阴影大部分分布在建筑和树木周围，如图8所示；图9中街道受周围建筑影响，分类的效果要差于其它类别；从上面的分析可以看出，各个类别均较好的区分开；下列表2给出了利用本发明所述的方法进行分类的分类精度、混淆矩阵以及Kappa系数，总体分类正确率为92.51%,Kappa系数为0.903,分类效果良好，无明显错分类；所用类别中，建筑、屋顶和街道同属于人工建筑，因材料近似而出现了一定程度的错分。表2分类结果统计<table>tableseeoriginaldocumentpage13</column></row><table>权利要求1、一种基于信息量维数序列的高光谱数据监督分类方法，其特征在于它包含以下步骤(一)读入某地区高光谱图像数据；(二)根据先验信息从光谱库中选择该区域可能存在的地物光谱作为参考光谱，或从图像中选取训练样本进行波段平均获取参考光谱；(三)获得参考光谱后，逐个计算参考光谱及所有测试光谱的信息量维数序列；(四)参考光谱和测试光谱的信息量维数序列计算完毕，测试光谱的信息量维数序列与所有参考光谱的信息量维数序列逐个进行向量角匹配，采用最小距离分类器确定与参考光谱最接近的一个或多个像元；(五)所有测试像元匹配完成，二值化匹配结果，将每类地物的匹配结果用二值图像表示，每幅图像中只包含一类地物。2、根据权利要求1所述的一种基于信息量维数序列的高光谱数据监督分类方法，其特征在于步骤(一)中所述的"读入某地区高光谱图像数据，，，是指在遥感影像处理平台ENVI4.2上打开高光谱数据文件来读入某地区高光谱图像数据。3、根据权利要求1所述的一种基于信息量维数序列的高光语数据监督分类方法，其特征在于步骤(二)中所述的"波段平均获取参考光镨，，，其具体实现方式如下外)=匸2>"),/"，2,3,…,iV5(1)式(1)中，；c广XM为训练样本，M为训练样本数目，W5为波段数目，及为A^维参考光谱；i第/个波段的值由X,~&第,'个波段值平均得到。4、根据权利要求1所述的一种基于信息量维数序列的高光语数据监督分类方法，其特征在于步骤(三)中所述的"逐个计算参考光谱及所有测试光谱的信息量维数序列"，其具体实现方式如下先给定光谱曲线的分段数目w以决定信息量维数序列的维数，光谱曲线信息量维数序列采用如下方法计算将光谱曲线等分为W个小区域，每个区域包含的波段数目为顺=一《)(2)式(2)中y/oor(JC)表示取小于x的最大整数；当iVB不能被iV整除时，前7V-l个区域包含的波段数目为,最后一个区间包含波段数目为M-(AM>MS;每个小区域的线度(小区域长度)大小为A,光谱曲线在该小区域的生长几率为f，不同小区域生长几率不同，用不同标度指数",来表征f=Z,a',z=1,2,3,…，iV(3)若线度Z,的大小趋于零，用下式表示局部分维a,=lim^，z'=1,2,3,…，W(4)计算每个小区域上的局部分维，得到一个由不同^所组成的序列构成的谱，用/(")表示；局部分维计算方法如下用线度&将线小区域A划分成一系列更小的区间，改变《4的大小，得到一系列的生长几率凡，对風),ln(《2)，…風"和(ln(。,ln(。,...,,}U为不同的&的个数)用最小二乘拟合直线，斜率就为该小区域上的局部分维",.，该局部分维就是对应小区间上的信息量维数，所有小区域的信息量维数构成光i普曲线的信息量维数序列/(");光谱曲线在小区域生长的几率《釆用如下方法计算设波段数目为，的光谱曲线X-0c,A,...,^)，每个分量x,为每个波段分量，定义x的光谱信息度量"化r^，其中顺用线度为A的区域包含了光谱曲线从、到(7々0的y-"i个波段，则该区域中光谱曲线生长的几率为A为光谱曲线x的概率分布，w为线度丄,下小区域的个数。5、根据权利要求1所述的一种基于信息量维数序列的高光谱数据监督分类方法，其特征在于其步骤(四)中所述的"测试光镨的信息量维数序列与所有参考光谱的信息量维数序列逐个进行向量角匹配，采用最小距离分类器确定与参考光谱最接近的一个或多个像元"，其具体实现方式如下将信息量维数序列当作向量，测试光谱的信息量维数序列与所有参考光谱得信息量维数序列计算向量角，向量角计算公式如下A和A分别为参考光谱和测试光语的信息量维数序列，iV为其维数，代A,D,)为A_和Z),所成的向量角度；每一测试光谱的信息量维数序列与所有参考光谱的信息量维数序列计算得到一个向量角，与测试光语的信息量维数序列所成角度最小的参考光语代表的类别为测试光谱所属类别，即如果<formula>formulaseeoriginaldocumentpage4</formula>则测试光谱属于类别y;c"表示参考光语类别数目。6、根据权利要求1所述的一种基于信息量维数序列的高光语数据监督分类方法，其特征在于其步骤(五)中所述的"二值化匹配结果，将每类地物的匹配结果用二值图像表示"，其具体实现方式如下在匹配结果图像中，将每个感兴趣类别设为"1"，其他类别或背景设为"0"，得到最终分类结果图，每一幅分类结果中只包含一类地物，其它均为背景。f=(Ea)《a)，'=1,2,…,w全文摘要一种基于信息量维数序列的高光谱数据监督分类方法，包括如下步骤(一)读入某地区高光谱图像数据；(二)从光谱库中选择参考光谱，或从图像中选取训练样本进行波段平均获取参考光谱；(三)逐个计算参考光谱与所有测试光谱的信息量维数序列；(四)测试光谱与所有参考光谱的信息量维数序列逐个进行向量角匹配，采用最小距离分类器分类；(五)二值化匹配结果，将每类地物的匹配结果用二值图像表示，每幅图像中只包含一类地物。基于信息量维数序列的高光谱数据分类方法将信息量维数引入到光谱域分析，综合了全波段匹配和局部量化特征匹配的优点，能获得较高的分类效率和分类精度，在高光谱数据分类和目标识别中有重要价值。文档编号G06K9/64GK101299237SQ20081011443公开日2008年11月5日申请日期2008年6月5日优先权日2008年6月5日发明者刘小刚,娜李,赵慧洁申请人:北京航空航天大学