高光谱影像半监督弱信号解混方法

1.本发明属于图像处理领域，涉及高光谱影像半监督弱信号解混方法。


背景技术：

2.高光谱遥感技术，也称成像光谱，在每个像素点上以数十或数百个连续间隔的窄光谱带记录地表反射或发射的辐射，已很好地应用于重金属污染、精准农业。由于传感器有限空间分辨率，图像中易包含由多种端元线性组合的混合像素。然而，端元通常被认为是代表图像中主要的地物，而当地物中由于低混合成份和有限地物覆盖造成强度低的弱信号时，对弱信号的讨论和研究通常不在传统高光谱分析的范式中但这一类弱信号地物往往具有较大的探测价值。
3.大多数情况下，弱信号表示一些小而有用的信息，如污染或矿物。水质监测的一个典型例子是涉及间接测量叶绿素浓度的藻类水华检测。但低叶绿素浓度如小于2mg/m3时污染水体的光谱仅在蓝色波长(400
–
500nm)范围内具有较高的反射率。另一个例子是可用于估算地区的绿地密度的归一化差异植被指数(ndvi)，它是一个无量纲指数，描述植被覆盖的可见光和近红外反射率之间的差异。随着植被覆盖度的增加，ndvi趋于饱和。也就是说，如果植被覆盖以及低叶绿素浓度与背景覆盖相比不够饱和，则观测地面的解释精度无法满足要求。在这些情况下，由于实际应用中信号较弱，传统的分解策略往往会遇到瓶颈。目前，大多数研究者更多关注传统的高光谱解译方法，而对强度低的地物弱信号的处理则没有明确的分析范式。
4.基于非负矩阵分解(nmf)的方法在高光谱解混中通过使用乘法更新规则交替优化目标变量来将高光谱数据近似地分解为端元矩阵和丰度矩阵的乘积。然而，当地物中存在强度低的弱信号时，采用nmf方法进行图像解译时精度不够。由于某些目标具有细微的混合成份，解决这样的问题的难度在于初始化和优化很困难。其次是当引入先验知识来检测和分离相应的弱信号时，它可能会在优化任务期间退化。在这方面，基于稀疏解混的模型可能是一个不错的选择，因为这类方法引入了现有的光谱库并且只优化了丰度。但当光谱库并不完整或者数据过大时，它可能会导致不匹配、高相干性或高计算复杂度问题。
5.半监督nmf模型在某种程度上可以嫁接在高光谱弱信号分析中，原因在于这类模型剥离了端元同时优化模式而对不同的端元进行差异化的优化。但现有基于半监督nmf的方法仅认为一个区域具有的某些土地覆盖物光谱特征可以从标准库或地面调查中检索到。当一些纯端元在执行解混任务之前已知并可用时，它们将被视为端元库并随后替换初始化的端元。然而，如果谱库和初始化端元之间的谱相关性低，则引入的端元库可能不匹配已经初始化的端元。其次，大量先前的工作证明基于nmf的方法可以在相对低纯度的情况下产生所需的解混结果。


技术实现要素：

6.有鉴于此，本发明的目的在于提供一种高光谱影像半监督弱信号解混方法。
7.为达到上述目的，本发明提供如下技术方案：
8.高光谱影像半监督弱信号解混方法，该方法包括以下步骤：
9.s1：获取卫星对指定地面范围的遥感观测影像，完成影像辐射矫正和大气矫正预处理工作；
10.s2：采用影像光谱波峰波谷粗分析和实地勘测相结合的方式估测出卫星所观测的地区可能存在的弱信号地物，包括藻类污染的水体和工业废水污染的地面；通过实验室测量弱信号所代表的标准光谱；
11.s3：对经过预处理的遥感影像进行初始化工作来获得影像中典型地物的初始光谱和对应的丰度分布范围；
12.s4：采用包括无监督模块、有监督模块和稀疏正则项的模型来估测影像中地物的光谱和对应的地物分布；同时获得弱信号地物的分布范围和混合成份。
13.可选的，所述无监督模块估算端元和丰度，其中端元包含由robust-osp初始化的端元和先验知识；
14.设表示包含l个波段和n个像元的高光谱图像矩阵，表示含有l个波段和
p
个初始化端元，表示由观测者添加的q个先验光谱用于检测高光谱数据中存在的弱小信号；无监督模块表述如下：
[0015][0016][0017]
其中端元矩阵和丰度矩阵需要满足非负约束且丰度矩阵需要满足丰度和为一约束；对于先验的优化过程会使得其出现退化，考虑一个有监督模块去最小化先验退化。
[0018]
可选的，所述有监督模块用于约束先验信息的退化，监测高光谱图像中的弱小信号；
[0019]
有监督模块表述如下：
[0020][0021]
其中a＝[0
p
×q；iq×q]表示指示矩阵用于从端元矩阵中提取出对应先验知识的端元光谱，v表示先验知识；v和的区别在于前者在模型中视为被观测对象不会被优化，而后者在模型中视为检测弱小信号的先验因而在模型中会被更新优化。
[0022]
可选的，所述无监督模块和有监督模块在高光谱弱小信号检测中的协同作用，并同时考虑丰度矩阵的稀疏特性，通过联合无监督模型和有监督模型，得到半监督重加权稀疏模型：
[0023][0024]
约简后，得到模型：
[0025][0026]
其中第一个保真项是无监督模块，第二个保真项为有监督模块，第三项为稀疏正则项；r表示对丰度矩阵非零元素进行重加权约束的矩阵用于迫使丰度矩阵变得更加稀疏。
[0027]
可选的，所述监测高光谱图像中的弱小信号具体为：
[0028]
s11：对于观测得到的高光谱图像数据y，采用端元初始化方法robust-osp获取典型地物的端元
[0029]
s12：对于观测得到的高光谱图像数据，通过野外调查或地物光谱反射波峰分析，获取该观测地区范围存在的弱地物信号并通过实验室光谱成像获得弱小信号的成像光谱并将其作为先验知识；
[0030]
s13：通过优化模型(5)获得弱小信号光谱所对应的地物分布情况。
[0031]
可选的，所述采用端元初始化方法robust-osp获取典型地物的端元具体为：
[0032]
s21：采用svd算法将高维数据y转为低维数据子空间x；
[0033]
s22：计算具有最大l2范数的数据点e1为初始低维端元；
[0034]
s23：将剩余数据投影至已确定端元所张成的子空间中，获得第二个低维端元e2，并更新端元集；
[0035]
s24：重复s23直至获得指定数目的低维端元；
[0036]
s25：将低维端元逆变换至高维端元。
[0037]
本发明的有益效果在于：本发明提出的hyperweak方法基于半监督稀疏非负矩阵设计，通过考虑两个保真度项来解决弱信号场景。第一个保真项无监督地估计由初始化端元和先验增广的端元和丰度。第二个保真项有监督地最小化先验退化。为了提高hyperweak的鲁棒性，添加了一个重加权稀疏正则项和基于robust-osp的端元初始化方法，使得所提出的hyperweak算法能有效且准确的估测弱信号的地物丰度成份和分布。
[0038]
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
[0039]
为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作优选的详细描述，其中：
[0040]
图1为六种提取的地物光谱和标准光谱的可视化对比；(a)为tree；(b)为water；(c)为dirt；(d)为road；(e)为sagebrush；(f)为water+montmor；
[0041]
图2为六种提取的地物丰度和标准地物丰度分布的可视化对比；
[0042]
图3为弱小信号分析的半监督模型流程示意图；
[0043]
图4为robust-osp算法的端元初始化示意图。
具体实施方式
[0044]
以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。需要说明的是，以下实施例中所提供的图示仅以示意方式说明本发明的基本构想，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合。
[0045]
其中，附图仅用于示例性说明，表示的仅是示意图，而非实物图，不能理解为对本发明的限制；为了更好地说明本发明的实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；对本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。
[0046]
本发明实施例的附图中相同或相似的标号对应相同或相似的部件；在本发明的描述中，需要理解的是，若有术语“上”、“下”、“左”、“右”、“前”、“后”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此附图中描述位置关系的用语仅用于示例性说明，不能理解为对本发明的限制，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语的具体含义。
[0047]
本发明包括以下步骤：
[0048]
s1：获取卫星对指定地面范围的遥感观测影像，完成影像辐射矫正和大气矫正预处理工作；
[0049]
s2：采用影像光谱波峰波谷粗分析和实地勘测相结合的方式估测出卫星所观测的地区可能存在的弱信号地物，包括藻类污染的水体和工业废水污染的地面；通过实验室测量弱信号所代表的标准光谱；
[0050]
s3：对经过预处理的遥感影像进行初始化工作来获得影像中典型地物的初始光谱和对应的丰度分布范围。
[0051]
s4：采用模型(5)来估测影像中主要地物的光谱和对应的地物分布。同时获得弱信号地物的分布范围和混合成份。
[0052]
无监督优化：
[0053]
hyperweak的无监督部分侧重于估算端元和丰度，其中端元包含由robust-osp初始化的端元和先验知识。假定表示包含l个波段和n个像元的高光谱图像矩阵，表示含有l个波段和
p
个初始化端元，表示由观测者添加的q个先验光谱用于检测高光谱数据中存在的弱小信号。无监督优化模块表述如下：
[0054][0055][0056]
其中端元矩阵和丰度矩阵需要满足非负约束且丰度矩阵需要满足丰度和为一约束。然而上述优化过程中，对于先验的优化过程会使得其出现退化。因此考虑一个有监督模块去最小化先验退化。
[0057]
有监督优化：
[0058]
有监督模块用于约束先验信息的退化，以有效监测高光谱图像中的弱小信号。有监督模型可以表述如下：
[0059][0060]
其中a＝[0
p
×q；iq×q]表示指示矩阵用于从端元矩阵中提取出对应先验知识的端元光谱，v表示先验知识。v和的区别在于前者在模型中视为被观测对象不会被优化，而后者在模型中视为检测弱小信号的先验因而在模型中会被更新优化。
[0061]
总体模型：
[0062]
本方案考虑无监督模块和有监督模块在高光谱弱小信号检测中的协同作用，并同时考虑丰度矩阵的稀疏特性，通过联合无监督模型和有监督模型提出了如下的半监督重加权稀疏模型：
[0063][0064]
对上述模型进行约简，有如下模型：
[0065][0066]
其中第一个保真项是无监督模块，第二个保真项为有监督模块，第三项为稀疏正则项。r表示对丰度矩阵非零元素进行重加权约束的矩阵用于迫使丰度矩阵变得更加稀疏。
[0067]
本发明中的弱小信号检测方法和现有传统解混方法运行在标准jasper ridge实验数据中，其中jasper ridge包含4种正常地物(tree，water，dirt，road)和两种弱信号地物(sagebrush和water+montmor)。实验结果展现在表1中。
[0068]
表1算法在jasper ridge数据集上的实验结果
[0069][0070]
由表1可知，对比现有传统高光谱解混的方法，本发明在高光谱图像存在弱小信号的场景下具有更低的均方根误差结果和最低的信号退化水平。同时，对于正常地物，本发明的亦能获得较好的端元提取能力。图1和图2展现了本发明从图像中所估计的6种地物(短实线)和标准光谱库(长实线)的可视化对比。由图1和图2所示，本发明所估计的端元光谱基本拟合了真实的端元光谱，同时所估计出的地物丰度和成份也基本符合真实的地物分布信息。
[0071]
图2中六种提取的地物丰度(第二行)和标准地物丰度分布(第一行)的可视化对比，从左到右为tree，water，dirt，road，sagebrush和water+montmor。
[0072]
图3展示了本发明在高光谱图像中进行弱小信号分析的示意图。该示意图可以简要的描述为如下步骤：
[0073]
步骤1：对于观测得到的高光谱图像数据y，采用端元初始化方法robust-osp获取典型地物的端元
[0074]
步骤2：对于观测得到的高光谱图像数据，通过野外调查或地物光谱反射波峰分析，获取该观测地物范围存在的弱地物信号并通过实验室光谱成像获得弱小信号的成像光谱并将其作为先验知识。
[0075]
步骤3：通过优化模型(5)获得弱小信号光谱所对应的地物分布情况。
[0076]
图4展示了本发明中所涉及到的重要步骤，即robust-osp初始化算法的初始化过程。该过程可以表述为以下步骤：
[0077]
步骤1：采用svd算法将高维数据y转为低维数据子空间x。
[0078]
步骤2：计算具有最大l2范数的数据点e1为初始低维端元
[0079]
步骤3：将剩余数据投影至已确定端元所张成的子空间中，获得第二个低维端元e2，并更新端元集。
[0080]
步骤4：重复第三步直至获得指定数目的低维端元。
[0081]
步骤5：将低维端元逆变换至高维端元。
[0082]
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。