CN102707289A - 基于航向角旋转的多波束侧扫声纳图像实时拼接方法 - Google Patents

基于航向角旋转的多波束侧扫声纳图像实时拼接方法 Download PDF

Info

Publication number
CN102707289A
CN102707289A CN2012101977202A CN201210197720A CN102707289A CN 102707289 A CN102707289 A CN 102707289A CN 2012101977202 A CN2012101977202 A CN 2012101977202A CN 201210197720 A CN201210197720 A CN 201210197720A CN 102707289 A CN102707289 A CN 102707289A
Authority
CN
China
Prior art keywords
course angle
slit
data
ping
beam side
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
Application number
CN2012101977202A
Other languages
English (en)
Other versions
CN102707289B (zh
Inventor
叶秀芬
李鹏
邓莹莹
希吉尔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201210197720.2A priority Critical patent/CN102707289B/zh
Publication of CN102707289A publication Critical patent/CN102707289A/zh
Application granted granted Critical
Publication of CN102707289B publication Critical patent/CN102707289B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明提供的是一种基于航向角旋转的多波束侧扫声纳图像实时拼接方法。采集多波束侧扫声纳所有通道中的原始数据,并设定拼接图像初始位置;移除因为第一个回波点返回之前换能器接收水下噪声产生的水柱区;将移除水柱区后的原始数据中的第一ping的N个通道数据写入画布,并记录当前的ping数和当前航向角;读取下一ping N个通道的数据和对应的GPS及航向角信息,计算前后两ping数据航向角的差值,若航向角无变化,则按GPS信息顺次拼接,若航向角有变化,将数据进行旋转后拼接,旋转角度为航向角差值;修补拼接后图像中出现的缝隙。本发明是一种精确度高、显示效果好、计算量低的多波束侧扫声纳图像实时拼接方法。

Description

基于航向角旋转的多波束侧扫声纳图像实时拼接方法
技术领域
本发明涉及的是一种图像拼接方法,具体地说是一种多波束水声图像实时拼接的新方法。
背景技术
广袤的海洋空间对于人类的生存和发展具有重要意义,随着人类对海洋的开发和利用,对海地地形地貌进行探测,绘制海底地形图,借以掌握详尽的海洋地理信息显得至关重要。当前海洋探测中,主要以声学探测系统为主要手段。但是船体在拖拽侧扫声纳探测过程中,不可避免地会受到各种因素的影响,由于声纳斜距成像的特点,GPS定位不准及异常跳变等不稳定因素存在,导致了侧扫声纳图像无法真实的显示海底地貌,从而无法获得精确的海底地形图像。
声纳图像拼接的目的,就是读取原始声纳数据中的信息,进行相应的处理,消除斜距的影响,并根据GPS数据将声纳文件中记录的信息还原成真实的海底地貌。当前,国外一些著名的声纳设备制造商都研发了相应的声纳成像软件,但价格昂贵,且其分析技术为商业机密无法公开,对于当前广泛使用的多波束侧扫声纳数据的拼图研究更是罕见,而现有的单波束声纳图像处理方法,逐一计算每个回波点坐标,将其应用于多波束声纳图像处理,积累误差大,计算时间长,无法进行精准的实时拼接,因此,对多波束侧扫声纳图像的实时拼接方法的研究具有重要的意义。
发明内容
本发明的目的在于提供一种精确度高、显示效果好、计算量低的基于航向角旋转的多波束侧扫声纳图像实时拼接方法。
本发明的目的是这样实现的:
(1)采集多波束侧扫声纳所有通道中的原始数据,并设定拼接图像初始位置;
(2)移除因为第一个回波点返回之前换能器接收水下噪声产生的水柱区;
(3)将移除水柱区后的原始数据中的第一ping的N个通道数据写入画布,并记录当前的ping数和当前航向角;
(4)读取下一ping N个通道的数据和对应的GPS及航向角信息,计算前后两ping数据航向角的差值,若航向角无变化,则按GPS信息顺次拼接,若航向角有变化,将数据进行旋转后拼接,旋转角度为航向角差值;
(5)修补拼接后图像中出现的缝隙。
本发明的方法可以进一步说明为:
(1)读取声纳文件中左右两侧N个通道的原始数据,并设定拼接图像的初始位置;
(2)移除原始数据中因为第一个回波点返回之前换能器接收水下噪声产生的水柱区;
(3)以拖鱼接收到的GPS坐标为中心,将移除水柱区后的原始数据中的第一ping的N个通道数据写入画布,并记录当前的ping数x和当前航向角αx
(4)记录[x+1]ping的N个通道数据,读取其航向角值α[x+1],计算前后两ping数据航向角差值Δα,若Δα=0,则以当前GPS坐标为中心顺次拼接,若Δα≠0,则以当前GPS坐标为中心进行旋转后拼接,其旋转矩阵为:
C ′ = cos Δα - sin Δα sin Δα cos Δα ;
(5)前后两ping数据拼接后,计算两个GPS中心点之间的距离,设两GPS中心分别为trace[x]、trace[x+1],则距离
Figure BDA00001771024700022
当航向角无变化时,设Δd=d-w,其中w为一ping数据宽度,若Δd≠0,则两ping数据间存在无效像素即缝隙,取阈值ε,若Δd<ε,则缝隙较小,采用最邻近插值,若Δd>ε,或航向角有变化时,则缝隙较大或存在三角形缝隙,采用双线性插值,填补缝隙中的无效像素。
本发明的特点是:步骤(1)中同时读入左右两侧共N个通道的回波数据,并在步骤(2)、(3)中将N个通道的数据作为整体进行处理,避免了逐一进行回波点坐标计算的积累误差,同时提高了计算速度。声纳文件中航向角记录有两个,一个来自拖鱼,一个来自船体,所述步骤(4)中航向角信息采用来自船体的记录,是由于拖鱼是在水中,受到潮汐、海流、鱼群等外界干扰比较大,因此其信息存在严重的误差及延迟,采用来自船体的航向角数据,使得拼接精度更高。
本发明与现有技术相比的优势在于:针对多波束侧扫声纳工作特点,同时对每ping的N个波束数据进行处理,无需逐一计算回波点坐标,减少了积累误差,降低了计算量,同时采用来自船体的航向角数据,该数据与拖鱼航向角数据相比,受到潮汐及洋流影响较小,因此使得拼图更加精确,针对拼图后的缝隙,结合最邻近插值和双线性插值填补拼接数据间的无效像素,提高了显示效果。
附图说明
图1是本发明方法的基本流程图;
图2拼接后未进行缝隙填补的图像;
图3拼接后利用插值修补缝隙后的图像。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
结合图1,本发明基于航向角旋转的多波束侧扫声纳图像实时拼接方法的具体步骤如下:
第一步是读取声纳文件中左右两侧N个通道的的原始信息,并设定拼接图像的初始位置;
第二步是移除原始数据中因为第一个回波点返回之前换能器接收水下噪声产生的水柱区;
第三步以拖鱼接收到的GPS坐标为中心,将移除水柱区后的原始声纳文件中的第一ping的N个通道数据作为整体写入画布,避免了逐一计算回波点所带来的积累误差,并记录当前的ping数x和当前航向角αx
第四步记录[x+1]ping的N个通道数据,读取其航向角值α[x+1],计算前后两ping数据航向角差值Δα,判断航向角是否变化,若Δα=0,则航向角不变,若Δα≠0,则航向角发生改变;
第五步根据航向角的变化情况对数据进行处理,若Δα≠0,则以当前GPS坐标为中心进行旋转,其旋转矩阵为:
C ′ = cos Δα - sin Δα sin Δα cos Δα ;
第六步是将旋转后的数据进行拼接,若Δα=0,则以当前GPS坐标为中心顺次拼接。其中,声纳文件中航向角记录有两个,一个来自拖鱼,一个来自船体,由于拖鱼是在水中,受到潮汐、海流、鱼群等外界干扰比较大,因此其信息存在严重的误差及延迟,因此,用于旋转的航向角信息采用来自船体的记录,使得拼接精度更高;
第七步是前后两ping数据拼接后,计算两个GPS中心点之间的距离,设两GPS中心分别为trace[x]、trace[x+1],则距离
Figure BDA00001771024700032
当航向角无变化时,设Δd=d-w,其中w为一ping数据宽度,若Δd≠0,则两ping数据间存在无效像素即缝隙,若Δd=0,则两ping数据间不存在缝隙;
第八步是插值处理,取阈值ε,若Δd<ε,则缝隙较小,采用最邻近插值,若Δd>ε,或航向角有变化时,则缝隙较大或存在三角形缝隙,采用双线性插值,填补缝隙中的无效像素。

Claims (9)

1.一种基于航向角旋转的多波束侧扫声纳图像实时拼接方法,其特征是:
(1)采集多波束侧扫声纳所有通道中的原始数据,并设定拼接图像初始位置;
(2)移除因为第一个回波点返回之前换能器接收水下噪声产生的水柱区;
(3)将移除水柱区后的原始数据中的第一ping的N个通道数据写入画布,并记录当前的ping数和当前航向角;
(4)读取下一ping N个通道的数据和对应的GPS及航向角信息,计算前后两ping数据航向角的差值,若航向角无变化,则按GPS信息顺次拼接,若航向角有变化,将数据进行旋转后拼接,旋转角度为航向角差值;
(5)修补拼接后图像中出现的缝隙。
2.根据权利要求1所述的基于航向角旋转的多波束侧扫声纳图像实时拼接方法,其特征是:所述采集多波束侧扫声纳所有通道中的原始数据是读取声纳文件中左右两侧N个通道的原始数据。
3.根据权利要求1或2所述的基于航向角旋转的多波束侧扫声纳图像实时拼接方法,其特征是:所述将移除水柱区后的原始数据中的第一ping的N个通道数据写入画布是以拖鱼接收到的GPS坐标为中心。
4.根据权利要求1或2所述的基于航向角旋转的多波束侧扫声纳图像实时拼接方法,其特征是:所述将数据进行旋转后拼接的旋转矩阵为:
C ′ = cos Δα - sin Δα sin Δα cos Δα ,
其中Δα为前后两ping数据航向角差值。
5.根据权利要求3所述的基于航向角旋转的多波束侧扫声纳图像实时拼接方法,其特征是:所述将数据进行旋转后拼接的旋转矩阵为:
C ′ = cos Δα - sin Δα sin Δα cos Δα ,
其中Δα为前后两ping数据航向角差值。
6.根据权利要求1或2所述的基于航向角旋转的多波束侧扫声纳图像实时拼接方法,其特征是所述修补拼接后图像中出现的缝隙的方法为:前后两ping数据拼接后,计算两个GPS中心点之间的距离,设两GPS中心分别为trace[x]、trace[x+1],则距离
Figure FDA00001771024600021
当航向角无变化时,设Δd=d-w,其中w为一ping数据宽度,若Δd≠0,则两ping数据间存在无效像素即缝隙,取阈值ε,若Δd<ε,则缝隙较小,采用最邻近插值,若Δd>ε,或航向角有变化时,则缝隙较大或存在三角形缝隙,采用双线性插值,填补缝隙中的无效像素。
7.根据权利要求3所述的基于航向角旋转的多波束侧扫声纳图像实时拼接方法,其特征是所述修补拼接后图像中出现的缝隙的方法为:前后两ping数据拼接后,计算两个GPS中心点之间的距离,设两GPS中心分别为trace[x]、trace[x+1],则距离
Figure FDA00001771024600022
当航向角无变化时,设Δd=d-w,其中w为一ping数据宽度,若Δd≠0,则两ping数据间存在无效像素即缝隙,取阈值ε,若Δd<ε,则缝隙较小,采用最邻近插值,若Δd>ε,或航向角有变化时,则缝隙较大或存在三角形缝隙,采用双线性插值,填补缝隙中的无效像素。
8.根据权利要求4所述的基于航向角旋转的多波束侧扫声纳图像实时拼接方法,其特征是所述修补拼接后图像中出现的缝隙的方法为:前后两ping数据拼接后,计算两个GPS中心点之间的距离,设两GPS中心分别为trace[x]、trace[x+1],则距离
Figure FDA00001771024600023
当航向角无变化时,设Δd=d-w,其中w为一ping数据宽度,若Δd≠0,则两ping数据间存在无效像素即缝隙,取阈值ε,若Δd<ε,则缝隙较小,采用最邻近插值,若Δd>ε,或航向角有变化时,则缝隙较大或存在三角形缝隙,采用双线性插值,填补缝隙中的无效像素。
9.根据权利要求5所述的基于航向角旋转的多波束侧扫声纳图像实时拼接方法,其特征是所述修补拼接后图像中出现的缝隙的方法为:前后两ping数据拼接后,计算两个GPS中心点之间的距离,设两GPS中心分别为trace[x]、trace[x+1],则距离
Figure FDA00001771024600024
当航向角无变化时,设Δd=d-w,其中w为一ping数据宽度,若Δd≠0,则两ping数据间存在无效像素即缝隙,取阈值ε,若Δd<ε,则缝隙较小,采用最邻近插值,若Δd>ε,或航向角有变化时,则缝隙较大或存在三角形缝隙,采用双线性插值,填补缝隙中的无效像素。
CN201210197720.2A 2012-06-15 2012-06-15 基于航向角旋转的多波束侧扫声纳图像实时拼接方法 Expired - Fee Related CN102707289B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210197720.2A CN102707289B (zh) 2012-06-15 2012-06-15 基于航向角旋转的多波束侧扫声纳图像实时拼接方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210197720.2A CN102707289B (zh) 2012-06-15 2012-06-15 基于航向角旋转的多波束侧扫声纳图像实时拼接方法

Publications (2)

Publication Number Publication Date
CN102707289A true CN102707289A (zh) 2012-10-03
CN102707289B CN102707289B (zh) 2016-03-09

Family

ID=46900219

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210197720.2A Expired - Fee Related CN102707289B (zh) 2012-06-15 2012-06-15 基于航向角旋转的多波束侧扫声纳图像实时拼接方法

Country Status (1)

Country Link
CN (1) CN102707289B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103292792A (zh) * 2013-04-26 2013-09-11 国家海洋局第二海洋研究所 一种适用海底探测与假地形处理的实测svp重构方法
CN103364787A (zh) * 2013-03-28 2013-10-23 哈尔滨工程大学 一种多波束侧扫声纳图像镶嵌缝隙消除方法
CN105372663A (zh) * 2015-12-01 2016-03-02 宁波工程学院 一种面向侧扫声纳图像的重采样方法
CN105389777A (zh) * 2015-10-23 2016-03-09 首都师范大学 无人机序列影像快速无缝拼接系统
CN105869193A (zh) * 2016-03-31 2016-08-17 哈尔滨工程大学 基于uuv的侧扫声纳图像辅助判读方法
CN104536005B (zh) * 2014-11-17 2017-04-12 哈尔滨工程大学 一种基于盲区校正的多波束侧扫声呐斜距失真消除方法
CN115469315A (zh) * 2022-11-02 2022-12-13 广东智能无人系统研究院 一种采用倾斜侧扫模式的声纳成像方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5699318A (en) * 1996-07-25 1997-12-16 Northrop Grumman Corporation Topographic composition sonar map
CN101650426A (zh) * 2008-12-26 2010-02-17 中国科学院声学研究所 一种合成孔径声纳图像数据拼接系统及方法
CN102012513A (zh) * 2010-06-29 2011-04-13 上海大学 前视声纳图像的实时拼接方法和系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5699318A (en) * 1996-07-25 1997-12-16 Northrop Grumman Corporation Topographic composition sonar map
CN101650426A (zh) * 2008-12-26 2010-02-17 中国科学院声学研究所 一种合成孔径声纳图像数据拼接系统及方法
CN102012513A (zh) * 2010-06-29 2011-04-13 上海大学 前视声纳图像的实时拼接方法和系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
H. PAUL JOHNSON ET AL.: "The geological interpretation of side-scan sonar", 《REVIEWS OF GEOPHYSICS》, vol. 28, no. 4, 30 November 1990 (1990-11-30), pages 357 - 380 *
YET-CHUNG CHANG ET AL.: "SIDESCAN SONAR IMAGE PROCESSING: CORRECTING BRIGHTNESS VARIATION AND PATCHING GAPS", 《JOURNAL OF MARINE SCIENCE AND TECHNOLOGY》, vol. 18, no. 6, 31 December 2010 (2010-12-31), pages 785 - 789 *
李廷桓: "论旁侧声纳图像镶嵌技术", 《海洋与湖沼》, vol. 22, no. 6, 30 November 1991 (1991-11-30) *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103364787A (zh) * 2013-03-28 2013-10-23 哈尔滨工程大学 一种多波束侧扫声纳图像镶嵌缝隙消除方法
CN103364787B (zh) * 2013-03-28 2017-11-28 哈尔滨工程大学 一种多波束侧扫声纳图像镶嵌缝隙消除方法
CN103292792A (zh) * 2013-04-26 2013-09-11 国家海洋局第二海洋研究所 一种适用海底探测与假地形处理的实测svp重构方法
CN104536005B (zh) * 2014-11-17 2017-04-12 哈尔滨工程大学 一种基于盲区校正的多波束侧扫声呐斜距失真消除方法
CN105389777A (zh) * 2015-10-23 2016-03-09 首都师范大学 无人机序列影像快速无缝拼接系统
CN105372663A (zh) * 2015-12-01 2016-03-02 宁波工程学院 一种面向侧扫声纳图像的重采样方法
CN105869193A (zh) * 2016-03-31 2016-08-17 哈尔滨工程大学 基于uuv的侧扫声纳图像辅助判读方法
CN105869193B (zh) * 2016-03-31 2018-09-28 哈尔滨工程大学 基于uuv的侧扫声纳图像辅助判读方法
CN115469315A (zh) * 2022-11-02 2022-12-13 广东智能无人系统研究院 一种采用倾斜侧扫模式的声纳成像方法
CN115469315B (zh) * 2022-11-02 2023-03-10 广东智能无人系统研究院 一种采用倾斜侧扫模式的声纳成像方法

Also Published As

Publication number Publication date
CN102707289B (zh) 2016-03-09

Similar Documents

Publication Publication Date Title
CN102707289A (zh) 基于航向角旋转的多波束侧扫声纳图像实时拼接方法
US9651698B2 (en) Multi-beam bathymetric chart construction method based on submarine digital depth model feature extraction
US6084827A (en) Dual-head multibeam sonar apparatus and method for tracking objects underwater
CN106199562B (zh) 基于机载激光雷达测量海底地形的海面误差校正方法
CN101334278B (zh) 一种数码遥感地质测绘的方法及装置
Firing et al. Equatorial subthermocline currents across the Pacific
Menna et al. A photogrammetric approach to survey floating and semi-submerged objects
CN105352476A (zh) 船载水岸线水上水下一体化测量系统集成方法
CN105258684A (zh) 基于激光点云为约束的多波束低掠射角波束归位方法
CN108089155A (zh) 一种深海环境下单水听器声源被动定位方法
Pydyn et al. Exploration and reconstruction of a medieval harbour using hydroacoustics, 3‐D shallow seismic and underwater photogrammetry: A case study from Puck, southern Baltic Sea
CN103729846A (zh) 基于不规则三角网的LiDAR点云数据边缘检测方法
CN103592650A (zh) 基于图形处理器的三维声纳成像系统及其三维成像方法
CN113640808B (zh) 浅水区海底电缆埋深探测方法及装置
CN110926434A (zh) 一种浅地层剖面测量系统软件架构及实现方法
CN101592731B (zh) 一种基于航迹线文件的侧扫声纳拖鱼航迹处理方法
CN104515499A (zh) 新型内河航道断面测量系统及其测量方法
CN101290223A (zh) 海底地形测量波束角效应改正方法
Leng et al. Satellite derived active-passive fusion bathymetry based on GRU model
White Utilization of LIDAR and NOAA's vertical datum transformation tool (VDatum) for shoreline delineation
Parrott et al. Integration of multibeam bathymetry and LiDAR surveys of the Bay of Fundy, Canada
Gaglianone et al. Investigating submerged morphologies by means of the low-budget “GeoDive” method (high resolution for detailed 3D reconstruction and related measurements)
CN107064942A (zh) 一种高精度多普勒信号模拟方法
Ludvigsen et al. Photogrammetric models for marine archaeology
CN103364787B (zh) 一种多波束侧扫声纳图像镶嵌缝隙消除方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160309

Termination date: 20210615