一种高光谱遥感影像端元提取方法
【专利摘要】本发明公开了一种高光谱遥感影像端元提取方法，针对纯净像元假设不一定成立的情况，该方法将线性混合模型与粒子群优化算法相结合，配合高维粒子编码方式，推导出新的有效的粒子进化方程。同时，实验证明该方法针对端元提取问题更加有效。
【专利说明】
一种高光谱遥感影像端元提取方法
技术领域
[0001] 本发明属于高光谱遥感影像处理技术领域，涉及一种高光谱遥感影像端元提取方 法，尤其涉及一种线性混合模型约束的粒子群优化的高光谱遥感影像端元提取方法。
【背景技术】
[0002] 高光谱遥感影像具有光谱分辨率高，图谱合一的特性。每个像元可以提供几乎连 续的地物光谱曲线，使利用高光谱遥感反演陆地细节成为可能。高光谱遥感技术已经成为 当前遥感领域的前沿技术之一。然而，由于高光谱传感器对每一个像元提供很窄的成像波 段，这一波段具有较低的能量，使得只有扩大瞬时视场面积，传感器才能获取到有效信息。 因此，高光谱传感器一般具有较低空间分辨率，在进行对地成像时，受自然界地物分布复杂 多样性的影响，像元很难由单一均匀的地表覆盖类型组成，因此高光谱遥感影像中存在大 量的混合像元。
[0003] 混合像元分解技术是解决混合像元问题的有效途径，该方法主要包括端元提取和 丰度估计两个核心任务。其中，端元提取是丰度估计的前提条件。针对混合像元端元提取问 题，国内外学者提出了很多方法，传统的端元提取算法基于纯净像元存在假设。该类方法认 为高光谱图像中每类地物至少存在一个纯净的像元。经典的端元提取算法包括:纯像元指 数、N-FINDR、正交子空间投影、顶点成分分析等。
[0004] 然而，由于成像条件的影响以及地物的复杂性，纯净像元存在的假设往往不成立。 那么传统的基于纯净像元假设的算法就失效了，现有的一些端元生成算法发展还不成熟。 因此，需要发展一种有效的端元提取算法，适用于纯净像元不存在的情况。
【发明内容】
[0005] 为了解决上述技术问题，本发明提供了一种线性混合模型约束的粒子群优化算法 用于高光谱遥感影像端元提取方法。
[0006] 本发明所采用的技术方案为:一种高光谱遥感影像端元提取方法，其特征在于:将 高光谱遥感影像分别读入大小为I Xη的矩阵R= {ri,...，rj,...，rn}中，矩阵中每列向量rj = (rij，. . .，rij，. . .，Γυ)τ为像素的光谱辐射值，1为高光谱遥感影像的波段数，η为遥感影 像的像素数;对高光谱遥感影像矩阵R进行如下操作：
[0007] 步骤 1:初始化种群，Xi(t) = (ei,i(t)，ei,2(t)，· · ·，ei,D(t))i = l，2, · · ·，Μ;其中，t 表示迭代的次数，初始化过程中令I = O5X1U)表示种群中的一个粒子;M表示该种群中粒子 的数目；eij(0)eR，表示初始化过程中第i个粒子中的第j个向量;基于纯净像元存在假设的 顶点成分分析VCA(Vertex Component Analysis)算法提取的结果作为初始的种群；
[0008] 步骤2:计算每初始粒子的个体最优位置P1(O)及所有粒子的全局最优位置G(O);
[0009] Pi(O)=Xi(O)
[0010] (?：(0) = ai-gmiri/(^(〇))
[0011]
[0012] 其中，f( ·)表示适应度函数;VOl(X1U))表示求解X1U)所构成的单形体的体积； 丰度aid (t)满，
f表示两个约束之间的比例系数。
[0013] 步骤3:粒子进化;粒子进化方程为：
[0014] Xi,j(t+l) = (pi,j(t)-aj(t) · (pi(t)*Pj(t)))/(l-aj(t))
[0015] ρ,Χ?^φΧ?γΡ,,Χ? + ^-φ^-Ο^?) 9j{t)~U{Q,\)
[0016] 其中，aj(t)为一个随机生成数，表示图像的混合程度;&(t)是一个随机的丰度向 量，需要满足的条件是:该向量中的第j个元素是〇,且满足和为一条件;A⑴表示的是一个〇 ~1的随机生成数；
[0017] 步骤4:更新粒子的个体最优位置及所有粒子的全局最优位置；
[0018] 粒子个体最优位置：
[0019]
[0020] 其中，{xKthx^U+iM表示的是X1U)粒子的第j个列向量已经更新;f( x)表示适 应度函数；
[0021] 所有粒子的全局最优位置：
[0022] G(t+1) = argmm/(/^(ti-l)).
[0023] 步骤5:判断是否达到循环终止条件；
[0024]如果否，重新执行第3步；
[0025] 如果是，则循环终止，输出最终的结果G(t)。
[0026] 作为优选，步骤1中所述基于纯净像元存在假设的VCA算法，其具体实现过程是： (1)选择适当的超平面，将影像中的像元数据投影到选择的超平面上得到投影后的数据，后 续的步骤在投影后的像元数据上操作；（2)将像元数据投影到随机方向，并将具有最大投影 的像元作为第一个端元；（3)通过迭代地将数据投影到与由已提取端元构成的子空间正交 的方向上的方法来提取其余端元，被极限投影所对应的像元作为新的端元，直到所有的像 元被提取出来。
[0027] 作为优选，步骤3中所述aj(t)选择aj(t)~U(-0.01，0.01)。
[0028]作为优选，步骤5中所述循环终止条件是X(t)收敛或达到最大迭代次数。
[0029]作为优选，所述方法还包括对最后的输出结果进行校正，校正公式为：
[0030]
[0031]
[0032] 其中，SVCAeRlwSVCA方法提取的端元结果，认为是混合像元;Slmmc- PSQeRM为线性 混合模型约束的粒子群优化算法端元提取的结果，认为是纯净像元;funmixing是丰度估计的 函数，选择无约束的最小二乘方法;Ανω是采用纯净像元Slmmopsq求解混合像元Svca得到的丰 度;校正之后的端元光谱曲线& 11131为最终的结果。
[0033] 本发明具有以下优点和积极效果：
[0034] (1)创造性的提出线性混合模型约束的粒子群优化算法，该方法没有速度更新公 式，根据端元提取问题，利用线性混合模型更新粒子位置。
[0035] (2)保持了高光谱图像光谱特征物理意义。
[0036] (3)本发明模型简单，参数少，易于实现。
【附图说明】
[0037]图1是本发明实施例的流程不意图。
[0038] 图2是本发明实施例的实验结果示意图。
【具体实施方式】
[0039] 为了便于本领域普通技术人员理解和实施本发明，下面结合附图及实施例对本发 明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不 用于限定本发明。
[0040] 本发明需要解决的问题包括：（1)在纯净像元不存在的情况下，粒子群优化算法不 能直接应用于端元提取问题；（2)高光谱图像光谱特征物理意义的保持。
[0041] 本发明关键发明点为创造性提出线性混合模型约束的粒子群优化算法，并由于解 决纯净像元不存在情况下的端元提取问题。
[0042] 请见图1，本发明提供的一种高光谱遥感影像端元提取方法，将高光谱遥感影像分 别读入大小为1\11的矩阵1?={>1，...，；1^，...，1' 11}中，矩阵中每列向量；1^=(1'1」，...， ...，ru)TS像素的光谱辐射值，1为高光谱遥感影像的波段数，η为遥感影像的像素数； 对高光谱遥感影像矩阵R进行如下操作：
[0043] 步骤 1:初始化种群，Xi(t) = (ei,i(t)，ei,2(t)，· · ·，ei,D(t))i = l，2, · · ·，Μ;其中，t 表示迭代的次数，初始化过程中令I = O5X1U)表示种群中的一个粒子;M表示该种群中粒子 的数目；eij(0)eR，表示初始化过程中第i个粒子中的第j个向量;基于纯净像元存在假设的 顶点成分分析VCA(Vertex Component Analysis)算法提取的结果作为初始的种群；
[0044] VCA算法是比较成熟的、效果较好的端元提取算法，简单的操作步骤为：（1)选择适 当的超平面，将影像中的像元数据投影到选择的超平面上得到投影后的数据，后续的步骤 在投影后的像元数据上操作；（2)将像元数据投影到随机方向，并将具有最大投影的像元作 为第一个端元；（3)通过迭代地将数据投影到与由已提取端元构成的子空间正交的方向上 的方法来提取其余端元，被极限投影所对应的像元作为新的端元，直到所有的像元被提取 出来。
[0045] 步骤2:计算每初始粒子的个体最优位置P1(O)及所有粒子的全局最优位置G(O);
[0049] 共屮·；衣不;Ui胆皮圆双；衣不水胛AUt斤平形体的体积；
[0046]
[0047]
[0048] 丰度ai,j( $系数。
[0050]
[0051]
[0052]
[0053] 其中，aj(t)为一个随机生成数，表示图像的混合程度，在算法中选择q(t) ~U(_ 0.01，0.01) 是一个随机的丰度向量，需要满足的条件是:该向量中的第j个元素是0， 且满足和为一条件;灼⑴表示的是一个〇~1的随机生成数；
[0054] 步骤4:更新粒子的个体最优位置及所有粒子的全局最优位置；
[0055] 由于本发明所用的高维编码方式，粒子的每个列向量都是一个端元光谱曲线，每 个列向量的进化都会引起适应度函数的变化。为了避免每个列向量同时变化带来的叠加结 果的影响，在本发明中，采用每个粒子的每个列向量更新的同时更新粒子的个体最优以及 所有粒子的全局最优位置。
[0056] 蛉早个优晶优仿罟.
[0057]
[0058]其中，{XdthX^U+l)}表示的是X1⑴粒子的第j个列向量已经更新;f(x)表示适 应度函数，步骤2已经给出定义；
[0059]所有粒子的全局最优位置：
[0060] G(t+1) = arg^min /(f (1+1))
[0061] 步骤5:判断X(t)是否收敛或达到最大迭代次数；
[0062]如果否，重新执行第3步；
[0063] 如果是，则循环终止，输出最终的结果G(t)。
[0064] 在上述过程中，无法保证高光谱图像光谱特征物理意义，因此会出现光谱曲线光 滑性很差的现象。后处理主要是是根据高光谱成像的物理机制，也就是线性混合模型对上 述的结果进行校正，校正公式为：
[0065]
[0066]
[0067] 其中，SvcaEIT^VCA方法提取的端元结果，认为是混合像元;Slmmc-psoER1?为线性 混合模型约束的粒子群优化算法端元提取的结果，认为是纯净像元;funmixing是丰度估计的 函数，选择无约束的最小二乘方法;Ανω是采用纯净像元Slmmopsq求解混合像元Svca得到的丰 度;校正之后的端元光谱曲线& 11131为最终的结果。
[0068]校正之后的端元光谱曲线就是最终的结果。
[0069] 基于矩阵R，对高光谱遥感影像进行端元提取，在具体实现的时候还有以下注意事 项：
[0070] (1)粒子的定义:采用高维编码方式，每个粒子为一个矩阵而不是一个向量，粒子 中的每个列向量代表一条待定的端元光谱曲线，这样的定义方式更符合问题的物理意义。
[0071] (2)自适应度函数f的选择。自适应度函数包括两个部分:一个部分是体积，负责粒 子向着单形体的内内部移动，避免单形体的体积过大，另外一个部分是丰度约束，负责粒子 向着单形体的外部移动，避免单形体的体积过小。
[0072]
[0073]其中，au为解混丰度，γ为比例系数，两个约束条件的比例。
[0074]以下通过对比试验来验证本发明的有益效果。
[0075] 本试验采用的数据为仿真数据，共224个波段，光谱覆盖范围为:0.38-2.5μπι，影像 尺寸64像素\64像素。分别采用￥04、1^-匪？、1^4、1^3、31341^和本发明方法进行端元提 取。
[0076] 试验采用的定量评价指标为光谱角距离，比较所提取的端元光谱与真实参考光谱 的光谱角距离，光谱角距离越小，说明该方法的效果越好。每个方法运行10次，取统计结果， 实验结果见附图2。从试验结果上看，我们的方法效果明显优于对比算法。
[0077] 应当理解的是，本说明书未详细阐述的部分均属于现有技术。
[0078] 应当理解的是，上述针对较佳实施例的描述较为详细，并不能因此而认为是对本 发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不脱离本发明权 利要求所保护的范围情况下，还可以做出替换或变形，均落入本发明的保护范围之内，本发 明的请求保护范围应以所附权利要求为准。
【主权项】
1. 一种高光谱遥感影像端元提取方法，其特征在于:将高光谱遥感影像分别读入大小 为1><11的矩阵1?={1'1,...,;^,...,1'。}中，矩阵中每列向重:^ = (1'川...，1'山...，1'^)了为像 素的光谱福射值，1为高光谱遥感影像的波段数，η为遥感影像的像素数;对高光谱遥感影像 矩阵R进行如下操作： 步骤 1:初始化种群,Xi(t) = (ei,i(t) ,ei,2(t), . . . ,ei,D(t))i = l ,2, . . . ,Μ;其中，t表不 迭代的次数，初始化过程中令t = 0;Xi(t)表示种群中的一个粒子;Μ表示该种群中粒子的数 目；eu(0)eR，表示初始化过程中第i个粒子中的第j个向量;基于纯净像元存在假设的顶点 成分分析VCA算法提取的结果作为初始的种群； 步骤2:计算每初始粒子的个体最优位置Pi(0)及所有粒子的全局最优位置G(0);其中，f( ·)表示适应度函数;vol(Xi(t))表示求解Xi(t)所构成的单形体的体积;丰度 ai, J (t)满巧丫表示两个约束之间的比例系数； 步骤3:粒子进化;粒子进化方程为： Xi,j(t+l) = (pi,j(t)-aj(t) · (pi(t)*0j(t)))/(l-aj(t)) iU')=巧(Ο'巧(Ο + [1-A (坤G,' W W ~1/化4 其中，a^t)为一个随机生成数，表示图像的混合程度;是一个随机的丰度向量，需 要满足的条件是:该向量中的第j个元素是0,且满足和为一条件;巧(0表示的是一个0~1的 随机生成数； 步骤4:更新粒子的个体最优位置及所有粒子的全局最优位置； 粒子个体最优位置：其中，|Xi(t)，Xi,j(t+l)}表示 的是Xi(t)粒子的第j个列向量已经更新;f(x)表示适应度函数； 所有粒子的全局最优位置：步骤5:判断是否达到循环终止条件； 如果否，重新执行第3步； 如果是，则循环终止，输出最终的结果G(t)。2. 根据权利要求1所述的高光谱遥感影像端元提取方法，其特征在于:步骤1中所述基 于纯净像元存在假设的VCA算法，其具体实现过程是：（1)选择适当的超平面，将影像中的像 元数据投影到选择的超平面上得到投影后的数据，后续的步骤在投影后的像元数据上操 作；（2)将像元数据投影到随机方向，并将具有最大投影的像元作为第一个端元；（3)通过迭 代地将数据投影到与由已提取端元构成的子空间正交的方向上的方法来提取其余端元，被 极限投影所对应的像元作为新的端元，直到所有的像元被提取出来。3. 根据权利要求1所述的高光谱遥感影像端元提取方法，其特征在于：步骤3中所述aj (t)选择 aj(t)~υ(-Ο.ΟΙ,Ο.ΟΙ)。4. 根据权利要求1所述的高光谱遥感影像端元提取方法，其特征在于:步骤5中所述循 环终止条件是Χ( t)收敛或达到最大迭代次数。5. 根据权利要求1所述的高光谱遥感影像端元提取方法，其特征在于:所述方法还包括 对最后的输出结果进行校正，校正公式为：其中，SvcaER1*e为VCA方法提取的端元结果，认为是混合像元;Slmmc-psoERW为线性混合 模型约束的粒子群优化算法端元提取的结果，认为是纯净像元;fun^ixing是丰度估计的函数， 选择无约束的最小二乘方法;Avca是采用纯净像元Slmmc-ps冰解混合像元SVCA得到的丰度;校 正之后的端元光谱曲线Sfinal为最终的结果。
【文档编号】G06T7/00GK105844619SQ201610156222
【公开日】2016年8月10日
【申请日】2016年3月18日
【发明人】杜博, 许明明, 张良培, 张乐飞
【申请人】武汉大学