高光谱遥感图像的线性全约束丰度估计方法
【专利摘要】本发明公开了一种高光谱遥感图像的线性全约束丰度估计方法，构建原光谱解混问题的对偶问题，通过全约束最小二乘法对丰度进行初始化，结合偏离度ρ因素利用原始对偶内点法优化丰度结果，具体包括以下步骤：S1:建立图像的线性光谱混合模型：S2:对线性光谱混合模型施加约束条件；S3:采用原始对偶内点法和全约束最小二乘丰度估计算法对线性光谱混合模型进行混合像元的丰度估计。该发明有效地克服了全约束最小二乘法精度不高的问题，同时也解决了原始对偶内点法速度慢的问题，可以作为一种解决高光谱遥感图像全约束丰度估计的有效手段。在基于高光谱遥感图像的高精度的解混以及地面目标的检测和识别方面具有重要的应用价值。
【专利说明】
高光谱遥感图像的线性全约束丰度估计方法
技术领域
[0001] 本发明涉及遥感图像处理技术领域，尤其涉及高光谱遥感图像的线性全约束丰度 估计方法。
【背景技术】
[0002] 高光谱遥感技术能够实现地物空间信息、辐射信息、光谱信息的同步获取，具有波 段信息丰富、分辨率高、图谱和一三大特点，因而包含了更多精细的地物信息。高光谱遥感 数据最终形成的是一个高光谱数据立方体，这个数据立方体就是一个包含了地物信息的数 据宝库，只有对这个数据宝库进行发掘、分析，才能得到有用的信息。
[0003] 由于受到遥感器分辨率的限制和地物信息复杂性的影响，高光谱遥感数据立方体 中包含了大量的混合像元。丰度估计就是对混合像元中端元所占的比例进行估计，所使用 的光谱混合模型主要有线性混合模型和非线性混合模型。线性混合模型假设混合像元由几 种纯地物按照一定的比例线性混合而成，具有明确的物理意义且容易建模，得到了广泛的 应用。
[0004] 根据实际的物理意义，混合像元的分解需要满足两项约束条件，一是非负约束 (Abundance Non-negative Constraint ANC)，即端元丰度不能出现负值;一是和为1约束 (Abundance Sum_t〇-〇ne Constraint，ASC)，即各端元丰度和为1。满足ANC和ASC约束的丰 度估计算法即为全约束丰度估计算法。现有技术中采用单一的算法如全约束最小二乘法进 行。
【发明内容】
[0005] 根据现有技术存在的问题，本发明公开了一种高光谱遥感图像的线性全约束丰度 估计方法，构建原光谱解混问题的对偶问题，通过全约束最小二乘法对丰度进行初始化，结 合偏离度P因素利用原始对偶内点法优化丰度结果，具体包括以下步骤：
[0006] S1:建立图像的线性光谱混合模型：
[0007] S2:对线性光谱混合模型施加约束条件；
[0008] S3:采用原始对偶内点法和全约束最小二乘丰度估计算法对线性光谱混合模型进 行混合像元的丰度估计:在该过程中融合偏离度P作为丰度估计的一个因素，用偏离度P表 征混合像元中不能被解析的部分，将端元数量、地形、遥感器、端元纯净性、噪声因素所引起 的影响使用解析偏差表示;从而获得高光谱图像中每个混合像元的最佳偏离度P Z，得到最 佳的丰度估计结果实现原图像的重构。
[0009] 所述线性光谱混合模型为：
[0010]
⑴
[0011]其中：r是混合像元中任意的1维光谱向量，p是端元个数，M=[mi，m2, . . . .mp]是大 小为l*p的端元光谱矩阵，α是p维矢量系数，其各分量元素是混合像元r中端元mi的丰度，e 为1维矢量，表示模型误差项。
[0012] 将约束条件施加于（1)式上，可以得到ASC约束为，ANC约束为:α0〇(? = i=l 1，2，. .p)，在最小二乘法算法的基础上加入ASC和ANC约束，则全约束最小二乘法丰度估计 算法为：
[0013] min|aeQ{ (r-Ma)T(r-Ma)}
[0014]〇 (2)
[0015] S3中将采用原始对偶内点法和全约束最小二乘法相结合时具体为:设公式(2)为 原始问题，根据原始对偶理论，引入参数人"^仏私…如广则公式⑵的对偶问题如下：
[0016]
< 3 )
[0017 ]根据KKT条件，公式⑵（3)的最优解等价于求解如下参数方程：
[0018]
(4)
[0019] 其中X = diag(a1，a2, · · ·，ap)，Z = diag(zi，Z2, · · ·ζΡ)，e = (li，li, · · · lP)T，u为障碍 因子，u-0。公式(4)中第一个式子为原可行性，第二个式子为对偶可行性，第三个式子为互 补松弛性；
[0020] 根据原始对偶理论，对偶间隙为J(a)-D(a)=aTz根据公式(4)中的第三个式子，对 偶间隙简化表示为J(a)-D(a)=pu，p为端元个数，当u-0时，J(a)-D( a)-0,设对偶问题解 为1=(〇八，4，则<1、1则分别收敛到原问题(2)和对偶问题(3)的最优解，则1( 〇)即原始对偶 内点法全约束丰度估计的结果；
[0021 ]设第k次迭代的丰度估计的结果为Wk= (ak，Ak，zk)，丰度估计的步长为Δ W= ( Δ a， Λ λ，Λ z)，A X = diag( A αι，Λ a2，· · ·，Α αρ) Λ Z = diag( Αζι，Αζ2，··· Αζρ)，Αζ=(Αζι， Δ Z2, · · · δ Ζρ)τΔα = ( Adi, Δα2, · · ·，Δαρ)τ，采用牛顿法对公式(4)求解，可得：
[0022]
[0023] 根据公式(5)可以计算得到 AW=( Δα, Δλ，ΔΖ)，Eak=l-lTak，dk=(MTM)a k-Uk-zk-MTr，ζζ = I I dk I 12〇
[0024] 在原始对偶内点法算法中加入全约束最小二乘丰度估计算法以及偏离度p具体采 用如下方式：
[0025] Stepl运行FCLS算法，以运行结果初始化aQ。
[0026] Step2初始化P，需要根据具体情况设置，例P = 0.05。
[0027] St印3初始化λ°，Λ由于α°的准确性很高，因此设置为较小的数值，例如le-8， 得到W*3 = (yd)。4
-，p为端元个数，同时给定较小的正数0〈δ〈 1;给定终止 误差ε>〇，置Ν=0，最大迭代次数Ns;
[0028] Step4当（（ak)Tz彡ε )&&(ak彡ε )&&(zz彡ε )，则停止丰度的查找，ak为丰度估计最优 值；否则转到Step5。
[0029] Step5$yk+l=yk(l_5)，根据公式(8)计算 Δ#=(Δα，Δλ，Δζ)
[0030] Step6ffk+1=ffk+〇Affk,N=N+l〇
[0031] Step7判断光谱重构误差是否递减，递减则转到Step 8，否则#为丰度估计最优 值，算法结束。
[0032] Step8如果N> = NS，算法结束，否则转到Step3。如果将Step8中Ns=l，即为FCLS+ ηπρ(ν=ι，ρ)算法。经过数据测试发现，该算法分解精度的提高主要集中在迭代的首次。因 此，可以通过控制N限制迭代次数，牺牲部分分解精度来较大的降低算法的运行时间。
[0033] 在混合像元的丰度估计中，将由于端元数量、地形、遥感器、端元纯净性、噪声等一 切因素所引起的影响都使用解析偏差统一表示，P即表示这种偏离度的大小。对每一个混合 像元而言，都有一个最佳的P Z，使得该混合像元的分解精度达到最佳。
[0034] 由于采用了上述技术方案，本发明提供的高光谱遥感图像的线性全约束丰度估计 方法，通过构造原解混问题的对偶问题，结合传统的全约束最小二乘丰度估计FCLS方法，将 高光谱图像全约束丰度估计问题转化为二次规划的优化问题。进而，该发明引入了丰度的 和为一偏离度，对于系统不能解析的噪声等因素进行表示，并将对偶间隙和解混误差同时 作为优化目标，提高了算法的解混精度。该发明从全约束线性高光谱混合模型出发，结合最 小二乘方法和原始对偶内点法，有效地克服了全约束最小二乘法精度不高的问题，同时也 解决了原始对偶内点法速度慢的问题，可以作为一种解决高光谱遥感图像全约束丰度估计 的有效手段。在基于高光谱遥感图像的高精度的解混以及地面目标的检测和识别方面具有 重要的应用价值。
【附图说明】
[0035] 为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下， 还可以根据这些附图获得其他的附图。
[0036]图1为本发明中Cuprite数据图的示意图；
[0037]图2为本发明中从Cuprite数据中提取的5种矿物质的端元光谱图；
[0038]图3为本发明中模拟高光谱数据合成区域图；
[0039] 图4是模拟高光谱数据彩色展示图；
[0040] 图5是FCLS、PDIP和优化FCLS+PDIP(N=1，P)算法的SRMSE和P值关系示意图；
[0041]图6是真实高光谱图像Hydice数据第120波段示意图。
【具体实施方式】
[0042]为使本发明的技术方案和优点更加清楚，下面结合本发明实施例中的附图，对本 发明实施例中的技术方案进行清楚完整的描述：
[0043] 如图1所示的一种高光谱遥感图像的线性全约束丰度估计方法，构建原光谱解混 问题的对偶问题，通过全约束最小二乘法对丰度进行初始化，结合偏离度P因素利用原始对 偶内点法优化丰度结果，具体包括以下步骤：
[0044] S1:建立图像的线性光谱混合模型：
[0045] S2:对线性光谱混合模型施加约束条件；
[0046] S3:采用原始对偶内点法和全约束最小二乘丰度估计算法对线性光谱混合模型进 行混合像元的丰度估计:在该过程中融合偏离度P作为丰度估计的一个因素，用偏离度P表 征混合像元中不能被解析的部分，将端元数量、地形、遥感器、端元纯净性、噪声因素所引起 的影响使用解析偏差表示;从而获得高光谱图像中每个混合像元的最佳偏离度P Z，得到最 佳的丰度估计结果实现原图像的重构。
[0047] 实施例:分别用仿真数据和实际遥感图像数据为例，说明本发明的【具体实施方式】。 我们用重构误差和运行时间来衡量全约束解混算法的优劣，运行环境为Wind 〇ws7系统下的 Matlab7 · 0 平台。
[0048] 重构误差米用光谱均方根误差（Spectral Root Mean Square Error，SRMSE)作为 丰度评价标准，如公式(6)所示。SRMSE值越小表明丰度估计的结果越准确，在混合像元分解 中，可以代表混合像元分解的精度，是一个绝对评价标准。具体计算公式如下：
[0049] (6)
[0050] 其中，p是波段个数，N是像素的个数，rij是第i个混合像元在第j波段的光谱亮度 值，Μ为端元矩阵， αι是第i个混合像元的丰度估计值。
[0051] 模拟高光谱数据Cuprite数据是一个公开的高光谱数据，因此以Cuprite图像数据 为基础，提取5种矿物质A、B、C、D、E，分别是Kaolinite，Calcite，Alunite，Buddingtonite， Muscovi te，如图1所示，它们的福射亮度值曲线如图2所示。
[0052]以5种矿物质为端元合成波段189个波段的32*32大小的模拟图像。将图像划分为 八，8，(：，0,6，11^，11，〇共9个区域，如图3所示。将0&1(3^6作为基础物质，在9个区域内都添加。 以在A区域内每个像元添加80 %的Alunite和20 %的Cal cite;在B区域内添加80 %的 1311(1(1;[1^1:011；^6和20%的0&1(3；^6;在(]区域内添加80%的1118(3〇￥；^6和20%的0&1(3；^6;在0 区域内添加80 %的Kaolinite和20 %的Calc ite;在e区域内添加40 %的Alunite，40 %的 1311(1(1;[]^1:〇11；^6以及20%的0&1(3；^6;在;1^区域内添加40%的厶11111；^6,40%的]\1118(3〇￥；^6以及 20% 的Cal cite;在g区域内添加40 % 的Muscovi te，40 % 的Kaolinite 以及 20 % 的Calc ite; 在11区域内添加40%的1(&〇1;[11；^6,40%的1311(1(1;[1^1:011；^6以及20%的0&1(3；^6;在0区域内添 加20%的厶11111；^6,20%的1311(1(1;[1^1:011；^6,20%的0&1(3；!^6,20%的1(&01;[11；^6以及20%的 Muscovite;同时为每个像元中加入30:1的信噪比，完成模拟混合像元数据。从189个波段中 提取3个波段合成彩色图4所示。
[0053] 仿真实验中的测试数据共有4组，本别依据FCLS、PDIP、FCLS+PDIP(P)、FCLS+roiP (P，N=1)产生。PDIP、FCLS+PDIP(p)、FCLS+PDIP(p，N=1)中的参数 ε 设置为 1θ-8;0〈δ〈1，取δ =〇. 1;更新步长σ = 0.9。
[0054]实验1测试重构误差与偏离度ρ之间的关系。在这个实验中，我们改变偏离度ρ的 值，分别测试各算法进行全约束解混时的重构误差。表1、2分别给出
[0055] 了 FCLS、PDIP的数据结果。表1FCLS算法在不同P值下的丰度估计结果
[0056] 表1 FCLS算法在不同P值下的丰度估计结果
[0057]
[0058] 表2 roip算法在不同P值下丰度估计结果
[0059]
[UUOU」 从表1、2口」以宥出，迠刍的増刀口P值口」以进一步的陴仳SKMSE值，也就是提咼J分解 精度，说明P预设值在向最佳PZ靠近。如果继续增大P值，SRMSE值反而会升高，也就是降低了 分解精度，说明P的预设值在偏离最佳P Z。因此，要想获得最佳的丰度估计结果，就要找到最 佳PZ，而且是找到每个混合像元的最佳P Z，才能实现最佳的丰度估计。
[0061 ]从表1、2中还可以看出，PD IP算法可以事先精确的控制每个混合像元的P值，而 FCLS算法则做不到。FCLS算法得到的P值范围在预设值上下波动，从表1中可以看出，P均值 小于预设值。从分解结果上看，FCLS算法的SRMSE值略小于Η)ΙΡ算法，也就是分解精度较高， 说明FCLS算法比Η)ΙΡ算法在噪声处理上较强，对P值有一定的自适应性，有自动寻找最佳p z 的趋势。
[0062] roip算法可以事先精确的控制p值，这是该算法的特点，但是在丰度估计中并不是 优点，因为每个混合像元都有自己最佳Pz，因此事先精确控制Pz并不是一个好的选择。但是， 如果以FCLS算法分解结果的P值作为标准，事先设置在PDIP算法上，经过数据测试可以看 到，roiP算法得到的丰度值、SRMSE值与FCLS算法一致。这也说明了最佳的丰度估计在于对 噪声的处理，在于对最佳Pz的寻找。
[0063]由于事前并不知道最佳PZ，因此P的设置上具有盲目性，如果设置的不好，不仅起 不到提高分解精度的目的，反而会使分解经过更差。例如根据数据的模拟过程，信噪比为 30:1，因此如果设置P= 1/30，不管是FCLS算法还是Η)ΙΡ算法，SRMSE值反而高于P = 0时的分 解结果。如果用NSCLS算法做对比，NSCLS算法是实现了ASC约束，但是在实现ANC约束时，采 用标准化的方式进行，打破了目标函数最小化的约束，因此是一个分解效果比较差的全约 束丰度估计算法。NSCLS算法分解结果的SRMSE值为126.09609195，FCLS、roiP算法的SRMSE 在P = 1/30时已经比较接近该值了，说明FCLS、PDIP算法在P = 1/30时已经退化成了一个比 较差的算法了。因此，在FCLS、PDIP算法上如果添加P值不当，会严重影响算法的分解精度， 也说明了在FCLS、PDIP算法上直接添加P值，可行性不高。
[0064]然后，通过结合FCLS与PDIP，加入解混误差标准，使得该算法对P值有了较强的自 适应性，提高了分解精度，而且大幅度降低了在P值设置不合理的情况下分解精度的误差， 避免了 P值设置不合理的危险性。表3是优化FCLS+PDIP(N=1，p)算法在不同P值下的丰度估 计结果。图5是优化FCLS+PDIP (N= 1，P)算法与FCLS、PDIP算法分解结果的比较。
[0065] 表3优化FCLS+H)IP(N=1，p)算法在不同P值下丰度估计结果 [0066]
[0067] 从图5以及表1、2、3可以看到，优化FCLS+PDIP(N= 1，P)算法在分解精度上有了进 一步的提高，而且在P值设置较大时，并不会降低分解的精度，相对于原始FCLS、PDIP算法， 分解精度仍有提高。因此，本发明算法对P值有了一定的自适应性，强化了P值的自动寻找最 佳P Z的能力。
[0068]实验2测试各算法的运行时间，每种算法的每组参数设置都在测试数据上运行20 次后取平均值。PDIP算法花费时间较长，FCLS和PDIP算法想结合之后，虽然相对于roip能够 降低算法的部分运行时间，但是仍然较长。由于该算法分解精度的提高主要集中在首次，因 此可以通过牺牲部分分解精度，控制迭代次数(N=l)来有效的降低算法的运行时间。以P = 0.003为例，迭代次数、SRMSE值、时间如表4所示。
[0069] 表4 FCLS+H)IP(p)算法在Ρ = 0·003下的丰度估计
[0070]
[0071] 从表4可以看出，优化FCLS+PDIP(p)在分解精度上主要集中在首次FCLS+PDIP(N = 1，P)，牺牲之后的迭代次数可以有效的降低分解时间，使得分解时间靠近FCLS算法。针对不 同的P值，3种算法的运行时间结果如表5所示。
[0072] 表5各算法在不同P值下的运行时间
[0073]
[0075]结合之前的分析，优化FCLS+PDIP(N=1，p)算法在分解精度上有了较大的提高，而 且为了兼顾分解时间，控制N=l，使得优化FCLS+PDIP(N=1，p)算法成为了一个比较高效的 算法。在模拟高光谱数据集上，令P = 〇.01，比较本发明算法和原有算法的分解精度，以及本 发明算法的分解优势，如表6所示。
[0076]表6优化FCLS+PDIP(N=1，P = 0.01)算法相对于其他算法在分解精度上的提高比 例。
[0077]
[0078] 2、实际遥感高光谱数据
[0079]真实高光谱数据采用Hydice数据进行测试，共169波段，光谱空间分辨率为1.5米。 图像大小为64*64，图6显示的是第120波段的图像。Hydice数据主要包含interferer， grass，tree，road及其他地物。本发明分别从中提取3、6、9、12、18、24、36个端元进行全约束 丰度估计测试。
[0080]实验1测试各算法的解混精度。在真实高光谱数据上，由于地物环境的复杂性，噪 声通常较大，认为5%的噪声是一个合理的预设值，因此设置p = 0.05。验证不同端元下的丰 度估计结果，其重构误差和相他算法的精度提高比例如表7、8所示。
[0081 ] 表7全约束丰度估计算法的SRMSE值(P = 0.05)
[0082]
[0083] 表8优化FCLS+PDIP(N=1，P = 0.05)算法相对于其他算法在分解精度上的提高比 例
[0084]
[0085]在表7中，从纵向数据比较来看，每种算法都随着端元个数的增多，SRMSE值逐渐减 小，也就是说分解精度都在提高。从横向数据比较来看，FCLS算法SRMSE值相对于PDIP算法 只有微弱的降低，而本发明算法较大幅度的降低了 SRMSE值。结合表8可以看到，本发明算法 相对于FCLS、PDIP算法在分解精度上，提高比例基本都在10%以上，提高的幅度较大。因此， 可以看到本发明算法在端元确定的情况下，不用通过增加端元个数，也能够较大幅度的提 高分解精度。
[0086] FCLS+PDIP(N=1，p)算法虽然能够一定程度的自动寻找最佳pz，但仍然不能找到 最佳p z，可以通过设置不同的P值，找到一个比较精确的丰度估计结果。表9是在不同P的预 设值下的丰度估计结果。
[0087] 表9优化FCLS+PDIP(N=1，p)算法丰度估计结果SRMSE值 [0088]
[0089] 可以看到在P的预设值在0.1时，丰度估计的精度比P = 0.05时更好。在端元个数为 3、6时，P = 0.15则SRMSE值进一步减小，因此对于Hydice数据，在不同的端元个数下，都有其 对应的最佳PZ，要针对每一种情况进行特定的分析，才能得到较好的丰度估计结果。如果能 做到对每个混合像元都做到个性化分析，则能够得到更为精确的丰度估计结果。
[0090] 实验2测试各算法的运行时间，结果如表10所示。
[0091] 表10全约束丰度估计算法的运行时间
[0092]
[0093]从表10中可以看出，优化FCLS+PDIP (N = 1，P)算法随着端元数目的增多，相对于N 不限时节省的时间也越来越多，端元数目从9个开始，FCLS+PDIP(N=1，p)算法的执行时间 逐步靠近FCLS算法。因此，本发明FCLS+PDIP (N= 1，P)算法是一个兼顾了分解精度与分解时 间的比较高效的一个算法。
[0094]综上所述，可以得出两点结论:一是本发明提出了混合像元分解精度和P值之间深 刻关系，只有找到高光谱图像中每个混合像元的最佳PZ，才能得到最佳的丰度估计结果;二 是优化FCLS+PDIP(N= 1，P)算法虽然不能完全实现最佳pz的自动查找功能，但是避免了 P值 设置不合理的危险性，而且在一定程度上实现P值的自动查找功能，较大程度的提高了混合 像元的分解精度，而且兼顾了算法的运行时间，算法运行也较快，因此本发明仍然是一个比 较高效的算法。
[0095]以上所述，仅为本发明较佳的【具体实施方式】，但本发明的保护范围并不局限于此， 任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其 发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
【主权项】
1. 一种高光谱遥感图像的线性全约束丰度估计方法，其特征在于：构建原光谱解混问 题的对偶问题，通过全约束最小二乘法对丰度进行初始化，结合偏离度P因素利用原始对偶 内点法优化丰度结果，具体包括W下步骤： S1:建立图像的线性光谱混合模型： S2:对线性光谱混合模型施加约束条件； S3:采用原始对偶内点法和全约束最小二乘丰度估计算法对线性光谱混合模型进行混 合像元的丰度估计:在该过程中融合偏离度P作为丰度估计的一个因素，用偏离度P表征混 合像元中不能被解析的部分，将端元数量、地形、遥感器、端元纯净性、噪声因素所引起的影 响使用解析偏差表示;从而获得高光谱图像中每个混合像元的最佳偏离度Pz，得到最佳的 丰度估计结果实现原图像的重构。2. 根据权利要求1所述的高光谱遥感图像的线性全约束丰度估计方法，其特征还在于： 所述线性光谱混合模型为：(1) 其中：r是混合像元中任意的1维光谱向量，P是端元个数，Μ =[mi，m2，--mp ]是大小为 l*p的端元光谱矩阵，α是P维矢量系数，其各分量元素是混合像元r中端元mi的丰度，e为1维 矢量，表示模型误差项。3. 根据权利要求2所述的高光谱遥感图像的线性全约束丰度估计方法，其特征还在于： 将约束条件施加于（1)式上，可W得到ASC约束为：、NC约束为:αι>〇α = 1，2，..ρ)， 在最小二乘法算法的基础上加入ASC和ANC约束，则全约束最小二乘法丰度估计算法为： min|aeQ{ (r-Ma)^(r-Ma)}。 （2)4. 根据权利要求1所述的高光谱遥感图像的线性全约束丰度估计方法，其特征还在于： S3中将采用原始对偶内点法和全约束最小二乘法相结合时具体为:设公式(2)为原始问题， 根据原始对偶理论，引入参数λ，Z = (Zl，Z2，...Zp)T，则公式(2)的对偶问题如下：m 根据KKT条件，公式(2) (3)的最优解等价于求解如下参数方程：(4) 其中X = diag(ai,a2, . . .，ap)，Z = diag(zi,Z2, . . .Zp)，e=(li,li, . . .lp)T，u为障碍因子， u^O;公式(4)中第一个式子为原可行性，第二个式子为对偶可行性，第Ξ个式子为互补松 弛性； 根据原始对偶理论，对偶间隙为J(c〇-D(a)=aTz根据公式(4)中的第立个式子，对偶间 隙简化表示为J(a)-D(a)=pu，p为端元个数，当u一 ο时，J(a)-D(a) 一 0,设对偶问题解为W = (α，λ，u)，则α、W则分别收敛到原问题(2)和对偶问题(3)的最优解，则W(α)即原始对偶内点 法全约束丰度估计的结果； 设第k次迭代的丰度估计的结果为Wk=(αk，λk，zk)，丰度估计的步长为ΔW=(Δα，Δλ， Δζ)，AX = dia 邑（Δ曰 1，Δ曰2,. , Δ曰ρ)Δ Z = dia 邑（Δζι, Δζ2,. . . Azp), Δζ=( Δζι, Δ Ζ2，. . . Δζρ)ΤΔα = ( Δαι，Δα2，. . .，Δαρ)τ，采用牛顿法对公式(4)求解，可得：(5) 根据公式（5)可W计算得到ΔΙ=(Δα，Δλ，Δζ)，记址=1-1了〇>%化=(1了1)〇^1人^2^ ΜΤ;γ，ζζ = I I dk I 12。5.根据权利要求4所述的高光谱遥感图像的线性全约束丰度估计方法，其特征还在于： 在原始对偶内点法算法中加入全约束最小二乘丰度估计算法W及偏离度P具体采用如下方 式： St邱l运行FCLS算法，^运行结果初始化口<^; St邱2初始化P，根据实际情况设置P的数值； 8*6口3初始化人<\2<\得到胖*^=川，人<\2<〇，令，P为端元个数，设〇<δ<1;给定 终止误差ε〉0，置Ν=0，最大迭代次数化； Step4当（（ak)Tz《e)&&(址《ε)&&(ζζ《ε)，则停止丰度的查找，ak为丰度估计最优值； 否则转到Step5; St邱5令yk" = yk(l-S)，根据公式(5)计算 AWk=( Δα, Δλ, Δζ); St邱6Wk+i=Wk+o Δ wk，N=化1; St邱7判断光谱重构误差是否递减，递减则转到Step 8,否则ak为丰度估计最优值，算法 结束； St邱8如果N〉=化，算法结束，否则转到Step3。
【文档编号】G06T7/40GK105976404SQ201610281811
【公开日】2016年9月28日
【申请日】2016年4月29日
【发明人】宋梅萍, 安居白, 张建祎, 徐行伟
【申请人】大连海事大学