一种基于量子均值漂移的抗遮挡多运动车辆追踪方法
【专利摘要】本发明涉及一种基于量子均值漂移的抗遮挡运动多车辆追踪方法，包括如下步骤：1)读取第一帧图像，并用鼠标依次圈出要追踪的车辆。提取每一辆车所处区域图像的像素值，以此为基础建立均值漂移算法的数学模型；2)读取下一帧图像，提取每一辆车候选区域图像的像素值；3)判断车辆是否发生遮挡；4)如果车辆没有被遮挡，则以初始帧和当前帧灰度图像的像素点坐标为基础建立量子进化算法模型，搜索出当前帧车辆的最优位置。将当前帧搜索得到的车辆中心位置与前一帧的车辆位置进行比较，如果满足条件，则当前帧的车辆中心位置即为所追踪的车辆的中心位置，从而实现对视频图像中多车辆的实时追踪。
【专利说明】
-种基于量子均值漂移的抗遮挡多运动车辆追踪方法
技术领域
[0001] 本发明设及一种图像处理技术领域的抗遮挡运动多车辆追踪方法，具体设及一种 基于量子均值漂移的抗遮挡多运动车辆追踪方法。
【背景技术】
[0002] 存在遮挡的运动多车辆追踪一直是计算机视觉中具有挑战性的难题。目标遮挡的 主要特征是目标信息的逐渐丢失，因此追踪算法的关键在于捜索到足够多的目标信息，从 而较为准确地判断目标所在位置。遮挡增加了目标追踪的难度，导致目标追踪的不稳定性 甚至目标丢失。一个有效的追踪方法应该能够准确的判断遮挡的发生，并且当目标发生遮 挡时仍可W利用目标剩余的信息继续追踪目标。
[0003] 为了解决目标遮挡问题，人们提出了很多算法。王展青等人（见：王展青，凡友福， 张桂林.跟踪遮挡目标的一种鲁棒算法[J].计算机工程与应用，2007, 43 (27) :50-53)提 出了根据卡尔曼滤波器残差的大小来判断遮挡的方法，发生完全遮挡情况时，结合目标的 运动方向使用六点捜索策略来找回目标，但是如果目标运动方向发生变化时，目标丢失；田 伟等人（见：田伟，罗予频，华成英.遮挡情况下多目标跟踪的一种新方法[J].计算机工 程与应用，2007, 43(13) :33-36)提出了一种方法：将遮挡情况作为目标的一个状态分量， 在对遮挡估计的基础上计算目标的观测置信度，该方法对遮挡不严重的目标具有很好的效 果。
[0004] 均值漂移是一种核密度估计的图像特征分析方法，基于最优梯度下降的方法，通 过迭代的方法捜索目标，即先算出当前点的漂移均值点，移动该点到其漂移位置，然后W此 为新的起始点，继续移动，直到飘移到密度函数的局部极大值点或满足一定的条件为止，实 现对运动目标的追踪，能够实时追踪非刚性目标，对目标的变形、旋转等运动具有较好的适 用性。但是由于均值漂移方法在目标追踪过程中没有利用目标在空间中的运动方向和运动 速度信息，当周围环境存在干扰时，仅使用均值漂移容易丢失目标；而且，由于受均值漂移 方法理论本身原理的局限性，它只能对特征点近距离的像素点的概率密度函数进行泰勒展 开，当下一帖图像对应特征点偏离当前帖对应特征点较远时，均值漂移方法就会出现很大 误差，造成追踪失败。量子进化算法是最近发展起来的一种概率进化算法，具有很大的研究 价值。它W量子计算概念和原理为基础，用量子码和量子口作为更新算子来完成进化捜索。 目前量子进化算法已经应用到数值优化、组合优化、图形图像处理、电路设计、通信、多目标 优化等领域。量子进化算法建立在量子的态矢量表述基础上，用量子比特的几率幅来表示 染色体的编码，使得一条染色体可W表达多个态的叠加，并利用量子旋转口和量子非口实 现染色体的更新操作，从而实现目标的优化求解。量子进化算法是解决优化问题的一种有 效方法，具有非常高效的信息并行处理能力，但是由于它本身具有很大的随机性，使其在实 际应用中存在收敛速度慢、退化等问题，影响了其应用效果。
[0005] 本发明的目的在于针对均值漂移和量子进化算法的不足，提供一种基于量子均值 漂移的抗遮挡多运动车辆追踪方法，采用基于颜色直方图均值漂移并结合量子进化算法对 图像特征点进行捜索，同时引入帖间差分法判断目标的遮挡情况并线性预测当目标发生遮 挡时目标在当前帖的位置，有效克服了当下一帖图像对应特征点偏离当前帖对应特征点较 远时均值漂移方法出现的较大误差和目标发生遮挡时目标信息丢失的难题，取得了较好的 追踪效果。
【发明内容】
[0006] 本发明针对现有均值漂移理论不足，采用均值漂移与量子进化算法相结合来处理 动态背景下多辆车的追踪问题，同时引入了差分算法，既提高了目标追踪的准确率和实时 性，又有效解决了目标追踪过程中由于遮挡而识别效率低的问题，具有更高的可行性。本发 明的具体步骤如下：
[0007] step 1 :完成对视频图像信息的采集，读取第一帖图像，并用鼠标依次圈出要追踪 的车辆。提取每一辆车所在区域图像的像素值。
[0008] Step 2 :设车辆所在区域像素的中屯、坐标为X。，第i个像素的坐标用Xi表示，则目 标模型的第U种颜色分布特征值的概率为：
[0009]
[0010] 其中U= l，2，...，m，k(x)为核函数，h表示核函数的带宽，函数b(Xi)和5 (X)的 作用是判断Xi的颜色值是否属于第U特征值，C是一个标准化的常量系数，使得所有特征值 的概率和为1。
[0011] step 3 :读取下一帖图像，车辆可能存在的图像区域即目标候选区域，提取每一辆 车在目标候选区域图像的像素值。
[001引St巧4 :假设候选区域图像模型中屯、坐标为y。，XiQ = 1，2,…，而）表示该区域内 每个像素点的坐标，则候选模型的第U种颜色分布特征值的概率为：
[0013]
[0014] Step 5 :相似性函数是用来表示候选模型和目标模型之间的关系，即目标模型与 候选模型之间的相似度。在基于均值漂移的多目标追踪算法中，广泛使用化attacharyya 系数作为相似性函数。定义两个离散分布之间距离：
[0015]
[0016]
[0017]
[0018] St巧6 :建立卡尔曼滤波器模型：
[0019] 离散状态方程：x化）=A化)x化-1) + ?化-1)
[0020] 式中X化）表示状态变量构成的多维状态矢量，A和B是由系统特性确定的矩阵。
[0021] 离散量测方程：yk= C A+Vk。
[00过其中Sk= C A表示信号的真值，噪声Vk表示量测时引入的误差的随机向量。
[0023] 由下式判定车辆是否被严重遮挡。当f(k) = 1，则车辆处于无遮挡或者遮挡不严 重的状态，当f(k) = 0,则表明车辆处于严重遮挡或完全遮挡状态。
[0024]
[00巧]St巧7 :如果发生严重遮挡，利用利用运动状态线性预测车辆在当前帖中的位置 y。，并采用卡尔曼滤波捜索策略确认当前帖中车辆位置yi。
[0026] step 8 :如果车辆没有被严重遮挡，W前初始帖和当前帖灰度图像的像素值为基 础建立量子进化算法模型：
[0027] 1) W当前帖图像中候选模型区域的灰度图像的像素值作为初始种群
为吉'即：[0029]
，n为种群大小；m为量子染色体的长度；为定义如下的染色体： / = 1么.,，n.，进化代数初始化T = 0 ;
[002引。初始化种群Q(t)，若将全部染色体中的如，《（W，2,..，班）W及均置
[0030] 则表示在t = 0时刻，所有可能的线性叠加态W相同的概率出现。
[0031] 3)对种群Q(t)进行观测，每一个量子比特都会由计算机对应的产生一个[0， 1]之间的随机数r，当Kl 0i|2时，该位的观测态是"1"，否则是"0"。生成观测态种群： 巧f)=站成''.，4}，其中4( / = 1，2，...，h )是染色体9; </=1，2,..:，《 )观测产生的一个二进制串， 其长度是m。
[0032] 4)对生成的观测态种群P(t)中的所有个体评价其适应度值，保持适应度较高的 个体。
[0033] 5)对种群P(t)采用量子旋转口 U( 0 1)进行量子口操作，更新P(t)中适应度低的 个体，保存最优解。
[0034] 6)令进化代数T = T+1，转至3)续执行算法，直到满足停止条件为止，输出最优解 東1。
[003引 Step 9 :由Step 4和5计算计算候选模型输说)沁=1.....崎,并估计相似居
计算候选目标知执)地=1..,崎，并估计相似J
[003引 Step 10 :如勇
，继续计算幻W'l)，引。如果 li'i - i'ul < E '输出知，算法结束，输出最优解；否则知叫',，转至St巧4。
【附图说明】
[0037] 图1为算法流程图；
[0038] 图2为仿真结果。
[0039] 具体实施步骤
[0040] 如图1所示，【具体实施方式】如下：
[0041] Step 1 :完成对视频图像信息的采集，读取第一帖图像，并用鼠标依次圈出要追踪 的车辆。提取每一辆车所在区域图像的像素值。
[0042] Step 2 :设车辆所在区域像素的中屯、坐标为X。，第i个像素的坐标用Xi表示，则目 标模型的第U个颜色分布特征值的概率为：
[0043]
[0044] 其中U= 1，2，...，m，k (X)为核函数，h表示核函数的带宽，函数b (Xi)和5 (X)的 作用是判断Xi的颜色值是否属于第U特征值，C是一个标准化的常量系数，使得所有特征值 的概率和为1。
[0045] Step 3 :读取下一帖图像，车辆可能存在的图像区域即目标候选区域，提取每一辆 车在目标候选区域图像的像素值。
[004引 St巧4 :假设候选区域图像模型中屯、坐标为y。，XiQ = 1，2,…，而）表示该区域内 每个像素点的坐标，则候选模型的第U个颜色分布特征值的概率为：
[0047]
[0048] Step 5 :相似性函数是用来表示候选模型和目标模型之间的关系，即目标模型与 候选模型之间的相似度。在基于均值漂移的多目标追踪算法中，广泛使用化attacharyya 系数作为祐柄化品吿切证^窗部分布之间距离：
[0049]
[0050]
[0051]
[0052] Step 6 :建立卡尔曼滤波器模型：
[0053] 离散状态方程：x化）=A化)X化-1) + ?化-1)
[0054] 式中X (k)表示状态变量构成的多维状态矢量，A和B是由系统特性确定的矩阵。
[0055] 离散量测方程：yk = C A+Vk。
[005引其中Sk= C A表示信号的真值，噪声Vk表示量测时引入的误差的随机向量。
[0057] 由下式判定车辆是否被严重遮挡。当f(k) = 1，则车辆处于无遮挡或者遮挡不严 重的状态，当f(k) = 0,则表明车辆处于严重遮挡或完全遮挡状态。
[0058]
[0059] Step 7 :如果发生严重遮挡，利用利用运动状态线性预测车辆在当前帖中的位置 y。，并采用卡尔曼滤波捜索策略确认当前帖中车辆位置yi。
[0060] Step 8:如果车辆没有被严重遮挡，W前初始帖和当前帖灰度图像的像素值为基 础建立量子进化算法模型：
[0061] 1) W当前帖图像中候选模型区域的灰度图像的像素值作为初始种群
去'即：[0063]
3 n为种群大小；m为量子染色体的长度；接为定义如下的染色体： ./ = 1，2，....，《::，:进化代数初始化了 = 0;
[006引。初始化种群Q (t)，若将全部染色体中的如處（H2,。.典;W及為U=巧，...:，々)均置为
[0064] 则表示在t = 0时刻，所有可能的线性叠加态W相同的概率出现。
[0065] 3)对种群Q(t)进行观测，每一个量子比特都会由计算机对应的产生一个[0， 1]之间的随机数r，当r<| 0i|2时，该位的观测态是"1"，否则是"0"。生成观测态种群： 巧〇 =林1'，如'''，.<):，其中4(^ =化'.'，巧)是染色体9;(7' =化''.?)观^1产生的一个二进审。 串，其长度是m。
[0066] 4)对生成的观测态种群P(t)中的所有个体评价其适应度值，保持适应度较高的 个体。
[0067] 5)对种群P(t)采用量子旋转口 U( 0 1)进行量子口操作，更新P(t)中适应度低的 个体，保存最优解。
[0068] 6)令进化代数T = T+1，转至3)继续执行算法，直到满足停止条件为止，输出最优 解玄1。
[006引 St巧9 :由St巧4和5计算计算候选模型巧,,GVH" =1.并估计相似度
。计算候选目标化,化)他=1..."!)，并估计相似度
[0070] St巧10:如果
继续计算^^.1)，41。如果 |家1' 爲，:输出A，.算法结束，输出最优解；否则.? 知，转至St巧4。
【主权项】
1. 一种基于量子均值漂移的抗遮挡多运动车辆追踪方法，其特征在于，包括以下步 骤： (1) 完成对视频图像信息的采集，读取第一帧图像，并用鼠标依次圈出要追踪的车辆； 提取每一辆车所处区域图像的像素值； (2) 设车辆所在区域像素的中心坐标为X。，第i个像素的坐标用X1表示，则目标模型的 第u种颜色分布特征值的概率为：其中u = 1，2, ...，m，k(X)为核函数，h表示核函数的带宽，函数Mxi)和δ (X)的作 用是判断^处的颜色值是否属于第u特征值，C是标准化常量系数，使得所有特征值的概率 和为1 ; (3) 读取下一帧图像，提取每一辆车候选区域（车辆可能存在的区域称为候选区域）图 像的像素值； (4) 假设候选模型中心坐标为y(]，Xl表示该区域内每个像素点的坐标，i = 1，2,. . .，nk， 则候选模型的第u种颜色分布特征值的概率为：(5) 相似性函数用来表示候选模型和目标模型之间的关系，即目标模型与候选模型之 间的相似度；在基于均值漂移的多目标追踪算法中，广泛使用Bhattacharyya系数作为相 似性函数定义两个离散分布之间距离：则相似度系数为(6) 卡尔曼滤波器模型： 离散状态方程：x(k) = A(k)x(k-l) + co (k-1) 式中x(k)表示状态变量构成的多维状态矢量，A和B是由系统特性确定的矩阵； 尚散M测方程：yk= Ckxk+vk; 其中Sk= C Α表示信号的真值，^表示测量时引入的随机噪声；通过上式判定车辆是宵骰严重趣Θ ;刍fUO = 1，则半辆处t尤趣挡或者遮挡不严重 的状态，当f(k) = 0,则表明车辆处于严重遮挡或完全遮挡状态； (7) 如果发生严重遮挡，利用运动状态线性预测车辆在当前帧中的位置y。，并采用卡尔 曼滤波搜索策略确认当前帧中车辆位置y1; (8) 如果车辆没有被严重遮挡，以前初始帧和当前帧灰度图像的像素值为基础建立量 子进化算法模型，量子进化算法步骤如下： 1) 以当前帧图像中候选模型区域的灰度像素值作为初始种群e(〇 = _!"…￥丨，其中 η为种群大小；m为量子染色体的长度；iy丨为定义如下的染色体：： 1，2··· n，进化代数初始化T = O ; 2) 初始化种群Q(t)，若将全部染色体中的( / = U...,W )以及f) ( )均!P ：则意味着在t = O时刻，所有可能的线性叠加态以相同的概率出现； 3) 对种群Q(t)进行观测，每一个量子比特都会由计算机对应的产生一个[0，1] 之间的随机数r，当r〈| P112时，该位的观测态是"1"，否则是"0"。生成观测态种群： = ,其中< (_/ = 1，2…《 )是染色体f). U = 1A:.."h:)观测产生的一个二进制 串，其长度是m; 4) 对生成的观测态种群P (t)中所有个体计算适应度值，保留适应度较高个体； 5) 对种群P(t)采用量子旋转门U( Θ J进行量子门操作，更新p(t)中适应度低的个体， 保存最优解； 6) 令进化代数T = T+1，转至步骤3)继续执行算法，直到满足停止条件为止，输出最优 解>1 ; (9)由⑷和（5)计算候选模型{Α,(??)Κ? = 1.』）_，并估计相似度4夂H十 算候选目标(?ι)丨(μ = 1....w)，并估计相似度U=I (?ο)如果p[Kh)，釘《4Kb)名，令尺+Α)，继续计算/必 ΙΙΛ - ??| < "输出?,，算法结束；否则％ ^，转至⑷。2.根据权利要求1所述的一种基于量子均值漂移的抗遮挡多运动车辆追踪方法，其特 征在于，卡尔曼搜索策略的具体步骤如下： 在第k帧，车辆运动到B处时被部分遮挡，在k+Ι帧时车辆则处于完全被遮挡状态，此 时进入搜索阶段；利用卡尔曼滤波和平均速度V预测出第k+Ι帧中车辆的位置A，以B为中 心点，r为半径： 式中h是核函数的带宽； 在右半圆周上，设置5个等间距的搜索块，与B -起组成6个候选块。为了避免重复搜 索，车辆运动速度较小时，取半径为h ;车辆运动速度大于h时，为了增强搜索效率，取半径 为M ;当搜索块中有满足取|辦幻，则认为遮挡结束，恢复对车辆进行正常追踪；当 搜索不到车辆时，依据运动状态估计，认为车辆位于A位置，搜索结束，进入下一帧。
【文档编号】G06T7/20GK106023242SQ201510167462
【公开日】2016年10月12日
【申请日】2015年4月9日
【发明人】汤雅连, 程平
【申请人】广东易富网络科技有限公司