联合低秩张量分解与稀疏表示的高光谱异常目标检测算法_成宝芝.pdf

引用格式：成宝芝，张丽丽，赵春晖 联合低秩张量分解与稀疏表示的高光谱异常目标检测算法 电光与控制，,（）：,（）：联合低秩张量分解与稀疏表示的高光谱异常目标检测算法成宝芝，张丽丽，赵春晖（常州工学院计算机信息工程学院，江苏 常州；大庆师范学院机电工程学院，黑龙江 大庆；哈尔滨工程大学信息与通信工程学院，哈尔滨）摘 要：异常目标检测是当前高光谱图像处理中的一个研究热点。针对当前异常目标检测算法存在的问题，从解决高光谱图像中含有的背景、异常目标和噪声等相关量出发，利用高光谱图像的空间谱和光谱特性，提出了联合低秩张量分解和稀疏表示的新的高光谱图像异常目标检测算法。该算法首先利用低秩张量分解模型对高光谱进行图像恢复，使图像质量得到提升，从而使得异常目标变得突出，易于进行目标检测；然后，再利用稀疏差异指数进行异常目标检测，得到需要的异常检测结果；最后，利用真实的高光谱图像进行仿真实验，结果表明，新的异常目标检测算法具有检测精度高、虚警率低和鲁棒性好的特点。关键词：高光谱图像；异常目标检测；张量分解；稀疏表示中图分类号：文献标志码：,（,；,；,）：,：；引言高光谱是由二维的空间谱和一维的光谱构成的收稿日期：修回日期：基金项目：黑龙江省自然科学基金（）；大庆师范学院科学研究基金（）作者简介：成宝芝（），男，吉林东丰人，博士后，教授。“图谱合一”的遥感图像，高光谱图像通过极高的光谱分辨率能够对存在光谱差异的不同地物进行识别。近年来，高光谱图像在光谱解混、分类和目标检测等方面得到大量应用，其中，对不需要任何先验信息的异常目标检测的研究及其应用成为一个研究热点。从背景区域将异常目标检测和识别出来，主要是根据两个原则：一是异常目标的光谱特性明显偏离其邻域光谱；二第 卷 第 期 年 月 电 光 与 控 制 是异常目标在整个图像中是一个小概率发生事件，即异常目标存在稀疏特性。高光谱异常目标检测在军事和民用等领域都有广泛的应用，比如海上溢油检测、战场目标的准确识别以及农产品的质量检验等。总的来说，高光谱图像异常目标检测算法的研究主要分为基于统计的方法和基于表示的方法。文献提出的 异常检测方法，是基于统计的异常检测算法中的经典方法，该方法简化了高斯多元分布的背景模型，通过异常目标不同于背景特征分布这一高光谱图像特性，利用马氏距离进行异常目标检测。由于高光谱图像分布的复杂性，背景分布往往偏离高斯分布，因此，传统的 异常检测算子方法检测精度偏低。其他研究者对经典的 异常检测方法进行改进，先后提出快速 检测、全局、加权，（）等，这些算法都从不同角度对经典的 异常检测方法进行改进，提高了基于 异常检测算子检测方法的有效性和鲁棒性。近年来，基于表示的高光谱异常检测方法得到越来越多的重视。文献将稀疏表示引入高光谱异常检测中，通过光谱稀疏差异指数和空间稀疏差异指数进行加权，得到异常目标判决公式进而判断检测对象是异常目标还是背景，该方法得到了较好的异常检测结果。联合低秩矩阵和稀疏矩阵分解（，）进行异常检测的方法成为近期的研究热点，该方法的主要思想是：高光谱图像是由背景、异常目标和噪声组成。一方面，由于高光谱图像是由光谱和空间谱组成的三维图像数据，背景包含了主要信息内容，而且数百个光谱带具有互相关和均匀区域，因此，背景可以假定为低秩矩阵。另一方面，与背景信息相比，异常目标样本较少，被认为是稀疏的。高光谱图像中的噪声基本上由条纹噪声、高斯白噪声和光电噪声组成，通常可以模型化为独立分布的高斯随机变量。文献提出的 异常检测方法，用混合高斯分布代替单一分布，从而能够更好地检测和识别异常与不同类型的噪声；文献提出的 异常检测方法，将稀疏分量投影到背景正交子空间，由低秩分量估计背景正交子空间，从而抑制背景干扰、突出异常目标。本文通过对低秩矩阵联合稀疏表示方法的分析和研究，提出了联合低秩张量分解与稀疏表示（，）的高光谱异常目标检测算法，该算法首先利用低秩张量分解方法进行高光谱图像数据恢复，使得高光谱图像质量得到提升，异常目标信息能够更容易被检测到。在此基础上，引入稀疏表示方法，利用稀疏差异指数进行异常检测处理，得到异常目标检测结果。方法 张量分解理论高光谱图像为空间谱和光谱构成的三维谱。因此，高光谱图像可以利用张量进行表示和分解，如图 所示，其中，和 表示空间维，表示光谱维。图 高光谱数据张量表示和分解 首先，假设 为一个潜在干净的高光谱数据，其中，为空间谱的长和宽，为光谱维的维数。然后，利用 进行张量分解，核张量表示为 （）式中：，为核心矩阵；，和 是分解之后得到的因子矩阵。作为 的元素，表示不同成分之间的互连程度，由此，式（）可以分解为（）式中：为 的元素；，和 分别为张量分解后，和 的成分数目；，分别为，和 的元素。为了得到，式（）分解优化为 ，（），因为，和 有正交列，可以从式（）中求得，为单位矩阵。式（）的最小化问题等于如下的最大化问题，即 ，（），。低秩矩阵恢复模型低秩矩阵恢复模型（，）被引入图像处理中，用于处理数据低秩性恢复问题。给定一组图像数据样本，的目标是联合寻找所有数据的最低秩表示，该模型在低维子空间数据处理方面非常有效。假设 为一个含有噪声的图像数据，可表示为第 卷电 光 与 控 制成宝芝等：联合低秩张量分解与稀疏表示的高光谱异常目标检测算法 （）式中：为低秩矩阵；为稀疏误差矩阵。为了得到，可以将式（）转化为松弛的凸优化问题，即 ，（）式中：表示核范数；为正则化参数，用于平衡核范数和 范数之间的相对比例。对于式（）的求解，核心目的是得到恢复之后的，采用的方法为增广拉格朗日乘子算法，即交替方向方法（，）。通过对高光谱图像数据进行分析，低秩矩阵恢复模型实现了高光谱波段的矢量化，但是，高光谱图像波段间的相关性受到破坏。因此，基于高光谱图像的三维数据结构，采用张量结合低秩矩阵方法进行高光谱数据处理，能够保持高光谱图像波段间的空间结构和光谱特性的细节，使得数据处理结果更加精确。通过低秩张量分解之后的高光谱图像数据，再利用稀疏差异指数方法进行异常目标检测。由于原始高光谱图像受高光谱传感器、噪声、水汽等影响，使得图像数据产生降质。采用合适的方法对原始高光谱图像数据质量进行恢复，再进行异常目标检测，是一类具有创新性的方法。因此，高光谱图像数据恢复技术的采用非常重要，对比现有的高光谱图像数据恢复方法，采用文献提出的加权群稀疏正则化低秩张量分解（）方法进行图像数据恢复。假定高光谱图像数据为 （）式中：为含有噪声的高光谱图像；为潜在干净的高光谱图像；为除了高斯噪声之外的其他噪声构成的稀疏噪声；为高斯噪声。由此，模型表示为，（）式中：和 为正则化参数；为非负加权变量；为空间维度上的微分算子；为高斯噪声密度方差。在 模型中，充分利用了潜在干净的高光谱图像（数据集）和稀疏噪声 的先验知识，同时，低秩张量 分解能够保持空间和光谱的相关性，并利用加权群稀疏约束来保持所有波段的共享局部平滑度。稀疏分解表示及差异指数经过低秩张量模型处理过的高光谱图像数据集利用稀疏分解表示方法进行处理，得到需要的异常目标检测结果。表示为 （）式中：为 （）矩阵，是由背景 和异常目标构成的字典空间（和 分别是背景子字典矩阵和目标子字典矩阵）；表示高光谱图像的波段数，为背景的像元数，为目标的像元数；，为稀疏系数权向量（和 分别为背景系数权向量和目标系数权向量）。通过 的结果判断被检测的像元是背景还是目标，因此，的最优化问题可以表示为 （）。该最优化问题利用稀疏差异指数完成异常检测。假设 为权向量，为权向量 的维数，则光谱稀疏差异指数为，（）式中，为第 个波段的权向量。第 个波段的空间稀疏差异指数为，（）式中：是第 个波段的协同表示的权向量；是 的维数。由此空间协同的，表示为，（）式中，为波段总数。通过将光谱稀疏差异指数和空间稀疏差异指数协同处理高光谱图像，能够从高光谱图像背景信息中得到不同于背景信息分布的异常目标像元，从而检测出异常目标，得到高光谱图像处理要求的结果。算法实现步骤 算法利用低秩张量进行高光谱图像恢复，经过低秩矩阵恢复模型 方法处理后的高光谱图像，再利用稀疏差异指数进行异常目标检测，从而得到异常检测结果。具体步骤如下：）输入待检测的高光谱图像数据集；）利用基于低秩张量分解的 方法模型对原始的高光谱图像进行图像数据恢复，得到恢复之后的高光谱图像；）引入稀疏分解理论，利用稀疏差异指数对恢复后的高光谱图像进行异常目标检测；）输出高光谱图像异常目标检测结果，利用二值图像表示背景和异常信息。第 期 仿真实验和结果分析 真实高光谱图像数据为了验证 算法的异常检测性能，利用真实的高光谱图像进行仿真实验并对检测结果进行分析。圣地亚哥海军机场高光谱遥感图像被广泛用于异常目标检测，该图像来自于机载可见光 红外成像光谱仪（，），最初的波段数为 个，去除掉污染严重的波段后剩余 个，空间分辨率为，大小为 像元 像元。仿真实验中截取的图像大小为 像元 像元 像元，含有 个异常目标。如图 所示。图 仿真用的高光谱图像 仿真过程及结果分析为了验证本文所提 算法的检测率和鲁棒性，同时仿真了 算法、算法、算法和经典的 算法。算法是先利用总变分正则化低秩张量方法对原始高光谱图像进行图像提质，然后利用稀疏表示方法进行异常检测；算法是一类基于稀疏表示的异常目标检测方法；算法是一类非线性核函数 异常目标检测方法。这几类算法都具有典型性，仿真过程中使用的仿真实验软件是 ，操作系统为，电脑主频为，内存为 。对采用的 个波段的 图像先利用低秩张量方法进行图像数据恢复，结果如图 所示（以第 波段图像为例），低秩恢复之后的图像再利用稀疏差异指数进行异常目标检测，结果如图 所示。图 图像数据恢复后的高光谱图像（第 波段）（）图 种算法的异常检测结果 从检测结果能看到，算法异常目标检测精度明显高于其他 种算法。算法虽然也先进行了图像恢复，再采用稀疏表示方法进行异常检测，但算法有效性差于 算法。其他的 算法、算法和 算法异常能力明显低于 算法。图 是设定像元总数为 时得到的检测结果。为了更好地验证 算法的有效性和精确度，选取从 个像元，间隔为，即可得到总像元数变化过程中算法的检测能力，图 从检测到的目标像元数占总像元数、虚警的像元数占总像元数这两个角度比较了 种算法的检测能力。图 检测到的目标像元和虚警像元占总像元数的比较 通过图 可得，算法无论是在检测到的目标像元数还是虚警像元数方面都优于其他几种算法，说明该算法具有较好的精确度和鲁棒性。虽然 算第 卷电 光 与 控 制成宝芝等：联合低秩张量分解与稀疏表示的高光谱异常目标检测算法法也先进行了图像质量恢复，再利用稀疏表示进行异常检测，但是算法检测精度不高，主要是图像质量恢复效果一般；算法作为经典的基于稀疏表示的异常目标检测算法，具有较好的检测性能，但是由于对高光谱图像的空间谱和光谱协同处理不够，算法检测能力有一定局限性；算法、算法作为经典的非线性核函数、线性两类异常目标检测算法，主要是从光谱角度对高光谱进行异常目标检测，算法检测性能差于其他几种算法。因此，算法具有较好的算法检测性能。特性分析曲线是高光谱图像进行异常目标检测算法性能评价的重要指标之一，体现了检测概率 和虚警概率 之间的变化关系，分别定义为（）（）其中：为算法能够检测到的真实目标像元数；为被检测的高光谱图像含有的真实目标像元数；为算法检测到的高光谱图像含有的虚警像元数；为被检测的高光谱图像的总像元数。图 所示为 种算法的 性能检测对比情况。图 种算法的 性能检测比较 由图 可知，比较 算法、算法、算法、算法和 算法的 曲线能看到，算法在相同的虚警概率下，检测概率远远高于其他几种算法，而且具有一定的鲁棒性。因此，综合前文提到的 算法仿真验证结果可知，提出的算法能够较好地从高光谱图像中检测出异常目标信息，算法的检测精确度、有效性和鲁棒性都非常好。结论为了获得更好的高光谱图像异常目标检测性能，充分利用高光谱图像的光谱特性和空间特性，提出了联合低秩张量和稀疏表示（）的高光谱异常目标检测新算法。该算法利用低秩张量方法实现高光谱图像的质量恢复和提升