高光谱图像目标异常检测方法、系统、终端及存储介质与流程

1.本发明属于遥感图像处理领域，具体涉及一种高光谱图像目标异常检测方法、系统、终端及存储介质。


背景技术：

2.高光谱图像拥有众多的光谱域，可以提供地球表面的大量光谱信息和空间信息，由于高光谱图像的良好检测性能，它被广泛的应用于场景分类、混合、聚类、变化检测和目标检测等领域。高光谱图像的应用之一—高光谱异常检测已经是当前研究的热点之一，目前高光谱异常检测被广泛应用于军事和民用领域，并已经取得了不错的效果。高光谱异常检测是一种无监督的二分类问题，不需要先验知识。在高光谱异常检测中，将光谱与其周围的光谱显著不同的对象视为异常，将整个图像中除异常以外的对象定义为背景。
3.当前已有的很多方法都可以完成高光谱异常检测。reed
‑
xiaoli(rx)算法基于像素与背景之间的mahalanobis距离，并假设使用多元高斯背景分布来识别异常，详细的又分为使用整个图像建模背景的grx和使用局部双窗口进行建模的lrx。后来又有学者提出了非线性版本的rx(krx)，krx通过非线性映射函数扩展到与原始输入空间相关的特征空间。由于krx存在目标污染背景的问题，又有学者基于此问题对算法进行了改进，提出了在内核空间中使用鲁棒回归分析的鲁棒性非线性异常检测算法(rnad)。基于协作表示的检测器(crd)在后来被提出并被广泛使用。crd假设背景中的每个像素都可以使用其邻域像素进行近似，而异常像素则不能。假设背景像素是相邻像素的线性组合，并且通过最小化范数来增强协作表示。当前已有基于rx的检测器都是假定背景是均匀的，并使用多元高斯分布进行建模。但是实际情况中往往无法满足此假设，这将导致准确率下降。
4.除此以外，由于高光谱图像在获取时并不能完全屏蔽噪声和异常的干扰，因此在高光谱图像中某个波段可能会存在噪声和异常像素点。crd以及其派生算法都是假设每个波段都具有相同的重要性，这样的假设对于噪声波段的鲁棒性将会很差，从而降低异常检测的效果。
5.相关文献：
6.【1】reed i s,yu x.adaptive multiple
‑
band cfar detection of an optical pattern with unknown spectral distribution[j].ieee transactions on acoustics,speech,and signal processing,1990,38(10):1760
‑
1770.
[0007]
【2】molero j m,garz
ó
n e m,garc
í
a i,et al.analysis and optimizations of global and local versions of the rx algorithm for anomaly detection in hyperspectral data[j].ieee journal of selected topics in applied earth observations and remote sensing,2013,6(2):801
‑
814.
[0008]
【3】kwon h,nasrabadi n m.kernel rx
‑
algorithm:anonlinear anomaly detector for hyperspectral imagery[j].ieee transactions on geoscience and remote sensing,2005,43(2):388
‑
397.
[0009]
【4】zhao r,du b,zhang l.arobust nonlinear hyperspectral anomaly detection approach[j].ieee journal of selected topics in applied earth observations and remote sensing,2014,7(4):1227
‑
1234.
[0010]
【5】li w,du q.collaborative representation for hyperspectral anomaly detection[j].ieee transactions on geoscience and remote sensing,2014,53(3):1463
‑
1474.


技术实现要素：

[0011]
本发明的目的在于针对上述现有技术中的问题，提供一种高光谱图像目标异常检测方法、系统、终端及存储介质，能够有效避免噪声和异常波段对于高光谱异常检测的影响。
[0012]
为了实现上述目的，本发明有如下的技术方案：
[0013]
一种高光谱图像目标异常检测方法，包括以下步骤：
[0014]
通过最大互相关熵协作表示赋予高光谱图像每个波段一个合适的权重；
[0015]
对于每个波段具有一个合适的权重的高光谱图像，通过最大互相关熵协作表示得出像素异常得分，通过设置合适阈值选择出异常像素。
[0016]
作为本发明高光谱图像目标异常检测方法的一种优选方案，所述的高光谱图像定义为表示具有n个像素的高光谱图像，维度为d维；使用滑动双窗口来进行局部近似，外窗和内窗的尺寸分别用w
out
和w
in
来表示，因此，所在邻域的像素用来表示，因此，所在邻域的像素用来表示，其中s＝w
out
×
w
out
‑
w
in
×
w
in
。
[0017]
作为本发明高光谱图像目标异常检测方法的一种优选方案，所述高光谱图像的最大互相关熵协作表示目标函数如下：
[0018][0019]
其中，表示所使用的每个邻域像素的权重，λ为正则化参数；
[0020]
式中，越大代表协作表示的差异越大，则越小。
[0021]
作为本发明高光谱图像目标异常检测方法的一种优选方案，将所述最大互相关熵协作表示目标函数通过半二次技术进行重构得到如下表达式：
[0022][0023]
式中，是一个对角矩阵，对角线上的每一个元素代表该波段的权重；
[0024]
w
k
的定义如下：
[0025][0026]
作为本发明高光谱图像目标异常检测方法的一种优选方案，与中心像素差异较大的像素点对应较小的系数，利用对角正则化矩阵求解
[0027]
其中γ
i
为tikhonov正则化矩阵，表达式如下：
[0028][0029]
作为本发明高光谱图像目标异常检测方法的一种优选方案，采用迭代优化的方式对提出的目标函数进行优化，完成无监督高光谱图像目标异常检测。
[0030]
本发明还提出一种高光谱图像目标异常检测系统，包括：
[0031]
光谱波段加权模块，用于通过最大互相关熵协作表示赋予高光谱图像每个波段一个合适的权重；
[0032]
异常像素筛选模块，用于对于每个波段具有一个合适的权重的高光谱图像，通过最大互相关熵协作表示得出像素异常得分，通过设置合适阈值选择出异常像素。
[0033]
本发明还提出一种终端设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述的处理器执行所述的计算机程序时实现所述的高光谱图像目标异常检测方法的步骤。
[0034]
本发明还提出一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述的计算机程序被处理器执行时实现所述的高光谱图像目标异常检测方法的步骤。
[0035]
相较于现有技术，本发明至少具有以下有益效果：
[0036]
在高光谱图像目标异常检测过程中，通过对高光谱图像的不同波段自适应赋予不同权重，对于每个波段具有一个合适的权重的高光谱图像，通过最大互相关熵协作表示得出像素异常得分，通过设置合适阈值选择出异常像素，能有效避免波段噪声与异常像素点段对于检测准确度的影响，提高算法的鲁棒性。本发明提出的方法与当前主流高光谱图像目标异常检测方法相比具有明显优势，在公开数据集上的准确率都大幅领先其它方法。
[0037]
进一步的，本发明目标函数构建的是联合函数，同时求解波段权重和邻域像素点权重，使用迭代求解算法高效求解，避免了单独学习二者其中一个可能带来的次优解问题。
附图说明
[0038]
图1本发明与其它方法在aviris
‑
i数据集上roc曲线对比图；
[0039]
图2本发明与其它方法在abu
‑
airport
‑
1数据集上roc曲线对比图；
[0040]
图3 aviris
‑
i数据集伪彩色图；
[0041]
图4 aviris
‑
i数据集标签图；
[0042]
图5 aviris
‑
i数据集crd检测结果图；
[0043]
图6 aviris
‑
i数据集本发明检测结果图；
[0044]
图7 aviris
‑
ii数据集伪彩色图；
[0045]
图8 aviris
‑
ii数据集标签图；
[0046]
图9 aviris
‑
ii数据集crd检测结果图；
[0047]
图10 aviris
‑
ii数据集本发明检测结果图。
具体实施方式
[0048]
下面结合附图及实施例对本发明做进一步的详细说明。
[0049]
假设高光谱图像定义为表示具有n个像素的高光谱图像，维度为d维；使用滑动双窗口来进行局部近似，外窗和内窗的尺寸分别用w
out
和w
in
来表示，因此，所在邻域的像素用来表示，其中s＝w
out
×
w
out
‑
w
in
×
w
in
。
[0050]
s1、通过最大互相关熵协作表示赋予高光谱图像每个波段一个合适的权重；
[0051]
s2、对于每个波段具有一个合适的权重的高光谱图像，通过最大互相关熵协作表示得出像素异常得分，通过设置合适阈值选择出异常像素。
[0052]
本发明提出的高光谱图像目标异常检测方法基于信息理论学习(itl)中提出的最大熵理论(mcc)来构建检测器。itl提供了一种统一的方法来提高针对异常值鲁棒性的机器学习方法。mcc可以从数据中获得更多信息以进行调整，因此在处理非高斯信号时可以具有更大的优势。熵引入了一个新的度量，当两点接近时，该度量等于范数。随着两点逐渐远离，度量将接近范数，最终接近范数。因此，mcc对于异常值有着强鲁棒性。
[0053]
所述高光谱图像的最大互相关熵协作表示目标函数如下：
[0054][0055]
其中，表示所使用的每个邻域像素的权重，λ为正则化参数；
[0056]
式中，越大代表协作表示的差异越大，则越小。
[0057]
将所述最大互相关熵协作表示目标函数通过半二次技术进行重构得到如下表达式：
[0058][0059]
式中，是一个对角矩阵，对角线上的每一个元素代表该波段的权重；
[0060]
w
k
的定义如下：
[0061][0062]
与中心像素差异较大的像素点对应较小的系数，利用对角正则化矩阵求解与中心像素差异较大的像素点对应较小的系数，利用对角正则化矩阵求解其中γ
i
为tikhonov正则化矩阵，表达式如下：
[0063][0064]
采用迭代优化的方式对提出的目标函数进行优化，完成无监督高光谱图像目标异常检测。
[0065]
本发明提出的高光谱图像目标异常检测方法，能够有效避免噪声和异常波段对于高光谱异常检测的影响。通过本发明的方法给予噪声波段足够小的权重以至于不会影响推理结果，而给予贡献大的波段更多的权重，以便它能够提供更多的有效信息。
[0066]
本发明还提出一种高光谱图像目标异常检测系统，包括：
[0067]
光谱波段加权模块，用于通过最大互相关熵协作表示赋予高光谱图像每个波段一个合适的权重；
[0068]
异常像素筛选模块，用于对于每个波段具有一个合适的权重的高光谱图像，通过最大互相关熵协作表示得出像素异常得分，通过设置合适阈值选择出异常像素。
[0069]
本发明还提出一种终端设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述的处理器执行所述的计算机程序时实现所述的高光谱图像目标异常检测方法的步骤。
[0070]
本发明还提出一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述的计算机程序被处理器执行时实现所述的高光谱图像目标异常检测方法的步骤。
[0071]
所述的计算机程序可以被分割成一个或多个模块/单元，所述一个或者多个模块/单元被存储在所述存储器中，并由所述处理器执行，以完成本发明高光谱图像目标异常检测方法。
[0072]
所述终端可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备，也可以是处理器、存储器。处理器可以是中央处理单元(centralprocessingunit，cpu)，还可以是其他通用处理器、数字信号处理器(digitalsignalprocessor，dsp)、专用集成电路(applicationspecificintegratedcircuit，asic)、现成可编程门阵列(field
‑
programmablegatearray，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。存储器可用于存储计算机程序和/或模块，处理器通过运行或执行存储在所述存储器内的计算机程序和/或模块，以及调用存储在存储器内的数据，实现高光谱图像目标异常检测系统的各种功能。
[0073]
实施例
[0074]
一种高光谱图像目标异常检测方法，包括以下步骤：
[0075]
步骤一：输入高光谱图像正则项参数λ，σ以及滑动窗口大小w
out
和w
in
。
[0076]
步骤二：对于目标函数进行求解，固定w
i
(赋予初始值i)更新α
i
，目标函数将转化为：
[0077]
对于上述公式令其对α
i
求导并置导数为0，可以得到：
[0078]
根据此公式可对α
i
进行更新。
[0079]
步骤三：固定α
i
，通过公式更新w
i
。
[0080]
步骤四：重复进行步骤二和步骤三，直至收敛。
[0081]
步骤五：移动滑动窗口，对下一个像素进行处理。对其进行步骤二、步骤三和步骤四的操作。
[0082]
步骤六：通过得到高光谱全像素异常分数，设置阈值得到异常像素。
[0083]
表1是本发明基于最大熵准则的协作表示方法进行高光谱异常检测在公开数据集上与当前主流方法的auc值对比，其中，auc值越大越好。
[0084]
表1
[0085]
d
a
tasetgrxlrxcrdlsmadswcrdercrdoursaviris
‑
i0.91110.85840.99110.97580.99360.96780.9985aviris
‑
ii0.94030.84660.96960.97030.97980.96140.9855aviris
‑
iii0.87710.84110.98280.93180.98370.93490.9904abu
‑
airport
‑
10.82210.88860.9230.88060.92730.8970.9631abu
‑
airport
‑
20.84040.97060.97350.93150.96920.96930.9879salinas0.88720.91720.96060.93510.9510.94470.9643
[0086]
图1和图2是本发明与当前主流方法在aviris
‑
i和abu
‑
airport
‑
1数据集上的roc曲线对比图，图中的曲线越靠近左上角效果越好。图3，图4，图5和图6分别为aviris
‑
i数据集的伪彩色图，标签图，crd的检测效果图和本发明的检测效果图。图6，图7，图8和图9分别为aviris
‑
ii数据集的伪彩色图，标签图，crd检测效果图和本发明的检测效果图，这些检测图中像素越亮说明该像素为异常点的可能性越高。对比可见，本发明提出的方法与当前主流高光谱图像目标异常检测方法相比具有明显优势，在公开数据集(aviris
‑
i数据集以及aviris
‑
ii数据集)上的准确率都大幅领先于其它方法。
[0087]
以上所述的仅仅是本发明的较佳实施例，并不用以对本发明的技术方案进行任何限制，本领域技术人员应当理解的是，在不脱离本发明精神和原则的前提下，该技术方案还可以进行若干简单的修改和替换，这些修改和替换也均属于权利要求书所涵盖的保护范围之内。