基于像素选择过程的高光谱图像异常目标检测方法
【专利摘要】本发明提供一种基于像素选择过程的高光谱图像异常目标检测方法，包括以下步骤：1】对高光谱图像进行归一化预处理，并设定检测窗口尺寸；2】采用滑动窗口技术，获得待检测高光谱图像块X，X∈RN×D；其中，R是实数空间，N是图像块中的像素点数目，D是高光谱数据的光谱波段数；3】计算每个像素点Xi的重构误差εi；4】构建顶点-边缘权重图；5】构建像素选择过程模型；6】重复步骤2】-步骤5】，直至完成整幅高光谱图像异常目标的检测，获得最终的异常概率图。该高光谱图像异常目标的检测方法，解决了现有的检测技术需要对背景分布进行假设，适应性差，背景与目标的差异性不足，检测能力低的技术问题。
【专利说明】
基于像素选择过程的高光谱图像异常目标检测方法
技术领域
[0001] 本发明属于图像处理技术领域，具体设及一种针对高光谱图像进行异常目标检测 的方法。
【背景技术】
[0002] 高光谱图像包含丰富的光谱信息，能够为地物目标微小差异的识别提供辨别性线 索。而高光谱图像的异常目标检测正是基于此特性的典型应用，其主要目的是识别图像场 景中显著偏离背景光谱特性的目标。运种检测实质是一种二分类问题，即将待观测的像素 分类为背景或目标。与有监督的目标检测问题不同，高光谱的异常目标检测没有任何关于 目标或背景的先验光谱信息，其仅仅依靠对所设定的参考背景进行建模，寻找不满足背景 模型分布的点，进而定义为异常。由于运样的技术特性相比于有监督的目标检测更符合实 际情况，因而已经被应用在很多领域，如精准农业，公共安全，智能防御等。
[0003] 目前，根据对背景分布不同的假设，高光谱图像的异常检测算法可W分为两类。
[0004] 一是基于均匀光谱分布假设的异常检测。运种方法假设整个背景服从均匀的 光谱分布，具有代表性的算法就是由Reed和化等在文献"I. Reed and X.化.Adaptive Multiple-band CFA民 Detection of An Optical Pattern with Unknown Spectral Distribution. IEEE Transactions on Acoustics，Speech and Signal Process! ng，38 (10) : 1760 - 1770, 1990"提出的RX算法。该方法假设背景服从相同的多元正态分布， 利用马氏距离衡量待检测像素与背景光谱分布的偏离程度。该方法存在两点不足：第一点 是高斯分布不能准确地描述真实的高光谱数据，其主要原因来源于由地表材料的多样性和 信号获取过程的干扰性所造成的高光谱图像复杂的分布特性；第二点是小样本问题，利用 较小数量的样本进行高维背景协方差矩阵的计算，会导致病态矩阵处理问题。研究者也相 继提出很多方法克服RX算法存在的问题。一类是RX算法的变种，但是由于其本质的假设 没有改变，所W固有的问题还是存在；另一类是基于核理论的算法，将数据映射到高维特征 空间增强目标与背景的辨别性，但是较高的计算负担依然是该类核算法难W突破的瓶颈。 阳0化]二是基于多类光谱分布假设的异常检测。运种方法假设复杂的背景包含多类不同 的分布。典型的技术路线是先分别评估不同类别的信息，然后综合分析，其代表性的算法 Carlotto 在文献"M. J. Carlotto. A Ccluster-based Approach for Detecting Man-made Objects and Changes in Imagery. IEEE Transactions on Geoscience and Remote Sensing, 43(2) :374 - 387, 2005"提出的基于聚类的检测方法。该算法首先利用聚类算法 将高光谱图分割成几类，然后在每一类利用原始的RX算法实现异常目标的检测。运类算法 的检测性能依赖于聚类数目的设定：当聚类数目低于真实数据聚类数目，会导致类间的差 异性降低，检测性能下降；当聚类数目高于真实数据聚类数目，会导致异常目标被误分为一 类，造成不可检测。
【发明内容】
[0006]本发明目的是提供一种高光谱图像异常目标的检测方法，解决了现有的检测技术 需要对背景分布进行假设，适应性差，背景与目标的差异性不足，检测能力低的技术问题。 阳007] 本发明的技术解决方案是：
[0008] 一种基于像素选择过程的高光谱图像异常目标检测方法，包括W下步骤：
[0009] 1】对高光谱图像进行归一化预处理，并设定检测窗口尺寸；
[0010] 2】采用滑动窗口技术，获得待检测高光谱图像块X，X e 其中，R是实数空间， Ν是图像块中的像素点数目，D是高光谱数据的光谱波段数；
[0011] 3】计算每个像素点Xi的重构误差ε 1:
[0012] 3. 1】在高维数据空间，寻找图像块X中每个像素点Xi的Κ个最近邻；
[0013] 3. 2】进行局部线性表示，最小化重构误差
计算得到线性 权重系数
[0014] 3. 3】保持局部嵌入关系不变，通过线性变换，将高维空间数据Xi映射到低维向量 空间YiG Rd;最小化损失函数
计算得到每个像素点的重构误差 ε 1; 阳01引 4】构建顶点-边缘权重图：
[0016] 将图像块Xi中的每一个像素点作为图模型中的顶点V，将重构误差ε 1作为顶点 权重，利用公式
计算得到图像块中任意两点之间的欧式距离a，，作 为连接它们的边缘权重，完成图模型的构建；
[0017] 5】构建像素选择过程模型：
[0018] 5. 1】利用顶点权重和边缘权重值计算相似矩阵i = n ?Π ，其中Π =diag( ε );构 建像素选择目标函数：
[0019]
[0020] S. t. Ρ e Δ
[0021] 其中Ρ是异常概率向量，它的每一个元素表示图像块每一中像素点的异常概率， 满足 A = {Ρ > 0,忡=U ;
[0022] 5. 2】利用复制动态算法求解像素选择目标函数；给定初始值
通过一阶 复制方程的离散时间形式，即公式
[0023]
[0024] 进行迭代计算，当迭代收敛时，近似的稳定解等于局部解；
[00巧]5. 3】求解得到向量P，其包含的每个元素对应所处理的图像块中每个像素点的异 常概率。
[00%] 6】重复步骤2】-步骤5】，直至完成整幅高光谱图像异常目标的检测，获得最终的 异常概率图。
[0027] 步骤1中所述的检测窗口包括内窗口和外窗口；所述内窗口为保护窗口，用于减 少异常目标的干扰；所述外窗口用于限定检测范围。
[0028] 本发明的优点在于：
[0029] (1)本发明将流形学习技术与图理论相结合，避免了对背景数据分布进行假设，提 局了算法应用性；
[0030] (2)本发明中构建的顶点-权重图，充分考虑了每个像素的特性W及与其他所有 像素的关系，能够有效抑制数据中噪声的干扰，提高异常目标的检测能力。
【附图说明】
[0031] 图1为本发明基于像素选择过程的高光谱图像异常目标检测方法的流程图。
【具体实施方式】
[0032] 参照图1，本发明的实现方法主要如下：
[003引 （1)采用滑窗策略，提取每次处理的图像块。
[0034] (la)首先对获得的高光谱图像进行归一化预处理。然后，采用内外窗口，（内部 窗口作为保护窗口，减少异常目标的干扰；外部窗口，限定检测范围），提取高光谱图像块 X e RWXD，其中N是图像块中的像素点数目，D是高光谱数据的光谱波段数。
[00对 似采用流形学习技术，计算每个像素点的重构误差。
[0036] 姑）在高维空间，寻找图像块X中每个像素 Xi的K个最近邻，并进行局部线性表 示，最小化重构误差
可W计算得到线性表示系数wf
[0037] (2b)保持局部嵌入关系不变，通过线性变换，将高维空间数据Xi映射到低维向量 空间Rd。最小化损失函I
可W计算得到每个像素点的重构 误差。
[0038] (3)构建顶点-边缘权重图。
[0039] 将步骤（2b)中计算的重构误差ε，作为顶点权重。计算图像块中所有点之间的欧 式距离，作为图中的边缘连接两个顶点。利用
得到最终的边缘对称权 重矩阵Α。 W40] (4)异常像素选择过程，获得异常概率图，计算AUC值。
[0041] (4a)利用顶点权重和边缘权重值计算相似矩阵/'?二rl·们（Π = diag(<'-'））。与子集 选择类似，构建像素选择目标函数为：
[0042]
[0043] s. t. P e Δ，
[0044] 其中P是异常概率向量，它的每一个元素表示图像块每一中像素点的异常概率， 其满足A =巧> 0, iTp = 1}。最大化目标函数，实际上对应的是为具有较大的顶点权重 和边缘权重的顶点分配较大的概率，即越是异常的像素点将W较大的概率从所构建的图关 系中挑选出来。
[0045] (4b)利用复制动态算法求解像素选择目标函数。给定初始值
通过一阶 复制方程的离散时间形式，局部解可W通过
[0046]
[0047] 得到。当迭代收敛时候，近似的稳定解等于局部解。最终从求解得到的向量P的 每个元素对应图像块中每个像素的异常概率，本发明从P中选出与滑窗中屯、待检测像素所 对应的元素，作为该高光谱像素的异常概率。
[0048] 本发明的效果可W通过W下仿真实验做进一步的说明。 W例 1.仿真条件 阳化0] 本发明是在中央处理器为Intel佩Core i3-21303.40GHZ、内存16G、WINDOWS 8 操作系统上，运用MATLAB软件进行的仿真。
[0051] 实验中使用的实验数据是从美国陆军工程研究和开发中屯、的网站下载的HYDICE urban数据。该数据的光谱范围为400-2500皿，包括210个波段。实验中去除了低信噪比 和水汽吸收波段，保留了 160个光谱波段。同时，从原始图像截取80X100的子图进行异常 目标的检测。 阳化引 2.仿真内容
[0053] 按如下步骤用本发明方法进行人高光谱图像的异常检测：
[0054] 首先，在HYDICE urban数据，完成本发明算法（基于图像素选择）的实验。为 了证明算法的有效性，综合考虑算法的流形性、崭新性和多样性，选择了 5个基准算 法 RX、GRX、RSAD、WSCF、SV孤进行比较。其中，GRX 在文献"Z.化an,比 Sun, K. Ji, Z. Li, and Η. Zou. Local Sparsity Divergence for Hyperspectral Anomaly Detection. IEEE Geoscience and Remote SensingLetters, 11 (10):1697 - 1701，2014."有详 细的介绍；RSAD 算法是文献"B. Du and L Zhang. Random-selection-based Anomaly Detector for Hyperspectral Imagery.IEEE Transactions on Geoscience and Remote Sensing, 49 巧）：1578 - 1589, 2011，，提出的；啡5〔尸是文献"^-]\16过11〇61，]\1 Guillaume,and S. Bourennane. Whitening Spatial Correlation Filtering for Hyperspectral Anomaly Detection. In Proceedings of IEEE International Conference on Acoustics, Speech, and Signal Processing, 2005:333 - 336."提出的；SVDD 是文 南犬"A. Banerjee, P. Burlina,and C. Diehl.A Support Vector Method for Anomaly Detection in Hyperspectral Imagery. IEEE Transactions on Geoscienceand 民emote Sensing, 44巧）：2282 - 2291，2006"提出的。
[0055] 其次，将实验中每种算法的检测结果与地面真值进行比较，通过计算检测率和虚 警率，绘制接受者操作特征曲线，进而计算接受者特征操作曲线下的面积值（AUC)，作为检 测性能的衡量指标，结果如表1所示。
[0056] 表1高光谱异常检测结果
[0057]
[005引从表1可见，本发明的检测AUC值显著地高于另外5种基准算法。本发明由于将 流形学习技术与图理论相结合，避免了对背景数据分布进行假设，提高了算法应用性；本发 明构建的顶点-权重图，充分考虑了每个像素的特性W及与其他所有像素的关系，能够有 效抑制数据中噪声的干扰，提高异常目标的检测能力，因此获得了更好的检测结果，进一步 验证了本发明的先进性。
【主权项】
1. 一种基于像素选择过程的高光谱图像异常目标检测方法，其特征在于：包括W下步 骤： 1】对高光谱图像进行归一化预处理，并设定检测窗口尺寸； 2】采用滑动窗口技术，获得待检测高光谱图像块X，X G 其中，R是实数空间，N是 图像块中的像素点数目，D是高光谱数据的光谱波段数； 3】计算每个像素点Xi的重构误差e 1: 3. 1】在高维数据空间，寻找图像块X中每个像素点Xi的K个最近邻； 3. 2】进行局部线性表示，最小化重构误！计算得到线性权重 系数W，； 3. 3】保持局部嵌入关系不变，通过线性变换，将高维空间数据Xi映射到低维向量空间YiG Rd;最小化损失函数 计算得到每个像素点的重构误差e 1; ， 4】构建顶点-边缘权重图： 将图像块X;中的每一个像素点作为图模型中的顶点V，将重构误差e 1作为顶点权重， 利用公式计算得到图像块中任意两点之间的欧式距离曰1，，作为连接 它们的边缘权重，完成图模型的构建； 5】构建像素选择过程模型： 5. 1】利用顶点权重和边缘权重值计算相似矩闻；构建 像素选择目标函数：S. t. P G A 其中P是异常概率向量，它的每一个元素表示图像块每一中像素点的异常概率，满足 A = {P > 0, fp = 1}; 5. 2】利用复制动态算法求解像素选择目标函数；给定初始通过一阶复制 方程的离散时间形式，即公式进行迭代计算，当迭代收敛时，近似的稳定解等于局部解； 5. 3】求解得到向量P，其包含的每个元素对应所处理的图像块中每个像素点的异常概 率。 6】重复步骤2】-步骤5】，直至完成整幅高光谱图像异常目标的检测，获得最终的异常 概率图。2. 根据权利要求1所述的基于像素选择过程的高光谱图像异常目标检测方法，其特征 在于：步骤I中所述的检测窗口包括内窗口和外窗口；所述内窗口为保护窗口，用于减少异 常目标的干扰；所述外窗口用于限定检测范围。
【文档编号】G06T7/00GK105989597SQ201510079650
【公开日】2016年10月5日
【申请日】2015年2月13日
【发明人】袁媛, 王 琦, 马单丹
【申请人】中国科学院西安光学精密机械研究所