发明名称 一种基于概率模型的运动摄像机自动目标跟踪方法
摘要 本发明提供一种基于概率模型的运动摄像机自动目标跟踪方法,具体实现步骤为:首先,初始化启动跟踪的参数和实时运动目标跟踪的参数;其次,利用速度差异函数定位初始运动区域,并采用改进的最优聚类函数分割运动目标;然后,对当前帧检测到的运动目标和背景搜索窗建立颜色概率直方图模型,利用后验概率函数、颜色直方图相似度函数调整目标的大小和位置;最后,使用连续自适应均值漂移法寻找下一帧运动目标,实现在动态背景下对运动目标进行实时跟踪。本发明提高了运动目标跟踪的可靠性,可有效解决运动摄像机下初始目标选取、多目标跟踪、目标和背景存在相同颜色以及目标自身颜色变化问题。
申请公布号 CN102332166B 申请公布日期 2014.04.16
申请号 CN201110287392.0 申请日期 2011.09.26
申请人 北京航空航天大学 发明人 赵龙;伊威;肖军波
分类号 G06T7/20(2006.01)I 主分类号 G06T7/20(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 李新华;贾玉忠
主权项 1.一种基于概率模型的运动摄像机自动目标跟踪方法,在通过速度差异和最优聚类获取运动目标的基础上,通过建立改进的颜色概率直方图模型,增加目标颜色分布的可信度,并根据运动目标的位置、大小和颜色采用均值偏移法,对实时运动目标进行可靠跟踪,其特征在于包括以下步骤:步骤(1)初始化启动跟踪的参数和实时运动目标跟踪的参数;步骤(2)利用速度差异函数定位初始运动区域,并利用改进的最优聚类函数分割目标;步骤(3)对当前帧检测到的运动目标和背景搜索窗建立颜色概率直方图模型,利用后验概率函数、颜色直方图相似度函数调整目标的大小和位置;步骤(4)使用连续自适应均值漂移法寻找下一帧运动目标,实现实时跟踪;其中,所述步骤(2)中利用速度差异函数定位初始运动区域的具体步骤为:步骤a1、计算训练帧中最小的角点数为<img file="FDA0000425971110000011.GIF" wi="1253" he="208" />式中,FC<sub>max</sub>代表最大的角点数;TC代表所训练的帧的长度;C<sub>t</sub>(i,j)代表t时刻坐标为(i,j)位置的角点;步骤a2、计算角点全局速度为G<sub>speed</sub>=MostFrequence{S<sub>Speed</sub>(i1,j1)…S<sub>Speed</sub>(in,jn)…S<sub>Speed</sub>(iN,jN)}    (2)式中,1≤n≤N,(in,jn)∈D,D代表整幅图像,S<sub>Speed</sub>(in,jn)对应坐标为(in,jn)位置角点的速度,其中(i1,j1)≠(in,jn);步骤a3、计算局部运动区域为|S<sup>t</sup><sub>Speed</sub>(in,jn)-S<sup>t-1</sup><sub>Speed</sub>(in,jn)|<Th0,t∈(1,TC-1)    (3)|L<sup>t</sup><sub>speed</sub>(in,jn)-L<sup>t-1</sup><sub>speed</sub>(in,jn)|<Th0,t∈(1,TC-1)    (4)|L<sub>speed</sub>(in,jn)-G<sub>speed</sub>|>Th1,1≤n≤N    (5)式中,S<sup>t</sup><sub>Speed</sub>(in,jn)代表t时刻第n个角点的速度;G<sub>speed</sub>代表一帧中所有角点的全局运动速度;L<sub>speed</sub>(in,jn)代表第n个角点像素的局部运动速度;TC代表所训练的帧的长度;N代表最小的角点数;所述步骤(2)中利用改进的最优聚类函数分割目标的步骤为:步骤b1、类间离散度函数为<maths num="0001"><![CDATA[<math><mrow><mrow><mi>L</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>&xi;</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mo>|</mo><msub><mi>m</mi><mi>&xi;</mi></msub><mo>-</mo><mover><mi>m</mi><mo>&OverBar;</mo></mover><mo>|</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中m<sub>ξ</sub>为任一簇所含样本的均值;<img file="FDA0000425971110000022.GIF" wi="57" he="61" />为所有样本的均值,k表示总的聚类簇数;步骤b2、类内离散度函数为<maths num="0002"><![CDATA[<math><mrow><mi>W</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>&xi;</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><munder><mi>&Sigma;</mi><mrow><mi>p</mi><mo>&Element;</mo><msub><mi>C</mi><mi>&xi;</mi></msub></mrow></munder><mo>|</mo><mi>p</mi><mo>-</mo><msub><mi>m</mi><mi>&xi;</mi></msub><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,m<sub>ξ</sub>为任一簇所含样本的均值;p为对应簇中的任一样本;k表示总的聚类簇数;步骤b3、两离散度距离函数为<maths num="0003"><![CDATA[<math><mrow><mi>F</mi><mrow><mo>(</mo><mi>S</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>L</mi><mo>/</mo><mi>W</mi><mo>=</mo><mfenced open='' close=''><mtable><mtr><mtd><mn>1</mn><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>&xi;</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><munder><mi>&Sigma;</mi><mrow><mi>p</mi><mo>&Element;</mo><msub><mi>C</mi><mi>&xi;</mi></msub></mrow></munder><mo>|</mo><mi>p</mi><mo>-</mo><msub><mi>m</mi><mi>&xi;</mi></msub><mo>|</mo><mo>,</mo><mi>if</mi><mrow><mo>(</mo><mi>L</mi><mo>=</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><munderover><mi>&Sigma;</mi><mrow><mi>&zeta;</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mo>|</mo><msub><mi>m</mi><mi>&xi;</mi></msub><mo>-</mo><mover><mi>m</mi><mo>&OverBar;</mo></mover><mo>|</mo><mo>/</mo><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>&xi;</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><munder><mi>&Sigma;</mi><mrow><mi>p</mi><mo>&Element;</mo><msub><mi>C</mi><mi>&xi;</mi></msub></mrow></munder><mo>|</mo><mi>p</mi><mo>-</mo><msub><mi>m</mi><mi>&xi;</mi></msub><mo>)</mo></mrow><mo>,</mo><mi>if</mi><mrow><mo>(</mo><mi>L</mi><mo>&NotEqual;</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,S代表整个索引空间。
地址 100191 北京市海淀区学院路37号