基于图像处理技术的红外小目标的检测与跟踪(共11页)

1、对于(duy)低慢小目标的侦测和处理(chl)对策 -基于图像处理技术的红外小目标的检测(jin c)与跟踪摘 要：验证了一种基于红外小目标视频图像序列的跟踪算法，主要研究了基于形心计算的跟踪方法和基于kalman滤波器多帧数据关联方法的跟踪法。分别仿真验证，并从实现结果出发得出了两种算法的适用范围和各自存在的不足。0 引言 随着近十几年信息技术的飞速发展，计算机硬件的处理能力不断提高，存储成本大幅下跌，一些研究人员开始重点研究计算机视觉中有关运动的问题。与处理单幅图像相比，图像序列引入了新的时间维以及时间相关性约束，这一额外的约束激发了人们对视频理解的研究.视频序列目标跟踪是指对传感器摄取到

2、的图像序列进行处理与分析，一旦目标被确定，就可获得目标的特征参数选择。 。 由于视频跟踪具有广泛的应用范围，因而引起了世界范围内广大研究者的兴趣。 在1996年至1999年间，美国国防高级研究项目署（DARPA）资助卡内基梅隆大学、戴维SARNOFF研究中心等著名大学和公司合作，联合研制视频监视与监控系统VSAM，主要研究目的是开发用于战场及普通民用场景的自动视频理解技术。DARPA在2000年又资助了重大项目HID计划，其任务是开发多模式的监控技术以实现远距离情况下人的检测、分类和识别，以增强国防、民用等场合免受恐怖袭击的保护能力。在2008年，DARPA资助了一项研究实时流视频监视的技术。

3、 国内许多研究所和大学也投入了大量的精力致力于图像跟踪的研究。如清华大学人机交互与媒体集成研究所在人脸跟踪、视觉监控、等方面取得了许多科研成果；北京自动化研究所模式识别国家重点实验室在交通场景监控、人体跟踪、智能轮椅手势导航等领域进行了深入研究。 1 设计(shj)背景视频序列目标跟踪是指对传感器摄取到的图像序列进行处理与分析，充分利用传感器采集得到信息来对目标进行稳定跟踪的过程。一旦目标被确定，就可获得目标的位置、速度、加速度等运动(yndng)参数，进而获得目标的特征参数。 运动分析的目的是根据运动线索揭示关于环境有价值的信息，这些信息可以用来执行更高级的视觉任务，如活动识别、视觉检索等。

4、由于采集到的图像序列通常(tngchng)是二维图像，许多研究者开始使用简单的二维表示处理运动分析问题，避免了恢复三维世界的结构、属性等信息。通过视频目标跟踪系统通过对图像的去噪、对图像进行灰度变换、图像融合、图像增强等处理，可以快速的筛选并提取有用的信息。在军事上，视频序列目标跟踪技术广泛应用于精确制导、战场机器人自主导航、无人机助降，靶场光电跟踪等领域。在民用上，该技术主要应用在智能视频监控、智能交通管制、医疗影像诊断等方面。 2 总体方案2.1 常用的跟踪算法2.1.1基于特征的目标跟踪方法基于特征的方法利用了特征位置的变化信息,首先从图像序列中抽取显著特征,然后在序列图像上寻找特征点的

5、对应关系,特征匹配算法大都引入了刚体约束条件,已有的技术包括结构匹配、树匹配等。之后利用序列帧图像中目标特征点的对应关系,解算当前帧图像中目标位置。2.1.2基于(jy)相关的目标跟踪方法基于相关的目标跟踪方法是把一个预先存储(cn ch)的目标样板作为识别和确定目标位置的依据,然后用目标样板与实时图像中的各个(gg)子区域图像进行比较,找出和目标模板最相似的一个子图像位置,作为当前目标的位置,这就是相关跟踪的基本思想,这种方法也叫做图像匹配。2.1.3基于对比度的目标跟踪方法基于对比度的目标跟踪方法又称波门跟踪方法,该方法适用于目标和背景具有明显对比度的目标跟踪。该方法需设计一个波门,波门的

6、尺寸略大于目标尺寸,使目标不受波门外的背景和噪声干扰的影响。跟踪波门可分为固定式和自适应式两种。前者在跟踪过程中波门的大小始终不变;后者则是在跟踪目标的过程中波门随目标的大小变化而变化。波门跟踪算法可分为矩心跟踪算法、边缘跟踪算法、双边缘跟踪算法和区域平衡跟踪算法等。2.2 本方案实现的技术难点要解决图像序列的去除噪声问题，由于目标过小故要选择适当的去除噪声方法，对融合了多种噪声的图像要选用综合的去除噪声方法。最终的跟踪系统要实现自适应滤波。要解决小目标丢失或被遮挡后再次出现的问题，要保证目标再次出现后的准确识别和跟踪。3 原理(yunl)描述3.1 形心跟踪(gnzng) 视频图像序列 滤波

7、去噪 背景提取 二值化处理 设置波门 计算形心 目标物体跟踪 兴趣区获取 灰度转换摄像头标定物体保证在镜头内 图(1) 形心跟踪(gnzng)系统流程图如图（1）所示，形心跟踪首先要对图像进行灰度化处理，使彩色图像转换为8位灰度图像以方便以后的运输和处理过程。3.1.1滤波(lb)去噪本次所用图像序列的噪声课近似认为是椒盐噪声，在去噪声时可以用中值滤波法，本方法可以有效的消除(xioch)椒盐噪声并保持图像中的边缘和细节。中值滤波法是把领域内的所有像素按顺序排列，然后用中间值代替中心元素的灰度值作为输出。二维中值滤波(lb)定义为： （1）其中；f(x,y) 为原始图像阵列； g(x,y) 为

8、中值滤波后图像阵列；Median为中值滤波算子，取中值；A 为滤波窗口，大小为kl。中值滤波的缺点在于运算速度随窗口的扩大而变慢，故为了提高检测速度本次验证时采用的滤波窗口为窗口。3.1.2二值化处理二值化处理是利用同一区域具有某种共同的灰度特性进行图像的分割，其基本原理就是选取一个适当的灰度阈值设置一个阈值，然后对每个像素进行检测，当检测的像素的灰度值小于或等于阈值时，将此像素的灰度值变为0，当检测的像素的灰度值大于阈值时，将此像素灰度值变为255。 （2） 其中是设定的阈值；是像素原有的灰度值；是处理后的灰度值。本次设计我们所读取的目标跟踪序列，其背景与目标的灰度基本保持不变。对此类视频序

9、列图像的阈值确定可以采用直方图分布的办法。即初步了解到所要跟踪的目标区域和目标所在的背景区域各自处在一个灰度值的范围，在直方图的表现上为具有双峰分布，通过选取波分之间的低谷所对应的灰度值作为分割阈值就可以将目标区域的背景区域分割开来。3.1.3形心计算(j sun)对于二值化图像在假设物体的面密度均匀分布的境况下，形心就是要探测物体的几何中心。同过形心的计算可以快速的判断目标物体的具体位置(wi zhi)。当目标姿态改变时，形心的位置变动较小。因此形心跟踪比较平稳，同时由于系统效应，算法的抗干扰能力也较强。对于离散的数字图像，其形心公式可以(ky)定义为： （3）其中， 、是目标形心坐标；是图

10、像函数（即图像上处像素点的灰度）；是横坐标范围；是纵坐标范围。3.1.4添加波门算法在计算出目标形心的基础上，以形心为对称中心分别在上下，左右画线，连接在一起形成一个矩形框，即跟踪波门。波门的大小有目标大小而定，主要目的是可以标出目标的位置。过大的波门有较高的容错性，可以减小跟踪算法的精度不够而产生的跟踪偏差，但由于波门较大所以跟踪精度随之降低。3.2 基于(jy)kalman滤波器多帧数据关联(gunlin)方法 图(2) kalman滤波(lb)流程图Kalman滤波器是一个线性递归滤波器,用于对动态系统的状态序列进行线性最小方差误差估计，从而对下一个状态的运动做最优估计。预测时具有无偏、

11、稳定和最优的特点。3.2 .1 kalman滤波首先要引入一个离散控制过程的系统。该系统可用一个线性随机微分方程来描述： X(k)=A X(k-1)+B U(k)+W(k) （4）再加上系统的测量值：Z(k)=H X(k)+V(k) （5）上两式子中，X(k)是k时刻的系统状态，U(k)是k时刻对系统的控制量。A和B是系统参数，对于多模型系统，他们为矩阵。Z(k)是k时刻的测量值，H是测量系统的参数，对于多测量系统，H为矩阵。W(k)和V(k)分别表示过程和测量的噪声。他们被假设成高斯白噪声，他们的协方差分别是Q，R。之后来用他们结合他们的协方差来估算系统的最优化输出。利用(lyng)系统的过

12、程(guchng)模型，来预测下一状态的系统(xtng)。假设现在的系统状态是k，根据系统的模型，可以基于系统的上一状态而预测出现在状态：X(k|k-1)=A X(k-1|k-1)+B U(k) (6)其中，X(k|k-1)是利用上一状态预测的结果，X(k-1|k-1)是上一状态最优的结果，U(k)为现在状态的控制量，如果没有控制量，它可以为0。P(k|k-1)=A P(k-1|k-1) A+Q (7)其中，P表示协方差，P(k|k-1)是X(k|k-1)对应的协方差，P(k-1|k-1)是X(k-1|k-1)对应的协方差，A表示A的转置矩阵，Q是系统过程的协方差。式子1，2就是卡尔曼滤波器5

13、个公式当中的前两个，也就是对系统的预测。结合预测值和测量值，可以得到现在状态(k)的最优化估算值X(k|k)：X(k|k)= X(k|k-1)+Kg(k) (Z(k)-H X(k|k-1) (8)其中Kg为卡尔曼增益：Kg(k)= P(k|k-1) H / (H P(k|k-1) H + R) (9)为了要另卡尔曼滤波器不断的运行下去直到系统过程结束，还要更新k状态下X(k|k)的c协方差：P(k|k)=（I-Kg(k) H）P(k|k-1) (10)其中I为1的矩阵，对于单模型单测量，I=1。当系统进入k+1状态时，P(k|k)就是式子(2)的P(k-1|k-1)。3.2 .1 联合(lin

14、h)数据关联在第1次到第k次扫描所获得的全部(qunb)有效测量已知的情况下，第k次扫描(somio)时，第i（）个观测时间点内均为正确的测量的概率，称之为正确关联概率用表示。 （11）公式（4）中，第k次扫描第i个测量为正确测量的事件；第1次到第k次扫描所获得的全部有效测量的集合； 第k次测量所获得的测量数目。 （12）公式（5）中， ， 是有效测量皆来自目标的条件下的目标状态估计值。3.2 .1 确认矩阵引入了确认（聚）矩阵，用来反映聚集内扫描的所有目标和测量的关系 ，设m和n为某一具体聚集内点迹（测量）和目标数，则m（n+1）聚矩阵定义 （13）其中是是、一个二值变量（取0或1），用以表

15、示量测j是否落入目标t的确认门之内。t=0表示没有目标，测量源于杂波或虚警，。 图（3）两目标(mbio)示意图如图（3）2个目标(mbio)（t1、t2）和3个测量(cling)（z1、z2、z3），其确认矩阵为： 4分析1 对比度矩心跟踪法通过对原图像进行图像阈值分割、形心计算，最终实现目标的跟踪。此种算法具有简单、快速和实时性高的特点。但是，此种算法的实用性范围极其有限，一般来说只适用于目标图像与背景图像的灰度值相差较大的情形下才能实现，否则将得不到目标的形心。除此之外，该算法受图像成像过程中引入的噪声的影响比较大，在环境较为恶劣的情况下将会失去对目标的有效跟踪。因而，对比度矩心跟踪法面

16、临背景环境的限制以及成像系统噪声的限制。2基于kalman滤波器多帧数据关联方法不必对原图像进行预处理，基本保留了图像的全部信息，这对于目标所在的背景区域的了解至关重要；其次，此跟踪算法非常适合于在环境复杂多变的情况下进行目标的跟踪。但不足之处是由于在相关匹配的过程中涉及到大量的数学运算，因而系统跟踪的实时性将大大降低，为进一步提高系统的实时响应，相应的处理器、存储器等硬件设备的性能要求较高。即使在匹配的过程中进一步降低搜索区域，但那样有可能得到的是局部最优值，其最终结果是真正的目标没有被跟踪到。因而，基于kalman滤波器多帧数据关联方法面临实时性的提高以及快速搜索区域算法的进一步优化。主要

17、(zhyo)参考文献1 张然，吕高杰，张国华 光电目标图像(t xin)自动跟踪技术研究J， HYPERLINK /kns50/Navi/Bridge.aspx?DBCode=CJFD&LinkType=BaseLink&Field=BaseID&TableName=CJFDBASEINFO&NaviLink=%e7%94%b5%e5%85%89%e4%b8%8e%e6%8e%a7%e5%88%b6&Value=DGKQ t _blank 电光(dingung)与控制, HYPERLINK /kns50/Navi/Bridge.aspx?DBCode=CJFD&LinkType=IssueLi

18、nk&Field=BaseID*year*issue&TableName=CJFDYEARINFO&Value=DGKQ*2008*09&NaviLink=%e7%94%b5%e5%85%89%e4%b8%8e%e6%8e%a7%e5%88%b6 t _blank 2008,15(9)：65-682 王唯合，王宏志，赵兰英 运动目标检测跟踪方法研究，长春工业大学学报, HYPERLINK /kns50/Navi/Bridge.aspx?DBCode=CJFD&LinkType=IssueLink&Field=BaseID*year*issue&TableName=CJFDYEARINFO&Value=DGKQ*2008*09&NaviLink=%e7%94%b5%e5%85%89%e4%b8%8e%e6%8e%a7%e5%88%b6 t _blank 2007,28(4)：411-4143 张岩; 崔智社; 龙腾; 图象序列中机动目标的形心跟踪， HYPERLINK /KNS50/Navi/Bridge.aspx?LinkType=BaseLink&DBCode=cjfd&TableName=cjfdbaseinfo&Field=BaseID&Value=HKXB&NaviLink