发明名称 基于概率主成分分析和压缩感知的在线目标跟踪方法
摘要 本发明涉及一种基于概率主成分分析和压缩感知的在线目标跟踪方法,将压缩感知理论与PPCA相结合,得到了一个非常简洁的目标特征子空间表示模型,并且使用增量PCA算法在线更新该目标子空间表示模型,使得跟踪算法能够达到很好的实时性。同时我们将PPCA理论运用到计算候选目标与目标子空间表示模型的视觉相似性上,得到了两个距离DFFS和DIFS,而不仅仅是PCA中的重构误差,使得视觉相似性的计算更为准确,这样保证了跟踪算法的准确性。
申请公布号 CN103310463B 申请公布日期 2016.01.20
申请号 CN201310241336.2 申请日期 2013.06.18
申请人 西北工业大学 发明人 李映;宋旭;冉晨
分类号 G06T7/20(2006.01)I 主分类号 G06T7/20(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种基于概率主成分分析和压缩感知的在线目标跟踪方法,其特征在于步骤如下:步骤1:在第一帧中标记目标x<sub>1</sub>,初始化N个粒子及其权值<img file="FDA0000808465040000011.GIF" wi="319" he="147" />所述x<sub>1</sub>为第一帧中目标图像块的仿射变换参数;步骤2:使用经典粒子滤波算法对前T帧图像跟踪目标,得到初始目标样本集合A=[y<sub>1</sub>,y<sub>2</sub>,...,y<sub>T</sub>],A的每一列是这T帧图像中目标图像块的特征表示向量,y<sub>i</sub>是高维度的图像特征向量在测量矩阵P上的投影,是压缩了的低维度特征向量;步骤3:对A‑I<sub>A</sub>进行SVD分解,得到目标特征子空间的表示模型;步骤4:对当前帧是第t帧的前一帧中粒子<img file="FDA0000808465040000012.GIF" wi="190" he="85" />按其权值<img file="FDA0000808465040000013.GIF" wi="158" he="84" />的正比关系进行重采样,再运用高斯运动模型,得到当前帧中粒子<img file="FDA0000808465040000014.GIF" wi="154" he="84" />步骤5:求出当前帧中粒子<img file="FDA0000808465040000015.GIF" wi="124" he="81" />对应图像块的特征表示向量<img file="FDA0000808465040000016.GIF" wi="160" he="85" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mrow><mo>{</mo><msubsup><mi>y</mi><mi>t</mi><mi>i</mi></msubsup><mo>}</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><mo>=</mo><mi>P</mi><msubsup><mrow><mo>{</mo><msubsup><mi>x</mi><mi>t</mi><mi>i</mi></msubsup><mo>}</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup></mrow>]]></math><img file="FDA0000808465040000017.GIF" wi="343" he="80" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mi>w</mi><mi>i</mi><mi>t</mi><mi>h</mi></mrow></mtd><mtd><mrow><mi>p</mi><mi>r</mi><mi>o</mi><mi>b</mi><mi>a</mi><mi>b</mi><mi>i</mi><mi>l</mi><mi>l</mi><mi>i</mi><mi>t</mi><mi>y</mi></mrow></mtd><mtd><mfrac><mn>1</mn><mrow><mn>2</mn><mi>s</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>w</mi><mi>i</mi><mi>t</mi><mi>h</mi></mrow></mtd><mtd><mrow><mi>p</mi><mi>r</mi><mi>o</mi><mi>b</mi><mi>a</mi><mi>b</mi><mi>i</mi><mi>l</mi><mi>l</mi><mi>i</mi><mi>t</mi><mi>y</mi></mrow></mtd><mtd><mrow><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><mi>s</mi></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mrow><mi>w</mi><mi>i</mi><mi>t</mi><mi>h</mi></mrow></mtd><mtd><mrow><mi>p</mi><mi>r</mi><mi>o</mi><mi>b</mi><mi>a</mi><mi>b</mi><mi>i</mi><mi>l</mi><mi>l</mi><mi>i</mi><mi>t</mi><mi>y</mi></mrow></mtd><mtd><mfrac><mn>1</mn><mrow><mn>2</mn><mi>s</mi></mrow></mfrac></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000808465040000018.GIF" wi="821" he="407" /></maths>s=2or s=3and i=1,2,...,L.j=1,2,...,H其中,取L=50,H是高维特征向量的维数;步骤6:计算出当前帧中第i个粒子<img file="FDA0000808465040000019.GIF" wi="42" he="77" />对应图像块的特征表示向量<img file="FDA00008084650400000114.GIF" wi="55" he="73" />与目标特征子空间的视觉相似性的概率表述<img file="FDA00008084650400000110.GIF" wi="220" he="79" />将该值作为粒子<img file="FDA00008084650400000111.GIF" wi="48" he="76" />新的权值<img file="FDA00008084650400000112.GIF" wi="321" he="76" />然后运用最大后验概率准则MAP准则,得到当前帧中具有最大权值的粒子作为该帧目标的状态估计值<img file="FDA00008084650400000113.GIF" wi="85" he="82" />以此为对当前帧的跟踪结果;若当前帧是最后一帧,则结束,否则继续执行;步骤7:当已经跟踪的帧数等于m帧时,执行步骤8,否则转到步骤4;其中:m为更新频率,2&lt;m&lt;10;步骤8:以m个跟踪结果得到一个增量矩阵B=[y<sub>n+1</sub>,y<sub>n+2</sub>,...,y<sub>n+m</sub>],其中的每一列是这m帧中跟踪到的目标图像块对应的特征表示向量;步骤9:使用增量PCA算法更新得到C=[C,B],构成目标特征空间的子空间表示模型;然后重复转到步骤4;其中第一次更新时令C=[A,B]:增量PCA算法流程如下:输入:当前目标子空间模型U<sub>d</sub>,Σ<sub>d</sub>,u,增量矩阵B,B的均值I<sub>B</sub>,遗忘因子f;输出:新的目标子空间模型<img file="FDA0000808465040000021.GIF" wi="231" he="78" />a.更新均值:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>I</mi><mi>C</mi></msub><mo>=</mo><mfrac><mrow><mi>f</mi><mo>&CenterDot;</mo><mi>n</mi></mrow><mrow><mi>f</mi><mo>&CenterDot;</mo><mi>n</mi><mo>+</mo><mi>m</mi></mrow></mfrac><mi>u</mi><mo>+</mo><mfrac><mi>m</mi><mrow><mi>f</mi><mo>&CenterDot;</mo><mi>n</mi><mo>+</mo><mi>m</mi></mrow></mfrac><msub><mi>I</mi><mi>B</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000808465040000022.GIF" wi="619" he="146" /></maths>I<sub>C</sub>是C的均值;b.令<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><mi>B</mi><mo>&OverBar;</mo></mover><mo>=</mo><mo>&lsqb;</mo><msub><mi>y</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>I</mi><mi>B</mi></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>y</mi><mrow><mi>m</mi><mo>+</mo><mi>n</mi></mrow></msub><mo>-</mo><msub><mi>I</mi><mi>B</mi></msub><mo>,</mo><msqrt><mfrac><mrow><mi>n</mi><mo>&CenterDot;</mo><mi>m</mi></mrow><mrow><mi>m</mi><mo>+</mo><mi>n</mi></mrow></mfrac></msqrt><mrow><mo>(</mo><msub><mi>I</mi><mi>B</mi></msub><mo>-</mo><mi>u</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0000808465040000023.GIF" wi="892" he="152" /></maths>c.求解<img file="FDA0000808465040000024.GIF" wi="555" he="81" />其中QR(·)表示QR分解运算;d.令<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mrow><mi>f</mi><mo>&CenterDot;</mo><msub><mi>&Sigma;</mi><mi>d</mi></msub></mrow></mtd><mtd><mrow><msubsup><mi>U</mi><mi>d</mi><mi>T</mi></msubsup><mover><mi>B</mi><mo>&OverBar;</mo></mover></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><msup><mi>B</mi><mrow><mo>*</mo><mi>T</mi></mrow></msup><mrow><mo>(</mo><mover><mi>B</mi><mo>&OverBar;</mo></mover><mo>-</mo><msub><mi>U</mi><mi>d</mi></msub><msubsup><mi>U</mi><mi>d</mi><mi>T</mi></msubsup><mover><mi>B</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000808465040000025.GIF" wi="579" he="221" /></maths>对R进行SVD分解,有<img file="FDA0000808465040000028.GIF" wi="255" he="99" />e.最后得到<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>U</mi><mi>d</mi><mo>*</mo></msubsup><mo>=</mo><mo>&lsqb;</mo><msub><mi>U</mi><mi>d</mi></msub><mo>,</mo><msup><mi>B</mi><mo>*</mo></msup><mo>&rsqb;</mo><mover><mi>U</mi><mo>&OverBar;</mo></mover><mo>,</mo></mrow>]]></math><img file="FDA0000808465040000026.GIF" wi="348" he="87" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>&Sigma;</mi><mi>d</mi><mo>*</mo></msubsup><mo>=</mo><mover><mi>&Sigma;</mi><mo>&OverBar;</mo></mover><mo>,</mo><msup><mi>u</mi><mo>*</mo></msup><mo>=</mo><msub><mi>I</mi><mi>C</mi></msub><mo>.</mo></mrow>]]></math><img file="FDA0000808465040000027.GIF" wi="397" he="88" /></maths>
地址 710072 陕西省西安市友谊西路127号