一种基于红外图像的对空目标检测方法与流程

本发明涉及计算机视觉领域，特别是指一种基于红外图像的对空目标检测方法。
背景技术：
：红外目标的检测方法大致可分为两类：先跟踪后检测(tbd，trackbeforedetect)、先检测后跟踪(dbt，detectbeforetrack)。tbd方法一般需要一些先验信息，使用相邻的多张图像，根据计算出的可能的路径的后验概率估计出目标位置。这些方法对于机动性强、移动速度快的目标不太适用。而且tbd方法一般计算量较大，实时性较差。与tbd方法相比，dbt方法实现简单、计算速度快、计算资源需求少、需要的先验知识少，一直是红外目标检测方法研究中备受关注的方向。现有的dbt方法有基于传统图像处理方法实现的方法，有基于边缘检测方法实现的方法，这些方法对噪声较为敏感，易受到杂波干扰。而基于学习的方法依赖于训练样本以及特征的设计，由于对空红外目标样本较少，难以覆盖全面，因此基于学习的方法不易于实现。近年来有不少基于人类视觉机理的方法被提出来，但这些方法也或多或少存在由于噪声的干扰而导致虚警率较高的问题。技术实现要素：本发明要解决的技术问题是提供一种基于红外图像的对空目标检测方法，以解决现有技术所存在的由于噪声干扰和目标亮度不均匀导致的虚警率高的问题，并提高检出率。为解决上述技术问题，本发明实施例提供一种新的基于红外图像的对空目标检测方法，包括：提取红外图像各像素点邻域的对比度特征，生成对比度图像；提取红外图像各像素点的邻域亮度差分特征，生成邻域亮度差分图像；融合对比度图像和邻域亮度差分图像，生成融合后的图像；对融合后的图像进行置零操作的双阈值处理，提取目标位置。进一步地，所述提取红外图像各像素点的邻域亮度差分特征，生成邻域亮度差分图像包括：定义红外图像中点(x,y)处的邻域亮度差分指标d(x,y)为：其中，dk(x,y)为点(x,y)处第k个尺度的亮度差分；所述点(x,y)处第k个尺度的亮度差分dk(x,y)定义为：其中，ak(x,y)为点(x,y)处第k个尺度的邻域内的像素值的平均值；所述点(x,y)处第k个尺度的邻域内的像素值的平均值ak(x,y)定义为：其中，f(x,y)为图像中点(x,y)处的像素值，ωk为点(x,y)处第k个尺度的邻域，nk为ωk内的像素点的个数；所述点(x,y)处第k个尺度的邻域ωk定义为：ωk＝{(u,v)|max{|u-x|,|v-y|}≤k},k＝1,2,...,l其中，l为所考虑的最大尺度；所述生成的邻域亮度差分图像中点(x,y)处的值即为红外图像中点(x,y)处的邻域亮度差分d(x,y)。进一步地，所述融合对比度图像和邻域亮度差分图像包括：定义图像中点(x,y)处的最终的检测指标r(x,y)为：r(x,y)＝c(x,y)×d(x,y)其中，c(x,y)为点(x,y)处邻域的对比度，d(x,y)为点(x,y)处的邻域亮度差分；点(x,y)处的最终检测指标r(x,y)即为点(x,y)处邻域的对比度和点(x,y)处的邻域亮度差分的乘积；所述融合对比度图像和邻域亮度差分图像生成的融合后的图像中，点(x,y)处的值即为点(x,y)处的最终的检测指标r(x,y)。进一步地，所述对融合后的图像进行置零操作的双阈值处理，提取目标位置包括：利用t1、t2对融合后的图像做置零的阈值处理，得到图像i1、i2；其中图像i2中，点(x,y)处的值i2(x,y)为：图像i1中，点(x,y)处的值i1(x,y)为：其中，r(x,y)为融合后的图像中，点(x,y)处的值，也即红外图像中，点(x,y)处的最终的检测指标；进一步地，对于图像i2中的每个连通域，在i1中寻找与之有交集的连通域，得到图像i0；其中图像i0在点(x,y)处的值i0(x,y)为：进一步地，对于每个i0中的连通域λ，提取目标位置[x0y0]t：其中，s为连通域λ中所有像素的像素值之和，即：通过上述方法便可提取出红外图像中的对空目标位置。本发明的上述技术方案的有益效果如下：上述方案中，提取红外图像各像素点邻域的对比度特征，生成对比度图像；提取红外图像各像素点的邻域亮度差分特征，生成邻域亮度差分图像；融合对比度图像和邻域亮度差分图像，生成融合后的图像；对融合后的图像进行置零操作的双阈值处理，提取目标位置。这样，可以充分利用局部对比度与邻域亮度差分两种指标的特点，并通过置零操作的双阈值法，削弱噪声的干扰和目标亮度不均匀的影响，降低虚警率，同时提高检出率。附图说明图1为本发明实施例提供的一种基于红外图像的对空目标检测方法的流程示意图；图2为本发明实施例提供的一种基于红外图像的对空目标检测方法的详细流程示意图；图3为本发明实施例提供的在目标和背景像素的邻域亮度差分图像；图4为本发明实施例提供的用于测试实施例性能的测试样例图像；图5为本发明实施例提供的本发明与不同现有方法的roc曲线。具体实施方式为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。本发明针对现有技术所存在的由于噪声干扰导致的虚警率高的问题，提供一种新的基于红外图像的对空目标检测方法。如图1所示，本发明实施例提供的基于红外图像的对空目标检测方法，包括：s101，提取红外图像各像素点邻域的对比度特征，生成对比度图像；s102，提取红外图像各像素点的邻域亮度差分特征，生成邻域亮度差分图像；s103，融合对比度图像和邻域亮度差分图像，生成融合后的图像；s104，对融合后的图像进行置零操作的双阈值处理，提取目标位置。本发明实施例所述的基于红外图像的对空目标检测方法，提取红外图像各像素点邻域的对比度特征，生成对比度图像；提取红外图像各像素点的邻域亮度差分特征，生成邻域亮度差分图像；融合对比度图像和邻域亮度差分图像，生成融合后的图像；对融合后的图像进行置零操作的双阈值处理，提取目标位置。这样，可以充分利用局部对比度与邻域亮度差分两种指标的特点，并通过置零操作的双阈值法，削弱噪声的干扰和目标亮度不均匀的影响，降低虚警率，同时提高检出率。为了更好地理解本发明实施例所述的基于红外图像的对空目标检测方法，对所述方法进行详细说明，如图2所示，所述方法具体可以包括：a11，提取红外图像各像素点邻域的对比度特征，生成对比度图像考虑到对空目标距离红外传感器一般较远，在红外图像上所占据的尺寸较小，且外观通常表现为一个高亮的光斑。而背景则表现为相对亮度较暗的大面积连续区域。根据这样的特点，并考虑目标在红外图像中可能存在一定尺寸的变化，本实施例采取如下方式提取图像各点处的局部对比度。首先，定义红外图像中，点(x,y)处第k个尺度的邻域ωk为：ωk＝{(u,v)|max{|u-x|,|v-y|}≤k},k＝1,2,...,l(1)其中，l为所考虑的最大尺度。对红外图像中的每个点(x,y)以及所考虑的所有尺度k＝1,2,...,l，在点(x,y)处的第k个尺度的邻域ωk内，计算点(x,y)处第k个尺度的邻域内的像素值的平均值ak(x,y)：其中，f(x,y)为图像中点(x,y)处的像素值，nk为ωk内的像素点的个数。计算点(x,y)处第k个尺度的局部对比度ck(x,y)：ck(x,y)＝(ak(x,y)-al(x,y))2(3)即点(x,y)处第k个尺度的邻域内的像素值的平均值与所考虑的最大尺度(第l个尺度)的邻域内的像素值的平均值之差的平方。计算点(x,y)处的局部对比度c(x,y)：即从不同尺度的局部对比度中取最大值作为点(x,y)处的局部对比度指标。最后，所述提取红外图像各像素点邻域的对比度特征，生成的对比度图像中，点(x,y)处的值即为红外图像点(x,y)处的局部对比度c(x,y)。通过这种考虑多尺度的方式，计算不同尺度邻域与考虑的最大尺度邻域内平均灰度之差的平方，利用了目标与背景间的高对比度的特点，并考虑到了目标在红外图像中的尺度变化，可突出红外图像中不同尺度的目标。a12，提取红外图像各像素点的邻域亮度差分特征，生成邻域亮度差分图像红外图像中的对空目标与背景之间的高对比度还体现在目标与背景之间有明显的不连续性。将a11中所述的点(x,y)处第k个尺度的邻域内的像素值的平均值ak(x,y)与第k个尺度的邻域内的像素值的平均值ak+1(x,y)相减可得到类似于图3所示的相邻尺度邻域内平均灰度的差分，横轴为公式(1)中的尺度参数k，纵轴为公式(5)中的dk(x,y)。图3中(a)图为红外图像中目标所在像素点的邻域亮度差分，(b)图为红外图像中背景所在像素点的邻域亮度差分。从图3中可以看出，目标所在像素点的邻域亮度差分有一个明显的凸出，而背景所在像素点的邻域亮度差分则较为平滑。这也是本发明中利用邻域亮度差分作为检测红外图像中对空目标的依据。首先，对红外图像中的每个点(x,y)和所有考虑的尺度k＝1,2,...,l，计算点(x,y)处第k个尺度的邻域内的像素值的平均值。其定义与a11中的相同，这里不再赘述。计算点(x,y)处第k个尺度的邻域内的亮度差分dk(x,y)：dk(x,y)＝ak(x,y)-ak+1(x,y)(5)即相邻尺度邻域内像素值的平均值的差分。计算点(x,y)处的邻域亮度差分指标d(x,y)：即从不同尺度的邻域亮度差分中取最大值作为点(x,y)处的邻域亮度差分指标。最后，所述提取红外图像各像素点的邻域亮度差分特征，生成的邻域亮度差分图像中，点(x,y)处的值即为红外图像中点(x,y)处的邻域亮度差分d(x,y)。这样，便可以找到点(x,y)不同尺度的邻域亮度差分中的最大值。当点(x,y)位目标所在像素点时，该指标将得到较大的值，且目标与背景之间的不连续性，也就是亮度的突变越明显，该指标得到的值也就越大。而相反，当点(x,y)为背景像素时，该指标得到的值则较小。a13，融合对比度图像和邻域亮度差分图像，生成融合后的图像a11和a12中生成的对比度图像和邻域亮度差分图像是的目标更为突出，但均存在不同程度的杂波干扰，这也是导致最终虚警率较高的主要原因之一。因此本发明中融合以上两种特征，进一步对目标区域进行增强，并消除杂波的干扰。融合a11和a12中生成的对比度图像和邻域亮度差分图像，即将对比度图像和邻域亮度差分图像进行元素间乘法，得到最终的融合后的图像。融合后的图像中点(x,y)处的值通过以下公式计算：r(x,y)＝c(x,y)×d(x,y)(7)其中，c(x,y)为点(x,y)处邻域的对比度，d(x,y)为点(x,y)处的邻域亮度差分。点(x,y)处的最终检测指标r(x,y)即为点(x,y)处邻域的对比度和点(x,y)处的邻域亮度差分的乘积。最后，所述融合对比度图像和邻域亮度差分图像，生成的融合后的图像中，点(x,y)处的值即为点(x,y)处的最终的检测指标r(x,y)。这样通过融合两种特征，进一步突出红外图像中目标所在的像素，避免由单一特征中存在的不足导致的杂波的残留，并使目标更为明显。a14，对融合后的图像进行置零操作的双阈值处理，提取目标位置在融合了两个特征之后，在融合后的图像中，绝大部分背景杂波被滤除，而目标区域也变得更为突出。但融合后的图像中，仍难免存在低亮度的杂波。同时，对于不均匀的目标，单一阈值也容易将亮度不均匀的目标错误地提取为两个。因此本发明中采用双阈值的方法，即通过一大一小两个阈值，对目标的位置进行提取。本实施例中阈值t1、t2通过以下公式计算得到：t1＝μ+k1×σ(8)t2＝μ+k2×σ(9)其中，μ为所述融合后的图像的均值，即σ为所述融合后的图像的标准差，即n为所述融合后的图像总的像素数目；k1、k2为自选参数；r(x,y)为融合后的图像中，点(x,y)处的值，也即红外图像中，点(x,y)处的最终的检测指标。在计算得到一大一小两个阈值之后，利用两个阈值分别对融合后的图像应用置零操作的阈值处理，得到图像i1、i2。其中图像i2中，点(x,y)处的值i2(x,y)为：图像i1中，点(x,y)处的值i1(x,y)为：进一步地，对于应用大阈值获得的图像中的每个连通域，在应用小阈值获得的图像中寻找与之有交集的连通域，得到图像i0。在本实施例中，t1＜t2，因此对于图像i2中的每个连通域，在i1中寻找与之有交集的连通域，得到图像i0。其中图像i0在点(x,y)处的值i0(x,y)为：进一步地，在得到图像i0之后，对于每个i0中的连通域λ，提取目标位置[x0y0]t：其中，s为连通域λ中所有像素的像素值之和，即：为了验证本发明的效果，本实施例使用如下评价指标：(1)检出率(pd，probabilityofdetection)：检测出的正确的目标的个数与数据集中总的目标的个数的比值，衡量检测方法的检测能力，数值越高越好。(2)虚警率(fa，falsealarmratio)：存在虚警的帧数与数据集中总的帧数的比值，数值越小越好。(3)受试者工作特征(roc，receiveroperationcharacteristics)曲线：检出率与虚警率的曲线，可以直观地看出检出率随虚警率的变换趋势。本实施例所述的基于红外图像的对空目标检测方法具有以下优点：(1)针对应用单一指标对红外图像中的对空目标进行检测的方法容易残留过多的杂波导致虚警率较高的问题，本发明提出了一种新的检测指标，并融合了复数个检测指标对红外图像中的对空目标进行检测。该方法能够克服单一指标导致的杂波残留过多的问题，有效降低虚警率。(2)本发明中的检测指标采用了多尺度的计算方法，对红外图像中对空目标尺寸的变化具有一定的适应性。(3)针对融合后的图像可能残留的较暗的杂波，和红外图像中目标亮度可能不均匀导致应用单一阈值提取时可能出现错误提取的问题，本发明在提取目标的过程中，应用了双阈值的提取方法，即使用较大的阈值确定目标的大致位置，消去融合后的图像中可能残留的杂波，应用较小的阈值确定目标的形状，防止目标由于亮度不均匀而被错误地提取成两个，有效降低检测方法的虚警率。本实施例中，应用了两组共70张对空目标的红外图像序列来讨论所述的基于红外图像的对空目标检测方法的性能。如图4所示，图4中(a)、(b)分别为两个序列中的一幅图像。(a)所在的序列中，图像云层较多，目标较小，且明显与厚重的云层混杂在了一起。(b)所在的序列中，图像云层较为稀薄，但目标的尺寸较大，且存在一定目标亮度不均的现象。在实施过程中，首先对输入图像提取红外图像各像素点邻域的对比度特征，生成对比度图像；再对输入图像提取红外图像各像素点的邻域亮度差分特征，生成邻域亮度差分图像。然后，将生成的对比度图像和邻域亮度差分图像融合，生成融合后的图像。最后，对融合后的图像进行置零操作的双阈值处理，提取目标位置。本发明应用的基于红外图像的对空目标检测方法的流程图如图1、图2所示。在本实施例中，对于所述提取红外图像各像素点邻域的对比度特征和所述提取红外图像各像素点的邻域亮度差分特征时所考虑的最大的尺度l设置为l＝4，目标位置提取时的参数k1设置为k1＝1。为了验证本发明所述方法的性能，本文通过实验将其与其他方法进行了比较，测得了不同方法的roc曲线，实验结果如图5所示。并几个方法中选择性能优异的方法测试在不同阈值的条件下，检测方法对输入图像的检出率与虚警率，实验结果如表1至表3所示。图5中，基于加权局部差分度量(wldm，weightedlocaldifferencemeasure)的方法、基于局部突变加权信息熵(lmwie，localmutationweightedinformationentropy)的方法、基于局部对比度量(lcm，localcontrastmeasure)的方法和基于平均灰度绝对差分最大值(agadm，averagegrayabsolutedifferencemaximum)的方法分别为现有的不同的红外目标检测方法。而在对比检出率与虚警率时，实验选择了性能优异的wldm方法。实验中对不同的方法通过实验选择了检测效果最优异的参数来与本发明所述的方法进行对比。从图5中可以看出，本发明所述的方法可以在较高的检出率的条件下，保持较低的虚警率，具有较好的检测性能。表1中显示了本实施例中应用不同方法在图4(a)所在序列测得的检出率与虚警率的数据。从中可以看出，本发明在随着阈值的变化，检出率一直保持100％，而虚警率在阈值达到一定是降为零；而其他方法在随着阈值的变化过程中，虚警率虽有所下降，但检出率亦有所降低。从表2及表3中可以看出，本文的方法在随着阈值的提高，虚警率在不断下降；而wldm方法随着阈值的不断增大，虚警率呈现出先下降后上升的走向。wldm虚警率的升高，是由于其采用的单阈值的目标提取策略，导致在阈值增大到一定程度后，由于目标体积较大，亮度不均，导致目标区域错误地被分割为两个区域，从而导致虚警率的升高。与其他几种方法相比，本发明所述的检测方法具有明显优势。表1不同方法在图4(a)序列中阈值参数k取不同值时的检出率与虚警率表2本发明在图4(b)序列中阈值参数k2取不同值时的检出率与虚警率k21112131415161718191101pd100％100％100％100％100％100％100％95％82％55％37％fa67％0％0％0％0％0％0％0％0％0％0％表3wldm在图4(b)序列中阈值参数k取不同值时的检出率与虚警率k0.00.10.20.30.40.50.60.70.80.91.0pd100％100％100％100％97％97％95％95％92％90％0％fa100％30％5％17％35％57％72％72％65％50％0％需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。以上所述是本发明的优选实施方式，应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。当前第1页12