发明名称 基于回归分析的粒子滤波目标跟踪方法
摘要 本发明公开了一种基于回归分析的粒子滤波目标跟踪方法,主要解决现有技术在粒子滤波中,由于粒子匮乏导致跟踪失败和计算效率低的问题。其实现步骤是:1.在粒子滤波框架下,根据前一时刻回归分析的结果预测目标状态得到候选目标,提取候选目标特征;2.通过计算候选目标特征与模板之间的距离,得到对应粒子的权值;3.根据粒子的权值对粒子进行重采样,然后估计目标状态;4.计算回归分析的预测状态和目标状态的偏差,进一步更新粒子集,实现对目标的可靠跟踪。本发明对目标遮挡问题具有较好的鲁棒性,同时计算效率较高,可用于智能视频监控、机器人导航、道路交通管制等系统。
申请公布号 CN106408590A 申请公布日期 2017.02.15
申请号 CN201610917508.7 申请日期 2016.10.21
申请人 西安电子科技大学 发明人 苏镇镇;张文博;姬红兵;徐步;王磊
分类号 G06T7/20(2017.01)I;G06T7/292(2017.01)I 主分类号 G06T7/20(2017.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于回归分析的粒子滤波目标跟踪方法,包括:(1)初始化步骤:(1a)读入k‑1时刻的图像I<sub>k‑1</sub>,根据k‑1时刻的目标状态X<sub>k‑1</sub>,初始粒子数为N<sub>0</sub>,初始化k‑1时刻的粒子集<img file="FDA0001135659830000011.GIF" wi="195" he="87" />其中,<img file="FDA0001135659830000012.GIF" wi="85" he="63" />表示k‑1时刻第i个粒子的状态区间,i为粒子序号,取值为1,2,…,N,N表示粒子总数,k表示时刻;(1b)初始化目标跟踪窗:B<sub>k‑1</sub>=(r<sub>k‑1</sub>,c<sub>k‑1</sub>)<sup>T</sup>,其中r<sub>k‑1</sub>和c<sub>k‑1</sub>分别表示k‑1时刻目标跟踪窗的长度和宽度值,T表示向量转置运算;(1c)根据目标初始状态与目标跟踪窗B<sub>k‑1</sub>,初始化目标的特征模板M;(2)目标状态预测步骤:(2a)读入k时刻的图像I<sub>k</sub>,根据前n个时刻的速度,通过最小二乘法预测k‑1时刻的速度<img file="FDA0001135659830000013.GIF" wi="99" he="62" />并按照该速度<img file="FDA0001135659830000014.GIF" wi="70" he="63" />对k‑1时刻粒子集<img file="FDA0001135659830000015.GIF" wi="163" he="95" />进行预测,得到k时刻预测粒子集<img file="FDA0001135659830000016.GIF" wi="155" he="95" />其中<img file="FDA0001135659830000017.GIF" wi="43" he="62" />为k时刻第i个粒子的状态预测值;(2b)根据k时刻预测粒子集<img file="FDA0001135659830000018.GIF" wi="126" he="91" />和目标跟踪窗B<sub>k‑1</sub>,确定k时刻候选目标集<img file="FDA0001135659830000019.GIF" wi="158" he="95" />并求解候选目标集<img file="FDA00011356598300000110.GIF" wi="132" he="95" />中与候选目标对应的特征集<img file="FDA00011356598300000111.GIF" wi="150" he="83" />其中<img file="FDA00011356598300000112.GIF" wi="51" he="70" />为k时刻第i个候选目标,表示以<img file="FDA00011356598300000113.GIF" wi="46" he="62" />为中心、B<sub>k‑1</sub>为长宽所界定出的矩形区域;(3)计算权值集步骤:(3a)求取特征集<img file="FDA00011356598300000114.GIF" wi="120" he="83" />与特征模板M之间的距离集<img file="FDA00011356598300000115.GIF" wi="147" he="87" />其中d<sub>i</sub>表示第i个候选目标特征V<sub>i</sub>与特征模板M之间的距离;(3b)根据距离集<img file="FDA00011356598300000116.GIF" wi="121" he="83" />计算候选目标的权值集<img file="FDA00011356598300000117.GIF" wi="158" he="95" />其中ω<sub>i</sub>表示第i个候选目标的权值;(4)重采样步骤:(4a)对权值集<img file="FDA00011356598300000118.GIF" wi="139" he="91" />归一化得到<img file="FDA00011356598300000119.GIF" wi="149" he="95" />(4b)根据归一化后的权值集<img file="FDA00011356598300000120.GIF" wi="155" he="91" />利用系统重采样算法对k‑1时刻粒子集<img file="FDA00011356598300000121.GIF" wi="163" he="87" />进行重采样,得到k时刻的重采样粒子集<img file="FDA00011356598300000122.GIF" wi="160" he="95" />其中<img file="FDA00011356598300000123.GIF" wi="59" he="63" />表示k时刻第i个重采样粒子的状态区间;(5)估计目标状态步骤:(5a)根据k时刻的更新粒子集<img file="FDA0001135659830000021.GIF" wi="163" he="95" />估计k时刻的目标状态X<sub>k</sub>;(5b)根据k时刻的目标状态X<sub>k</sub>和目标跟踪窗B<sub>k‑1</sub>,确定k时刻的目标估计范围T:<img file="FDA0001135659830000022.GIF" wi="814" he="143" />其中,x和y分别表示的横、纵坐标,m<sub>k</sub>和n<sub>k</sub>分别表示k时刻目标状态X<sub>k</sub>的横坐标和纵坐标,r<sub>k</sub>和c<sub>k</sub>分别表示k时刻目标跟踪窗的长度和宽度值;(6)判断是否迭代:检查下一时刻的信息是否到达,若没有信息到达,则目标跟踪过程结束,输出估计目标状态;否则,执行步骤(7);(7)回归分析及预测步骤:(7a)在k时刻的目标状态X<sub>k</sub>中,取出目标的位置信息P<sub>k</sub>=[m<sub>k</sub> n<sub>k</sub>]<sup>T</sup>;(7b)根据k‑1时刻的预测速度<img file="FDA0001135659830000023.GIF" wi="95" he="62" />进一步预测目标在k时刻的位置<img file="FDA0001135659830000024.GIF" wi="67" he="73" />(7c)计算k时刻目标状态与预测状态之间的偏差:<img file="FDA0001135659830000025.GIF" wi="251" he="70" />(7d)根据偏差E选取粒子数,同时更新粒子集<img file="FDA0001135659830000026.GIF" wi="163" he="95" />并令k=k+1,返回步骤(2)进行迭代。
地址 710071 陕西省西安市雁塔区太白南路2号