基于多特征融合的红外弱小目标检测方法

1.本发明涉及红外图像目标检测领域，具体为基于多特征融合的红外弱小目标检测方法。


背景技术：

2.红外目标检测是远程预警、精确制导等系统的核心技术之一，其中红外弱小目标检测是经典的难题。主要难点在于：1)成像距离很远，目标在图像中不仅缺乏特征信息而且极易与噪点混淆；2)高强度的杂波边缘很容易被误检为目标，使得虚警增加；3)既弱又小的目标往往会淹没在复杂背景中，从而造成严重的漏检。因此，复杂背景下的红外弱小目标检测一直是一项极具挑战性的任务。
3.目前，关于红外弱小目标检测的研究还很少。与之相关的红外小目标检测方法大体上可分为以深度学习为代表的数据驱动法和基于数学、物理知识建模的模型驱动法两类。其中，由于缺乏红外弱小目标图像集、训练开销大且实时性差等问题，导致数据驱动法在实际工程中应用较少。模型驱动的红外小目标检测方法可进一步分为单帧方法和多帧方法。不过，由于无法获取目标的精确运动模型，使得多帧方法在红外弱小目标检测中效果受限。因此，高性能的单帧方法在红外小目标检测中更受重视。
4.自处理思路方面分析，红外小目标的单帧检测方法包括三类：一类可以称之为恢复背景法，其思想是先基于稀疏和低秩矩阵恢复背景，再用原图像减去背景图像，从而达到凸显目标的目的，但是存在对像素值敏感、耗时较长等问题。第二类是背景抑制法，主要是采用滤波技术抑制均匀背景和杂波。但是，该类方法对背景边缘很敏感，制约了其在空天目标检测等复杂场景中的使用。第三类是新兴的抑制背景增强目标并举法，主要是利用人类视觉的局部对比度敏感特性，通过提取目标与背景的局部对比度增强目标。该类多数方法在背景变化平稳的场景中检测效果都较好，但用于复杂场景的红外弱小目标检测时，无法有效抑制复杂背景中的高强度杂波边缘，致使虚警率高。因此，复杂背景下的红外弱小目标检测方法还有较大的发展空间。


技术实现要素：

5.本发明针对复杂背景下红外弱小目标普遍存在检测率低、虚警率高的问题，提出一种基于多特征融合的红外弱小目标检测方法，用于复杂背景下红外图像弱小目标的检测，有效提高检测的准确率及鲁棒性。
6.本发明采用如下技术方案：基于多特征融合的红外弱小目标检测方法，包括以下步骤：首先，基于红外弱小目标局部灰度值大的特性，利用目标与邻域背景的灰度对比度对真实目标进行增强、部分复杂背景进行抑制；其次，利用红外弱小目标自身灰度信息符合二维高斯分布的特性，通过计算红外弱小目标边缘像素点与中心像素点之间的距离和灰度差的协方差来检测目标，得到显着图一；再次，利用红外弱小目标与邻域相似度低的特性，通过计算相似性因子来检测目标，得到显着图二；最后，将显着图一与显着图二进行点乘，通
过融合红外弱小目标多种特性得到最终显着图，对最终显着图计算阈值进行分割得到最终检测结果。
7.上述的基于多特征融合的红外弱小目标检测方法，利用目标与邻域背景的灰度对比度对真实目标进行增强、部分复杂背景进行抑制具体包括以下步骤：将待检测红外图像i
in
转换为灰度图像i；以灰度图像i中的一个像素点s为中心构建滑动窗口，中心单元为s0，其周围的单元s
i
是局部背景区域；取局部背景区域中灰度均值最大的若干个单元，计算这几个单元的灰度均值μ
max
；比较中心单元的灰度均值和μ
max
的大小，若对中心单元像素进行增强，灰度值设置为i(i,j)*(μ0/μ
max
)2，i(i,j)表示中心单元中(i,j)位置的像素灰度值，否则进行抑制，灰度值设置为0；移动滑动窗口，直至灰度图像i所有像素都经过上述计算，得到背景抑制图像t。
8.上述的基于多特征融合的红外弱小目标检测方法，通过计算红外弱小目标边缘像素点与中心像素点之间的距离和灰度差的协方差来检测目标，得到显着图一具体包括以下步骤：以背景抑制图像t中的一个像素点t为中心构建滑动窗口，如果该滑动窗口内像素值不为0，计算该滑动窗口内每个像素到中心像素的距离与灰度差，从而计算整个滑动窗口的协方差系数；移动滑动窗口，直至背景抑制图像t所有像素都经过上述计算，得到由每个滑动窗口的协方差系数组成的协方差显着图cov。
9.上述的基于多特征融合的红外弱小目标检测方法，利用红外弱小目标与邻域相似度低的特性，通过计算相似性因子来检测目标，得到显着图二具体包括以下步骤：以背景抑制图像t中的一个像素点p为中心构建滑动窗口，中心单元称为sm0，其周围的s个背景单元sm
ij
是局部背景区域；分别计算背景单元sm
ij
与中心单元sm0之间的相似性因子sf
i
；移动滑动窗口，直至背景抑制图像t所有像素都经过上述计算，得到s个相似性矩阵，其中第i个相似性矩阵由中心单元与第i个背景单元的相似性因子组成；s个相似性矩阵中相同位置的元素最小值 min(sf
i
)作为相似性矩阵sf对应位置的元素的值，相似性矩阵sf计算方式如下：其中，min(sf
i
)为s个相似性矩阵中相同位置的元素最小值，max(sf
i
)为s个相似性矩阵中相同位置的元素最大值；相似性矩阵sf即为相似性显着图sf。
10.上述的基于多特征融合的红外弱小目标检测方法，将显着图一与显着图二进行点乘，通过融合红外弱小目标多种特性得到最终显着图，对最终显着图计算阈值进行分割得到最终检测结果，具体包括以下步骤：对协方差显着图和相似性显着图进行点乘操作，得到最终检测显着图sm：sm＝cov
⊙
sf，计算出分割阈值th，对最终检测显着图sm进行阈值分割，得到最终的检测结果图i
out
。
11.总体而言，本发明所提出的技术方案充分利用红外弱小目标的局部先验信息，具有以下技术特征及有益效果：
12.(1)由于红外小目标的形状是各向同性的，而背景边缘则是各向异性的。它们在特征上具有差异性，因此需要对窗口细化分块即形成单元。
13.(2)目标与邻域背景存在灰度差异，因此需要利用目标与背景之间的灰度差异增强目标，抑制背景。
14.(3)为分离目标，不仅需要考虑目标与背景的灰度对比，还需利用滑动窗口中心像素与边缘像素的灰度差与距离差的协方差值、相似性因子等先验知识进行分割目标，抑制背景。
15.(4)显著性图仍可能存在一定的杂波或噪声，但目标已明显增强，因此利用恒虚警率法分割，去除剩余虚警。
16.本发明一种基于多特征融合的复杂背景红外弱小目标检测新方法，有效抑制了红外图像中的复杂杂波，同时提升了弱小目标检测的精度。
附图说明
17.图1为本发明的总体框架图；
18.图2为本发明的基本流程图；
19.图3为本发明计算灰度对比度的滑动窗口分布图；
20.图4为本发明计算相似性因子的滑动窗口分布图；
21.图5(a)为本发明实施例中红外原始图像；
22.图5(b)为本发明实施例采用本发明提出的算法后得到的最终检测显着图；
23.图5(c)为本发明实施例的检测结果图。
具体实施方式
24.下面将结合附图及具体实施方式对本发明进行进一步阐述：
25.参照图1和图2，本实施例中的一种复杂背景下红外弱小目标检测方法包括以下步骤：
26.步骤1：输入一幅大小为h
×
w的待检测红外图像i
in
；
27.步骤2：将图像i
in
转换为灰度图像i；
28.步骤3.1：以灰度图像i中的一个像素点s为中心构建一个19
×
19的滑动窗口，分为9个单元，中心单元s0大小为5
×
5，背景区域单元s1、s3、s6、s8大小均为7
×
7，s2、s7大小为5
×
7，s4、s5大小为7
×
5，如图3所示；
29.步骤3.2：计算所有背景区域单元的池化亮度，即单元内所有像素灰度的均值，具体求法为：
[0030][0031]
其中，h
ij
表示局部背景中第i个背景区域单元的第j个像素的灰度值，则μ
i
表示第i个背景区域单元的灰度值均值，n表示区域块内像素点的个数；
[0032]
步骤3.3：取背景区域8个单元中灰度均值最大的4个单元，计算4个单元的灰度均值μ
max
，具体求法为：
[0033][0034]
其中，max4(
·
)表示取
·
中最大的四个灰度值，μ
max
代表其均值；
[0035]
步骤3.4：比较中心单元的灰度均值和μ
max
的大小，若对中心单元像
素进行增强，灰度值设置为i(i,j)*(μ0/μ
max
)2(i(i,j)表示中心单元中(i,j)位置的像素灰度值)，否则进行抑制，灰度值设置为0，具体求法为：
[0036][0037]
步骤3.5：从左至右，从上之下移动滑动窗口，步长为5，直至灰度图像i所有像素都经过上述步骤3.1
‑
步骤3.4计算，得到背景抑制图像t；
[0038]
步骤4.1：以背景抑制图像t中的一个像素点t为中心构建一个5
×
5的滑动窗口，如果该滑动窗口内像素值不为0，计算该滑动窗口内每个像素到中心像素的距离与灰度差，从而计算整个滑动窗口的协方差系数，该系数作为滑动窗口中心像素的值，滑动窗口内每个像素到中心像素的距离均值与灰度差均值具体求法为：
[0039][0040][0041]
其中，i、j分别表示滑动窗口内像素点的x、y方向坐标，为其均值，本发明取3，i
ij
为坐标为(i,j)点的像素灰度值，为中心像素点的像素灰度值，本发明取i
3,3
，n为滑动窗口内像素的个数，本发明取25；
[0042]
步骤4.2：滑动窗口的协方差系数具体求法为：
[0043][0044][0045][0046]
其中，cov(dis,gv)表示滑动窗口的协方差系数，gv表示滑动窗口中心点与边缘像素点的灰度值差，dis表示滑动窗口中心点与边缘像素点的距离差，sl 为滑动窗口的边长，本发明取5；
[0047]
步骤4.3：从左至右，从上之下移动滑动窗口，步长为1，直至背景抑制图像t所有像素都经过上述步骤4.1
‑
4.2计算，得到由每个滑动窗口的协方差系数组成的协方差显着图cov；
[0048]
步骤5.1：以背景抑制图像t中的一个像素点p为中心构建一个25
×
25的滑动窗口，
该窗口均为9单元，中心单元称为sm0，其周围的单元sm
ij
是局部背景区域，如图4；
[0049]
步骤5.2：分别计算8个背景单元与中心单元之间的相似性因子sf
i
，具体求法为：
[0050][0051]
式中，t代表背景抑制图像，(s,t)表示像素点坐标,t(s,t)表示该点的像素灰度值；
[0052]
步骤5.3：从左至右，从上之下移动滑动窗口，步长为1，直至图像t所有像素都经过上述步骤5.2计算，得到8个相似性矩阵，其中第i(i＝1,2,
…
,8)个相似性矩阵由中心单元与第i个背景单元的相似性因子组成；
[0053]
步骤5.4：取8个相似性矩阵中相同位置的元素最小值min(sf
i
)作为相似性矩阵sf对应位置的元素的值，相似性显着图sf具体求法为：
[0054][0055]
步骤6：对协方差显着图cov和相似性显着图sf进行点乘操作，得到最终检测显着图：
[0056]
sm＝cov
⊙
sf
[0057]
步骤7：利用公式计算出分割阈值th，对显着图sm进行阈值分割，得到最终的检测结果图i
out
:
[0058]
th＝μ+λ
×
σ
[0059]
其中，μ和σ为最终显着图sm的均值和标准差，λ为固定参数，本发明取 4。