高光谱图像异常检测方法技术

本发明专利技术公开了高光谱图像异常检测方法

【技术实现步骤摘要】
高光谱图像异常检测方法、装置、设备及介质


[0001]本专利技术涉及高光谱异常检测领域，更具体地说，它涉及高光谱图像异常检测方法
、
装置
、
设备及介质
。

技术介绍

[0002]高光谱图像
(HSI)
利用数百个从红外到紫外连续变化的窄带光谱实现了优异的光谱分辨率，可以准确描述物体的反射光谱，从而显着提高物体的检测和分类能力
。HSI
的独特优势引领了各种应用，例如目标检测和分类
。
具体来说，根据是含有目标的确切先验信息，高光谱目标检测可大致分为两类：匹配目标和异常检测，异常检测可用于各个领域，例如精准农业
、
矿物识别和军事监视
。
[0003]目前，基于自编码器
(AE)
的检测器被提出用于高光谱异常检测
。
这类方法都假设网络获得的背景重建误差较小
。
相比之下，网络获得的异常的重建误差则较大
。
因此，重建误差被视为异常
。
为了实现端到端检测，一种基于全卷积
AE
的自动高光谱异常检测网络被提出，称为
Auto
‑
AD。
该网络采用自适应权重来减少背景重建的误差
。
然而，该自编码器无法充分学习到背景的一些潜在特性，如低秩性和局部光滑性
。
同时，不同于基于模型的检测方法，
Auto
‑
AD
>为黑盒模型，其学习到的特征是不可解释的
。
为了结合模型驱动的先验和基于
AE
网络的优点，一种结合模型驱动和数据驱动的方法被提出用于高光谱异常检测，称为
DeepLR。
它采用自适应权重的
AE
网络来重建背景，并设计了低秩约束网络的输出以利用背景的低秩属性
。
然而，
DeepLR
仅利用了背景的低秩性，忽略了背景的局部光滑性
。

技术实现思路

[0004]本申请的目的是提供一种高光谱图像异常检测方法
、
装置
、
设备及介质，本专利技术首先构建了卷积自编码器，利用卷积自编码器强大的特征学习能力提取高光谱预处理数据的特征，在卷积自编码器的输出上构建包含背景潜在特性约束的目标函数，即低秩和稀疏约束，以利用背景中的低秩和稀疏特性，进一步的，本专利技术采用高效的交替方向乘子法，对目标函数进行求解，以调整卷积自编码器的网络参数，从而使得卷积自编码器能够更加充分学习到高光谱预处理数据的背景特性，从而降低背景重建时的误差，进而提高异常检测性能
。
[0005]本申请的上述技术目的是通过以下技术方案得以实现的：
[0006]本申请的第一方面，提供了一种高光谱图像异常检测方法，方法包括：
[0007]对高光谱图像光谱维度的每个切片进行归一化处理，得到高光谱预处理数据；
[0008]将所述高光预处理谱数据输入预先构建好的具有跨连接层的卷积自编码器中，以对高光谱的背景进行重构，输出高光谱预处理数据的背景重构张量；
[0009]沿着光谱维度将所述背景重构张量拉伸为矩阵，对所述矩阵施加使卷积自编码器的网络输出最优的目标函数，其中目标函数由自编码器残差约束项和背景潜在特性的低秩约束项和稀疏约束项构成；
[0010]采用交替方向乘子法迭代更新所述目标函数的低秩约束项
、
稀疏约束项和残差约束项，以调整所述卷积自编码器的网络参数，在迭代次数达到收敛阈值时，确定当前所述卷积自编码器输出的背景重构张量；
[0011]将所述高光谱预处理数据减去当前所述卷积自编码器输出的背景重构张量，得到背景重构后的误差张量；
[0012]对误差张量每个光谱维度的切片求平方后累加开方，得到检测图；
[0013]对所述检测图进行异常检测，得到异常检测结果
。
[0014]在一种实现方案中，所述预先构建好的具有跨连接层的卷积自编码器包括编码器和解码器，编码器和解码器均由卷积层构成，编码器的卷积层输出通过跨连接层与对应的解码器的卷积层输出拼接
。
[0015]在一种实现方案中，所述目标函数的表达式为：其中，表示高光谱预处理数据张量，
w
表示立方体数据的宽，
h
表示立方体数据的高，
b
表示立方体数据的光谱带宽维度，
||
·
||
F
表示
F
范数，
||
·
||
*
表示矩阵的核范数，表示差分算子，
||
·
||1表示
l1范数，
λ
表示用于平衡低秩约束的非负权重参数，
β
表示用于平衡稀疏约束的非负权重参数，表示自编码器的输出，
θ
表示自编码器的网络参数，表示将自编码器输出的背景重构张量沿着光谱维度拉伸成矩阵，表示对执行一阶差分运算
。
[0016]在一种实现方案中，采用交替方向乘子法迭代更新所述目标函数含网络参数项和所述目标函数的低秩约束项和稀疏约束项，包括：
[0017]引入第一松驰变量和第二松驰变量对所述目标函数的低秩约束项和稀疏约束项进行替换，得到替换后的目标函数，替换后的目标函数的表达式为：其中，
C1表示第一松驰变量，
C2表示第二松驰变量；
[0018]利用拉格朗日乘子法构建变量替换后的目标函数，构建的目标函数的表达式为：其中
,P1和
P2均为拉格朗日乘子，
μ
为拉格朗日权重系数；
[0019]依据重构的目标函数分别更新第一松弛变量和第二松弛变量，以分别计算出更新后第一松弛变量和第二松弛变量的解；
[0020]采用
Adam
算法，对代入更新后第一松弛变量和第二松弛变量的解的含有网络参数项的代价函数项进行求解，得到使卷积自编码器代价函数最小化的网络参数
。
[0021]在一种实现方案中，在迭代次数达到收敛阈值时，确定当前所述卷积自编码器输出的背景重构张量，具体为：
[0022]预设收敛阈值，计算该迭代次数下含有网络参数项的代价函数值；
[0023]在所述代价函数值小于收敛阈值时，停止迭代，确定当前所述卷积自编码器输出的背景重构张量
。
[0024]在一种实现方案中，所述代价函数值小于收敛阈值的表达式为：
k
表示迭代次数，
θ
k+1
表示第
k+1
次迭代跟新的网络参数，
θ
k
表示第
k
次迭代跟新的网络参数，表示自编码器代价函数在第
k+1
次的值，表示自编码器代价函数在第
k
次本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种高光谱图像异常检测方法，其特征在于，方法包括：对高光谱图像每个光谱维度的切片进行归一化处理，得到高光谱预处理数据；将所述高光谱预处理数据输入预先构建好的具有跨连接层的卷积自编码器中，以对高光谱的背景进行重构，输出高光谱预处理数据的背景重构张量；沿着光谱维度将所述背景重构张量拉伸为矩阵，对所述矩阵施加使卷积自编码器的网络输出最优的目标函数，其中目标函数由背景潜在特性的低秩约束项和稀疏约束项及自编码器的残差项构成；采用交替方向乘子法迭代更新所述目标函数的低秩约束项和稀疏约束项，以调整所述卷积自编码器的网络参数，在迭代次数达到收敛阈值时，确定当前所述卷积自编码器输出的背景重构张量；将所述高光谱预处理数据减去当前所述卷积自编码器输出的背景重构张量，得到背景重构后的误差张量；对所述误差张量每个光谱维度的切片求平方后累加开方，得到检测图；对所述检测图进行异常检测，得到异常检测结果
。2.
根据权利要求1所述的高光谱图像异常检测方法，其特征在于，所述预先构建好的具有跨连接层的卷积自编码器包括编码器和解码器，编码器和解码器均由卷积层构成，编码器的卷积层输出通过跨连接层与对应的解码器的卷积层输出拼接
。3.
根据权利要求1所述的高光谱图像异常检测方法，其特征在于，所述目标函数的表达式为：其中，表示高光谱预处理数据张量，
w
表示立方体数据的宽，
h
表示立方体数据的高，
b
表示立方体数据的光谱维度，｜｜
·
｜｜
F
表示
F
范数，｜｜
·
｜｜
*
表示矩阵的核范数，表示差分算子，
||
·
||1表示
l1范数，
λ
表示用于平衡低秩约束的非负权重参数，
β
表示用于平衡稀疏约束的非负权重参数，表示自编码器的输出，
θ
表示自编码器的网络参数，表示将自编码器输出的背景重构张量沿着光谱维度拉伸成矩阵，表示对执行一阶差分运算
。4.
根据权利要求3所述的高光谱图像异常检测方法，其特征在于，采用交替方向乘子法迭代更新所述网络参数和所述目标函数的低秩约束项和稀疏约束项，包括：引入第一松驰变量和第二松驰变量对所述目标函数的低秩约束项和稀疏约束项进行替换，得到替换后的目标函数，替换后的目标函数的表达式为：其中，
C1表示第一松驰变量，
C2表示第二松驰变量；利用拉格朗日乘子法对替换后的目标函数进行重构，得到重构后的目标函数，重构的目标函数的表达式为：其中
,P1和
P2均为拉格朗日乘子，
μ
为拉格朗日权重系数；依据重构后的目标函数分别更新第一松弛变量和第二松弛变量，以分别计算出更新后第一松弛变量和第二松弛变量的解；采用
Adam
算法，对代入...

【专利技术属性】
技术研发人员：舒勤，冯茂原，王畅，崔忠马，
申请(专利权)人：四川大学，
类型：发明
国别省市：