发明名称 一种基于时空处理块的连续小波变换目标跟踪方法
摘要 一种基于时空处理块的连续小波变换目标跟踪方法,它有四大步骤:步骤一、获得时空三维处理块B<sub>t</sub>和时空连续小波函数<img file="DDA0000055672870000011.tif" wi="198" he="47" />步骤二、采用Morlet小波作为小波基函数,通过函数变换得到用于对处理块B<sub>t</sub>进行时空连续小波变换的时空三维小波函数<img file="DDA0000055672870000012.tif" wi="937" he="88" />步骤三、根据处理块B<sub>t</sub>的第Z<sub>f</sub>+1帧图像窗与图像序列S的第t-1帧图像的对应关系,得到运动目标在处理块B<sub>t</sub>的第Z<sub>f</sub>+1帧图像窗中的运动参数<img file="DDA0000055672870000013.tif" wi="771" he="60" />步骤四:根据处理块B<sub>t</sub>的第Z<sub>f</sub>+1帧图像窗与图像序列S的第t-1帧图像的对应关系,以及处理块B<sub>t</sub>的第Z<sub>f</sub>+2帧图像窗与图像序列S的第t帧图像的对应关系,得到运动目标在图像序列S的第t帧图像中的运动参数<img file="DDA0000055672870000014.tif" wi="439" he="52" />本发明构思科学合理,程序简约,在机器视觉与模式识别技术领域内具有较好的实用价值和广阔地应用前景。
申请公布号 CN102156993B 申请公布日期 2014.08.20
申请号 CN201110095674.0 申请日期 2011.04.15
申请人 北京航空航天大学 发明人 王睿;赵岩
分类号 G06T7/20(2006.01)I;G06K9/00(2006.01)I 主分类号 G06T7/20(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 一种基于时空处理块的连续小波变换目标跟踪方法,其特征在于:该方法具体步骤如下:步骤一:读入待处理图像序列S中前Z帧图像作为图像序列S<sub>0</sub>,将图像序列S<sub>0</sub>的初始帧图像复制Z<sub>f</sub>次并补充到S<sub>0</sub>前面构成初始过渡图像序列S<sub>1</sub>;设定待处理图像序列S中的帧数为t,判断待处理图像序列S是否有第t‑1+Z帧图像存在即预先获得的图像序列、输入即实时拍摄的图像序列,根据判断结果,采用两种方式由前一过渡图像序列S<sub>t‑1</sub>构造过渡图像序列S<sub>t</sub>;根据运动目标在待处理图像序列S的第t‑1帧图像中的运动参数构造针对待处理图像序列S中第t帧图像的时空三维处理块B<sub>t</sub>,假设目标的初始运动参数已知;步骤二:采用Morlet小波作为小波基函数,通过函数变换得到用于对处理块B<sub>t</sub>进行时空连续小波变换的时空三维小波函数<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>y</mi><mi>w</mi></msub><mrow><mo>(</mo><mover><mi>&xi;</mi><mo>&RightArrow;</mo></mover><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>a</mi><mrow><mo>-</mo><mfrac><mn>3</mn><mn>2</mn></mfrac></mrow></msup><mi>y</mi><mrow><mo>(</mo><msup><mi>a</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>c</mi><mrow><mo>-</mo><mfrac><mn>1</mn><mn>3</mn></mfrac></mrow></msup><msup><mi>r</mi><mrow><mo>-</mo><mi>&theta;</mi></mrow></msup><mrow><mo>(</mo><mover><mi>&xi;</mi><mo>&RightArrow;</mo></mover><mo>-</mo><mover><mi>b</mi><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mo>,</mo><msup><mi>a</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>c</mi><mfrac><mn>2</mn><mn>3</mn></mfrac></msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000492118960000011.GIF" wi="1232" he="136" /></maths>简写为<img file="FDA0000492118960000012.GIF" wi="303" he="94" />式中各符号含义如下:<img file="FDA0000492118960000013.GIF" wi="282" he="91" />为空间坐标,t为时间坐标,a为尺度参数,c为速率参数,θ为方向参数且<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>r</mi><mrow><mo>-</mo><mi>&theta;</mi></mrow></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mi>&theta;</mi></mtd><mtd><mi>sin</mi><mi>&theta;</mi></mtd></mtr><mtr><mtd><mo>-</mo><mi>sin</mi><mi>&theta;</mi></mtd><mtd><mi>cos</mi><mi>&theta;</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000492118960000014.GIF" wi="532" he="173" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>b</mi><mo>&RightArrow;</mo></mover><mo>=</mo><mrow><mo>(</mo><msub><mi>b</mi><mi>x</mi></msub><mo>,</mo><msub><mi>b</mi><mi>y</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000492118960000015.GIF" wi="269" he="100" /></maths>为位移参数,τ为时移参数;步骤三:根据处理块B<sub>t</sub>的第Z<sub>f</sub>+1帧图像窗与待处理图像序列S的第t‑1帧图像的对应关系,得到运动目标在处理块B<sub>t</sub>的第Z<sub>f</sub>+1帧图像窗中的运动参数<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>w</mi><mi>t</mi><mrow><msub><mi>Z</mi><mi>f</mi></msub><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mrow><mo>(</mo><msup><mi>x</mi><mrow><msub><mi>Z</mi><mi>f</mi></msub><mo>+</mo><mn>1</mn></mrow></msup><mo>,</mo><msup><mi>y</mi><mrow><msub><mi>Z</mi><mi>f</mi></msub><mo>+</mo><mn>1</mn></mrow></msup><mo>,</mo><msubsup><mi>c</mi><mi>t</mi><mrow><msub><mi>Z</mi><mi>f</mi></msub><mo>+</mo><mn>1</mn></mrow></msubsup><mo>,</mo><msubsup><mi>&theta;</mi><mi>t</mi><mrow><msub><mi>Z</mi><mi>f</mi></msub><mo>+</mo><mn>1</mn></mrow></msubsup><mo>,</mo><msubsup><mi>a</mi><mi>t</mi><mrow><msub><mi>Z</mi><mi>f</mi></msub><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000492118960000016.GIF" wi="1025" he="102" /></maths>对处理块B<sub>t</sub>进行时空连续小波变换,估计运动目标在处理块第Z<sub>f</sub>+2帧图像窗中的运动参数;式中各符号含义如下:<img file="FDA0000492118960000017.GIF" wi="124" he="96" />表示在该帧图像窗中目标的纵坐标,<img file="FDA0000492118960000018.GIF" wi="132" he="92" />表示在该帧图像窗中目标的横坐标,<img file="FDA0000492118960000019.GIF" wi="125" he="96" />表示在该帧图像窗中目标运动的速率,<img file="FDA00004921189600000110.GIF" wi="130" he="98" />表示在该帧图像窗中目标运动的方向,<img file="FDA00004921189600000111.GIF" wi="132" he="102" />表征在该帧图像窗中目标相对图像窗的运动尺度;步骤四:根据处理块B<sub>t</sub>的第Z<sub>f</sub>+1帧图像窗与待处理图像序列S的第t‑1帧图像的对应关系,以及处理块B<sub>t</sub>的第Z<sub>f</sub>+2帧图像窗与待处理图像序列S的第t帧图像的对应关系,得到运动目标在待处理图像序列S的第t帧图像中的运动参数<img file="FDA0000492118960000021.GIF" wi="618" he="94" />判断待处理图像序列S是否已读取完,根据判断结果,决定跟踪过程是否结束;式中各符号含义如下:<img file="FDA0000492118960000022.GIF" wi="70" he="99" />表示在该帧图像中目标的纵坐标,<img file="FDA0000492118960000023.GIF" wi="73" he="97" />表示在该帧图像中目标的横坐标,<img file="FDA0000492118960000024.GIF" wi="60" he="86" />表示在该帧图像中目标运动的速率,<img file="FDA0000492118960000025.GIF" wi="65" he="93" />表示在该帧图像中目标运动的方向,<img file="FDA0000492118960000026.GIF" wi="71" he="102" />表征在该帧图像中目标相对于处理块B<sub>t</sub>的第Z<sub>f</sub>+2帧图像窗的运动尺度。
地址 100191 北京市海淀区学院路37号