发明名称 基于激光点云为约束的多波束低掠射角波束归位方法
摘要 本发明公开了一种基于激光点云为约束的多波束低掠射角波束归位方法,属于海洋测绘技术领域,其首先根据沿岸区域内的若干个声速剖面和沿航迹线附近内的表层声速值,利用经验正交函数分析法构建区域内浅层修正后的三维声速场,完成测深数据初始声速改正;然后以同目标激光点云提取的特征点为最外侧声波束位置约束,获取初始声速改正后特征点的坐标偏差,据此计算低掠射层内初始旋转角和初始放大系数,进而按角度线性加权内插60度以上剩余各波束对应的旋转角和放大系数,最终实现波束几何改正。本发明通过这种办法,实现了浅层海水中声波束的准确归位,解决了水岸线上下点云无缝拼接的问题。
申请公布号 CN105258684A 申请公布日期 2016.01.20
申请号 CN201510719540.X 申请日期 2015.10.30
申请人 山东科技大学 发明人 阳凡林;卜宪海;王明伟;马跃;石波;宿殿鹏;卢秀山;景冬
分类号 G01C13/00(2006.01)I;G01S19/42(2010.01)I;G01S7/52(2006.01)I;G01S15/88(2006.01)I 主分类号 G01C13/00(2006.01)I
代理机构 济南舜源专利事务所有限公司 37205 代理人 陈海滨
主权项 一种基于激光点云为约束的多波束低掠射角波束归位方法,采用船载水岸线一体化测量系统,其水上部分包括三维激光扫描仪和组合导航系统,水下部分包括多波束测深仪,多波束测深仪带有具有一定偏转角度的多波束换能器,水上部分与水下部分通过一个平台支架固连;其中的三维激光扫描仪用于测量水面以上地物地貌,多波束测深仪用于测量水下地物地貌,组合导航系统用于为三维激光扫描仪和多波束测深仪提供定位信息、时间信息、姿态信息和航向信息;船载水岸线一体化测量系统还包括表层声速仪和声速剖面仪两个辅助设备;其特征在于:所述的多波束低掠射角波束归位方法按照如下步骤进行:步骤1:对船载水岸线一体化测量系统进行标定,确定多波束测深仪、三维激光扫描仪和组合导航系统之间的空间相对位置关系;步骤2:选取一段具有水上下公共目标的测区,并利用船载水岸线一体化测量系统沿岸匀速进行测量,获取水上下点云同一坐标系下坐标;步骤3:对存储的多波束数据按时间分区;步骤4:根据表层声速仪和组合导航系统中的GPS数据存储的时间标签,获取每一个分区内行船附近的平均表层声速值c(A<sub>ts</sub>)和船附近内一点的地理坐标(x<sub>At</sub>,y<sub>At</sub>);步骤5:利用声速剖面仪均匀获取沿岸测量区域内多个声速剖面数据svp<sub>n</sub>,并记录每一个采样点的地理位置(x<sub>N</sub>,y<sub>N</sub>);步骤6:对步骤5中得到的声速剖面数据svp<sub>n</sub>进行层化处理,获得不同位置点处同一水层上的声速c(z<sub>i</sub>)和每一层声速的平均值<img file="FDA0000834266210000011.GIF" wi="135" he="79" />步骤7:通过步骤6层化处理后的声速剖面基于公式(1)得到测区内每一个点处的声速值:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>c</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>c</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>z</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>k</mi></munderover><msub><mi>&alpha;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>f</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>z</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000834266210000012.GIF" wi="1332" he="132" /></maths>其中加权系数α<sub>k</sub>(x,y)为地理坐标,f<sub>k</sub>(z<sub>i</sub>)为正交函数,<img file="FDA0000834266210000013.GIF" wi="115" he="85" />为每一层声速的平均值;步骤8:根据各分区内平均表层声速值c(A<sub>ts</sub>),对该区域内(x<sub>At</sub>,y<sub>At</sub>)处拟合的声速剖面表层部分进行修正;步骤9:根据各区域修正后的声速剖面,进行声速改正,获取多波束测深仪中的多波束换能器坐标系下水上下同目标特征点处激光点与多波束点的平均坐标差,计算低掠射层内初始几何旋转角<img file="FDA0000834266210000014.GIF" wi="72" he="79" />和初始放大系数<img file="FDA0000834266210000015.GIF" wi="93" he="78" />进而按角度线性加权内插计算低掠射层内60度以上剩余波束角对应的旋转角<img file="FDA0000834266210000016.GIF" wi="77" he="71" />和放大系数<img file="FDA0000834266210000017.GIF" wi="93" he="71" />步骤10:对低掠射层的固定层部分内的波束进行初始旋转角为<img file="FDA0000834266210000021.GIF" wi="74" he="78" />和初始放大系数为<img file="FDA0000834266210000022.GIF" wi="69" he="77" />的几何改正,对低掠射层的自由层部分各波束进行旋转角为<img file="FDA0000834266210000023.GIF" wi="76" he="76" />放大系数为<img file="FDA0000834266210000024.GIF" wi="70" he="70" />的几何改正。
地址 266590 山东省青岛市经济技术开发区前湾港路579号