发明名称 一种动态平台姿态和高度测量方法
摘要 本发明提供一种动态平台的姿态和高度测量方法。技术方案是利用两台正交安装、竖直向下扫描的二维激光扫描仪获得路面扫描线;使用RANSAC算法估计出扫描线的表达式方程,利用两扫描线的表达式方程得到道路平面方程;以道路平面方程为参考解算平台的姿态和高度。本发明使用的设备简单,成本低。测量结果不受运动平台自身加速度的影响,没有漂移,测量精度高。
申请公布号 CN103090863A 申请公布日期 2013.05.08
申请号 CN201310038521.1 申请日期 2013.01.31
申请人 中国人民解放军国防科学技术大学 发明人 冯莹;曹毓;雷兵;魏立安;程金龙;马相路;赵立双
分类号 G01C21/00(2006.01)I;G01C21/10(2006.01)I;G01C3/00(2006.01)I 主分类号 G01C21/00(2006.01)I
代理机构 国防科技大学专利服务中心 43202 代理人 王文惠
主权项 1.一种动态平台姿态和高度测量方法,其特征在于,将两台二维激光扫描仪固定在姿态和高度待测的平台上,两台扫描仪均以竖直向下的方式同步扫描路面,并且两台扫描仪的扫描平面相互正交;建立扫描仪系统坐标系X<sub>S-scanner</sub>,原点为两个扫描仪的扫描平面与平台所在平面三个平面的交点,X轴平行于平台所在平面但垂直于平台前进方向,Y轴平行于平台所在平面且与平台前进方向相反,Z轴垂直于平台所在平面且向上为正向,X轴、Y轴和Z轴构成右手坐标系;实施下述步骤:第一步,采集数据:平台在路面上移动时,两台扫描仪S<sub>1</sub>和S<sub>2</sub>同步扫描,采集数据;设S<sub>1</sub>和S<sub>2</sub>在某一时刻T分别扫描获得一条扫描线,令其分别为L<sub>1</sub>和L<sub>2</sub>;每条扫描线由n个扫描点组成,L<sub>1</sub>在扫描仪系统坐标系X<sub>S-scanner</sub>中第i个扫描点的坐标为<img file="FDA00002804187300011.GIF" wi="67" he="59" />(i=1,2,…,n),L<sub>2</sub>在扫描仪系统坐标系X<sub>S-scanner</sub>中第i个扫描点的坐标为<img file="FDA00002804187300012.GIF" wi="72" he="59" />(i=1,2,…,n);第二步,计算T时刻扫描线L<sub>1</sub>和L<sub>2</sub>的表达式方程:利用扫描线L<sub>1</sub>和L<sub>2</sub>所包含的扫描点坐标,使用RANSAC算法估计T时刻L<sub>1</sub>和L<sub>2</sub>的表达式方程,分别为:<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>L</mi><mn>1</mn></msub><mo>:</mo></mtd><mtd><msub><mi>a</mi><mn>1</mn></msub><mi>x</mi><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mi>z</mi><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><mo>=</mo><mn>0</mn><mo>,</mo></mtd><mtd><mi>y</mi><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>L</mi><mn>2</mn></msub><mo>:</mo></mtd><mtd><msub><mi>a</mi><mn>2</mn></msub><mi>x</mi><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mi>z</mi><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mo>=</mo><mn>0</mn><mo>,</mo></mtd><mtd><mi>x</mi><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced></math>]]></maths>公式一上式中,(x,y,z)表示在扫描仪系统坐标系X<sub>S-scanner</sub>中任意一点的坐标;第三步,计算道路平面方程:由扫描线L<sub>1</sub>和L<sub>2</sub>的直线方程,得到T时刻道路平面方程,如公式二所示:<maths num="0002"><![CDATA[<math><mrow><mi>A</mi><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mfrac><mrow><msub><mi>b</mi><mn>1</mn></msub><msub><mi>t</mi><mn>1</mn></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mrow><mo>(</mo><mi>y</mi><mo>-</mo><mfrac><mrow><msub><mi>b</mi><mn>2</mn></msub><msub><mi>t</mi><mn>2</mn></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>+</mo><mi>C</mi><mrow><mo>(</mo><mi>z</mi><mo>+</mo><mfrac><mrow><msub><mi>c</mi><mn>1</mn></msub><mo>/</mo><msub><mi>b</mi><mn>1</mn></msub><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mo>/</mo><msub><mi>b</mi><mn>2</mn></msub><mo>+</mo><msub><mi>a</mi><mn>1</mn></msub><msub><mi>t</mi><mn>1</mn></msub><mo>+</mo><msub><mi>a</mi><mn>2</mn></msub><msub><mi>t</mi><mn>2</mn></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></math>]]></maths>公式二其中:<maths num="0003"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>A</mi><mo>=</mo><msub><mi>a</mi><mn>1</mn></msub><msub><mi>b</mi><mn>2</mn></msub></mtd><mtd></mtd><mtd></mtd></mtr><mtr><mtd><mi>B</mi><mo>=</mo><msub><mi>a</mi><mn>2</mn></msub><msub><mi>b</mi><mn>1</mn></msub><mo>,</mo></mtd><mtd><msub><mi>t</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mfrac><mrow><msub><mi>a</mi><mn>1</mn></msub><msub><mi>b</mi><mn>2</mn></msub></mrow><msub><mi>b</mi><mn>1</mn></msub></mfrac><mrow><mo>(</mo><msub><mi>c</mi><mn>2</mn></msub><msub><mi>b</mi><mn>1</mn></msub><mo>-</mo><msub><mi>c</mi><mn>1</mn></msub><msub><mi>b</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow><mrow><msubsup><mi>a</mi><mn>1</mn><mn>2</mn></msubsup><msubsup><mi>b</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>a</mi><mn>2</mn><mn>2</mn></msubsup><msubsup><mi>b</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>b</mi><mn>1</mn><mn>2</mn></msubsup><msubsup><mi>b</mi><mn>2</mn><mn>2</mn></msubsup></mrow></mfrac><mo>,</mo></mtd><mtd><msub><mi>t</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mfrac><mrow><msub><mi>a</mi><mn>2</mn></msub><msub><mi>b</mi><mn>1</mn></msub></mrow><msub><mi>b</mi><mn>2</mn></msub></mfrac><mrow><mo>(</mo><msub><mi>c</mi><mn>2</mn></msub><msub><mi>b</mi><mn>1</mn></msub><mo>-</mo><msub><mi>c</mi><mn>1</mn></msub><msub><mi>b</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow><mrow><msubsup><mi>a</mi><mn>1</mn><mn>2</mn></msubsup><msubsup><mi>b</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>a</mi><mn>2</mn><mn>2</mn></msubsup><msubsup><mi>b</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>b</mi><mn>1</mn><mn>2</mn></msubsup><msubsup><mi>b</mi><mn>2</mn><mn>2</mn></msubsup></mrow></mfrac></mtd></mtr><mtr><mtd><mi>C</mi><mo>=</mo><msub><mi>b</mi><mn>1</mn></msub><msub><mi>b</mi><mn>2</mn></msub></mtd><mtd></mtd><mtd></mtd></mtr></mtable></mfenced></math>]]></maths>公式三第四步,解算平台姿态角和高度:T时刻平台姿态角由如下公式计算出:<maths num="0004"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>&phi;</mi><mo>=</mo><mi>a</mi><mi>tan</mi><mrow><mo>(</mo><mi>B</mi><mo>/</mo><mi>C</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&theta;</mi><mo>=</mo><mi>a</mi><mi>sin</mi><mo>[</mo><mi>sign</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mi>A</mi><mi>C</mi></mfrac><mo>)</mo></mrow><msqrt><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mrow><mo>(</mo><mi>B</mi><mo>/</mo><mi>A</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>C</mi><mo>/</mo><mi>A</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></msqrt><mo>]</mo></mtd></mtr></mtable></mfenced></math>]]></maths>公式四其中φ和θ分别为平台的俯仰角和横滚角;T时刻平台上待测点M的高度值H由如下公式计算出:<maths num="0005"><![CDATA[<math><mrow><mi>H</mi><mo>=</mo><mfrac><mrow><mo>|</mo><mi>A</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>M</mi></msub><mo>-</mo><mfrac><mrow><msub><mi>b</mi><mn>1</mn></msub><msub><mi>t</mi><mn>1</mn></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>M</mi></msub><mo>-</mo><mfrac><mrow><msub><mi>b</mi><mn>2</mn></msub><msub><mi>t</mi><mn>2</mn></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>+</mo><mi>C</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>M</mi></msub><mo>+</mo><mfrac><mrow><msub><mi>c</mi><mn>1</mn></msub><mo>/</mo><msub><mi>b</mi><mn>1</mn></msub><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mo>/</mo><msub><mi>b</mi><mn>2</mn></msub><mo>+</mo><msub><mi>a</mi><mn>1</mn></msub><msub><mi>t</mi><mn>1</mn></msub><mo>+</mo><msub><mi>a</mi><mn>2</mn></msub><msub><mi>t</mi><mn>2</mn></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>|</mo></mrow><msqrt><msup><mi>A</mi><mn>2</mn></msup><mo>+</mo><msup><mi>B</mi><mn>2</mn></msup><mo>+</mo><msup><mi>C</mi><mn>2</mn></msup></msqrt></mfrac></mrow></math>]]></maths>公式五其中,(x<sub>M</sub>,y<sub>M</sub>,z<sub>M</sub>)为待测点M在扫描仪系统坐标系中的坐标。
地址 410073 湖南省长沙市开福区德雅路109号