基于变分光流实现车辆跟踪的方法

基于变分光流实现车辆跟踪的方法
【技术领域】
[0001] 本发明涉及信息处理技术领域，尤其涉及智能交通中车辆跟踪的技术领域，具体 是指一种基于变分光流实现车辆跟踪的方法。
【背景技术】
[0002] 智能交通视频监控算法主要包含车辆检测、车辆跟踪和车辆行为分析三个步骤。 车辆检测实现初始交通场景中车辆目标的提取；车辆跟踪对提取的车辆目标进行实时的运 动描述，并生成各自的运动轨迹；行为分析基于产生的运动轨迹的规律来判断车辆是否发 生违章行为或处于异常状态。
[0003] 车辆跟踪技术具有重要的应用价值和发展前景，近年来发展迅速。车辆跟踪试图 在各帧图像之间确定目标相关信息参数的相互关系，通过前、后帧之间的对应匹配来获取 目标的轨迹信息。常见跟踪算法可归为以下四类。
[0004] (1)点跟踪法：包括单点跟踪和多点跟踪，主要利用目标历史位置、速度等信息实 现邻帧之间的目标关联，如kalman滤波法，光流法，SIFT匹配法。
[0005] (2)核跟踪法：核跟踪法通过对车辆外观模型在连续帧之间进行匹配来计算目标 的运动，包括基于模版的方法、基于概率模型的方法和多视觉模型的方法。这类算法涉及 的三个基本要素是目标外观模型、搜索策略和相似性度量。常采用的搜索策略有穷举法、 目标函数优化法和统计性方法等；而相似性度量标准一般采用平方差和（Sum of Squared Difference，SSD)、互相关系数、正则化相关系数等。基于模版的方法实质上就是进行模版 匹配。基于概率模型的方法将目标区域的特征用概率分布模型表示，通过模型匹配实施跟 踪，如基于颜色直方图的目标跟踪、基于LBP纹理特征的目标跟踪、基于GM(1，1)模型的跟 踪、基于马尔可夫随机场|旲型的跟踪、基于均值漂移的跟踪、基于质心迭代的跟踪、基于粒 子滤波器的跟踪和基于多特征自适应融合的目标跟踪。其中基于多特征自适应融合的目标 跟踪常采用的特征包括点、线、角点、颜色等。基于多视觉模型的方法能够适应目标姿态变 化引起的外观变化，常见的有基于主成份分析（PCA)和基于分类器的方法，如Adaboost分 类器。
[0006] (3)结构模型跟踪法：根据先验知识建立目标的几何模型，跟踪时先根据上一帧结 果预测当前帧目标姿态，将处于预测姿态的目标模型投影到图像平面，计算匹配误差；然后 通过优化预测姿态、最小化目标函数获得当前跟踪结果。常用的结构模型为3D线框模型。
[0007] (4)剪影跟踪法：首要任务是获取目标的轮廓，精确得到目标占据的区域；如基于 主动轮廓模型的跟踪法。
[0008] 此外，车辆跟踪算法根据算法驱动方式不同，还可分为自下而上的数据驱动方法， 典型代表为Mean Shift算法，和自上而下的模型驱动方法，典型代表为粒子滤波算法。
[0009] 上述的车辆跟踪算法存在以下缺点：
[0010] (1)目标的模型建立、轮廓获取难度大，且目标函数定义或优化较复杂；
[0011] (2)搜索策略计算量较大，很难达到实时跟踪要求
[0012] (3)目标特征不能及时得到更新，姿态变化易导致跟丢；
[0013] (4)跟踪过程与检测结果相对独立，无交互；
[0014] (5)对光线变化、相邻目标遮挡等因素敏感；
[0015] (6)对由于遮挡、错误消去等原因造成的短暂丢失的目标无法继续跟踪。

【发明内容】

[0016] 本发明的目的是克服了上述现有技术的缺点，提供了一种能够实现准确实时跟踪 多个车辆、添加抗光照因子弥补光流场在光照突变条件下鲁棒性不高的不足、适用于对算 法实时性要求较高的电子警察产品或电子卡口产品中、具有更广泛应用范围的基于变分光 流实现车辆跟踪的方法。
[0017] 为了实现上述目的，本发明的基于变分光流实现车辆跟踪的方法具有如下构成：
[0018] 该基于变分光流实现车辆跟踪的方法，其主要特点是，所述的方法包括以下步 骤：
[0019] (1)从视频电子警察中获取车辆影像并确定目标跟踪区域；
[0020] (2)在所述的目标跟踪区域中提取特征点；
[0021] (3)采用变分光流方程计算各个特征点的运动矢量；
[0022] (4)根据所述的各个特征点的运动矢量确定车辆的位移。
[0023] 较佳地，所述的从视频电子警察中获取车辆影像并确定目标跟踪区域，包括以下 步骤：
[0024] (11)从视频电子警察中获取拍摄的图片并检测到车辆位置；
[0025] (12)采用图片中能覆盖车尾的图像区域作为目标跟踪区域。
[0026] 较佳地，所述的在所述的目标跟踪区域中提取特征点，具体为：
[0027] 在所述的目标跟踪区域中采用Harris角点检测算法提取特征点。
[0028] 更佳地，所述的在所述的目标跟踪区域中采用Harris角点检测算法提取特征点， 包括以下步骤：
[0029] (21)对所述的目标跟踪区域中的各个像素点按照以下公式计算相关矩阵：
【主权项】
1. 一种基于变分光流实现车辆跟踪的方法，其特征在于，所述的方法包括以下步骤： (1) 从视频电子警察中获取车辆影像并确定目标跟踪区域； (2) 在所述的目标跟踪区域中提取特征点； (3) 采用变分光流方程计算各个特征点的运动矢量； (4) 根据所述的各个特征点的运动矢量确定车辆的位移。
2. 根据权利要求1所述的基于变分光流实现车辆跟踪的方法，其特征在于，所述的从 视频电子警察中获取车辆影像并确定目标跟踪区域，包括以下步骤： (11) 从视频电子警察中获取拍摄的图片并检测到车辆位置； (12) 采用图片中能覆盖车尾的图像区域作为目标跟踪区域。
3. 根据权利要求1所述的基于变分光流实现车辆跟踪的方法，其特征在于，所述的在 所述的目标跟踪区域中提取特征点，具体为： 在所述的目标跟踪区域中采用Harris角点检测算法提取特征点。
4. 根据权利要求3所述的基于变分光流实现车辆跟踪的方法，其特征在于，所述的在 所述的目标跟踪区域中采用Harris角点检测算法提取特征点，包括以下步骤： (21) 对所述的目标跟踪区域中的各个像素点按照以下公式计算相关矩阵： A = \ν(.\\ ν)??;. B = u(.v, v)?/,2. C = D = u'(.v,v)? (/,/,). ' ， 其中，A、B、C、D为各个像素点的相关矩阵，w (x, y)为各个像素点对应的高斯系数，Ix、 Iy为各个像素点对应的偏导； (22) 按照如下公式计算各个像素点的Harris角点的角点响应函数R : R = (AB-CD) 2-k (A+B)2 ; 其中，k 取（0· 04 ~0· 06); (23) 在所述的目标跟踪区域内提取各个像素点的角点响应函数R的极大值，将极大值 大于系统预设阈值的像素点确定为特征点。
5. 根据权利要求1所述的基于变分光流实现车辆跟踪的方法，其特征在于，所述的采 用变分光流方程计算各个特征点的运动矢量，包括以下步骤： (31) 对各个特征点计算数据项和平滑项，并按照如下公式建立基于变光分流算法的能 量泛函方程： E (u，V) = Edata (u，V) + a E smooth (u, v); 其中，Edata(u，v)为数据项，数据项为各个特征点的灰度值守恒约束，Es_th(u，v)为平滑 项，平滑项为各个特征点的光流平滑约束，u和V分别为各个特征点在x、y方向上的速度分 量，α为影响数据项与平滑项所占比重的权值。 (32) 对所述的能量泛函方程进行求解得到各个特征点的运动矢量（u，V)。
6. 根据权利要求5所述的基于变分光流实现车辆跟踪的方法，其特征在于，所述的对 各个特征点计算数据项和平滑项，包括以下步骤： (311)按照如下公式计算各个特征点的数据项Edata(u，V):
(312)按照如下公式计算各个特征点的平滑项Es_th(u，v):
7. 根据权利要求6所述的基于变分光流实现车辆跟踪的方法，其特征在于，所述的对 所述的能量泛函方程进行求解得到各个特征点的运动矢量，包括以下步骤： (321) 初始化各个特征点的运动矢量的u和V,使得（Uci, Vci) = (0, 0); (322) 按照如下公式建立稳态方程并求出u和V的变化值Λ u和Λ V :
其中，Ix为点（x+u^y+Vi)周围四邻域的图像X偏导的平均值，Iy为点（x+i^y+Vi)周 围四邻域的图像y偏导的平均值，It为点（x+Ui，y+Vi)周围四邻域的此帧与下一帧灰度值的 差值； (323) 按照如下公式确定各个特征点下一次迭代的运动矢量值： ui+1 = Ui+Au vi+1 = Vi+ Δ v ; (324) 判断u和ν的变化值Au和Λν是否趋向于0,如果是，则继续步骤（325)，否则 继续步骤（322); (325) 将当前ui+1和vi+1的值确定为对应特征点的运动矢量。
8. 根据权利要求5所述的基于变分光流实现车辆跟踪的方法，其特征在于，所述的根 据所述的各个特征点的运动矢量确定车辆的位移，包括以下步骤： (41) 根据所述的各个特征点的运动矢量计算（u，V)的平均值（); (42) 对运动矢量（Ui，Vi)在（M ±10)范围内的光流计算平均值并将该平均值确 定为车辆的位移。
【专利摘要】本发明涉及一种基于变分光流实现车辆跟踪的方法，其中包括从视频电子警察中获取车辆影像并确定目标跟踪区域，即从视频电子警察中获取拍摄的图片并检测到车辆位置后采用图片中能覆盖车尾的图像区域作为目标跟踪区域；在所述的目标跟踪区域中提取特征点；采用变分光流方程计算各个特征点的运动矢量；根据所述的各个特征点的运动矢量确定车辆的位移。采用该种基于变分光流实现车辆跟踪的方法，可应用于对算法实时性要求较高的电子警察产品或电子卡口产品中，可以实现准确、实时地跟踪多个车辆，并添加抗光照因子弥补光流场在光照突变条件下鲁棒性不高的不足，抗光照干扰，避免环境光巨变引起的跟丢情况，方法应用简便，具有更广泛的应用范围。
【IPC分类】G06T7-20
【公开号】CN104658006
【申请号】CN201310598033
【发明人】蒋慧涛
【申请人】上海宝康电子控制工程有限公司
【公开日】2015年5月27日
【申请日】2013年11月22日