一种基于时空域背景抑制的红外弱小目标检测方法技术

本发明专利技术属于红外图像处理领域，主要涉及一种基于时空域背景抑制的红外弱小目标检测方法。本发明专利技术用以解决复杂背景下的红外运动弱小目标检测问题。首先使用指导滤波抑制空域中较为平稳的背景杂波；其次，利用红外图像序列中目标的运动信息，在时域上采用梯度权重滤波方法抑制时域中变化较为缓慢的背景；然后将时空域背景抑制结果相融合，得到背景抑制后的弱小目标图像；最后利用自适应阈值分割图像，检测出弱小目标。该发明专利技术在目标检测时不仅利用了红外弱小目标的空间灰度信息，还充分利用了目标的时域运动信息，分别从时域和空域抑制背景杂波，因此大大提高了复杂背景下运动弱小目标的检测性能。

【技术实现步骤摘要】
一种基于时空域背景抑制的红外弱小目标检测方法
本专利技术属于红外图像处理领域，主要涉及一种基于时空域背景抑制的红外弱小目标检测方法。
技术介绍
在红外搜索跟踪系统(IRST)中，来袭目标距离探测器较远，在红外图像中通常显示为复杂背景中的几个像素点，同时，由于大气衰减及干扰，红外图像中目标与背景的对比度和信噪比较低，这就为后续的目标检测带来了困难。如何从低对比度、低信噪比的红外图像中准确、稳定的检测出目标，就成为IRST中的一项关键技术。近年来，由于红外探测技术在军事上的重大意义，许多研究者对红外弱小检测进行了深入的研究，提出了许多检测方法。主要有时域、空域、频域、小波变换、偏微分方程等滤波方法。这些方法分别从不同的角度解决红外弱小目标的检测问题，收到了一定的效果，但针对低信噪比的复杂背景红外图像，这些算法就显示出背景抑制效果差，检测虚警率高，算法复杂度高等缺陷。针对弱小目标在空域上只有灰度信息，检测难度大的问题，许多学者将空域检测方法与目标时域运动信息相结合，提出了时空域融合的目标检测算法，收到了一定的效果。在《基于时空域融合的红外弱小目标检测算法》(见弹箭与制导学报。2011年31卷(2期)：P225-227，作者：胡涛涛，樊翔，马东辉)中所述的一种时空域融合的红外弱小目标检测方法，首先用tophat变换对图像在空域上进行背景抑制，然后对背景抑制后的序列图像使用三帧差分法进行运动目标的检测，并对检测的结果进行或运算来累计能量，再利用形态学闭运算连接断裂轨迹，最后通过阈值检测出目标轨迹。该方法对背景变化缓慢、信噪比较高，目标运行速度较快的红外运动弱小目标序列的检测具有一定效果，但该方法同时存在明显的不足：1、tophat滤波方法对于信噪比较低的红外弱小目标图像背景抑制效果较差，并且结果与结构元的选择有非常大的关系，结构元选取不当有可能无法检测出弱小目标。2、三帧差分法可能导致目标强度变弱，且对低速运动的目标检测效果并不好，对背景变化敏感。3、算法仅能检测出目标运行轨迹，无法给出目标当前帧位置，不具有实时性。在《Smalltargetdetectionusingbilateralfilterandtemporalcrossproductininfraredimages》(见InfraredPhysics&Technology。54(2011)：P403-411，作者：Tae-WukBae等)中，作者提出了一种时空域融合的红外弱小目标检测方法。该方法的步骤如下：1、在时域上对n帧图像每一像素点求时域向量积来提取目标运行轨迹。2、根据时域目标轨迹图像的灰度值来生成参数索引图，使不同的灰度值对应不用的σd和σr。3、根据选定的σd和σr，对图像在空域上利用双边滤波进行背景抑制，得到空域背景抑制后的图像。4、将1的结果与3的结果进行点乘。5、选取阈值进行分割，得到目标检测结果。算法存在如下问题：1、时域上利用时域向量积只能检测到目标轨迹，且云层边缘起伏较大时，虚警率会增加。2、空域滤波受时域滤波影响较大，当时域结果存在较多虚警及杂波时，会使空域滤波结果变差。3、时空域结果融合时采用点乘，在目标轨迹上容易形成虚警点。
技术实现思路
针对低信噪比复杂背景下序列红外运动弱小目标检测中虚警率高的问题，本专利技术提出了一种时空域结合的红外序列图像运动目标检测方法。该方法分别在时域和空域上进行背景抑制。在时域上，充分利用目标运动信息，使用梯度权重滤波抑制背景，得到时域背景抑制后的图像；在空域上利用目标灰度信息，对单帧图像使用指导滤波进行背景预测，进而抑制背景，最后将时域滤波与空域滤波的结果相融合，使用自适应阈值分割图像，检测出目标。本方法能够实时地检测出目标位置，大大降低虚警概率，并且简单有效。本专利技术是通过以下技术方案实现的，本专利技术包括以下步骤：A、时域背景抑制：(1)取N帧图像，绘制图像中每一像素点在N帧图像中灰度值的变化的曲线(时域曲线)f(m，n，k)＝xk(m，n)k＝1，2…N，(m，n)为像素点的位置坐标，k为图像的帧数，x为灰度值；(2)对于每个像素点的时域曲线，计算曲线上每点的梯度值g：gk(m，n)＝|[xk(m，n)-xk-1(m，n)]+[xk(m，n)-xk+1(m，n)]|(3)使用高斯核计算时域曲线上每一点的权重W：其中，ε为调节参数；(4)对时域曲线进行梯度权重滤波，滤波器在第k帧的输出Pok为：其中，R为归一化参数，(5)用原时域曲线减去梯度权重滤波的结果，得到第N帧时域背景抑制后的图像：xN′＝xN-PoNB、在空域上对图像进行背景抑制：(1)计算指导滤波中图像每一像素点的权值：μk和σ2为指导图像I在滤波窗中的均值和方差，ωk为滤波窗口，ε为调节参数，调整滤波器的平滑程度，|ω|为窗口ωk中像素的个数；(2)计算指导滤波在每一点(m，n)处的输出：其中，L为滤波窗口的半径；(3)第N帧图像空域背景抑制后的结果如下：IsoutN＝PN-QNC、将A的结果与B的结果做与运算，得到第n帧背景初步抑制后的图像：ItempN＝IsoutN·xND、将A的结果作为原图像，C的结果作为指导图像，进行指导滤波，得到背景抑制结果：其中E、采用自适应阈值分割将图像二值化，得到最终目标检测结果：其中Th为阈值：Th＝μ+10σ2，μ、σ2分别为图像的均值和方差。本专利技术采用指导滤波来实现空域背景预测，由于其在平滑图像的同时具有良好的边缘保持特性，因此可以对图像中的边缘细节进行更准确地预测，进而有效抑制空域中较为平稳的背景及边缘杂波；在时域中，本专利技术利用目标的运动信息，提出了一种梯度权重的滤波方法进行时域背景预测，可以有效抑制时域中变化较为缓慢的背景杂波；时域与空域背景抑制独立进行，可以充分利用红外图像中弱小目标的灰度信息与运动轨迹信息；采用指导滤波的方法融合时域与空域背景抑制的结果，巧妙的将时空域结果相融合，并且能够进一步抑除背景杂波，为后续的目标检测提供更好的基础。附图说明通过参照附图更详细地描述本专利技术的示例性实施例，本专利技术的以上和其它方面及优点将变得更加易于清楚，在附图中：图1为本专利技术的一种基于时空域背景抑制的红外弱小目标检测方法的流程示意图；图2为不同的参数下对同一幅图像利用指导滤波平滑后的结果；(a)原图(b)ε＝0.1，ωk＝2(c)ε＝0.1，ωk＝4(d)ε＝0.4，ωk＝2；图3为红外图像中几个特征像素点的时域曲线；图4为两组包含弱小运动目标的复杂天空背景红外图像检测结果。(a)为原图；(b)为空域背景抑制后的结果；(c)为时域背景抑制后的结果；(d)为时空域融合后的背景抑制结果；(e)为使用自适应阈值分割后的检测结果。具体实施方式在下文中，现在将参照附图更充分地描述本专利技术，在附图中示出了各种实施例。然而，本专利技术可以以许多不同的形式来实施，且不应该解释为局限于在此阐述的实施例。相反，提供这些实施例使得本公开将是彻底和完全的，并将本专利技术的范围充分地传达给本领域技术人员。在下文中，将参照附图更详细地描述本专利技术的示例性实施例。参照图1，本专利技术的具体实现步骤如下：本方法具体实现步骤如下：步骤1、对序列红外图像的前N帧进行时域滤波，抑制背景杂波，图4(a)为序列红外图像中的一帧。1.1输入1～N本文档来自技高网...

【技术保护点】
一种基于时空域背景抑制的红外弱小目标检测方法，其特征在于：该方法包含以下步骤：A、时域背景抑制：(1)取N帧图像，绘制图像中每一像素点在N帧图像中灰度值的变化曲线(时域曲线)f(m，n，k)＝xk(m，n) k＝1，2…N，(m，n)为像素点的位置坐标，k为图像的帧数，x为灰度值；(2)对于每个像素点的时域曲线，计算曲线上每点的梯度值g：gk(m，n)＝|[xk(m，n)‑xk‑1(m，n)]+[xk(m，n)‑xk+1(m，n)]|(3)使用高斯核计算时域曲线上每一点的权重W：wk(m,n)=e-gk(m,n)2/ϵ2]]>其中，ε为调节参数；(4)对时域曲线进行梯度权重滤波，滤波器在第k帧的输出Pok为：Pok(m,n)=1RΣl=k-2k+2wl(m,n)xl(m,n)]]>其中，R为归一化参数，(5)用原时域曲线减去梯度权重滤波的结果，得到第N帧时域背景抑制后的图像；xN′＝xN‑PoNB、空域背景抑制：对输入第N帧图像进行指导滤波，得到背景的估计图像(即背景预测)，用原图像减去经指导滤波后的背景图像，得到空域背景抑制后的图像；C、将A的结果与B的结果做与运算，得到第N帧背景初步抑制后的图像；D、将A的结果作为原图像，C的结果作为指导图像，进行指导滤波，得到背景抑制结果；E、采用自适应阈值法将D的结果二值化，得到最终目标检测结果。...

【技术特征摘要】
1.一种基于时空域背景抑制的红外弱小目标检测方法，其特征在于：该方法包含以下步骤：A、时域背景抑制：(1)取N帧图像，绘制图像中每一像素点在N帧图像中灰度值的时域曲线f(m，n，k)＝xk(m，n)k＝1，2…N，(m，n)为像素点的位置坐标，k为图像的帧数，x为灰度值；(2)对于每个像素点的时域曲线，计算曲线上每点的梯度值g：gk(m，n)＝|[xk(m，n)-xk-1(m，n)]+[xk(m，n)-xk+1(m，n)]|(3)使用高斯核计算时域曲线上每一点的权重W：其中，ε为调节参数；(4)对时域曲线进行梯度权重滤波，滤波器在第k帧的输出Pok为：其中，R为归一化参数，(5)用原时域曲线减去梯度权重滤波的结果，得到第N帧时域背景抑制后的图像；xN′＝xN-PoNB、空域背景抑制：对输入第N帧图像进行指导滤波，得到背景的估计图像，用原图像减去经指导滤波后的背景图像，得到空域背景抑制后的图像；C、将A的结果与B的结果做与运算，得到第...

【专利技术属性】
技术研发人员：秦翰林，李佳，延翔，周慧鑫，牟媛，宗靖国，韩姣姣，曾庆杰，郝静雅，倪曼，刘上乾，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61