基于鲁棒背景估计方法的高光谱图像局部目标检测方法
【专利摘要】本发明公开了一种基于鲁棒背景估计方法的高光谱图像局部目标检测方法，用于解决现有高光谱图像局部目标检测方法虚警率高的技术问题。技术方案是首先利用基于光谱角的聚类方法，在白化空间中对输入图像进行聚类。在检测过程中，引入MCD估计方法估计待测像元所属聚类的背景参数，从而提高检测性能。在RIT提供的数据集上目标测试结果显示，反映虚警数的平均得分为2.8，较Halper的改进方法降低了4.4。卫星拍摄的AVIRIS数据集上的测试结果表明，在100％的检测率下，虚警率为0.11％，而Halper的方法虚警率为0.29％左右，全局方法的虚警率为0.82％；可见本发明方法明显降低了虚警率。
【专利说明】基于鲁棒背景估计方法的高光谱图像局部目标检测方法
【技术领域】
[0001] 本发明涉及一种高光谱图像局部目标检测方法，特别涉及一种基于鲁棒背景估计 方法的高光谱图像局部目标检测方法。
【背景技术】
[0002] 高光谱图像是利用成像光谱仪将视场中观测到的各种地物的光谱信息记录下来 得到影像数据，它结合了地面物质的空间和光谱信息，以精细的光谱分辨率为地物的分类 和检测提供了可靠依据。
[0003] 大部分的检测方法通过计算待测像元与背景均值的距离，或者通过计算待测像元 与已知的目标特征光谱之间的相似程度，将目标从背景中区分出来。自适应余弦一致性评 估器算法（Adaptive Cosine/Coherence Estimator, ACE)等传统的非结构化目标检测算法 是一类全局检测算法。其通过直接对整幅高光谱图像计算背景均值和协方差矩阵，实现对 图像背景参数的估计。当场景复杂时，图像中大量存在的异常点对背景造成污染，导致全局 背景均一性的下降，从而破坏了全局检测算法的检测性能。
[0004] 近年来，通过对图像进行分类或利用滑动窗口构造相对较均匀的局部背景，再对 局部背景使用传统检测算法进行检测的局部方法越来越受到人们的关注。但现有的多数局 部检测算法，仍然沿用全局算法用来估计背景参数的经验估计方法。该方法容易受到异常 数据的干扰，因此影响局部检测算法的鲁棒性。
[0005] 针对这一问题，近来也有相关的研究提出了一些改进方法，例如Marin S. Halper 在文献"Global, local, and stochastic background modeling for target detection in mixed pixels, 2010,（7695) :769527-769527-13"中，利用全局RXD(Reed-X Detector)和全 局ACE算法估计数据中的异常点和感兴趣目标，形成掩膜，在估计局部背景参数时利用掩 膜将异常点和疑似目标排除在外。然而，该方法受全局算法对异常点和感兴趣目标的估计 结果影响很大，缺乏稳定性，不能很好的应对复杂场景。
【发明内容】
[0006] 为了克服现有高光谱图像局部目标检测方法虚警率高的不足，本发明提供一种基 于鲁棒背景估计方法的高光谱图像局部目标检测方法。该方法首先利用基于光谱角的聚类 方法，在白化空间中对输入图像进行聚类，形成具有较高均一性的具备背景。在检测过程 中，引入MCD估计方法估计待测像元所属聚类的背景参数，以获得更可靠的估计结果，从而 提高算法的检测性能。在RIT (Rochester Institute of Technology)提供的数据集上目 标测试结果显示，在各测试目标上的平均得分（反映虚警数）为2. 8,较Halper的改进方法 降低了 4. 4。卫星拍摄的AVIRIS数据集上的测试结果表明，在100%的检测率下，虚警率为 0. 11%，而Halper的方法虚警率为0. 29%左右，全局方法的虚警率为0. 82% ;可见本发明 方法明显降低了虚警率。
[0007] 本发明解决其技术问题所采用的技术方案是：一种基于鲁棒背景估计方法的高光 谱图像局部目标检测方法，其特点是包括以下步骤：
[0008] 步骤一、对于输入高光谱图像，使用MCD估计方法估计出全局背景参数并进行白 化，得到白化数据Dw。
[0009] (1)对于一幅具有p个波段、n个像元的高光谱图像Xn = {Xp x2, ? ? ?，xj，用Xn 中下标在集合Ht中的成员计算其协方差矩阵E t和均值ut，Ht为1到n内的整数集合，即 //, C彳1，2，...，/7丨，|Ht| =h, t为当前迭代代数。
[0010] ⑵计算每一像元Xi到当前均值向量ut的马氏距离，
【权利要求】
1. 一种基于鲁棒背景估计方法的高光谱图像局部目标检测方法，其特征在于包括以下 步骤： 步骤一、对于输入高光谱图像，使用MCD估计方法估计出全局背景参数并进行白化，得 到白化数据Dw ; (1) 对于一幅具有P个波段、n个像元的高光谱图像Xn = {xp x2,. . .，xn}，用Xn中 下标在集合Ht中的成员计算其协方差矩阵E t和均值ut，Ht为1到n内的整数集合，即 //, c彳1，2,"丨，|Ht| =h，t为当前迭代代数； (2) 计算每一像元Xi到当前均值向量+的马氏距离，
i = 1，? ? ?，n ; ⑶将得到的各距离排序dta)，得到序列dta)彡dta)彡…彡d ta); (4)由上述序列构造下次迭代所需的集合Ht+1= ，…，jh}; 重复上述步骤⑴到步骤⑷产生序列det( E J彡det( E 2)彡det( E 3)彡…，该 序列收敛于det ( E m) = 0或det ( E m) = det ( E 时，将此时的协方差矩阵和均值向 量作为估计结果；在产生初始集合％时，先从Xn中随机抽出p+1个成员组成集合V，计算V 上的协方差矩阵[v和均值uv ;若det ( E v) = 0,从Xn中每次增选一个成员来扩充V，直至 det ( E v)尹0 ;对V运用一次所述估计背景参数的步骤（1)到（4)，将步骤（4)产生的集合 H 作为 HQ ;h 的取值为[(n+p+1) /2, n]，且 h>p ; 步骤二、对白化数据Dw进行聚类， (1) 将第一个像元标记为第一类； (2) 对于每一像元Xi计算其与现有每一聚类中心Cj, j = 1，2,. . .，1的夹角，1为当前 聚类中心数目； (3) 若存在一聚类中心Ck，使得\和(；的夹角小于阈值e，则将Xi归入Ck所代表的类 另IJ ;否则，将Xi作为新增聚类中心C1+1，同时1 = 1+1 ; (4) 对于步骤（3)中得到的聚类结果，用每类的均值更新其聚类中心；若各聚类中心较 上次迭代有变化，则转步骤（2);否则，转步骤（5); (5) 取消类成员数量小于A ^的聚类，对其中成员重新使用步骤（2)和步骤（3)进行 标记过程；其中，入为给定阈值，n为图像中像元总数； 获得m个聚类LBs = {Clp Cl2,. . .，Clm}和少量无标签的像元； 步骤三、对于每一聚类Clk，k = 1，2,. . . .，m，使用MCD估计方法，估计其背景参数：协方 差矩阵Ek和均值向量uk; 步骤四、对于每一待测像元，使用其所属类别的局部背景参数作为ACE目标检测算法 的输入，按如下公式得到各像元对应的输出：
其中，s为目标特征光谱，E k和uk分别为第k聚类的协方差矩阵和均值；对于聚类后 仍无标签的像元，使用MCD估计方法估计得到的全局背景参数作为ACE目标检测算法的输 入，由各像元对应的算法输出值组成基于MCD估计方法的高光谱图像局部ACE目标检测的 结果。
【文档编号】G06T7/00GK104346812SQ201410538062
【公开日】2015年2月11日 申请日期:2014年10月13日 优先权日:2014年10月13日
【发明者】张艳宁, 魏巍, 严杭琦, 张磊, 李飞, 王波波 申请人:西北工业大学