专利名称：一种高光谱图像的端元提取方法
技术领域：
本发明涉及高光谱图像分析领域，特别是涉及一种高光谱图像的 端元提取方法。
背景技术：
高光谱技术实现了 20世纪遥感技术的一次飞跃，是当今国际遥 感领域关注的焦点之一。高光谱图像是对某个特定目标进行多谱段成 像所得到的一组二维图像，其中每个谱段对应一幅二维图像，并在二 维图像的基础上多了一维光谱信息，具有图谱合一的特性，能够反映 被测对象的物理特性，用于复杂环境中的目标检测。高光谱图像在空 间遥感、军事侦察、农业及环境检测等领域，如环境成像、地质研究、 职务及矿物的鉴定以及农作物的分析领域都有着广泛应用。高光谱数 据的端元提取是理解高光谱数据，继而对数据进行进一步分析(如混 合像元分解、目标检测)等的前提条件，如何从高光谱图像中提取端 元成为一个热点问题。
基于顶点成分分析(Vertex Component Analysis,简称VCA)的 端元提取方法因为其计算复杂度低而得到广泛应用。VCA端元提取 方法的原理如下
假设像元的线性混合模型表达为r = wM/ + "，其中r代表/xl 维像元的光谱向量，/为波段数，7k^[附,，...，w，...，附J是端元矩 阵，m,表示第/个端元向量，《为所覆盖区域中的端元数；《是一个 尺度因子，用来建模由于表面地形起伏引起的亮度变化(o^O); A， A，…，A]t是羊度向量；"是加性噪声。
如图l所示，线性混合模型下，观测光谱向量集合组成一个凸锥 Cr = f e: r = ft^/ ,[l]T/ = 2 0,w2 o}。当w = 1时，该凸锥简化为一
3个单体，表示为^"xei ':x二MA[l]Xl],-20卜该单体的顶点对应 端元向量。凸锥O投影到恰当选择的超平面上的投影组成一个单体 &，该单体的顶点做适当投影就可以得到单一地物的光谱向量。
将高光谱图像转换为二维观测光谱矩阵/ =[ ...， ]，其
中，r,代表第/个像元的光谱向量，N为图像中包含的像元总数，则 端元的具体求解步骤如下
1、首先用奇异值分解(SVD)对观测光谱矩阵降维，变换到《 维子空间X^"V ，其中，％是由SVD的左变换矩阵的前《个向量组 成的矩阵。
2、 将X投影到超平面上得到单体&: m.v=[ZJ:,,/{[XT:,,},其 中，w = mean(X)， w是lx^的向量。
3、 初始给定一个方向/: /=[ (/-^L4#) H/[H (/-^4#)wl|]， 其中，J = [e |0|...|0]， ^是应为qgp矩阵，用于存储估计端元信号的 投影，e =
v.ll 。将上式存储起来即存储像元索引[/"^"],=匕
5、 求取下一个投影方向M]:,-[7b,将此投影方向带入/=[ (/
II (/-wl门，/为投影方向，爿为顶点矩阵，j# 为顶点矩阵的伪逆，w为与/同维度的向量。重复执行^=/￥和& = argmax^, [v]:JI的运算，每重复执行一次，随机产生一个向量
/正交于矩阵J的列张成的空间，并且y被投影到/上，就可以得到 对应于极值的像元位置。
6、 利用M=仏[XU血e得到端元的光谱曲线。
因此，上述VCA算法沿任意投影方向对样本集进行迭代投影便可得到全部端元，但须满足以下前提条件高光谱图像所构成的样本 空间中一定存在纯点，即端元； 一定存在某个空间，使得该样本空间
在其上的投影一定是一个凸集(凸集的顶点近似为纯点，即端元)；
投影的极值 一 定对应凸集的 一 个顶点(其佘投影顶点的获得都是在此 基础上迭代投影得到的)。
可见，利用VCA方法提取端元的核心是，首先选择初始投影方
向，确定初始顶点，再依据"新的投影方向必须与已知顶点垂直"的原
则迭代投影确定全部顶点。因此，VCA方法在端元提取中的精度主 要取决于投影方向。传统VCA算法中，当J已知时，给定w即可确 定新的投影方向。若令5 = /-^<#,则^ = 0,也就是说，S中的每一
个行向量都与j中的列向量(即顶点)垂直。于是，由万中的行向
量为基即可组成一个与已有顶点垂直的超平面。,二S^w即为W在该
超平面上的投影。根据立体几何知识， 一个斜向量对平面的垂直投影 大于它向平面上任何其他方向的投影。考虑一种特殊情况，当w的方 向与某样本的方向重合时，依据上述原则，该样本就具有了特殊优先 权，得到较大的投影值。这种情况下得到了一个顶点，而事实上，对 后面的迭代过程来说它可能不是最优的选择(比如它作为 一个顶点并 不显著，周围点的投影值与其相差甚微)，这将影响精度。
在利用VCA方法提取端元的过程中，通过产生一个零均值高斯 随机向量来生成初始投影向量，由于高斯向量选取的随机性使得最后 生成的端元集合不稳定而且精度不高，不具有普适性。
发明内容
为了克服现有VCA方法在提取端元的过程中所存在的缺陷，本 发明的目的是提供一种基于顶点成分分析的高光谱图像的端元提取 方法，该方法通过选取具有一般性的投影方向生成高稳定性和高精度 的端元。
为达到上述目的，本发明的技术方案提供一种高光谱图像的端元提取方法，所述方法包括如下步骤
Sl，获取高光谱图像并将其转换为二维光谱矩阵，对所述二维光
谱矩阵降维，并将降维后的矩阵投影到超平面得到单形体； S2，确定投影方向；
S3，将所述单形体投影到所述投影方向上得到投影结果；
S4，根据上述投影结果获取投影极值及其对应的像元位置，并建
立像元位置索引；
S5，利用上述投影极值及位置索引得到端元的光谱曲线；
S6，判断是否达到预定的端元个数，如果是，则结東；如果否，
转下一步；
S7，利用S4中得到的位置索引求取下一个投影方向，并重复步 骤S3-S5;
其特征在于，所述投影方向为观测光谱向量集合的中心方向。 其中，所述投影方向为/=[ (/-^4#) (V,-v2) ]/[ II (v,-v2) l门，/为投影方向，J为顶点矩阵，^为顶点矩阵的伪逆，
vh=max{4; sGS}, v"=min{4; ^ S}， S为观测光谱向量集合，观测光 谱向量^为^维，^表示5在第/个维度上的数值。
本发明所提供的高光谱图像的端元提取方法，在现有VCA方法 的基础之上作出改进。由于该方法中投影方向的选取和确定具有 一般 性，任何样本对于投影方向的选择都没有特殊优先权，因此保证了改 进后的方法能够达到更好的提取精度以及更好的稳定性。
图l是传统VCA端元提取方法的原理示意图； 图2是本发明高光谱图像的端元提取方法的流程示意图； 图3是本发明实施例高光谱图像的端元提取方法中作为端元生成 具有局部相关性的渐近图像；
图4是利用本发明实施例高光谱图像的端元提取方法进行端元提丰度值矩阵重建的图像；
图5a-5d是设定的端元曲线；
图6a-6d是利用本发明实施例高光谱图像的端元提取方法所提取 的端元曲线。
具体实施例方式
下面结合附图和实施例，对本发明的具体实施方式
作进一步详细 描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
为了保证VCA算法的稳定性和普适性，就要使投影方向的选取 和确定具有一般性，即任何样本对于w的选择都没有特殊优先权。为 此，选择凸集的中心方向为投影方向w。由于凸集的中心方向不容易 确定，因此釆用其它方法予以近似。考虑一个由样本集构成的凸集， 用一个超立方体将其包含在体内部。设样本集为S,样本s为d维， 以4表示每个维度上的数值，则超立方体可以用两个向量vl和v2的 组合来表示，其中vl产maxW.; seS}， v2z=min{4.; ^曰^}。于是，凸 集的中心方向w就可以表示为
vv = vl _ v2
将上式代入公式/=[ wll]，得到本发明
中选取的投影方向
<formula>formula see original document page 7</formula>如图2所示为本发明提供的 一种高光谱图像的端元提取方法，所 述方法包括如下步骤获取高光谱图像并将其转换为二维光谱矩阵， 对所述二维光谱矩阵降维，并将降维后的矩阵投影到超平面得到单形
体；确定投影方向；将所述单形体投影到所述投影方向上得到投影结 果；根据上述投影结果获取投影极值及其对应的像元位置，并建立像 元位置索引；利用上述投影极值及位置索引得到端元的光谱曲线；利 用位置索引求取下一个投影方向并进行在新的投影方向上投影、求取投影极值及其对应的像元位置以及端元的光谱曲线，直至所获得的端 元达到预定的端元数量。本发明中的投影方向选择为观测光谱向量集
合的中心方向。优选地，所述投影方向选择为/=[ (w v2)]/[ II (/-^4#) (v广v2) l门，/为投影方向,^为顶点矩阵，爿# 为顶点矩阵的伪逆，v!产maxW; ^化v2,.=min{《；^￡5"}， S为观测光 谱向量集合，观测光谱向量^为"维，d,表示s在第/个维度上的数 值。
以下用模拟实验数据和真实实验数据对本发明的端元提取方法 进行了实验测试与分析，并与传统VCA算法及改进算法ATGP-VCA 进行了比较。
由于所得的顶点(端元)主要应用于图像重建，所以我们釆用原 始图像与重建图像间的SED值作为评判标准。对于一幅图像来说， 重建图像与原图像的SED值越小，则重建效果越好。
1)正确性验证
用四种物质的反射系数的部分波段，作为端元生成具有局部相关 性的渐进图像，如图3所示。通过端元提取后，应用端元及丰度值矩 阵重建的图像如图4所示。设定的端元曲线如图5a-5c所示，利用本 发明提取的端元曲线如图6a-6c所示。由上面的图像对比结果可以看 出，该方法能够正确地提取端元并重建图像。
以下定量分析验证本发明端元提取方法的正确性。计算提取的端 元与原始设定的端元之间的SAM参数，可得到SAM参数矩阵如下
0.6280 0 0.7167 0.9703
0 + 0細0i 0.6280 0.2803 0.5562
0.2803 0.7167 0 0.3745
0.5562 0.9703 0.3745 0
SAM参数的意义是两个向量之间的偏角，当两个向量越接近， 偏角越小。由以上上矩阵可以看出，图形上相似的一对端元的SAM 值接近O。
82)精确度验证
对同样的模拟数据，以SED为评价标准对其进行评价。如表1 所示，可以看出，原始的VCA每次运行都会得到不同的结果，并且 平均SED数量级在l(T16左右。而本发明所提供的方法可以具有运行 稳定的特质，并且将数量级降低到10-21 。使用ATGP-VCA能够产生 更好的结果，而将本发明中对w的选择方法应用到ATGP-VCA中产 生了更好的结果，SED值与ATGP-VCA同级并且略小。而且相比 ATGP-VCA的方法，本发明所提供的方法计算量更小，更易实现。
表1模拟图像的SED评价表
方法原始VCA选定w为中心方 向的改进ATGP-VCA选定w为中心方向的 ATGP-VCA
SED2.7560e-0161.0375e-0216.6158e-0233.2323e-023
值1.8498e-017L0375e匿0216.6158e-0233.2323e-023
4.2285e-0161.0375e-0216.6158e-0233.2323e-023
对模拟数据的实验说明了问题，但不能说明对真实的数据也能产 生同样的效果，因此本发明还对真实数据也进行了精度的检测。真实
数据采用ENVI软件自带的Cuprite地区图像文件名为Cup95eff,对 其提取ll个端元。依据SED的评价结果如表2所示。
表2真实图像的SED评价表
方法原始VCAw为中心方向的VCAATGP-VCAw为中心方向的ATGP-VCA
SED值1.2296e+0081.1766e+0081.5867e+0081.3874e+008
1.7141e+008U766e+0081.5867e+0081.3874e+008
1.3055e+008U766e+0081.5867e+0081.3874e+008
可以看到选择w为中心方向的改进对精度还是有略微提高的，但
后两种情况的结果并不是很理想。当然仍然可以看到对ATGP-VCA的 改进相对于ATGP-VCA还是有提高的。 一组数据说明，当存在噪声或 者实际端元并非纯线性组合时，改进方法并没有起到大的效果，这种情况下的算法值得进一步研究。
由以上实施例可以看出，本发明实施例通过釆用选取和确定具有 一般性的投影方向，并对该实施例中的方法与传统的高光谱图像基于 顶点成分分析的端元提取方法进行正确性和精确度的定量和定性比 较，证实了本发明提供的端元提取方法能够达到更好的提取精度以及 更好的稳定性。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领 域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以 做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
权利要求
1、一种高光谱图像的端元提取方法，所述方法包括如下步骤S1，获取高光谱图像并将其转换为二维光谱矩阵，对所述二维光谱矩阵降维，并将降维后的矩阵投影到超平面得到单形体；S2，选取投影方向；S3，将所述单形体投影到所述投影方向上得到投影结果；S4，根据上述投影结果获取投影极值及其对应的像元位置，并建立像元位置索引；S5，利用上述投影极值以及位置索引得到端元的光谱曲线；S6，判断是否达到预定的端元个数，如果是，则结束；如果否，转下一步；S7，利用S4中得到的位置索引求取下一个投影方向，并重复步骤S3-S5；其特征在于，所述投影方向为观测光谱向量集合的中心方向。
2、 如权利要求1所述的高光谱图像的端元提取方法，其特征在 于，所述投影方向为/= [(/_ A4# X " - v2)]/[ II ( /- A4弁X Vl - v2) II ]， /为投影方向，^为顶点矩阵，#为顶点矩阵的伪逆，"尸maxW.; ^ e 5},V2,:min^;^e^,S为观测光谱向量集合，观测光谱向量^为d维，《 表示^在第/个维度上的数值。
全文摘要
本发明公开了一种高光谱图像的端元提取方法，所述方法包括将降维后的高光谱图像光谱矩阵投影到超平面得到单形体；选择观测光谱向量集合的中心方向为投影方向；将所述单形体投影到所述投影方向上；根据投影结果获取投影极值及其对应的像元位置，并建立像元位置索引；利用上述投影极值及位置索引得到端元的光谱曲线；利用位置索引求取下一个投影方向并进行在新的投影方向上投影、求取投影极值及其对应的像元位置以及端元的光谱曲线，直至所获得的端元达到预定的端元数量。本发明所提供的方法能够达到更好的提取精度以及更好的稳定性。
文档编号G06T7/00GK101540049SQ200910082780
公开日2009年9月23日 申请日期2009年4月29日 优先权日2009年4月29日
发明者珺 何, 余乐军, 波 孙, 张熙予, 融 肖, 陶文靖 申请人:北京师范大学