发明名称 一种多激光跟踪仪快速互瞄定向方法
摘要 本发明涉及一种多激光跟踪仪快速互瞄定向方法,对于某两台激光跟踪仪,包括:在测量空间布设两台激光跟踪仪的站位,分别为第一站位和第二站位,分别将一个球座置于各自激光跟踪仪的俯仰转座上,且偏离原点一定距离;用两台激光跟踪仪互相扫瞄对方球座上的标记点,为姿态1,分别记录两台激光跟踪仪的测量数据;将两台激光跟踪仪的水平角和俯仰角都转180度左右,为姿态2,相互瞄准对方球座上的标记点,分别记录两台激光跟踪仪的测量数据;在两个站位接近中间的连线上放置一个公共点,两个跟踪仪分别测量该点并记录数据;求解两台激光跟踪仪站位之间的位姿关系;步骤六、以上述过程得到的位姿关系为初值,经过优化得到最终的位姿关系数据。
申请公布号 CN106524904A 申请公布日期 2017.03.22
申请号 CN201610832917.7 申请日期 2016.09.19
申请人 天津大学 发明人 林嘉睿;邾继贵;杨凌辉;任永杰;高扬
分类号 G01B11/00(2006.01)I 主分类号 G01B11/00(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 程毓英
主权项 一种多激光跟踪仪快速互瞄定向方法,对于某两台激光跟踪仪,包括以下步骤:步骤一、在测量空间布设两台激光跟踪仪的站位,分别为第一站位和第二站位,分别将一个球座置于各自激光跟踪仪的俯仰转座上,且偏离原点一定距离;步骤二、用两台激光跟踪仪互相扫瞄对方球座上的标记点,为姿态1,分别记录两台激光跟踪仪的测量数据;步骤三、将两台激光跟踪仪的水平角和俯仰角都转180度左右,为姿态2,相互瞄准对方球座上的标记点,分别记录两台激光跟踪仪的测量数据;步骤四、在两个站位接近中间的连线上放置一个公共点,两个跟踪仪分别测量该点并记录数据;步骤五、求解两台激光跟踪仪站位之间的位姿关系,方法如下:步骤5‑1)设两台激光跟踪仪一个为第一跟踪仪,另一个为第二跟踪仪,第二跟踪仪的原点在第一跟踪仪坐标系下的矢量方向为p<sub>12</sub>,相应第一跟踪仪的原点在第二跟踪仪坐标系下的矢量方向为p<sub>21</sub>,由于Z轴平行得到姿态角<img file="FDA0001116683340000011.GIF" wi="455" he="71" />步骤5‑2)平移向量的求取需要求得两站原点之间的距离,设两台跟踪仪的中心连线与公共点之间的夹角分别为φ<sub>1</sub>和φ<sub>2</sub>,则得到两个站位原点之间的距离为:<maths num="0001"><math><![CDATA[<mrow><mi>d</mi><mo>=</mo><msqrt><mrow><msup><msub><mi>d</mi><mrow><mn>1</mn><mi>m</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>d</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub><mn>2</mn></msup><mo>-</mo><mn>2</mn><msub><mi>d</mi><mrow><mn>1</mn><mi>m</mi></mrow></msub><msub><mi>d</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><mi>&pi;</mi><mo>-</mo><msub><mi>&phi;</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&phi;</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow></msqrt></mrow>]]></math><img file="FDA0001116683340000012.GIF" wi="984" he="87" /></maths>其中,d<sub>1m</sub>和d<sub>2m</sub>分别为第一跟踪仪和第二跟踪仪的原点到公共点之间的距离,则第二跟踪仪相对于第一跟踪仪的平移向量T=d*p<sub>21</sub>,由此可得到位姿关系为姿态角θ和平移向量T;步骤六、以上述过程得到的位姿关系为初值,经过优化得到最终的位姿关系数据:步骤6‑1)对于两个标记点,设第一跟踪仪的标记点为P1,第二跟踪仪上的标记点为P2;在姿态1中,标记点P2相对与第一跟踪仪的坐标系L1,球坐标为(d<sub>1</sub>,α<sub>1</sub>,β<sub>1</sub>),d<sub>1</sub>表示标记点到跟踪仪坐标系原点之间的距离,α<sub>1</sub>是水平角,β<sub>1</sub>是垂直角,标记点P2相对于第一跟踪仪坐标系L1的关系表示为三维坐标<img file="FDA0001116683340000013.GIF" wi="123" he="71" />同理在L1坐标系下标记点P1变换后的坐标为<img file="FDA0001116683340000014.GIF" wi="97" he="67" />变化矩阵为一个旋转矩阵<img file="FDA0001116683340000015.GIF" wi="99" he="70" />在L2坐标系下点P2变换后的坐标为<img file="FDA0001116683340000016.GIF" wi="107" he="71" />变化矩阵为一个旋转矩阵<img file="FDA0001116683340000017.GIF" wi="98" he="71" />在姿态2中,L1坐标系下点P1变换后的坐标为<img file="FDA0001116683340000018.GIF" wi="98" he="71" />变化矩阵为一个旋转矩阵<img file="FDA0001116683340000019.GIF" wi="99" he="70" />在L2坐标系下点P2变换后的坐标为<img file="FDA00011166833400000110.GIF" wi="107" he="70" />变化矩阵为一个旋转矩阵<img file="FDA00011166833400000111.GIF" wi="99" he="70" />在这里变换矩阵<img file="FDA00011166833400000112.GIF" wi="75" he="70" />和<img file="FDA00011166833400000113.GIF" wi="74" he="71" />与姿态1中的标记点球坐标中的角度值相关,同理<img file="FDA00011166833400000114.GIF" wi="73" he="70" />和<img file="FDA00011166833400000115.GIF" wi="78" he="70" />与姿态2中标记点球坐标中的角度值相关;在L1坐标系和L2坐标系下公共点的坐标分别为<sub>L2</sub>P<sub>3</sub>和<sub>L1</sub>P<sub>3</sub>;步骤6‑2)标记点坐标之间的关系和两台激光跟踪仪间的位姿关系参数表示为:<maths num="0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>2</mn></mrow><mrow><mi>f</mi><mn>1</mn></mrow></mmultiscripts><mn>1</mn></msub><mo>-</mo><mi>R</mi><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>1</mn></mrow><mrow><mi>f</mi><mn>1</mn></mrow></mmultiscripts><mn>1</mn></msub><mo>-</mo><mi>T</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>2</mn></mrow><mrow><mi>f</mi><mn>1</mn></mrow></mmultiscripts><mn>2</mn></msub><mo>-</mo><mi>R</mi><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>1</mn></mrow><mrow><mi>f</mi><mn>1</mn></mrow></mmultiscripts><mn>2</mn></msub><mo>-</mo><mi>T</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>2</mn></mrow><mrow><mi>f</mi><mn>2</mn></mrow></mmultiscripts><mn>1</mn></msub><mo>-</mo><mi>R</mi><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>1</mn></mrow><mrow><mi>f</mi><mn>2</mn></mrow></mmultiscripts><mn>2</mn></msub><mo>-</mo><mi>T</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>2</mn></mrow><mrow><mi>f</mi><mn>2</mn></mrow></mmultiscripts><mn>2</mn></msub><mo>-</mo><mi>R</mi><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>1</mn></mrow><mrow><mi>f</mi><mn>2</mn></mrow></mmultiscripts><mn>2</mn></msub><mo>-</mo><mi>T</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mmultiscripts><mi>R</mi><mprescripts/><mrow><mi>L</mi><mn>2</mn></mrow><none/></mmultiscripts><mn>3</mn></msub><mo>-</mo><msub><mi>R</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><msub><mi>P</mi><mn>3</mn></msub><mo>-</mo><mi>T</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001116683340000021.GIF" wi="494" he="415" /></maths>位姿关系中的姿态角θ可以表示为旋转矩阵R:<maths num="0003"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mrow><mi>sin</mi><mi>&theta;</mi></mrow></mtd><mtd><mrow><mi>cos</mi><mi>&theta;</mi></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001116683340000022.GIF" wi="549" he="228" /></maths>消去公式中的<img file="FDA0001116683340000023.GIF" wi="339" he="71" />和<img file="FDA0001116683340000024.GIF" wi="107" he="71" />重新整理方程得:<maths num="0004"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msup><mmultiscripts><mi>R</mi><mprescripts/><mrow><mi>L</mi><mn>1</mn></mrow><mrow><mi>f</mi><mn>1</mn></mrow></mmultiscripts><mi>T</mi></msup><mo>&CenterDot;</mo><msup><mi>R</mi><mi>T</mi></msup><mrow><mo>(</mo><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>2</mn></mrow><mrow><mi>f</mi><mn>1</mn></mrow></mmultiscripts><mn>1</mn></msub><mo>-</mo><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msup><mmultiscripts><mi>R</mi><mprescripts/><mrow><mi>L</mi><mn>1</mn></mrow><mrow><mi>f</mi><mn>2</mn></mrow></mmultiscripts><mi>T</mi></msup><mo>&CenterDot;</mo><msup><mi>R</mi><mi>T</mi></msup><mrow><mo>(</mo><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>2</mn></mrow><mrow><mi>f</mi><mn>2</mn></mrow></mmultiscripts><mn>1</mn></msub><mo>-</mo><mi>T</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msup><mmultiscripts><mi>R</mi><mprescripts/><mrow><mi>L</mi><mn>2</mn></mrow><mrow><mi>f</mi><mn>1</mn></mrow></mmultiscripts><mi>T</mi></msup><mrow><mo>(</mo><mi>R</mi><mo>&CenterDot;</mo><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>1</mn></mrow><mrow><mi>f</mi><mn>1</mn></mrow></mmultiscripts><mn>2</mn></msub><mo>+</mo><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msup><mmultiscripts><mi>R</mi><mprescripts/><mrow><mi>L</mi><mn>2</mn></mrow><mrow><mi>f</mi><mn>2</mn></mrow></mmultiscripts><mi>T</mi></msup><mrow><mo>(</mo><mi>R</mi><mo>&CenterDot;</mo><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>1</mn></mrow><mrow><mi>f</mi><mn>2</mn></mrow></mmultiscripts><mn>2</mn></msub><mo>+</mo><mi>T</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>2</mn></mrow><none/></mmultiscripts><mn>3</mn></msub><mo>-</mo><mi>R</mi><mo>&CenterDot;</mo><msub><mmultiscripts><mi>P</mi><mprescripts/><mrow><mi>L</mi><mn>1</mn></mrow><none/></mmultiscripts><mn>3</mn></msub><mo>-</mo><mi>T</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001116683340000025.GIF" wi="998" he="255" /></maths>步骤6‑3)用非线性最优化方法LM算法求解上述方程组,以步骤5得到的位姿关系为迭代初值,得到精确优化的位姿关系。
地址 300072 天津市南开区卫津路92号