基于时空显著性和四元余弦变换的红外弱小目标检测方法
【技术领域】
[0001] 本发明属于红外图像处理领域，涉及红外弱小目标显著性检测。
【背景技术】
[0002] 红外探测系统在红外警戒领域中有着广泛的应用，其中弱小目标检测技术在预警 系统、精确打击武器和防空系统等领域中发挥非常重要的作用。由于红外目标有着信噪比 低，对比度差，目标像素点少等特点，所以红外弱小目标检测是一个具有很大挑战性的关键 课题。
[0003] 在计算视觉领域，视觉显著性已经得到了学者们越来越多的关注，在这个过程中， 也形成了一些视觉显著性计算模型。视觉显著性已经应用在图像压缩、卫星遥感图像、图像 与场景分析、视频处理等方面。通过显著性检测模型的计算，得到的显著性检测图包含了可 能引起人眼关注的区域，抑制了部分背景，极大的简化了通常对整个图像区域处理的过程， 提高了计算速度和计算精度。
[0004] 随着红外检测技术的不断发展，已经提出了许多检测方法，主要有时域、空域、小 波变换等方法，随着视觉显著性提出和发展，显著性检测也用在了红外弱小目标检测中。基 于显著性的红外弱小目标检测方法目前也有一些，其中基于相位谱傅里叶变换(PFT)和相 位谱四元组傅里叶变换(PQFT)这两种算法都有一定的检测效果，PQFT是对PFT的扩展，在目 标检测中取到了更好的效果。这是因为PQFT具有最多支持四个数据通道的能力，可以为傅 里叶变换提供更多的信息，而传统傅里叶变换只能使用灰度信息。基于四元组傅里叶变换 的红外弱小目标检测方法，首先，利用S0DD滤波器得到四张二阶方向导数图，组成一个四元 组;其次，对四元组进行傅里叶变换;然后，只利用傅里叶变换的相位谱进行反傅里叶变换， 得到空域的图像;最后，对空域变换的结果取幅度谱，并用高斯低通滤波器对图像滤波，得 到最终的红外弱小目标检测图。这种方法在目标对比度比较大的时候取得了比较好的效 果，但是对于目标昏暗对比度比较低的红外图像，在检测结果中，目标亮度比较低，检测结 果不理想，而且只能检测单帧图片，没有加入运动信息。
【发明内容】
[0005] 结合现有技术，针对在目标亮度低和对比度低的情况下，得到的检测结果不理想 的问题，本发明提出了一种基于视觉显著性和四元余弦变换(QDCT)的红外弱小目标检测的 方法。
[0006] 本发明方法利用图像的时域特征和空域特征，得到两个特征通道图；结合类高斯 分布的红外弱小目标具有各向同性的特征，而背景具有局部方向的特征，S0DD滤波器作为 方向滤波器对输入视频的每一帧滤波，得到另外两个特征通道图，四个特征通道图构成一 个四元组;最终的红外弱小目标检测结果是先利用四元组余弦变换，然后用符号函数处理， 再进行反四元组余弦变换得到。
[0007] 本发明是通过以下技术方案实现的，基于时空显著性和四元余弦变换的红外弱小 目标检测方法，包括以下步骤：
[0008]步骤1:采集2D视频图像，视频图像为I(x，y，t)，其中t = l · · · N，N为视频总帧 数，（x，y)为像素点的位置；
[0009] 步骤2:计算二阶方向导数图，并构建四元组；
[0010] 步骤2.1:对当前帧图像计算0°，45°，90°，135°这4个方向的二阶方向导数图；
[0011]步骤2.2:将0°，90°的二阶方向导数图进行融合，将45°，135°的二阶方向导数图进 行融合，获得两个特征通道图；
[0012] 步骤2.3:计算当前帧的强度特征通道图b(X，y，t)和运动特征通道图m(X，y，t);
[0013] 步骤2.4:将步骤2.2和步骤2.3共获得的4个特征通道图加权相加，构成四元组；
[0014] 步骤3:对步骤2获得四元组进行左四元余弦变换或右四元余弦变换，获得新的四 元组；
[0015]步骤4:计算检测结果图；
[0016] 步骤4.1:定义四元组的符号函数sgn(q)为：
[0018]其中X 〇，X i，X 2，X 3是四元组的四个分量，| q |是四元组的模，
[0019] 步骤4.2:对四元余弦变换的符号函数做四元反余弦变换，得到空域的图像，然后 用高斯低通滤波器对空域图像的幅度谱滤波，获得检测结果图。
[0020] 进一步的，所述步骤2的四元组的具体构建方法为：
[0021] (1)二阶方向导数图的计算公式；
[0022] S0DD滤波器是通过facet model设计得到的，通过facet model可以得到三个卷积 核：
[0024]二阶方向导数图的表达式为：
[0026 ] K4，K5，K6是W4，W5，W6分别对待处理帧I (X，y，t)卷积的结果，1为所需要的方向，α为1 与X轴的夹角，β为1与y轴的夹角；
[0027] (2)计算二阶方向导数图；
[0028]本发明中1分别为0°，45°，90° ,135°，得到四幅二阶方向导数图dQ，d45，d9Q，d135:
[0029] do(x，y，t) = 2K4(x，y，t)
[0030] d45(x，y，t) = 2K4(x，y，t)+2K5(x，y，t)+2K6(x，y，t)
[0031] d9〇(x,y,t) = 2K6(x,y,t) (3)
[0032] di35(x，y，t)=2K4(x，y，t)-2K5(x，y，t)+2K6(x，y，t)
[0033]
[0034] (3)融合四个二阶方向导数图，得到两个特征通道图
[0035] di(x,y,t) = do(x,y,t) · dg〇(x,y,t)
[0036] d2(x，y，t) = d45(x，y，t)di35(x，y，t) (4)
[0037]
[0038] 这里的"·"是两个相互垂直的二阶方向导数图对应像素点的乘积；
[0039] (4)强度特征通道b(x，y，t)和运动特征通道m(x，y，t)的计算公式：
[0040] b(x,y,t)= | I(x,y,t) (5)
[0041] m(x,y,t)= |b(x,y,t)-b(x,y,t-T) (6)
[0042] 其中τ = 3,运动特征通道是用帧差法得到的。
[0043] (5)待处理帧的四个特征通道共同构成四元组的，四元组表示为：
[0044] q(x，y，t)=m(x，y，t)+di(x，y，t)i+d2(x，y，t)j+b(x，y，t)k (7)
[0045] 其中i，j，k满足jk = i2 = j2 = k2 = _l，i 丄 j 丄k。
[0046] 进一步的，所述步骤3中四元余弦变换的方法为：
[0047] 四元余弦变换有左四元余弦变换和右四元余弦变换，这两种变换的结果是两个不 同的四元组，这是因为i j = k而j i = -k，四元组q (X，y，t)的分辨率为Μ X N，两种四元余弦变 换的计算公式为：
[0050]
[0051] ，yQ是单位纯虚四元组，和cf都是变换的系数，表 达式相似，而/C和/C都是余弦函数，表达式也相似，分别表示为：
[0054]
[0057]
[0058] 进一步的，所述步骤4中最终检测结果图的计算方法为：
[0059] (1)四元组q的符号函数sgn(q)定义为：
[0061 ]其中叉〇，叉1，叉2，13是四元组的四个分量，|9|是四元组的模，
[0062] (2)左反四元余弦变换和右反四元余弦变换的计算公式为：
[0065]
[0066] (3)最终的红外弱小目标检测结果是通过对四元余弦变换的符号函数做四元反余 弦变换，得到空域的图像，然后用高斯低通滤波器对空域图像的幅度谱滤波，这里用左余弦 变换，右余弦变换是类似的，计算方法如下：
[0067] q,(x,y,t) = IQDCTL(sgn(QDCTL(q(x,y,t))) (13)
[0068] DM(x,y,t)=g*|q，(x,y,t) |2 (14)
[0069] 其中，δ是高斯低通滤波器的方差，取δ = 1.5。
[0070] 本发明采用S0DD滤波器把具有局部方向性的杂波滤波为类条状分布的背景杂波， 得到四幅二阶方向导数图，通过两组相互垂直的通道分别相乘，得到两个特征通道图，加入 强度特征通道图和运动特征通道图，共同组成一个四元组，用四元余弦变换的符号函数的 反四元余弦变换得到最终的红外弱小目标检测结果，对四元余弦变换取符号函数的目的是 抑制背景，检测具有类高斯分布的显著的红外弱小目标。
【附图说明】
[0071] 图1为构建四元组的基本流程；
[0072] 图2为红外弱小目标检测的基本流程；
[0073] 图3为红外弱小目标检测的结果（只取了每个视频的一帧，上面是原图，下面是对 应的检测结果）。
【具体实施方式】
[0074] 1.计算二阶方向导数图
[0075] 二阶方向导数图是通过二阶方向导数滤波器在四个不同的方向滤波