基于差分搜索的高光谱图像解混方法
【技术领域】
[0001] 本发明设及高光谱遥感领域，更具体的说，是设及一种基于差分捜索的高光谱图 像解混方法。
【背景技术】
[0002] 遥感图像经历了全色（黑白）、彩色、多光谱成像阶段，到20世纪80年代光谱成像仪 出现后，进入高光谱遥感阶段。成像光谱仪通过对物质表面反射或发射不同波长的电磁波 来获取影像，将地物中所含物质的光谱与反映地物存在格局的空间影响结合在一起，即把 遥感影像和光谱合为一体，实现谱像一体化。同时记录二者信息，W像元为单位，在各个波 段分别成像，在二维空间的基础上增加一维光谱信息，把运些图像叠合在一起，形成高光谱 图像立方体，每一个像元对应着一条光谱曲线。高光谱图像中每个像元中包含地物类型的 种类是不同的，若像元中只包含一种地物类型，则称为纯像元;若像元中包含多种地物类 型，则称为混合像元。混合像元在高光谱图像中更广泛的存在，成为高光谱图像研究的主要 障碍。由于混合像元的存在，传统的基于像元级的方法难W达到实际应用的精度。所W，对 混合像元的分析，必须从像元级深入到亚像元级，进入到像元内部，从传统的像元解混发展 到亚像元解混。基于亚像元级的高光谱图像解混，能够提高解混的精度。线性混合模型是假 设物体之间没有相互作用，认为每个像元是图像中各个端元的线性混合。线性混合模型具 有建模简单、效率高、物理含义明确的优点，是目前使用广泛的解混模型。
[0003] 在线性混合模型中，高光谱图像的丰度必须满足"非负性"和"和为一性"。目前许 多方法都是先从高光谱图像中提取端元信息，然后反演求得丰度信息。如VCA(Vertex Component Analysis顶点成分分析）、N-FINDR(Fast Autonomous Spectral Endmember Determination In Hyperspectral Data高光谱数据中的快速自主谱端元确定）和MVSA (Minimum Volume Simplex Analysis最小单形体分析)方法。它们基于凸面几何理论，把高 光谱数据映射到一个高维平面，在该平面中通过寻找最大单形体来实现端元提取。运些方 法都假设高光谱图像中含有各个地物的"纯像元"，但运样的假设并不符合实际情况。
【发明内容】
[0004] 本发明是在盲信号分离框架下解决高光谱图像解混问题。在盲源分离理论中，可 W根据源信号之间的独立性实现分离，称为独立成分分析（Independent Component Analysis, ICA)。但已有方法证明不能将ICA直接应用于高光谱图像解混，因为高光谱图像 的丰度具有ANC和ASC特性，意味着丰度之间不相互独立。
[000引本发明在ICA处理过程中添加相应的约束项，使其能够应用于高光谱图像解混。一 般情况下为了简化估计过程，在ICA中都假设独立成分的数量与观测到混合信号的数量相 等。在高光谱问题中得到观测图像数量往往都大于端元数量，所W用ICA处理高光谱图像解 混问题的首要任务是对观测图像进行降维，即把X由L维降到P维。目前针对高光谱图像降维 大都采用主成分分析(PCA)方法，通过特征值比较来选取P个成分。降维后的r，即为ICA的 观测信号，丰度矩阵的估计Υ计算方法如下
[0006] Y=WX' (1)
[0007] 其中，W是解混矩阵。
[0008] 传统的方法采用梯度下降法求解混矩阵W，但梯度下降法容易陷入局部极值，而不 能得到有效的解混矩阵W。针对梯度下降法的不足，近年来许多学者采用仿生智能优化方法 来进行改进。该种方法能对解空间进行大范围的捜索，通过不断进行位置的更新来获得全 局最优解。本发明采用DS优化方法来得到高光谱图像中的解混矩阵W。
[0009] 本发明的具体过程为ICA独立性测度的选取、丰度约束项的构造和用DS(差分捜 索，Differential Search)方法优化求解。
[0010] 1 ICA独立性测度
[0011] 衡量ICA的测度有多种，如峭度、负赌、高阶累积量和互信息等。本发明采用互信息 作为独立性的度量，其计算方法为
[001 引
（2)
[0013]其中，KY)取值范围为[0，w)，当为0时表示变量之间相互独立。高光谱图像的丰 度在凸集合中通常是分散的，并且它们之间具有弱独立性，因此互信息可W很好地表示丰 度之间的特性。本发明通过最小化互信息得到弱独立性的丰度信息，在方法过程中Ι(γ)始 终大于0,运证明所估计的丰度不是完全独立的。yi的概率密度函数未知，所W不能直接得 ?jH(yi)。本发明采用Gram-Charlier展开式来逼近H(yi)，其展开式为
[0016] 计算式独立性之前需要对高光谱图像进行白化，使期望值为0,协方差为1。白化是 一个预处理过程，在ICA中至关重要，能够对观测信号进行去相关。如果没有限制，可能所有 的变量都能使互信息的值最小，只有在满足E{yyT}二I条件下，最小化互信息才能保证分离 结果的独立性。
[0017] 2丰度约束项
[0018] ANC和ASC特性使丰度之间不相互独立，仅仅依靠互信息是不能从高光谱图像中提 取出有效的丰度信息。本发明根据丰度的特性，在ICA中添加相应的约束项来使最后提取出 的丰度满足要求。
[0019] 各端元在像元中所占的比例值都是非负。基于此，在线性混合模型中，定义丰度非 负约束的目标函数为
[0020]
(4)
[0021] 其中，η为正整数，取n=l; 乂;;表示丰度估计矩阵y中所有数值为负的元素。当y中 所有元素非负，则fANC(Y)=〇。所W通过最小化fANC(Y)，实现丰度非负特性。
[0022] 高光谱图像中，各端元在每一像素点的丰度值之和必须为一。基于此，在线性混合 模型中，定义丰度和为一约束的目标函数为
[0023]
(5)
[0024] 其中，m为正整数，取m= 1。当丰度满足ASC约束，则fAsc(Y) =0，否则fAsc(Y)〉0。通过 最小化fAse(Y)，实现丰度和为一特性。
[002引结合式(2)~(5)得到高光谱图像解混的目标函数为
[0026] min f(Y) = I(Y)+yifANc(Y)+y2fAsc(Y) (6)
[0027] 其中，μι、μ2为两个约束所占的比例系数。
[002引3差分捜索方法
[0029] 本发明用基于DS优化方法的ICA(DS-ICA)来实现高光谱图像解混。差分捜索方法 是±耳其学者Pinar Civicioglu提出的一种仿生智能优化方法，与W往的智能优化方法不 同，它是模仿布朗提出的生物体在迁徙过程中的随机运动而提出的。有些生物体随周期性 气候变化会发生相应的迁徙行为，它们从一个容量和食物匿乏的栖息地迁徙到其它有利的 栖息地，当大量个体同时进行迁徙的时候，运些个体就形成了一个超级生物体，超级生物体 不断更新自己当前的位置，使其向着更加肥沃的区域迁徙。当超级生物体中某个个体探索 到新的肥沃区域，它就停留在运一段时间，简称停歇地(个体寻找各自停歇地是根据布朗随 机运动来实现），然后向着更加肥沃的区域继续探索。其它个体也会得知该停歇地的位置， 并且向着该方向进行迁徙。在运个过程中超级生物体总是向着全局最优位置进行迁徙，从 而能够获得全局最优解。DS方法的具体过程如下。
[0030] 1)生物个体位置的初始化定义为 [0031 ] hij = rand(upj-lowj)+lowj (7)
[0032] 其中，}1。表示个体的初始位置，^' = {1，2,3，-，，0}，0表示所优化问题的维数，1 = α，2,3，···，Τ}，Τ表示个体的总数;啡神日lowj分别表示第j维个体位置的上限值和下限值。将 初始化后的个体代入到相应的适应度函数，得到个体初始的适应度值。
[0033] 2)迁徙过程中依据下式寻找超级生物体的停歇地，根据布朗随机运动来寻找停歇 地，能有效避免陷入局部最优，使超级生物体向全局最优的区域迁徙。
[0034] Site = F*os+RXmapX (dinor-Pos) (8)
[003