一种基于非负独立成分分析的布谷鸟搜索高光谱解混方法
【技术领域】
[0001] 本发明涉及高光谱图像的预处理手段领域，更具体的说，是涉及一种结合非负独 立成分分析方法和布谷鸟搜索优化方法的无监督的高光谱图像解混方法。
【背景技术】
[0002] 光谱解混作为高光谱图像的一种预处理手段，不仅是实现地物精确分类和识别的 重要前提，而且是遥感技术向定量化深入发展的重要条件，它对发挥高光谱图像，在社会经 济和国防建设方面的效益有着重要的现实意义。高光谱图像是由高光谱遥感成像仪捕获的 包含空间和光谱信息的图像，能够帮助研究者们更好地提取地物信息。尽管高光谱图像的 光谱分辨率很高，但因自然界地物复杂多样以及遥感成像仪空间分辨率有限，单像元点往 往包含多种不同类别的地物，意味着"混合像元"广泛存在，使得传统的基于像元级光谱统 计特性的分类方法无法直接使用，这不仅制约着地物识别的精度，而且成为高光谱遥感向 定量化方向深入发展的主要障碍。因此，如何有效地解决混合像元问题，使遥感应用由像元 级达到亚像元级，实现"光谱解混"，是高光谱图像处理技术面临的难题之一。"光谱解混"过 程就是将混合像元分解为一些主要成分(端元)和每个成分所占比例(丰度)的过程。
[0003] 线性混合模型和非线性混合模型是高光谱解混方法中的两个基本模型，相比于非 线性混合模型，线性混合模型因其模型简单且物理意义明确而应用地更为广泛。高光谱线 性混合模型忽略不同端元间的多次散射，把每个像元近似为不同端元光谱以其丰度为权 重的线型组合。考虑到，高光谱数据为L个波段上的N个不同的像元，则某单像元点的观测值 X RLX1可以表示为
[0004] x=As+e (1)
[0005] 式中，义=|^1，12丄，11^]护><1儿是波段数;4=[31，32丄，3卩]护 ><[5是端元光谱矩阵，已1 Κ?χ1代表第i个端元光谱，P是端元数;s = [Si，S2，L，sp]T Rpxl是每个像元点的丰度向量;e R1 X1是误差项。其中，丰度向量s对应于每个端元所占比例，丰度必须满足丰度非负约束与丰 度和为一约束，即
[0007] 对于高光谱图像的N个不同像元所有观测向量组成观测矩阵X RLXP，可以表示为
[0008] X=AS+E (3)
[0009] 式中，S RPXN是N个像元的丰度向量组成的丰度矩阵，E对应得误差矩阵。现存的很 多方法解混方法都是基于线性混合模型的，主要有三大类:基于几何原理的方法、基于稀疏 回归的方法和基于统计方法的方法。基于几何方法的解混方法，如N-FINDR，顶点成分分析 (vertex component analysis，VCA)，最小体积单形体分析（Minimum Volume Simplex Analysis，MVSA)等，以观测数据中存在"纯端元"为前提，把高光谱数据集构成单形体的顶 点作为端元，往往实际地物分布并不能保证所有端元都有"纯"像元存在，故处理实际遥感 数据时效果不理想。
[0010] 基于稀疏回归的方法，假设观测像元能够用一组已知的"纯"的端元光谱线性表示 (光谱库作为已知端元），与整个光谱库的维数先比，每个像元点中的端元数很小，因此可以 通过寻找该光谱库中能够表示混合像元的最佳光谱子集实现解混。然而，该解混结果很大 程度上依赖合适光谱库的选取，很难获取与观测数据集相同条件下的光谱库，这是一个耗 时且高代价的过程。
[0011] 基于统计方法的解混方法是一种典型的无监督的方法，在不需要先验知识和人工 干预的条件下，实现端元光谱的自动提取。实际地物的分布情况往往非常复杂，导致端元的 无监督自动提取成为当前研究的一个难点。值得一提的是，随着盲信号分离技术的兴起和 发展，为如何用无监督的方法更快更精确地实现高光谱盲解混提供了另外一条可能途径， 己经成为当今遥感领域里的一个研究热点。独立分分析方法作为一种典型性的盲分离算 法，未用到光谱的先验信息，对观测像元的纯度不做要求，具有很高的研究价值。
【发明内容】
[0012] 本发明的目的在于提供基一种基于非负独立成分分析的布谷鸟搜索高光谱解混 方法。
[0013] 本发明一种结合非负独立成分分析方法和布谷鸟搜索优化方法的高光谱图像解 混方法，具体包括以下步骤：
[0014] 步骤一、对观测数据进行白化预处理，去除二阶相关性，对于高光谱图像解混，采 用利用主成分分析方法(Princ ipa 1 Components Ana 1 ys i s，PCA)对观测数据进行白化降维 预处理，具体步骤如下：
[0015] (1)对原始观测数据不进行去均值处理，保证估计出的丰度源是非负的；（2)先求 原始数据的自相关矩阵Corr = XXT/N，确保丰度向量间使相关的；（3)采用PCA计算出Corr的 前P个较大特征值组成的对角矩阵D p和对应的特征向量矩阵E p，利用公式
求出白化矩阵X'=VX RPXN，从而使得白化后的数据同时得到降维；
[0016] 步骤二：
[0017] 在非负独立成分分析方法的基础上，引入丰度和为一约束项构建目标函数为
[0019] 式中，炉为权重；
，为非负独立成分分析方法的目标 函数，其中，Y_=min(Y，0)是矩阵Y的负的部分；
，为丰度和为 一约束的目标函数，其中，λ是一个较大的正数，1RPX1的全一列矢量；
[0020] 利用基于Gives矩阵的QR分解理论，降低解混矩阵维数，将解混矩阵W转化为Gives 矩阵乘积的形式如下：
[0021] W=TpiLT2Ti
[0022] (2)
[0023] Ti = TipTi,piLTi2,T2 = T2pT2,piLT23,Tpi = Tpi,p
[0024] 式中，Tij是P阶旋转矩阵，旋转角度θ[0，231];
[0025] 采用布谷鸟搜索算法进行优化的迭代公式为
[0026]
[0027] 式中，Μ是参数维数，jiwsf表示第k个巢在第t代的位置，表示点对点乘法，α〇取值 0.01，是全局最优巢位置，Levy(step)~I step 11β，0〈β2表示随机步长step服从 Levy分布；
[0028] 采用布谷鸟搜索算法对目标函数进行优化具体步骤如下：
[0029] (1)初始化种群，设置巢个数为M，搜索空间维数为i/ = Cf，随机初始化巢位置
计算每个巢的目标函数的适度值，记录当前的最优适度 值fmin和最优位置 ，best{1，2，L，M};
[0030] (2)保留前一代最优巢位其中t[l，Tmax]，并按位置更新公式（10)对鸟 巢位置进行更新，得到当前巢位置
，与上一代巢的位置
'进行比较，将适应值较好的巢位置作为当前最好位置，
[0031] (3)每个巢分配一个服从均匀分布的0-1的随机数rd与卵被宿主发现的概率pJi 行比较。若满足rd pa，则将当前最好位置gt保存;如果rd〈pa成立，则随机改变满足该条件的 部分巢位置得到一组新巢，并与当前最优位置gt比较，将适应值较好的位置取代较差的位 置，得到一组新的最优位置
[0032] (4)计算（3)中得到的ggt中最优巢位置和相应的最优适度值fmin，并进行 判断。若满足迭代终止条件，则输出全局最优值和对应的全局最优位置，反之，则返回步骤 二继续循环更新;一般情况下，迭代终止条件为:达到指定的最大迭代次数或规定的精度要 求;本发明中采前者作为终止条件；
[0033]步骤三、利用步骤二中得到的参数最优值计算解混矩阵，估计丰度矩阵；
[0034] 步骤四、利用非负最小二乘方法估计出端元光谱。
[0035] 本发明针对高光谱图像中"混合像元"广泛存在、线性混合模型中丰度非负约束和 丰度和为一约束的特点，在传统的非负独立成分分析方法的基础上，加入光谱线性混合模 型中的丰度和为一约束项，构成有约束的非负独立成分分析方法，引入用布谷鸟群智能搜 索优化方法，利用其全局随机搜索性对目标函数进行优化，避免陷入局部极值，且对初始化 选择没有严格要求;在优化过程中，为了减少参数维数并缩小布谷鸟搜索方法群体搜索范 围，利用矩阵的QR分解理论，将对分离矩阵的搜索转化为对一系列Gives矩阵的识别，从而 减小了计算量。仿真数据和真实高光谱数据的实验结果表明，本发明能有效地克服传统独 立成分分析方法在解决高光谱解混问题时的局限性，获取良好的精度。
【附图说明】
[0036] 图1是本发明的基于非负独立成分分析的布谷鸟搜索高光谱解混方法的流程图。
【具体实施方式】
[0037] 下面结合附图对本发明做进一步说明。
[0038] 有约束的非负独立成分分析方法
[0039] 独立成分分析方法是一个无监督的盲源分离方法，其目标是在仅给出观测数据X 的情况下，根据独立性测度，通过寻找一种线性变换分离出源，即
[0040] Y=ffX = US (4)
[0041] 式中，Y=[yi，y2，L，yL]T RLXN是源S的估计，W是通过最大化Y成分间的独立性得到 的正交分离矩阵。通常情况下，独立性测度采用峭度、负熵和互信息等。
[0042] 非负独立成分分析方法将传统独立成分分析方法对独立性要求与源的非负性结 合起来进行分离，使分离结果尽可能的独立且是非负，通过解相关和非负旋转结合就可以 实现，为了保证估计信号非负，通过最小化如下目标函数实现：
[0044] 式中，Y_ = min(Y，0)是矩阵Y的负的部分，若Y是非负的，Jnka(Y)接近于零；否则， Jnica(Y)是一个较大的正数。
[0045] 线性混合模型中的丰度和为一约束条件要求，对于每一个观测像元，其各个端元 丰度之和必须为一，其目标函数可以定义为
[0047]式中，λ是一个较大的正数，1RPX1的全一列矢量，当丰度矩阵中每个像元的丰度都 满足和为一约束时，JASe(Y)接近于零，否则，JASe(Y)是一个较大的正数。<