发明名称 基于跟踪时间预测的特征选择方法
摘要 本发明公开了一种基于跟踪时间预测的特征选择算法,包括:随机选择一帧二维图像的n个特征进行特征跟踪,获得特征匹配的结果,并利用该结果估计二维图像的运动信息;对二维图像进行分割,得到一组互相不重叠的图像子区域,将每个子区域的中心像素作为特征点,根据二维图像的运动信息和特征点的当前所在位置,通过前向迭代算法来预测特征点的跟踪时间,并以此作为特征点所在子区域内任一特征的预测跟踪时间,即作为每个子区域的预测跟踪时间;比较所有子区域的预测跟踪时间,在具有最大预测跟踪时间的子区域内进行特征提取。本发明提供的特征选择方法,在用于机器人从环境中选取特征时可以有效减少其自身定位估计的不确定性。
申请公布号 CN101916374B 申请公布日期 2012.10.03
申请号 CN201010262039.2 申请日期 2010.08.20
申请人 浙江大学 发明人 陈耀武;孟旭炯;史勇强;欧进利
分类号 G06K9/46(2006.01)I;G06T7/20(2006.01)I 主分类号 G06K9/46(2006.01)I
代理机构 杭州天勤知识产权代理有限公司 33224 代理人 胡红娟
主权项 1.一种基于跟踪时间预测的特征选择方法,其特征在于,包括以下步骤:(1)随机选择一帧二维图像的n个特征进行特征跟踪,获得特征匹配的结果;并利用所述的特征匹配的结果估计二维图像的运动信息,其过程如下:采用如式(I)所示的六参数仿射模型表示二维图像的运动:u<sub>2</sub>=a<sub>1</sub>×u<sub>1</sub>+a<sub>2</sub>×v<sub>1</sub>+a<sub>3</sub>                     (I)v<sub>2</sub>=a<sub>4</sub>×u<sub>1</sub>+a<sub>5</sub>×v<sub>1</sub>+a<sub>6</sub>式(I),[a<sub>1</sub>,a<sub>2</sub>,a<sub>3</sub>,a<sub>4</sub>,a<sub>5</sub>,a<sub>6</sub>]为所述的六参数仿射模型的参数,a<sub>1</sub>,a<sub>2</sub>,a<sub>4</sub>,a<sub>5</sub>表示图像的缩放和旋转,a<sub>3</sub>和a<sub>6</sub>分别对应上下和左右的平移运动;根据所述的特征匹配的结果,特征点p<sub>i</sub>在k时刻的坐标为<img file="FDA0000137501140000011.GIF" wi="174" he="60" />在k+1时刻的坐标为<img file="FDA0000137501140000012.GIF" wi="228" he="60" />由式(I)得到特征点p<sub>i</sub>的运动方程如式(II)所示:<maths num="0001"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>u</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>v</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>a</mi><mn>1</mn></msub></mtd><mtd><msub><mi>a</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>4</mn></msub></mtd><mtd><msub><mi>a</mi><mn>5</mn></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>u</mi><mi>i</mi><mi>k</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>v</mi><mi>i</mi><mi>k</mi></msubsup></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>a</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>6</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>II</mi><mo>)</mo></mrow></mrow></math>]]></maths>则,n个特征点对应的运动方程组如式(III)所示:<maths num="0002"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>u</mi><mn>1</mn><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>v</mi><mn>1</mn><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msubsup><mi>u</mi><mi>n</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>v</mi><mi>n</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>a</mi><mn>1</mn></msub></mtd><mtd><msub><mi>a</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>4</mn></msub></mtd><mtd><msub><mi>a</mi><mn>5</mn></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>u</mi><mn>1</mn><mi>k</mi></msubsup></mtd><mtd><msubsup><mi>v</mi><mn>1</mn><mi>k</mi></msubsup></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msubsup><mi>u</mi><mn>1</mn><mi>k</mi></msubsup></mtd><mtd><msubsup><mi>v</mi><mn>1</mn><mi>k</mi></msubsup></mtd></mtr><mtr><mtd></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd></mtd><mtd></mtd></mtr><mtr><mtd><msubsup><mi>u</mi><mi>n</mi><mi>k</mi></msubsup></mtd><mtd><msubsup><mi>v</mi><mi>n</mi><mi>k</mi></msubsup></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msubsup><mi>u</mi><mi>n</mi><mi>k</mi></msubsup></mtd><mtd><msubsup><mi>v</mi><mi>n</mi><mi>k</mi></msubsup></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>a</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>6</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>III</mi><mo>)</mo></mrow></mrow></math>]]></maths>求解式(III)所示的运动方程组,确定参数[a<sub>1</sub>,a<sub>2</sub>,a<sub>3</sub>,a<sub>4</sub>,a<sub>5</sub>,a<sub>6</sub>],从而确定表示二维图像运动的六参数仿射模型,以估计二维图像的运动信息;(2)对二维图像进行分割,得到一组互相不重叠的图像子区域,将每个子区域的中心像素作为特征点,根据二维图像的运动信息和特征点的当前所在位置,通过前向迭代算法来预测特征点的跟踪时间,并以此作为特征点所在子区域内任一特征的预测跟踪时间,即作为每个子区域的预测跟踪时间;其中,所述通过前向迭代算法来预测特征点的跟踪时间的方法如下:假设特征在当前时刻位于(u<sub>1</sub>,v<sub>1</sub>)处,则基于式(III)所示的运动方程组的k次前向迭代方程如式(IV)所示:<maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>u</mi><mo>^</mo></mover><mn>2</mn></msub><mo>=</mo><msub><mi>a</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>u</mi><mn>1</mn></msub><mo>+</mo><msub><mi>a</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>v</mi><mn>1</mn></msub><mo>+</mo><msub><mi>a</mi><mn>3</mn></msub></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>v</mi><mo>^</mo></mover><mn>2</mn></msub><mo>=</mo><msub><mi>a</mi><mn>4</mn></msub><mo>&times;</mo><msub><mi>u</mi><mn>1</mn></msub><mo>+</mo><msub><mi>a</mi><mn>5</mn></msub><mo>&times;</mo><msub><mi>v</mi><mn>1</mn></msub><mo>+</mo><msub><mi>a</mi><mn>6</mn></msub></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mo>.</mo><mo>.</mo><mo>.</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>IV</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mover><mi>u</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><msub><mi>a</mi><mn>1</mn></msub><mo>&times;</mo><msub><mover><mi>u</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>a</mi><mn>2</mn></msub><mo>&times;</mo><msub><mover><mi>v</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>a</mi><mn>3</mn></msub></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msub><mover><mi>v</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><msub><mi>a</mi><mn>4</mn></msub><mo>&times;</mo><msub><mover><mi>u</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>a</mi><mn>5</mn></msub><mo>&times;</mo><msub><mover><mi>v</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>a</mi><mn>6</mn></msub></mrow></math>]]></maths>该迭代模型在遇到以下两种情况时自动停止:1、<img file="FDA0000137501140000026.GIF" wi="467" he="54" />其中L和W为图像的长度和宽度;2、迭代次数满k次;因此,预测跟踪时间<img file="FDA0000137501140000027.GIF" wi="30" he="61" />近似表示为式(V)所示<maths num="0008"><![CDATA[<math><mrow><msub><mover><mi>t</mi><mo>^</mo></mover><mi>v</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mi>i</mi><mo>,</mo><mi>i</mi><mo>&lt;</mo><mi>k</mi></mtd></mtr><mtr><mtd><mi>k</mi><mo>,</mo><mi>i</mi><mo>&GreaterEqual;</mo><mi>k</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>(3)比较所有子区域的预测跟踪时间,在具有最大预测跟踪时间的子区域内进行特征提取。
地址 310027 浙江省杭州市西湖区浙大路38号