基于DSP的序列图像弱小运动目标检测

摘　要：针对强杂波背景中的弱小运动目标检测问题，在分析弱小目标的运动特性在时频变换中表现出显着差异的基础上，采用两级滤波方法检测运动目标，提出了一种基于DSP的序列图像弱小运动目标检测方法。该方法首先对序列图像进行预处理、设置门限过滤噪声，然后根据目标和杂波的时频特征差异识别出弱小运动目标。实验结果表明，该方法可以很好地抑制噪声的干扰，并能检测出复杂背景下的弱小运动目标。
关键词：时频分析；DSP；弱小目标检测；序列图像
       弱小运动目标检测的应用涉及到国防、交通和安全等各个领域。早期的先检测后跟踪的处理技术DBT(Detect before Track)已经不适应弱小目标探测与提取的要求。近年来，边检测边跟踪技术TBD(Track before Detect)逐渐成为弱小目标探测跟踪技术发展的主流。该类方法主要包括三维匹配滤波、多级假设检验序贯处理、动态规划、高阶相关处理方法等。此外，目前还有利用模糊推理、神经网络[1]、小波变换、分维分析等手段用于检测与跟踪弱小目标。但几乎所有这些边检测边跟踪算法都很少分析目标与噪声在空间、时间、频率及其联合多维度上表现出的细微特征，使得只能单纯利用目标信号轨迹上的相关性，即轨迹的连续性和一致性，通过类似于穷尽搜索的方式来检测目标轨迹。
     事实上，目标具有多样化的特征(包括物理特征和数学特征),而强度特征仅仅是其诸多特征中的一部分，充分利用信号的多方面特征信息可有效提高系统对目标和杂波的分辨能力[2]。本文主要创新点在于：利用弱小运动目标在时频多维联合域中分析目标和杂波的特征差异，系统分析运动目标和杂波的时频特征，指出弱小目标经过处是长时非平稳而短时平稳的随机信号，而杂波所在处通常为长时平稳信号。利用运动目标经过处时频发生的微小变化来检测目标，而杂波处的时频基本保持平稳的差异，有效去除背景杂波的干扰，于此检测出弱小运动目标。最后通过DSP实现弱小运动目标检测算法。
1 检测算法分析
     弱小运动目标的序列图像f(x,y,k)可表述为：
     
其中，s(x,y,k)、b(x,y,k)和n(x,y,k)分别为第k帧经过像素(x,y)的目标、杂波和噪声的幅度值，k为图像序列的帧数。在图像大小、帧数等参数已经确定的条件下，首先调用目标检测模板对序列图像进行处理，根据目标经过处的像素点会引起灰度值的瞬时起伏,而噪声点是服从时间独立的高斯分布这一特点，设置相应判决门限并根据判决准则得到一系列疑似目标点，然后根据目标与杂波在时频上存在的明显差异识别运动目标。

其中，f(xi,yi)为模板上像素点的灰度值。检验统计量T大于设置门限Fm时，则认为此点为疑似目标点，反之则为噪声。
　　一般情况下，采用的检测模板大小为9个像素，如图1所示。考虑到目标点及其他干扰所占像素值有可能大于实际所选用模板，所以做了适当的改进,将疑似目标点属于同一疑似目标斑块的像素点进行合并，如图2所示。
   
       假设Fi为3×3模板上任意一点,Fim为3×3模板中心点,将满足以下条件的疑似像素点进行合并：
   
     式(3)表明两个3×3模板上任意两点八邻域距离小于L1的情况(L1取n1个像素),式(4)表明两个3×3摸板上两中心点八邻域距离小于L2的情况(L2取n2个像素)。其中n1、n2根据实际目标点、杂波大小确定。式(5)表明将满足式(3)、式(4)的两个疑似目标点合并为一个新的疑似目标点，重复进行上述合并操作，直至将所有疑似目标点合并完为止。
1.2 时频分析及其实现
     时频检测方法体现了信号在时域和频域的联合分布信息，可为分析及分辨目标和杂波提供细微的空间、时间、频率及其联合多维度的特征。
     假设序列图像(x,y)按采样时间先后顺序形成三维图像(x,y,k)，将所有待处理原始序列图像中同一位置空间坐标为(x,y)的灰度值读取出来并按照时间的先后顺序存储为f(m)，即形成f(x,y,k)→f(m)的映射。然后对f(m)信号进行如下短时加窗处理:

其中, t为时频变换后的时间，如果假设图像序列中相邻帧之间的时间间隔为?驻t，则时间t即为?驻t×k;m为序列图像的第k帧，w为短时傅里叶变换之后的频率，与采样点数有关;W(N1)为短时窗函数，一般为能量聚集性强的汉明窗或汉宁窗，N1为短时窗函数的窗长。
  
2 系统结构
     系统结构如图3所示。首先将序列图像进行预处理(可以采用FPGA等逻辑器件处理)，使得图像序列每帧大小一致、时间间隔恒定，然后通过系统核心芯片DSP进行算法实现(即滤除噪声干扰、识别出疑似目标点),最后根据运动目标的空时频特性的差异识别出运动目标，输出检测结果[3]。