发明名称 一种基于双目视觉的空间目标旋转状态辨识方法
摘要 本发明涉及一种基于双目视觉的空间目标旋转状态辨识方法,技术特征在于:将双目相机置于目标附近的任意位置P点,以目标上任意的两个标志点A和B为特征点,计算旋转轴的方向向量,计算标志点A和B在旋转平面上的投影,计算操作目标的旋转角度,计算旋转角速度,得到目标的旋转角速度矢量和目标的旋转轴方程辨识结果。有益效果:计算量小,辨识速度快,辨识精度高,易于实现的特点,弥补了现有方法无法辨识目标旋转轴的缺陷,能够被用来对自旋卫星、失控旋转卫星等空间目标的旋转状态进行辨识,为空间旋转目标的近距离空间操作奠定了基础,有着巨大的应用潜力。
申请公布号 CN102981011A 申请公布日期 2013.03.20
申请号 CN201210524887.5 申请日期 2012.12.10
申请人 西北工业大学 发明人 黄攀峰;胡仄虹;刘正雄;孟中杰
分类号 G01P3/38(2006.01)I;G01P13/00(2006.01)I;G01P13/04(2006.01)I 主分类号 G01P3/38(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种基于双目视觉的空间目标旋转状态的辨识方法,其特征在于步骤如下:步骤1:将双目相机置于目标附近的任意位置P点,以目标上任意的两个标志点A和B为特征点,以t<sub>0</sub>时刻两点在相机坐标系中的位置向量r<sub>A0</sub>和r<sub>B0</sub>及t<sub>0</sub>+ΔT时刻两点在相机坐标系中的位置矢量r<sub>A1</sub>和r<sub>B1</sub>,计算角度量θ:<maths num="0001"><![CDATA[<math><mrow><mi>&theta;</mi><mo>=</mo><mi>arccos</mi><mo>[</mo><mfrac><mrow><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>A</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><mi>A</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>B</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><mi>B</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow></mrow><mrow><mo>|</mo><msub><mi>r</mi><mrow><mi>A</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><mi>A</mi><mn>0</mn></mrow></msub><mo>|</mo><mo>|</mo><msub><mi>r</mi><mrow><mi>B</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><mi>B</mi><mn>0</mn></mrow></msub><mo>|</mo></mrow></mfrac><mo>]</mo></mrow></math>]]></maths>若计算的|θ|>5°且|180-θ|>5°,转入步骤2,否则需要重新选择标志点A和B;步骤2计算旋转轴的方向向量:若r<sub>flag</sub>·r<sub>s</sub>≥0,则旋转轴的方向向量S=r<sub>s</sub>;若r<sub>flag</sub>·r<sub>s</sub><0,则旋转轴的方向向量S=-r<sub>s</sub>;其中:r<sub>flag</sub>=(r<sub>B0</sub>-r<sub>A0</sub>)×(r<sub>B1</sub>-r<sub>A1</sub>),<maths num="0002"><![CDATA[<math><mrow><msub><mi>r</mi><mi>S</mi></msub><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>A</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><mi>A</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>B</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><mi>B</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow></mrow><mrow><mo>|</mo><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>A</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><mi>A</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>B</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><mi>B</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>|</mo></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>步骤3计算标志点A和B在旋转平面上的投影r<sub>PA0</sub>、r<sub>PA1</sub>、r<sub>PB0</sub>、r<sub>PB0</sub>:r<sub>PA0</sub>=r<sub>A0</sub>-(r<sub>A0</sub>·S)S,r<sub>PB0</sub>=r<sub>B0</sub>-(r<sub>B0</sub>·S)S,r<sub>PA1</sub>=r<sub>A1</sub>-(r<sub>A1</sub>·S)S,r<sub>PB0</sub>=r<sub>B0</sub>-(r<sub>B0</sub>·S)S;计算投影长度比γ:<img file="FDA00002552137900013.GIF" wi="317" he="137" />若γ≥cos(80°)=0.174,则转入步骤4,否则需要返回步骤1,重新选择标志点A和B后重复步骤1~3;步骤4计算操作目标的旋转角度<img file="FDA00002552137900014.GIF" wi="55" he="35" /><img file="FDA00002552137900015.GIF" wi="743" he="147" />若<img file="FDA00002552137900016.GIF" wi="175" he="62" />则转入步骤5,否则需要回到步骤1,增大观测时间ΔT后重复步骤1~4;步骤5差分计算旋转角速度ω:<img file="FDA00002552137900021.GIF" wi="178" he="96" />计算目标旋转运动四元素如下:<maths num="0003"><![CDATA[<math><mrow><msub><mi>q</mi><mn>1</mn></msub><mo>=</mo><msub><mi>S</mi><mi>x</mi></msub><mi>sin</mi><mfrac><mi>&phi;</mi><mn>2</mn></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>q</mi><mn>2</mn></msub><mo>=</mo><msub><mi>S</mi><mi>y</mi></msub><mi>sin</mi><mfrac><mi>&phi;</mi><mn>2</mn></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>q</mi><mn>3</mn></msub><mo>=</mo><msub><mi>S</mi><mi>z</mi></msub><mi>sin</mi><mfrac><mi>&phi;</mi><mn>2</mn></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>q</mi><mn>4</mn></msub><mo>=</mo><mi>cos</mi><mfrac><mi>&phi;</mi><mn>2</mn></mfrac></mrow></math>]]></maths>式中,S<sub>x</sub>、S<sub>y</sub>和S<sub>z</sub>分别表示旋转轴方向向量S在三个坐标轴上的分量;得到目标旋转矩阵R为:<maths num="0007"><![CDATA[<math><mrow><mi>R</mi><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>q</mi><mn>1</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mn>2</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mn>3</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>q</mi><mn>4</mn><mn>2</mn></msubsup></mtd><mtd><mn>2</mn><mrow><mo>(</mo><msub><mi>q</mi><mn>1</mn></msub><msub><mi>q</mi><mn>2</mn></msub><mo>+</mo><msub><mi>q</mi><mn>3</mn></msub><msub><mi>q</mi><mn>4</mn></msub><mo>)</mo></mrow></mtd><mtd><mn>2</mn><mrow><mo>(</mo><msub><mi>q</mi><mn>1</mn></msub><msub><mi>q</mi><mn>3</mn></msub><mo>-</mo><msub><mi>q</mi><mn>2</mn></msub><msub><mi>q</mi><mn>4</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>2</mn><mrow><mo>(</mo><msub><mi>q</mi><mn>1</mn></msub><msub><mi>q</mi><mn>2</mn></msub><mo>-</mo><msub><mi>q</mi><mn>3</mn></msub><msub><mi>q</mi><mn>4</mn></msub><mo>)</mo></mrow></mtd><mtd><mo>-</mo><msubsup><mi>q</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>q</mi><mn>2</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mn>3</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>q</mi><mn>4</mn><mn>2</mn></msubsup></mtd><mtd><mn>2</mn><mrow><mo>(</mo><msub><mi>q</mi><mn>2</mn></msub><msub><mi>q</mi><mn>3</mn></msub><mo>+</mo><msub><mi>q</mi><mn>1</mn></msub><msub><mi>q</mi><mn>4</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>2</mn><mrow><mo>(</mo><msub><mi>q</mi><mn>1</mn></msub><msub><mi>q</mi><mn>3</mn></msub><mo>+</mo><msub><mi>q</mi><mn>2</mn></msub><msub><mi>q</mi><mn>4</mn></msub><mo>)</mo></mrow></mtd><mtd><mn>2</mn><mrow><mo>(</mo><msub><mi>q</mi><mn>2</mn></msub><msub><mi>q</mi><mn>3</mn></msub><mo>-</mo><msub><mi>q</mi><mn>1</mn></msub><msub><mi>q</mi><mn>4</mn></msub><mo>)</mo></mrow></mtd><mtd><mo>-</mo><msubsup><mi>q</mi><mn>1</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>q</mi><mn>3</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>q</mi><mn>4</mn><mn>2</mn></msubsup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>;</mo></mrow></math>]]></maths>步骤6计算旋转平面中心点C的坐标近似值<img file="FDA00002552137900027.GIF" wi="63" he="51" /><maths num="0008"><![CDATA[<math><mrow><msub><mover><mi>r</mi><mo>~</mo></mover><mi>C</mi></msub><mo>=</mo><mrow><mo>(</mo><msup><mi>A</mi><mi>T</mi></msup><mi>A</mi><mo>)</mo></mrow><msup><mi>A</mi><mi>T</mi></msup><mi>b</mi></mrow></math>]]></maths>式中:<maths num="0009"><![CDATA[<math><mrow><mi>A</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>I</mi><mo>-</mo><mi>R</mi></mtd></mtr><mtr><mtd><mi>I</mi><mo>-</mo><mi>R</mi></mtd></mtr><mtr><mtd><msup><mi>S</mi><mi>T</mi></msup></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><maths num="0010"><![CDATA[<math><mrow><mi>b</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>r</mi><mrow><mi>PA</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>Rr</mi><mrow><mi>PA</mi><mn>0</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>r</mi><mrow><mi>PB</mi><mn>1</mn></mrow></msub><mo>-</mo><mi>R</mi><msub><mi>r</mi><mrow><mi>PB</mi><mn>0</mn></mrow></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>其中,I表示3×3的单位矩阵;步骤7:目标的旋转角速度矢量的辨识结果为:ωS;目标的旋转轴方程的辨识结果为:<img file="FDA000025521379000211.GIF" wi="518" he="121" />其中:r<sub>Cx</sub>、r<sub>Cy</sub>和r<sub>Cz</sub>分别表示旋转中心坐标近似值<img file="FDA000025521379000212.GIF" wi="39" he="51" />的在三个坐标轴上的分量,即<img file="FDA000025521379000213.GIF" wi="380" he="80" />
地址 710072 陕西省西安市友谊西路127号