发明名称 自适应运动相关协作多目标跟踪方法
摘要 本发明公开了一种自适应运动相关协作多目标跟踪方法,利用多个目标之间的时间和空间信息,建立目标间支持度模型,在目标发生遮挡时,通过所述支持度模型估计出该目标的状态。本发明方法利用多个目标之间的时间和空间信息,通过建立目标间支持度模型,即目标间加速度的二维高斯加权模型,在目标发生遮挡时,通过这种支持度模型正确估计出该目标的状态。本发明方法可以有效的处理目标遮挡,在较为复杂的环境下,快速稳定的跟踪多个目标。
申请公布号 CN102214301A 申请公布日期 2011.10.12
申请号 CN201110209402.9 申请日期 2011.07.26
申请人 西南交通大学 发明人 权伟;陈锦雄;余南阳
分类号 G06K9/00(2006.01)I;G06T7/20(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 成都信博专利代理有限责任公司 51200 代理人 张澎
主权项 1.自适应运动相关协作多目标跟踪方法,利用多个目标之间的时间和空间信息,建立目标间支持度模型,在目标发生遮挡时,通过所述支持度模型估计出该目标的状态;包含如下步骤:(1)实时图像输入将实时视频图像作为要进行跟踪的输入图像;(2)确定跟踪目标及其区域预处理通过运动对象检测自动提取运动对象,将需要进行跟踪的对象作为跟踪目标,并加入到跟踪目标列表中;对每一个目标所在的区域及其邻近搜索区域进行局部标准差图像增强预处理;(3)计算每个目标的搜索区域相关度对每个目标:在目标邻近搜索区域计算相关度平面,相关度计算公式为:<maths num="0001"><![CDATA[<math><mrow><mi>c</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>s</mi><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mi>i</mi><mo>,</mo><mi>y</mi><mo>+</mo><mi>j</mi><mo>)</mo></mrow><mi>T</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>s</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mi>i</mi><mo>,</mo><mi>y</mi><mo>+</mo><mi>j</mi><mo>)</mo></mrow></msqrt><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>T</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></msqrt></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>目标相关度平面的最大相关度值记为c<sub>max</sub>,对应坐标记为(x<sub>m</sub>,y<sub>m</sub>);设相关度阈值为δ<sub>c</sub>,如果c<sub>max</sub>>δ<sub>c</sub>,则该目标标记为未遮挡;否则,标记为遮挡;(4)目标状态更新对每个标记为未遮挡的目标:根据(3)得到的目标当前坐标(x<sub>m</sub>,y<sub>m</sub>),更新目标状态,通过线性卡尔曼滤波方法进行状态更新,并作为本次跟踪过程中被遮挡目标的协作支持目标;(5)目标状态估计设t时刻,目标的状态向量为<img file="FDA0000078422940000021.GIF" wi="630" he="61" />其中,<img file="FDA0000078422940000022.GIF" wi="59" he="54" /><img file="FDA0000078422940000023.GIF" wi="391" he="59" />和<img file="FDA0000078422940000024.GIF" wi="60" he="59" />分别表示目标的横纵坐标,x方向和y方向的速度与加速度;则对每个标记为遮挡的目标,其支持度状态估计加权模型为:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>X</mi><mi>i</mi><mi>t</mi></msubsup><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow><mi>n</mi></munderover><msubsup><mi>W</mi><mi>ij</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><mi>C</mi><msubsup><mi>&Delta;X</mi><mi>ij</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>X</mi><mi>j</mi><mi>t</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>W</mi><mi>i</mi></msub><mi>C</mi><msubsup><mi>X</mi><mi>i</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&NotEqual;</mo></mrow><mi>n</mi></munderover><msubsup><mi>W</mi><mi>ij</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msub><mi>W</mi><mi>i</mi></msub></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>其中,n为未遮挡目标的个数,W<sub>ij</sub>为目标间支持度,即状态估计加权值,C为如下的矩阵:<maths num="0003"><![CDATA[<math><mrow><mi>C</mi><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn><mo>/</mo><mn>2</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn><mo>/</mo><mn>2</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>.</mo></mrow></math>]]></maths>如果全部目标均被遮挡,则直接通过目标前一时刻的状态,估计当前目标状态;(6)目标支持度更新目标之间的支持度计算公式为:<maths num="0004"><![CDATA[<math><mrow><mi>W</mi><mrow><mo>(</mo><msub><mi>&Delta;a</mi><mi>x</mi></msub><mo>,</mo><msub><mi>&Delta;a</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mrow><mn>2</mn><mi>&pi;&sigma;</mi></mrow><mn>1</mn></msub><msub><mi>&sigma;</mi><mn>2</mn></msub><msqrt><mn>1</mn><mo>-</mo><msup><mi>&rho;</mi><mn>2</mn></msup></msqrt></mrow></mfrac><mo>&times;</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mi>exp</mi><mo>{</mo><mo>-</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&rho;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><mo>[</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>&Delta;a</mi><mi>x</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><msup><msub><mi>&sigma;</mi><mn>1</mn></msub><mn>2</mn></msup></mfrac><mo>-</mo><mn>2</mn><mi>&rho;</mi><mfrac><mrow><mrow><mo>(</mo><msub><mi>&Delta;a</mi><mi>x</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>&Delta;a</mi><mi>y</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>&sigma;</mi><mn>1</mn></msub><msub><mi>&sigma;</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>&Delta;a</mi><mi>y</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><msup><msub><mi>&sigma;</mi><mn>2</mn></msub><mn>2</mn></msup></mfrac><mo>]</mo><mo>}</mo><mo>,</mo></mrow></math>]]></maths>其中,Δa<sub>x</sub>和Δa<sub>y</sub>分别为两个目标在x方向和y方向的加速度差;目标间支持度更新后,本次多目标跟踪任务完成,跳转到(2)。
地址 610031 四川省成都市二环路北一段111号西南交通大学科技处