CN112083429A - 一种基于l-m算法的多波束横摇残差改正方法 - Google Patents
一种基于l-m算法的多波束横摇残差改正方法 Download PDFInfo
- Publication number
- CN112083429A CN112083429A CN202010967899.XA CN202010967899A CN112083429A CN 112083429 A CN112083429 A CN 112083429A CN 202010967899 A CN202010967899 A CN 202010967899A CN 112083429 A CN112083429 A CN 112083429A
- Authority
- CN
- China
- Prior art keywords
- water depth
- ping
- error
- sampling
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明提供一种基于L‑M算法的多波束横摇残差改正方法,属于多波束测量数据误差处理方法研究领域。本发明主要通过建立水深与横摇残差之间的函数关系式,采用非线性最小二乘L‑M算法并结合剖面趋势线构建思想进行横摇残差提取,最终实现对海底地形的改正。其中,剖面趋势线构建的充分顾及了海底基本走势,综合考虑了横摇残差对水深的影响特点,从而使残差提取精度更高。该方法十分有望成为多波束数据处理中改正运动残差问题的关键方法之一。
Description
技术领域
本发明属于多波束测量数据误差处理方法研究领域,尤其涉及多波束横摇残差的改正方法。
背景技术
多波束测深系统(multi-beam echo sounder,MBES)已成为目前水下地形测量的主要手段之一,它的广泛应用使水下地形测量变得更为精细、高效的同时,也需要顾及更多误差对测量数据质量的影响,例如换能器安装偏差、姿态误差、声速误差、接收波束角测量误差等。如何减弱或消除这些误差的影响,提高多波束测量的数据质量一直是国内外专家学者研究的重点。多波束测量误差按照其性质可分为静态误差和动态误差,国内外许多研究人员对静态误差进行了较为深入的研究并提出了相应的改正方法,取得了较好的效果,而动态误差则因其影响机制复杂、影响因素多样而处理难度较大,但其对提高测量数据精度以及海底精细化探测同样产生较大影响,例如海底“起伏”现象。
多波束运动残差主要分为:姿态传感器固有误差(尺度误差和时延)、横/纵摇轴与参考系不完全对准的误差以及换能器固定杆晃动误差。实际情况下,由于舷侧安装的固定支点太高或者换能器固定杆太细而不足以抗衡测量船运动引起的换能器侧向阻力影响,为此在测量过程中不可避免地导致换能器固定杆发生晃动和轻微形变,那么换能器姿态不再是船体的姿态,而是附加了除测船姿态外的换能器抖动。多波束作业一般采用姿态传感器与换能器分开安装方式,因而姿态传感器测量的是测船的实时姿态,而不是实际换能器的实时姿态。运动残差的本质均是直接产生较大的横摇误差,使多波束数据姿态改正不彻底,从而间接导致数据条带边缘产生“起伏”现象。
当前对相关问题的研究多是针对某单一误差源的影响进行分析与改正,然而运动残差误差源并非单一存在多为共同影响。因此,亟需开展针对多误差源的运动残差改正方法研究,从而进一步提高多波束测深数据精度。
发明内容
本发明主要针对多波束测深系统由于存在运动残差而产生横摇误差,从而导致数据条带边缘呈现“波浪”状周期性起伏的问题,提出一种基于L-M(Levenberg-Marquardt)算法的多波束横摇残差改正方法。通过建立水深与横摇残差之间的函数关系式,采用非线性最小二乘L-M算法并结合剖面趋势线构建思想进行横摇残差提取,最终实现对海底地形的改正。
为了达到上述目的,本发明的技术方案为:
一种基于L-M算法的多波束横摇残差改正方法,具体包括以下步骤:
第一步,确定误差参数
残差改正的目的是获得一组横摇时序改正数,使利用改正横摇得到的水深计算值与真值达到最佳拟合,可以形象的认为是将海底“起伏”逐渐“熨平”的过程。结合横摇残差影响性质,采用非线性最小二乘L-M算法进行残差的提取,由于L-M算法具有较强的鲁棒性,因此更适合解决影响因素较为复杂的非线性运动残差问题。在构建误差方程之前,需要确定各个参数变量。首先,定义变量m为:
式中,Tping为波束采样周期;tpos为姿态数据采样采样周期;中括号表示取整。
由于多波束测深系统的水深数据和姿态数据的采样周期一般不同,并且姿态采样频率一般大于水深采样频率。因此,在姿态改正前需要对水深采样时刻的姿态信息进行内插求解,这里采用时域内反距离内插法。利用水深采样时刻附近若干姿态采样值及其采样时间进行定权内插。此时,设每ping用于内插的姿态采样值个数为s,将s个内插姿态采样值作为误差方程参数。
第二步,构建误差方程
在问题条带内沿航迹方向选取连续N个ping的水深数据,并进一步提取每ping中某一固定入射角度附近的若干波束水深数据作为观测值。选取同一条带中间偏外部40°~60°波束的水深数据,例如50°波束。为确保方程有解,必须满足观测方程个数大于参数个数,因此至少选取s+1个波束点水深数据。根据非线性最小二乘思想,建立各ping水深值与内插姿态采样值之间的误差方程关系式:
式中,F表示点位坐标归算算法;j为采样ping的序号,j=1,2,3…N;i为ping内的波束点序号,i=1,2,3…s+1;k为内插姿态采样的序号,k=1,2,3…s;θij、tij分别为第j号ping中选取的第i个波束点的回波角度和回波时间;rjk为第j号ping用于姿态内插的第k个横摇采样值;ajk为内插时rjk的权值,这里采用时域内反距离内插法,即基于时间间隔的反距离加权。
假设测量回波角度和时间无误差,根据L-M算法的线性化原则,将上式进行泰勒级数展开并取至一阶近似可得:
由此得到误差方程为:
多波束测线一般平行于等深线方向布设,因此沿测线方向水深变化相对稳定,并且每次选取的解算ping一般时长较短,从而可将海底描述为线性函数,即利用加权最小二乘法对沿航迹向左右舷50°附近波束测量水深进行线性拟合,构建左右舷两侧50°水深剖面趋势线。将两侧剖面中同一ping的测量水深与趋势线差值之和视为由横摇残差造成的剖面内该ping总水深误差,基于横摇残差对同一ping左右舷相同波束入射角处波束测量水深影响大小相同、方向相反的特点,对左右舷水深剖面趋势线进行改正,使同一ping在左右舷剖面中的测量水深与改正后趋势线水深之差,大小相同、符号相反,且均为总水深误差的一半。此时,改正后左右舷剖面趋势线即可认为是理想海底,将剖面趋势线中相应ping的水深值作为误差方程中常数项的近似值。
第三步,求解误差参数
式中,Lj为第j号采样ping的观测值与近似值的误差矩阵;权阵P要根据系数阵中所选波束入射角度余弦值大小而定,权大小由中央向边缘依次降低。研究发现,式(8)中法方程的系数方阵Aj TPAj接近奇异值成为病态矩阵。L-M算法采用岭估计的方法来改善方阵的病态性,即在系数方阵的主对角线上加上一个常数B:
式中,E为单位矩阵。选用不同的B值可以得到不同的岭估计值,当B取零时,即是最小二乘估计。
本发明的有益效果为;本发明提供的一种基于L-M算法的多波束横摇残差改正方法,主要针对多波束测深系统由于存在运动残差而产生横摇误差,从而导致水深数据姿态改正不彻底并且数据条带边缘呈现“波浪”状周期性起伏的问题,采用非线性最小二乘L-M算法并结合剖面趋势线构建思想进行横摇残差提取,最终实现对海底地形的改正。其中,剖面趋势线构建的充分顾及了海底基本走势,综合考虑了横摇残差对水深的影响特点,从而使残差提取精度更高。该方法十分有望成为多波束数据处理中改正运动残差问题的关键方法之一。
附图说明
图1为实验1数据处理前后的三维海底地形图,其中,(a)为处理前,(b)为处理后;
图2为实验1数据处理前后左右舷50°波束测量水深剖面的对比图;其中,(a)为右舷,(b)为左舷;
图3为实验2数据处理前后的三维海底地形图,其中,(a)为处理前,(b)为处理后;
图4为实验2数据处理前后左右舷50°波束测量水深剖面的对比图;其中,(a)为右舷,(b)为左舷;
图5为本发明方法的流程图。
具体实施方式
为使本发明解决的方法问题、采用的方法方案和达到的方法效果更加清楚,下面结合附图和实验对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实验仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部内容。
实验数据来源于2019年的多波束实测数据,采用的是R2sonic 2024浅水多波束测量系统,水深约为23米左右。利用Caris对水深数据进行处理时发现条带边缘出现了明显的规律性“起伏”,幅度约为一米左右,初步判定是运动残差导致。选取问题条带中长度200ping的数据段,提取该数据段内的姿态信息、水深信息、定位信息、声速信息、潮汐信息等相关数据内容。该次测量所采用的多波束测深系统的波束采样周期Tping为0.11s,姿态采样周期tpos为0.02s。利用剖面趋势线的构建思想对左右舷50°波束水深剖面进行计算,将剖面趋势线中相应ping的水深值作为误差方程中常数项的近似值。权阵P的确定原则为利用入射角余弦值作为权。在解算法方程时发现其系数方阵接近奇异值,利用岭估计法改善其病态性,最终确定的B值为0.01
观察图1可以发现,改正前边缘波束较为明显的规律性起伏问题在改正后得到很好的改善,沿航迹方向地形过渡更加平滑,更为接近真实海底变化。
表1实验1数据处理前后左右舷50°波束测量水深剖面各项统计数据。
表2实验2数据处理前后左右舷50°波束测量水深剖面各项统计数据。
图2中剖面比对数据见表1。通过表1可以看出,左右舷剖面水深变化范围明显缩小,右舷剖面水深变化区间由改正前的0.9米降至改正后0.5米,左舷则由0.96米降至0.44米。同时剖面水深标准差也明显降低,左舷剖面标准差降低约62%,右舷标准差降低约57%,左右舷起伏地形均得到一定程度的校正。改正前后左右舷剖面水深均值几乎不变,这也证明了该方法的可行性,即在对“起伏”进行抑制的同时也保证了海底的基本走势。
选取另一段问题数据进行实验2,利用实验1相同处理步骤对该数据进行处理,结果如图3和图4所示。
根据实验结果可以看出,边缘波束的规律性起伏问题同样得到了很好的削弱,右舷剖面水深变化区间由改正前的1.75米降至改正后0.76米,左舷则由1.66米降至0.79米,左右舷剖面标准差分别降低约58.2%和57.6%,并且改正前后剖面水深均值几乎不变。
最后应说明的是:以上各实验仅用以说明本发明的方法方案,而非对其限制;尽管已经对本发明进行了详细的说明,本领域的普通方法人员应当理解:其对前述所记载的方法方案进行修改,或者对其中部分或者全部方法特征进行等同替换,并不使相应方法方案的本质脱离本发明方法方案的范围。
Claims (1)
1.一种基于L-M算法的多波束横摇残差改正方法,其特征在于,具体包括以下步骤:
第一步,确定误差参数
采用非线性最小二乘L-M算法进行残差的提取,在构建误差方程之前,需要确定各个误差参数变量,定义误差参数变量m为:
式中,Tping为波束采样周期;tpos为姿态数据采样采样周期;中括号表示取整;
在姿态改正前需要采用时域内反距离内插法对水深采样时刻的姿态信息进行内插求解;利用水深采样时刻附近若干姿态采样值及其采样时间进行定权内插;设每ping用于内插的姿态采样值个数为s,将s个内插姿态采样值作为误差方程参数;
第二步,构建误差方程
在问题条带内沿航迹方向选取连续N个ping的水深数据,并进一步提取每ping中某一固定入射角度附近的若干波束水深数据作为观测值;选取同一条带中间偏外部40°~60°波束的水深数据作为观测值;必须满足观测方程个数大于参数个数,因此至少选取s+1个波束点水深数据;根据非线性最小二乘思想,建立各ping水深值与内插姿态采样值之间的误差方程关系式:
式中,F表示点位坐标归算算法;j为采样ping的序号,j=1,2,3…N;i为ping内的波束点序号,i=1,2,3…s+1;k为内插姿态采样的序号,k=1,2,3…s;θij、tij分别为第j号ping中选取的第i个波束点的回波角度和回波时间;rjk为第j号ping用于姿态内插的第k个横摇采样值;ajk为内插时rjk的权值,采用时域内反距离内插法,即基于时间间隔的反距离加权;
设测量回波角度和时间无误差,根据L-M算法的线性化原则,将上式进行泰勒级数展开并取至一阶近似获得:
由此得到误差方程为:
利用加权最小二乘法对沿航迹向左右舷50°附近波束测量水深进行线性拟合,构建左右舷两侧50°水深剖面趋势线;将两侧剖面中同一ping的测量水深与趋势线差值之和视为由横摇残差造成的剖面内该ping总水深误差,基于横摇残差对同一ping左右舷相同波束入射角处波束测量水深影响大小相同、方向相反的特点,对左右舷水深剖面趋势线进行改正,使同一ping在左右舷剖面中的测量水深与改正后趋势线水深之差,大小相同、符号相反,且均为总水深误差的一半;此时,改正后左右舷剖面趋势线即认为是理想海底,将剖面趋势线中相应ping的水深值作为误差方程中常数项的近似值;
第三步,求解误差参数
式中,Lj为第j号采样ping的观测值与近似值的误差矩阵;权阵P要根据系数阵中所选波束入射角度余弦值大小而定,权大小由中央向边缘依次降低;式(8)中法方程的系数方阵Aj TPAj接近奇异值成为病态矩阵;L-M算法采用岭估计的方法来改善方阵的病态性,即在系数方阵的主对角线上加上一个常数B:
式中,E为单位矩阵;选用不同的B值得到不同的岭估计值,当B取零时,即是最小二乘估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010967899.XA CN112083429B (zh) | 2020-09-15 | 2020-09-15 | 一种基于l-m算法的多波束横摇残差改正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010967899.XA CN112083429B (zh) | 2020-09-15 | 2020-09-15 | 一种基于l-m算法的多波束横摇残差改正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112083429A true CN112083429A (zh) | 2020-12-15 |
CN112083429B CN112083429B (zh) | 2023-08-15 |
Family
ID=73736794
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010967899.XA Active CN112083429B (zh) | 2020-09-15 | 2020-09-15 | 一种基于l-m算法的多波束横摇残差改正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112083429B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113093159A (zh) * | 2021-03-01 | 2021-07-09 | 中国人民解放军海军大连舰艇学院 | 多波束测深误差改进模型设计方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2010148968A (ru) * | 2010-11-30 | 2012-06-10 | Владимир Максимович Смолин (RU) | Способ съемки рельефа дна акватории с движущегося судна многолучевым эхолотом с вертикальным зондированием гидроакустическими сигналами поверхности дна и многолучевой эхолот с вертикальным зондированием гидроакустическими сигналами поверхности дна акватории для его осуществления |
CN104181523A (zh) * | 2013-05-21 | 2014-12-03 | 中国科学院声学研究所 | 一种基于横摇稳定策略的多波束测深方法及系统 |
CN104459678A (zh) * | 2014-12-09 | 2015-03-25 | 中国海洋石油总公司 | 多波束测深系统实时横摇补偿的方法 |
CN106500671A (zh) * | 2016-09-22 | 2017-03-15 | 天津大学 | 一种基于lm算法分解激光雷达波形确定海水深度的方法 |
CN108957462A (zh) * | 2018-05-22 | 2018-12-07 | 中国海洋大学 | 一种基于平坦海底的多波束水体数据处理方法 |
CN110081864A (zh) * | 2019-04-24 | 2019-08-02 | 长江水利委员会水文局长江上游水文水资源勘测局 | 一种顾及水深值的水深测量综合延迟改正方法 |
-
2020
- 2020-09-15 CN CN202010967899.XA patent/CN112083429B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2010148968A (ru) * | 2010-11-30 | 2012-06-10 | Владимир Максимович Смолин (RU) | Способ съемки рельефа дна акватории с движущегося судна многолучевым эхолотом с вертикальным зондированием гидроакустическими сигналами поверхности дна и многолучевой эхолот с вертикальным зондированием гидроакустическими сигналами поверхности дна акватории для его осуществления |
CN104181523A (zh) * | 2013-05-21 | 2014-12-03 | 中国科学院声学研究所 | 一种基于横摇稳定策略的多波束测深方法及系统 |
CN104459678A (zh) * | 2014-12-09 | 2015-03-25 | 中国海洋石油总公司 | 多波束测深系统实时横摇补偿的方法 |
CN106500671A (zh) * | 2016-09-22 | 2017-03-15 | 天津大学 | 一种基于lm算法分解激光雷达波形确定海水深度的方法 |
CN108957462A (zh) * | 2018-05-22 | 2018-12-07 | 中国海洋大学 | 一种基于平坦海底的多波束水体数据处理方法 |
CN110081864A (zh) * | 2019-04-24 | 2019-08-02 | 长江水利委员会水文局长江上游水文水资源勘测局 | 一种顾及水深值的水深测量综合延迟改正方法 |
Non-Patent Citations (4)
Title |
---|
张志伟 等: "多波束换能器安装偏差对海底地形测量的影响", 《海洋测绘》, vol. 36, no. 1, pages 51 - 54 * |
张志伟 等: "多波束测点位置归算横纵摇耦合效应分析", 《地理与地理信息科学》, vol. 32, no. 4, pages 38 - 41 * |
白福成: "多波束测深系统运动补偿新技术研究与硬件设计", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》, no. 04 * |
阳凡林 等: "多波束勘测运动传感器偏移的改正方法", 《武汉大学学报 信息科学版》, vol. 35, no. 7, pages 816 - 820 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113093159A (zh) * | 2021-03-01 | 2021-07-09 | 中国人民解放军海军大连舰艇学院 | 多波束测深误差改进模型设计方法 |
CN113093159B (zh) * | 2021-03-01 | 2023-12-22 | 中国人民解放军海军大连舰艇学院 | 多波束测深误差改进模型设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112083429B (zh) | 2023-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110081864B (zh) | 一种顾及水深值的水深测量综合延迟改正方法 | |
CN111025273B (zh) | 一种畸变拖曳阵线谱特征增强方法及系统 | |
US9500484B2 (en) | System and method for water column aided navigation | |
CN107179535A (zh) | 一种基于畸变拖曳阵的保真增强波束形成的方法 | |
CN101587187B (zh) | 一种测深声纳系统偏差校正方法 | |
CN110061716B (zh) | 一种基于最小二乘和多重渐消因子的改进kalman滤波方法 | |
EP2098882A2 (en) | Location measurement method based on predictive filter | |
CN109765523A (zh) | 基于自适应akf的单应答器斜距水声定位方法及系统 | |
CN108562287A (zh) | 一种基于自适应采样粒子滤波的水下地形辅助导航方法 | |
CN111537982B (zh) | 一种畸变拖曳阵线谱特征增强方法及系统 | |
CN104181523A (zh) | 一种基于横摇稳定策略的多波束测深方法及系统 | |
CN110673148A (zh) | 一种主动声纳目标实时航迹解算方法 | |
CN113156413B (zh) | 一种基于双程声路径的海底基准校准方法 | |
CN110221278B (zh) | 一种基于多传感器组合的合成孔径声呐运动补偿方法 | |
CN105277932B (zh) | 一种基于下变频波束形成中的多普勒频移校正方法 | |
Mohammadloo et al. | Correcting multibeam echosounder bathymetric measurements for errors induced by inaccurate water column sound speeds | |
CN115307643A (zh) | 一种双应答器辅助的sins/usbl组合导航方法 | |
CN111220146B (zh) | 一种基于高斯过程回归学习的水下地形匹配定位方法 | |
CN112083429B (zh) | 一种基于l-m算法的多波束横摇残差改正方法 | |
JP2533287B2 (ja) | 線形アレイ横方向運動補正方法 | |
CN116047494A (zh) | 通感一体背景下基于多距离划分的多扩展目标的跟踪方法 | |
CN117146830A (zh) | 一种自适应多信标航位推算和长基线的紧组合导航方法 | |
Collin et al. | On spatial uncertainty in a surface long baseline positioning system | |
CN115184910B (zh) | 一种河道断面的单波束测量波束角效应的改正方法 | |
Dewey et al. | Reynolds stresses and turbulent kinetic energy estimates from various ADCP beam configurations: Theory |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |