发明名称 一种基于高斯改进型粒子群粒子滤波的移动目标追踪方法
摘要 本发明属于信号处理领域,具体涉及用于快速跟踪外观变化的移动目标的基于高斯改进型粒子群粒子滤波的移动目标追踪方法。本发明包括:首先在视频第一帧中选取目标,初始化特征基;跟踪过程采集下一帧图像,获得跟踪目标的特征基;通过高斯改进型粒子群粒子滤波器滤波,存储对应于概率最大的粒子的图像窗口;当达到采集的图像达到要求的数目,进行观测值特征基的更新。本发明主要针对移动机器人快速跟踪外观变化的移动目标,提出一种能够提高复杂环境下快速跟踪移动目标的方法。本发明利用低维特征基空间表示描述跟踪的物体,采用高斯改进粒子群粒子滤波器进行跟踪目标位置的跟踪,从而提高了跟踪移动目标的速度。
申请公布号 CN104794736A 申请公布日期 2015.07.22
申请号 CN201510166890.8 申请日期 2015.04.09
申请人 哈尔滨工程大学 发明人 姚建均;余瀚;陈硕;肖蕊;王涛;牛庆涛
分类号 G06T7/20(2006.01)I 主分类号 G06T7/20(2006.01)I
代理机构 代理人
主权项 一种基于高斯改进型粒子群粒子滤波的移动目标追踪方法,其特征在于,包括如下步骤:(1)首先在视频第一帧中选取目标,初始化特征基;(2)跟踪过程采集下一帧图像,获得跟踪目标的特征基;(3)通过高斯改进型粒子群粒子滤波器滤波,存储对应于概率最大的粒子的图像窗口;(3.1)初始化粒子群飞行速度,设置学习因子,初始化粒子滤波器,首先由p(x<sub>0</sub>)抽取样本<img file="FDA0000696122430000011.GIF" wi="181" he="78" />其中N为采样粒子数;(3.2)使用状态转移方程f<sub>n</sub>(·)将样本在时间上向前传播,即<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mover><mi>x</mi><mo>~</mo></mover><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><msub><mi>f</mi><mi>n</mi></msub><mrow><mo>(</mo><msubsup><mover><mi>x</mi><mo>~</mo></mover><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>v</mi><mi>n</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000696122430000012.GIF" wi="403" he="85" /></maths>其中v<sub>n</sub>为噪声项,x<sub>n</sub>表示在t<sub>n</sub>时刻的状态变量;(3.3)计算<img file="FDA0000696122430000013.GIF" wi="60" he="66" />和<img file="FDA0000696122430000014.GIF" wi="99" he="84" />即<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>&mu;</mi><mo>&OverBar;</mo></mover><mi>n</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mover><mi>x</mi><mo>~</mo></mover><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup></mrow>]]></math><img file="FDA0000696122430000015.GIF" wi="303" he="149" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>&Sigma;</mi><mo>&OverBar;</mo></mover><mi>n</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mover><mi>x</mi><mo>~</mo></mover><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><msup><msubsup><mover><mi>x</mi><mo>~</mo></mover><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mi>T</mi></msup><mo>-</mo><msub><mover><mi>&mu;</mi><mo>&OverBar;</mo></mover><mi>n</mi></msub><msubsup><mover><mi>&mu;</mi><mo>&OverBar;</mo></mover><mi>n</mi><mi>T</mi></msubsup></mrow>]]></math><img file="FDA0000696122430000016.GIF" wi="573" he="156" /></maths>(3.4)计算重要性权值<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mover><mi>&omega;</mi><mo>~</mo></mover><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mi>p</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>n</mi></msub><mo>|</mo><msub><mi>x</mi><mi>n</mi></msub><mo>=</mo><msubsup><mi>x</mi><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000696122430000017.GIF" wi="445" he="85" /></maths>(3.5)将样本粒子的权值进行归一化,得<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>w</mi><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mfrac><msubsup><mover><mi>w</mi><mo>~</mo></mover><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mover><mi>w</mi><mo>~</mo></mover><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup></mrow></mfrac></mrow>]]></math><img file="FDA0000696122430000018.GIF" wi="318" he="228" /></maths>其中<img file="FDA0000696122430000019.GIF" wi="84" he="82" />为归一化前的权值,<img file="FDA00006961224300000110.GIF" wi="82" he="81" />为归一化后的权值,(3.6)求取对应最大的权值的粒子作为粒子群初始全局最优的粒子;(3.7)将粒子的权重作为粒子群的适应值;(3.8)根据粒子群算法利用下式来更新每个粒子的速度和位置:<img file="FDA00006961224300000111.GIF" wi="1415" he="95" /><img file="FDA00006961224300000112.GIF" wi="1103" he="99" /><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>x</mi><mi>ij</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>n</mi><mi>ij</mi><mi>n</mi></msubsup><mo>+</mo><msubsup><mi>v</mi><mi>ij</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mrow>]]></math><img file="FDA00006961224300000113.GIF" wi="311" he="87" /></maths>其中C<sub>1</sub>C<sub>2</sub>为常数,rand()为随机数;<img file="FDA00006961224300000114.GIF" wi="57" he="80" />为第i粒子在第j维的飞行速度,即粒子在一次迭代中移动的距离;<img file="FDA00006961224300000115.GIF" wi="61" he="78" />为第i个粒子在第j维位置;p<sub>ij</sub>为当前粒子最优位置,p<sub>gj</sub>为全局最优位置;(3.9)再一次计算各个样本的权值,<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mover><mi>&omega;</mi><mo>~</mo></mover><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mi>p</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>n</mi></msub><mo>|</mo><msub><mi>x</mi><mi>n</mi></msub><mo>=</mo><msubsup><mi>x</mi><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000696122430000021.GIF" wi="467" he="79" /></maths>(3.10)将各样本粒子的权值归一化,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>w</mi><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mfrac><msubsup><mover><mi>w</mi><mo>~</mo></mover><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mover><mi>w</mi><mo>~</mo></mover><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000696122430000022.GIF" wi="308" he="219" /></maths>(3.11)求取最大权值对应的粒子,如果粒子的权值小于设置的阈值,并且在设置的次数内跳转到步骤(3.6);(3.12)输出对应最大权值的粒子;(4)当达到采集的图像达到要求的数目,进行观测值特征基的更新。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室