专利名称：：一种高光谱遥感数据多类别监督分类方法
技术领域：
：本发明涉及一种高光谱遥感数据多类别监督分类方法，属于高光谱数据处理方法与应用
技术领域：
，适用于高光谱数据监督分类的理论方法和应用技术研咒。
背景技术：
：高光谱遥感数据监督分类的方法主要包括两大类基于光谱特征匹配的方法和基于统计分析模型的方法。由于高光谱数据获取过程中，受到大气、地形、光照等条件的影响，使得其获得的地物光谱特征变化较大，因此，基于光谱特征匹配的方法会出现不同地物之间混淆较大、不稳定等问题，使得分类结果精度降低；基于统计分析模型的方法主要是对高光谱数据样本的总体特征进行统计分析，根据样本采样点统计分布特征实现不同地物的分类。但是在统计分析模型建立过程中为了实现精确的分类，需要对样本进行分布假设以及大量的统计参数实现分类建模。对数(Logistic)两分类模型广泛的应用于医学和生物学，尤其在流行病学得到及其广泛的应用。1989年Hosmer和Lemeshow将Logistic回归模型的响应变量扩展到*(*>2),使得Logistic回归模型成功应用到*(*>2)类别的分类问题。近年来随着遥感技术的发展和数据处理的需求，Logisitc回归模型被引入到遥感图像分类应用中。与高斯分类模型相似，Logistic是监督分类模型，需要训练样本估计模型参数。但在Logistic回归模型中，仅仅需要对影响因变量的因素进行线性关系假设，不需要对图像中样本进行任何分布假设。通常情况下，需要对高光谱数据降维后才能实现Logistic建模，否则将带来庞大的计算负担，但是数据降维方法将在减少数据维数的同时丢失光谱的细节特征，从而导致相似类别的混淆。
发明内容本发明的目的在于克服现有技术的不足，提供一种不需要样本分布假设、参数估计少的高光谱遥感数据多类别监督分类方法。本发明的技术解决方案为一种利用分形理论实现较少训练样本条件下的高精度高光谱遥感数据多类别监督分类方法，该方法主要利用分形理论和统计分析模型，通过计算光谙的多重分形谱参数实现Logistic建模，最后利用最大似然估计模型参数估计，从而实现高光谱数据的监督分类。多重分形i普主要通过光谱概率测度计算和配分函数估计得到，并根据多重分形谱的特点选择a'、"mm、"max、/(j作为Logistic建模参数，从而在有效降低模型估计参数数目的同时增加类间可分性。本发明一种高光谱遥感数据多类别监督分类方法，其步骤如下(1)读入高光谱数据；(2)确定分类类别数目，并选择训练样本和测试样本；(3)计算多重分形谱参数；(4)建立基于多重分形谱参数的对数回归多类别分类模型；(5)利用极大似然估计进行模型求解；(6)利用概率最大原则进行分类，并计算分类精度。其中，步骤(1)中所述的读入高光i普数据为X-[x"x2,…，xJ,"为像元数。其中，步骤(2)中所述的分类类别数目为J，训练样本和测试样本选择根据参考图像和地面调查确定，且训练样本与测试样本不相同。其中，步骤(3)中所述的多重分形谱/(cO计算如下首先计算光谱概率测度S柳式中，iV(。表示光谱曲线划分为尺寸是5的一维小盒子的个数，《.(。表示盒子尺寸为5时第/个小盒子内所有光谱波段的辐亮度数值之和；然后估计配分函数，配分函数a(。为光谱概率测度f(。的g阶矩式中，《为权重因子，iV(。为尺度^下的盒子数目。再通过光谱概率测度计算与配分函数估计得到的尺度函数r(《)log2&(《)》,log2《+C(《)式中r(《)为尺度函数。最后，通过勒让德(Legendre)变换/()与r(《)联系起来，并得到如下关系式其中，《为权重因子。本发明建模过程中选择的多重分形谱参数为其中，《'为多重分形语/(")最大值点对应的Holder指数，"她和"^表示给定条件下概率奇异性的最小值和最大值，/("^)为光谱曲线上Holder指数取最小值的波段子集分形维数。其中，步骤(4)所述的基于参考类别NC建立的基于多重分形谱参数的对数回归多类别分类模型如下户尸WC其中，NC为参考类别，J为类别数目，j=l，2,...,J-l，S为输入变量属于第j个类别的概率，^、"min、"■、/K,J分别为多重分形谱参数，A,为模型待估计参数，j=l,2"."J-l,i=0,l,.."4。其中步骤(5)所述的利用极大似然方法进行模型参数估计过程如下对步骤(4)建立的模型进行极大似然参数估计，采用牛顿-拉臾逊(Newton-Raphson)迭代法解得M7为参考类别的Logistic模型参数的极大似然估计值为A=(U"，A2，U,4),7'=1,2,../-1。其中，步骤(6)所述的分类准则为概率最大原则，即z/f(x):maxCf;(x)，尸2(x),…A(;c))贝'JxeG,x属于第/个类别的概率最大，那么就将其归入类别G,;分类精度分析采用混淆矩阵和Kappa系数实现。本发明与现有技术相比的优点在于克服了传统的基于统计分析模型的高光语数据监督分类方法参数估计多、样本假设分布等局限，本方法利用多重分形理论和Logsitic模型，实现了高精度的地物类型分类。它具有以下的优点(1)采用统计分析模型，解决了由于光语特征匹配高光谱遥感数据监督分类方法受大气、光照、地形等外界影响较大带来的误差增多和可靠性下降等问题；(2)利用Logistic多类别模型进行高光谱数据分类，不需要对样本分布进行假设；(3)引入多重分形理论，建立基于多重分形参数的Logistic多类别分类模型，有效的利用了光谱曲线的差异和细节信息，并降低了模型估计参数数目、提高了类间可分性。图1为本发明中多重分形参数"'、"min、"max、/("，)的定义示意图；具体实施例方式为了更好的说明本发明涉及的高光谱遥感数据多类别监督分类方法，利用PHI航空高光i普数据进行江苏方麓茶场地区农作物精细分类。本发明一种高光谱遥感数据多类别监督分类方法，具体实现步骤如下(1)读入高光谱数据读入方麓茶场推扫高光谱成像仪(PHI)高光谱数据，数据大小为210x150x64,波段区间455~805nm;(2)确定分类类别数目，并选择训练样本根据参考图像，分类类别数目为」=6,训练样本和测试样本根据参考图像获得，具体的分类类别、训练样本与测试样本如下表所示；<table>tableseeoriginaldocumentpage7</column></row><table>(3)计算多重分形谱参数首先计算光谱概率测度<formula>formulaseeoriginaldocumentpage7</formula>式中，W(。表示光谱曲线划分为尺寸是S的一维小盒子的个数，S,(。表示盒子尺寸为5时第/个小盒子内所有光谗波段的辐亮度数值之和；然后估计配分函数，配分函数a(。为光谱概率测度《(。的^阶矩,=1式中，《为权重因子，A^)为尺度5下的盒子数目。再通过光谱概率测度计算与配分函数估计得到的尺度函数r(《)log2&(《)《,l。g2《+C(g)式中r(《)为尺度函数。最后，通过Legendre变换/(a)与r(《)联系起来，并得到如下关系式<formula>formulaseeoriginaldocumentpage7</formula>其中，《为权重因子。本发明建模过程中选择的多重分形谱参数为[">,，"，,/(",)]7，其中，"'为多重分形谱/(")最大值点对应的Holder指数，"随和"，表示给定条件下概率奇异性的最小值和最大值，/(,J为光谱曲线上Holder指数取最小值的波段子集分形维数，各参数定义如图l所示；(4)建立基于多重分形谱参数的Logistic多类别分类模型以水(W2)^作为参考类别，得到Logistic分类模型的判别函数如下lOg(^)=A。+A,"*+A2"min+A]"max+/^/(min)尸6尸6log(，)=A。++/52amm+A3"max+A4/("咖)(5)利用极大似然估计进行模型求解，得到的以水(W2)^作为参考类别，其余五个类别均具有一套独立的参数，得到Logistic分类模型的判别函数如下=譜".柳.-<5W.""m+S7.7脇_+"7.97氛.J《《/og^l")=-〃269+;7W9a'+277.2j。湖-+2W.卵/U(6)利用概率最大原则进行分类，并计算分类精度分类准则为概率最大原则，即,/S(x)-max(S(x),尸2(x),…A(x))则xeG,，x属于第/个类别的概率最大，那么就将其归入类别G,;分类精度分析采用混淆矩阵和Kappa系数实现，如下表所示分类精度。多重分形参数Logisitc分类模型对测试样本分类混淆矩阵8C4T6T7V2V13W2分类精度C47973000099.625%T616264000094.286%T72427301097.500%V200100240096.000%V130201970198.500%W200000譜100.00%Kappa系数0.981总体精度98.506%从利用本发明方法得到的分类结果可以看出总体分类精度达到98.506%,主要原因是经过多重分形参数的计算，类别间的差别较明显。水稻(C4)、茶(T7)、红薯(V2)、及水(W2)光谱可分性较高，基于多重分形参数的Logistic监督分类精度都在95。/。以上；对于较难区分的竹(T6)与香菜(V13)，通过提取光谱曲线的多重分形参数，有效的提高了地物间的可分性，建立的基于多重分形参数的Logistic监督分类方法对竹(T6)和香菜(V13)的分类精度分别达94%与98%;多重分形谱参数的kappa系数为0.981,分类结果的一致4交高。权利要求1、一种高光谱遥感数据多类别监督分类方法，其特征在于它包含以下步骤(1)读入高光谱数据；(2)确定分类类别数目，并选择训练样本和测试样本；(3)计算多重分形谱参数；(4)建立基于多重分形谱参数的对数回归多类别分类模型；(5)利用极大似然估计进行模型求解；(6)利用概率最大原则进行分类，并计算分类精度。2、根据权利要求1所述的一种高光谱遥感数据多类别监督分类方法，其中步骤(3)计算多重分形谱参数包括[">m,","_,/("m』7,其中，"'为多重分形i普/(")最大值点对应的Holder指数，"自和"^表示给定条件下概率奇异性的最小值和最大值，/("^)为光谱曲线上Holder指数取最小值的波段子集分形维数。3、根据权利要求1所述的一种高光谱遥感数据多类别监督分类方法，其中步骤(4)建立的基于多重分形谱参数的对数回归模型为<formula>formulaseeoriginaldocumentpage2</formula>其中，NC为参考类别，J为类别数目，j=l，2，...,J-l,巧为输入变量属于第j个类别的概率，a'、amin、《_、/(_)分别为多重分形谱参数，A,为模型待估计参数，j=l,2，...,J-l，i=0，l，...，4。全文摘要一种高光谱遥感数据多类别监督分类方法，包含以下步骤(1)读入高光谱数据；(2)确定类别数目，并选择训练样本和测试样本；(3)计算多重分形谱参数；(4)建立基于多重分形谱参数的对数回归多类别分类模型；(5)利用极大似然估计进行模型求解；(6)利用概率最大原则进行分类，并计算分类精度。该方法不需要对变量的概率分布进行任何假设，并且分类器模型中需要估计的参数数目较少，通过多重分形特征增加类内一致性的同时提高了类间的可分性。因而，该方法在训练样本的数量较小的情况下能够获得较高的分类精度。文档编号G06K9/66GK101667253SQ200910093539公开日2010年3月10日申请日期2009年9月25日优先权日2009年9月25日发明者娜李,牛志宇,贾国瑞,赵慧洁申请人:北京航空航天大学