基于自权重协同表示的高光谱异常目标检测方法技术

本发明专利技术提供了一种基于自权重协同表示的高光谱异常目标检测方法。首先，设定双窗结构的大小，并以待测像素点为中心，选取双窗结构内的所有像素点构建字典；然后，构建自权重协同表示的异常检测目标函数，将字典代入其中，并采用迭代更新方法求解目标函数，得到期望的稀疏向量；最后，利用稀疏向量计算待检测像素点与重建像素之间的误差，并通过与阈值比较判断该像素点是否属于异常目标。由于本发明专利技术将权重学习与协同表示结合起来构建目标函数，可以有效地提升异常目标检测效果。

Hyperspectral outlier detection based on self weight cooperative representation

【技术实现步骤摘要】
基于自权重协同表示的高光谱异常目标检测方法
本专利技术属高光谱异常目标检测
，具体涉及一种基于自权重协同表示的高光谱异常目标检测方法。
技术介绍
高光谱异常目标检测是近年来高光谱图像处理的热点研究领域之一，在民用和军事上具有广泛的应用。高光谱异常目标检测属于一种不需要任何先验信息的目标检测问题。传统的高光谱异常目标检测方法假设图像背景服从高斯分布的统计特性，然而由于图像背景的复杂性以及异常目标的干扰，通过该统计假设设计的异常检测算法的异常检测效果并不理想。近年来，协同表示被引入到目标检测，取得了较好的检测效果。文献“W.LiandQ.Du,Collaborativerepresentationforhyperspectralanomalydetection,IEEETrans.Geosci.RemoteSens.,vol.53,no.3,pp.1463–1474,Mar.2015.”公开了一种基于协同表示的高光谱异常目标检测方法。该方法以待检测像素点为中心构建双窗，以双窗之间的所有像素点作为字典，再通过稀疏表示模型求解稀疏向量，最后通过待检测像素点与重建像素之间的误差判断该像素点是否属于异常目标。虽然基于协同表示的高光谱异常目标检测方法取得了一定的检测效果，但该方法将高光谱图像的众多波段同等对待，忽视了冗余甚至噪声等低质量的波段对检测效果的不利影响。
技术实现思路
为了克服现有技术的不足，本专利技术提供一种基于自权重协同表示的高光谱异常目标检测方法。一种基于自权重协同表示的高光谱异常目标检测方法，其特征在于步骤如下：步骤1：对于高光谱图像的待检测像元i＝1,…,n，n为高光谱图像中的像元总数，d为高光谱图像的波段数，以yi为中心分别构建大小为ωout和win的双窗，双窗之间的所有像元构成字典Yc，即c表示字典包含的像元个数，c＝wout×wout-win×win，win＜wout＜min{h，w}，h和w分别为高光谱图像的行数和列数，n＝h×w，为字典中的第j个像元，j＝1,…,c；步骤2：构建高光谱异常检测目标函数如下：其中，αi为像元yi异常检测待求解的期望稀疏向量，∑为权值矩阵，是以自权重向量为对角元素的对角矩阵，σk表示自权重向量σ的第k个元素，λ为正则化参数，可取值范围为10-6至100；步骤3：采用迭代更新求解算法对目标函数进行求解，具体为：步骤3.1，初始化：按照初始化权重矩阵，I为单位阵；步骤3.2，固定∑，更新αi：固定αi，更新∑：∑＝diag(σ)(4)其中，mk表示向量m的第k个元素，向量m＝yi-Ycαi，σ为由其σk组成的向量，k＝1,…,d；步骤3.3，利用更新后的αi和∑计算目标函数值，并重复步骤3.2，直至此次计算的目标函数值和上一次计算的目标函数值之差小于10-6时，迭代终止，此时得到的αi即为像元yi异常检测的稀疏向量；步骤4：按下式计算像元yi的重建误差δi：δi＝||yi-Ycαi||.(5)之后对δi归一化，给定阈值A，如果δi>A，则判定像元yi为异常像素点，否则，判定像元yi为背景像素点；所述的阈值A的可取值范围为0至1之间。本专利技术的有益效果是：由于本专利技术将权重学习与协同表示结合起来，设计了统一的目标函数，高光谱图像中高质量的波段能被自动赋予较大的权重，而低质量的波段被自动赋予较小的权重，使得加权后的波段特征更具判别性，有效提升了异常目标检测的效果。附图说明图1是本专利技术的基于自权重协同表示的高光谱异常目标检测方法流程图；图2是采用传统方法进行高光谱异常目标检测结果图；图3是采用本专利技术方法进行高光谱异常目标检测结果图。具体实施方式下面结合附图和实施例对本专利技术进一步说明，本专利技术包括但不仅限于下述实施例。高光谱异常目标检测方法通过对每一个像素点计算其为异常目标的分值大小，当该分值大于一定的阈值时，则判断该像素点属于异常目标像素，否则属于背景像素点。假定待处理的高光谱图像是一个三维矩阵，包含图像的行数、列数和波段数，为方便计算，可以将其转换成一个二维矩阵其中，d表示高光谱图像的波段数，n为高光谱图像所有像素点的个数。如图1所示，本专利技术提供了一种基于自权重协同表示的高光谱异常检测方法，其基本实现过程如下：1、通过双窗结构构建字典假定待测的像素点为以yi为中心构建双窗结构。令wout表示外窗的大小，win表示内窗的大小。选择双窗之间的所有像素点构建字典其中，c表示字典包含的像元个数，c＝wout×wout-win×win，win＜wout＜wmin{h，w}，h和ω分别为高光谱图像的行数和列数，有n＝h×ω，为字典中的第j个像元，j＝1,…,c。2、基于自权重协同表示的高光谱异常检测方法传统的基于协同表示的高光谱异常目标检测方法的目标函数是：其中，λ为正则化参数，αi为希望得到的像元yi异常检测的稀疏向量。本专利技术提出的基于自权重协同表示的高光谱异常检测方法，在上述方法的基础之上考虑到不同质量的波段会对检测结果产生不同的影响，因此需要对不同的波段赋予合适的权重。然而，高光谱异常目标检测实际上是一个无监督的二分类问题，并没有任何先验信息来判断不同波段的质量。为解决该问题，本专利技术将权重学习与协同表示结合起来，构建了统一的目标函数，高光谱图像中高质量的波段能被自动赋予较大的权重，而低质量的波段被自动赋予较小的权重，使得加权后的特征更具判别性。本专利技术所解决的问题可用数学模型表示为：上式有两个未知变量∑和αi，其中∑是以自权重向量为对角元素的对角矩阵，σ的第k个元素σk表示第k个波段所对应的权重。本专利技术中正则化参数λ的可取值范围为10-6至100。3、模型求解为求解目标函数(7)，本专利技术设计了一种迭代更新算法，通过不断迭代更新权重矩阵∑以及稀疏向量αi直到目标函数收敛。具体为：(1)初始化：初始化权重矩阵，令其中I为单位阵。(2)固定∑，更新αi：当∑固定时，求解问题(7)等价于求解：将上式对αi求导然后令其等于0，可得：(3)固定αi，更新∑：当αi固定时，问题(7)可简化表示为：定义向量其第k个元素表示为mk，公式(10)可进一步转化为：上式对未知的权重向量σ有两个约束，由于约束σT1＝1，0≤σk≤1可简化为非负约束σk≥0。先只考虑第一个约束σT1＝1，问题(11)所对应拉格朗日函数可表示为：其中，γ为拉格朗日乘子。将上式对σk求导然后令其等于0，那么：由于约束σT1＝1，可得：结合公式(13)和公式(14)，有：对于任意k，σk均是非负的，正好符合问题(11)对σ的非负约束，因此公式(15)中得到的∑，就是问题(8)的最优解。<本文档来自技高网...

【技术保护点】
1.一种基于自权重协同表示的高光谱异常目标检测方法，其特征在于步骤如下：/n步骤1：对于高光谱图像的待检测像元

【技术特征摘要】
1.一种基于自权重协同表示的高光谱异常目标检测方法，其特征在于步骤如下：
步骤1：对于高光谱图像的待检测像元n为高光谱图像中的像元总数，d为高光谱图像的波段数，以yi为中心分别构建大小为wout和win的双窗，双窗之间的所有像元构成字典Yc，即c表示字典包含的像元个数，c＝wout×wout-win×win，win＜wout＜min{h，w}，h和w分别为高光谱图像的行数和列数，n＝h×w，为字典中的第j个像元，j＝1,…,c；
步骤2：构建高光谱异常检测目标函数如下：



其中，αi为像元yi异常检测待求解的期望稀疏向量，∑为权值矩阵，是以自权重向量为对角元素的对角矩阵，σk表示自权重向量σ的第k个元素，λ为正则化参数，可取值范围为10-6至100；
步骤3：采用迭代更新求解算法对目标函数进行求解，具体为...

【专利技术属性】
技术研发人员：王榕，胡豪杰，何芳，聂飞平，王震，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西;61