CN113376639B - 基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法 - Google Patents

基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法 Download PDF

Info

Publication number
CN113376639B
CN113376639B CN202110813425.4A CN202110813425A CN113376639B CN 113376639 B CN113376639 B CN 113376639B CN 202110813425 A CN202110813425 A CN 202110813425A CN 113376639 B CN113376639 B CN 113376639B
Authority
CN
China
Prior art keywords
terrain
sonar
point
imaging
dimensional
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.)
Active
Application number
CN202110813425.4A
Other languages
English (en)
Other versions
CN113376639A (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.)
Fuzhou University
Original Assignee
Fuzhou 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 Fuzhou University filed Critical Fuzhou University
Priority to CN202110813425.4A priority Critical patent/CN113376639B/zh
Publication of CN113376639A publication Critical patent/CN113376639A/zh
Priority to US17/700,293 priority patent/US11841426B2/en
Application granted granted Critical
Publication of CN113376639B publication Critical patent/CN113376639B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8993Three dimensional imaging systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C15/00Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • G01S7/4004Means for monitoring or calibrating of parts of a radar system
    • G01S7/4026Antenna boresight
    • G01S7/403Antenna boresight in azimuth, i.e. in the horizontal plane
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/523Details of pulse systems
    • G01S7/526Receivers
    • G01S7/53Means for transforming coordinates or for evaluating data, e.g. using computers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/539Details 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/56Display arrangements
    • G01S7/62Cathode-ray tube displays
    • G01S7/6245Stereoscopic displays; Three-dimensional displays; Pseudo-three dimensional displays
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2215/00Indexing scheme for image rendering
    • G06T2215/16Using real world measurements to influence rendering

Abstract

本发明提出基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法,包括以下步骤;步骤S1、利用声呐从不同方位对桥墩基础冲刷区地形进行完整扫描成像,采集n组基础冲刷地形声像图;步骤S2、对采集的声像图以相同间距截取多个分析截面;步骤S3、提取图中分析截面内地形成像条带上下边缘的关键参数,转换至三维空间中,并以扇弧表述声呐波束扇面;步骤S4、根据同一分析截面内一组测点扇弧相对位置关系,识别其中的冲刷地形轮廓形态;步骤S5、对各地形轮廓形态进行逐一识别,分别提取出空间散点三维坐标数据;步骤S6、采用克里金插值法对空间散点数据进行插值拟合,完成基础冲刷地形的三维重构;本发明能实现对桥墩基础冲刷区地形的三维还原的目的。

Description

基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法
技术领域
本发明涉及水下桥墩基础冲刷地形检测技术领域,尤其是基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法。
背景技术
桥墩基础长期处于复杂的水环境中,在冲刷作用下使河床不断下切,基础埋深减小,形成基础冲刷病害。当遇到洪水期间冲刷作用更加强烈,老旧桥梁受到长期局部冲刷作用下,存在重大安全隐患,在洪水剧烈冲刷下极易发生倾覆。因此,对桥墩基础冲刷进行定期检测,掌握桥墩基础冲刷状况,可以为桥梁工程的结构状态安全评估、冲刷范围河床的预防处理等提供信息和支持。
传统的基础冲刷检测装置安装费用及维护成本高,在洪水期间容易受到杂物影响而使装置失效,且测量的数据较为离散,检测效率低。声呐设备相比于传统的检测装置更易于安装,单次对基础冲刷扫描成像可以获取多个冲刷深度数据,检测效率高,且可以通过移动测点位置对桥墩不同方位的基础冲刷地形进行探测,从而获得不同方位基础冲刷地形的侧向轮廓,但是采集的声学图像只能反映波束覆盖某一区域地形冲刷轮廓,难以完整还原桥墩基础周围冲刷整体的地形地貌信息。
因此,需要一种基于机械扫描声呐成像的桥墩基础冲刷区地形三维重构方法。
发明内容
本发明提出基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法,能实现对桥墩基础冲刷区地形的三维还原,达到从三维空间中直接查看基础冲刷地形地貌的目的。
本发明采用以下技术方案。
基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法,所述方法包括以下步骤;
步骤S1、利用机械扫描声呐设备从不同方位对桥墩基础冲刷区地形进行完整扫描成像,采集得到不同方位下的n组基础冲刷地形声像图;
步骤S2、对采集的各组声像图以相同间距截取多个分析截面;
步骤S3、提取图中分析截面内地形成像条带上下边缘的关键参数,根据声学成像原理将提取参数转换至三维空间中,并以扇弧表述声呐波束扇面;
步骤S4、根据同一分析截面内一组测点扇弧相对位置关系,识别该分析截面内的冲刷地形轮廓形态;
步骤S5、对不同方位、不同分析截面的地形轮廓形态进行逐一识别,分别提取出轮廓曲线上关键点三维坐标,得到所需的空间散点三维坐标数据;
步骤S6、采用克里金插值法对空间散点数据进行插值拟合,完成基础冲刷地形的三维重构。
所述步骤S3中,提取声像图中分析截面内地形成像条带上下边缘的关键参数,根据声学成像原理将提取参数转换至三维空间的具体步骤如下:
步骤3a、在距离声呐中心为
Figure BDA0003169315390000021
的分析截面
Figure BDA0003169315390000022
内,标记成像条带的上下边缘点为M、N,原始声呐中心点
Figure BDA0003169315390000023
和上下边缘点M、N距离分别为
Figure BDA0003169315390000024
波束扇面的大开角角度θ;
步骤3b、以上边缘点M为圆心,
Figure BDA0003169315390000025
为半径作圆弧,和分析截面
Figure BDA0003169315390000026
交点即为等效声呐中心点
Figure BDA0003169315390000027
等效声呐中心点
Figure BDA0003169315390000028
超出原始声呐中心所在平面的高度
Figure BDA0003169315390000029
以超出声呐所在平面距离
Figure BDA00031693153900000210
的等效声呐中心点
Figure BDA00031693153900000211
为圆心,以
Figure BDA00031693153900000212
为半径、波束扇面大开角角度θ作扇弧,得到上边缘点M在三维空间对应的扇弧
Figure BDA00031693153900000213
以下边缘点N为圆心,
Figure BDA00031693153900000214
为半径作圆弧,和分析截面
Figure BDA00031693153900000215
交点即为等效声呐中心点
Figure BDA00031693153900000216
高于原始声呐中心所在平面距离
Figure BDA00031693153900000217
Figure BDA00031693153900000218
为圆心,以
Figure BDA00031693153900000219
为半径、角度θ作扇弧,得到下边缘点N在三维空间对应于扇弧
Figure BDA00031693153900000220
其中,上标i表示桩号,第一个下标j表示一组测点相对桥墩的所在方位,第二个下标k表示一组测点下的具体测点,k取0、1、2分别表示一组测点中的左、中、右测点;
步骤3c、对一组测点中另外两个测点也采用相同的参数转换方式,得到同一分析截面内一组测点成像图对应的三维空间中的扇弧相对位置关系。
步骤S4中,根据同一分析截面内参数转换得到的一组测点对应的扇弧相对位置关系识别冲刷区地形轮廓形态时,具体步骤如下:
步骤4a、根据一组扇弧相对位置关系判定地形类型属于凹陷地形、隆起地形、连续地形、不确定地形中的哪一种,其中不确定地形需要结合相邻方位下确定的地形轮廓曲线对地形类型进行进一步判断;
步骤4b、筛选可用于形成地形轮廓曲线的关键点,按以下条件进行筛选,1)地形类型属于凹陷地形、隆起地形或连续地形;2)拟合的地形轮廓曲线位于两扇弧范围之间,且相邻扇面内交集部分为地形轮廓曲线必经过区域;3)各测点还原的上下扇弧线上均至少有一个地形轮廓点;
步骤4c、在满足上述条件下确定出各个扇弧线上地形轮廓曲线各个关键点,通过平滑曲线将各个关键点进行连接,得到该分析截面下波束覆盖区的地形轮廓曲线。
在步骤4a中,若在比较两侧辅助测点波束上侧扇弧边缘点和声呐所在平面距离时,距离较大值小于另一侧辅助检测点下侧扇弧边缘点到声呐所在平面的距离,即用公式表达为L02>l00>l02;且两侧波束上侧扇弧中点连线a00a02,中间波束上侧扇弧中点a01位于连线a00a02的下方,即可确定地形类型为凹陷地形类型。
所述步骤S6中,将得到的大量的空间散点数据通过克里金插值法进行插值拟合处理,实现桥墩基础冲刷区地形三维重构,得到基础冲刷等高线图和三维可视化图,能在三维空间中直观查看桥墩基础冲刷区的地形地貌信息。
本发明的有益效果在于:本发明所述方案,能通过声呐扫描成像,直观反映出桥墩基础周围冲刷区地形地貌信息,便于在三维空间中直接查看桥墩基础冲刷状况。
附图说明
下面结合附图和具体实施方式对本发明进一步详细的说明:
附图1是本发明的流程示意图;
附图2是本发明所述方法中,声呐桥墩基础冲刷地形成像分析截面截取的示意图;
附图3是本发明所述方法中,声像图提取参数转换至三维空间的示意图;
附图4是按本发明所述方法,根据一组测点下的扇弧相对位置关系来识别凹陷地形轮廓的示意图。
具体实施方式
如图所示,基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法,所述方法包括以下步骤;
步骤S1、利用机械扫描声呐设备从不同方位对桥墩基础冲刷区地形进行完整扫描成像,采集得到不同方位下的n组基础冲刷地形声像图;
步骤S2、对采集的各组声像图以相同间距截取多个分析截面;
步骤S3、提取图中分析截面内地形成像条带上下边缘的关键参数,根据声学成像原理将提取参数转换至三维空间中,并以扇弧表述声呐波束扇面;
步骤S4、根据同一分析截面内一组测点扇弧相对位置关系,识别该分析截面内的冲刷地形轮廓形态;
步骤S5、对不同方位、不同分析截面的地形轮廓形态进行逐一识别,分别提取出轮廓曲线上关键点三维坐标,得到所需的空间散点三维坐标数据;
步骤S6、采用克里金插值法对空间散点数据进行插值拟合,完成基础冲刷地形的三维重构。
所述步骤S3中,提取声像图中分析截面内地形成像条带上下边缘的关键参数,根据声学成像原理将提取参数转换至三维空间的具体步骤如下:
步骤3a如图3中的a所示,、在距离声呐中心为
Figure BDA0003169315390000041
的分析截面
Figure BDA0003169315390000042
内,标记成像条带的上下边缘点为M、N,原始声呐中心点
Figure BDA0003169315390000043
和上下边缘点M、N距离分别为
Figure BDA0003169315390000044
Figure BDA0003169315390000045
波束扇面的大开角角度θ;
步骤3b、以上边缘点M为圆心,
Figure BDA0003169315390000046
为半径作圆弧,和分析截面
Figure BDA0003169315390000047
交点即为等效声呐中心点
Figure BDA0003169315390000048
等效声呐中心点
Figure BDA0003169315390000049
超出原始声呐中心所在平面的高度
Figure BDA00031693153900000410
以超出声呐所在平面距离
Figure BDA00031693153900000411
的等效声呐中心点
Figure BDA00031693153900000412
为圆心,以
Figure BDA00031693153900000413
为半径、波束扇面大开角角度θ作扇弧,得到上边缘点M在三维空间对应的扇弧
Figure BDA0003169315390000051
以下边缘点N为圆心,
Figure BDA0003169315390000052
为半径作圆弧,和分析截面
Figure BDA0003169315390000053
交点即为等效声呐中心点
Figure BDA0003169315390000054
高于原始声呐中心所在平面距离
Figure BDA0003169315390000055
Figure BDA0003169315390000056
为圆心,以
Figure BDA0003169315390000057
为半径、角度θ作扇弧,得到下边缘点N在三维空间对应于扇弧
Figure BDA0003169315390000058
其中,上标i表示桩号,第一个下标j表示一组测点所在方位,第二个下标k表示一组测点下的具体测点,k取0、1、2分别表示一组测点中的左、中、右测点;
以图为例:步骤3b以图3中的b、c展示为:原始声呐中心点
Figure BDA0003169315390000059
和该分析截面的上下边缘点M、N距离分别为
Figure BDA00031693153900000510
以上边缘点M为圆心,
Figure BDA00031693153900000511
为半径作圆弧,和分析截面
Figure BDA00031693153900000512
交点即为等效声呐中心点
Figure BDA00031693153900000513
等效声呐中心点
Figure BDA00031693153900000514
超出原始声呐中心所在平面的高度
Figure BDA00031693153900000515
以超出声呐所在平面距离
Figure BDA00031693153900000516
的等效声呐中心点
Figure BDA00031693153900000517
为圆心,以
Figure BDA00031693153900000518
为半径、波束扇面大开角角度θ作扇弧,得到图三(c)上边缘点M在三维空间对应于扇弧
Figure BDA00031693153900000519
同样地,以下边缘点N为圆心,
Figure BDA00031693153900000520
为半径作圆弧,和分析截面
Figure BDA00031693153900000521
交点即为等效声呐中心点
Figure BDA00031693153900000522
高于原始声呐中心所在平面距离
Figure BDA00031693153900000523
Figure BDA00031693153900000524
为圆心,以
Figure BDA00031693153900000525
为半径、角度θ作扇弧,得到下边缘点N在三维空间对应于扇弧
Figure BDA00031693153900000526
其中,上标i表示桩号,第一个下标j表示一组测点相对桥墩的方位,第二个下标k表示一组测点下的具体测点,k取0、1、2分别表示一组测点中的左、中、右测点。
步骤3c、对一组测点中另外两个测点也采用相同的参数转换方式,得到同一分析截面内一组测点成像图对应的三维空间中的扇弧相对位置关系。
步骤S4中,根据同一分析截面内参数转换得到的一组测点对应的扇弧相对位置关系识别冲刷区地形轮廓形态时,具体步骤如下:
步骤4a、根据一组扇弧相对位置关系判定地形类型属于凹陷地形、隆起地形、连续地形、不确定地形中的哪一种,其中不确定地形需要结合相邻方位下确定的地形轮廓曲线对地形类型进行进一步判断;
步骤4b、筛选可用于形成地形轮廓曲线的关键点,按以下条件进行筛选,1)地形类型属于凹陷地形、隆起地形或连续地形;2)拟合的地形轮廓曲线位于两扇弧范围之间,且相邻扇面内交集部分为地形轮廓曲线必经过区域;3)各测点还原的上下扇弧线上均至少有一个地形轮廓点;
步骤4c、在满足上述条件下确定出各个扇弧线上地形轮廓曲线各个关键点,通过平滑曲线将各个关键点进行连接,得到该分析截面下波束覆盖区的地形轮廓曲线。
如图3所示,在步骤4a中,若在比较两侧辅助测点波束上侧扇弧边缘点和声呐所在平面距离时,距离较大值小于另一侧辅助检测点下侧扇弧边缘点到声呐所在平面的距离,即用公式表达为L02>l00>l02;且两侧波束上侧扇弧中点连线a00a02,中间波束上侧扇弧中点a01位于连线a00a02的下方,即可确定地形类型为凹陷地形类型。
如图4所示,当满足地型类型确定为凹陷地形类型时,拟合的地形轮廓曲线位于两扇弧范围之间,且相邻扇弧面内交集部分,即图4中A01C01D00c00和c02C02D01B01区域为地形轮廓曲线必经过区域;则各测点还原的上下扇弧线上均至少有一个地形轮廓点。
所述步骤S6中,将得到的大量的空间散点数据通过克里金插值法进行插值拟合处理,实现桥墩基础冲刷区地形三维重构,得到基础冲刷等高线图和三维可视化图,能在三维空间中直观查看桥墩基础冲刷区的地形地貌信息。
以上所述,仅为本发明较佳的具体实施方法,这些具体实施方法是基于整体构思下的不同实现方式,而本发明的保护范围并不仅限于此。仅用于解释说明本发明的功能及使用方法,而并非对于本发明的限制,本技术领域的普通技术人员在本发明的实质范围内,做出各种变化或替代,也应属于本发明的保护范畴。

Claims (4)

1.基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法,其特征在于:所述方法包括以下步骤;
步骤S1、利用机械扫描声呐设备从不同方位对桥墩基础冲刷区地形进行完整扫描成像,采集得到不同方位下的n组基础冲刷地形声像图;
步骤S2、对采集的各组声像图以相同间距截取多个分析截面;
步骤S3、提取图中分析截面内地形成像条带上下边缘的关键参数,根据声学成像原理将提取参数转换至三维空间中,并以扇弧表述声呐波束扇面;
步骤S4、根据同一分析截面内一组测点扇弧相对位置关系,识别该分析截面内的冲刷地形轮廓形态;
步骤S5、对不同方位、不同分析截面的地形轮廓形态进行逐一识别,分别提取出轮廓曲线上关键点三维坐标,得到所需的空间散点三维坐标数据;
步骤S6、采用克里金插值法对空间散点数据进行插值拟合,完成基础冲刷地形的三维重构;
所述步骤S3中,提取声像图中分析截面内地形成像条带上下边缘的关键参数,根据声学成像原理将提取参数转换至三维空间的具体步骤如下:
步骤3a、在距离声呐中心为
Figure FDA0003641593730000011
的分析截面
Figure FDA0003641593730000012
内,标记成像条带的上下边缘点为M、N,原始声呐中心点
Figure FDA0003641593730000013
和上下边缘点M、N距离分别为
Figure FDA0003641593730000014
波束扇面的大开角角度θ;
步骤3b、以上边缘点M为圆心,
Figure FDA0003641593730000015
为半径作圆弧,和分析截面
Figure FDA0003641593730000016
交点即为等效声呐中心点
Figure FDA0003641593730000017
等效声呐中心点
Figure FDA0003641593730000018
超出原始声呐中心所在平面的高度
Figure FDA0003641593730000019
以超出声呐所在平面距离
Figure FDA00036415937300000110
的等效声呐中心点
Figure FDA00036415937300000111
为圆心,以
Figure FDA00036415937300000112
为半径、波束扇面大开角角度θ作扇弧,得到上边缘点M在三维空间对应的扇弧
Figure FDA00036415937300000113
以下边缘点N为圆心,
Figure FDA00036415937300000114
为半径作圆弧,和分析截面
Figure FDA00036415937300000115
交点即为等效声呐中心点
Figure FDA00036415937300000116
高于原始声呐中心所在平面距离
Figure FDA00036415937300000117
Figure FDA00036415937300000118
为圆心,以
Figure FDA00036415937300000119
为半径、角度θ作扇弧,得到下边缘点N在三维空间对应于扇弧
Figure FDA0003641593730000021
其中,上标i表示桩号,第一个下标j表示一组测点相对桥墩的所在方位,第二个下标k表示一组测点下的具体测点,k取0、1、2分别表示一组测点中的左、中、右测点;
步骤3c、对一组测点中另外两个测点也采用相同的参数转换方式,得到同一分析截面内一组测点成像图对应的三维空间中的扇弧相对位置关系。
2.根据权利要求1所述的基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法,其特征在于:步骤S4中,根据同一分析截面内参数转换得到的一组测点对应的扇弧相对位置关系识别冲刷区地形轮廓形态时,具体步骤如下:
步骤4a、根据一组扇弧相对位置关系判定地形类型属于凹陷地形、隆起地形、连续地形、不确定地形中的哪一种,其中不确定地形需要结合相邻方位下确定的地形轮廓曲线对地形类型进行进一步判断;
步骤4b、筛选可用于形成地形轮廓曲线的关键点,按以下条件进行筛选,1)地形类型属于凹陷地形、隆起地形或连续地形;2)拟合的地形轮廓曲线位于两扇弧范围之间,且相邻扇面内交集部分为地形轮廓曲线必经过区域;3)各测点还原的上下扇弧线上均至少有一个地形轮廓点;
步骤4c、在满足上述条件下确定出各个扇弧线上地形轮廓曲线各个关键点,通过平滑曲线将各个关键点进行连接,得到该分析截面下波束覆盖区的地形轮廓曲线。
3.根据权利要求2所述的基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法,其特征在于:在步骤4a中,若在比较两侧辅助测点波束上侧扇弧边缘点和声呐所在平面距离时,距离较大值小于另一侧辅助检测点下侧扇弧边缘点到声呐所在平面的距离,即用公式表达为L02>l00>l02;且两侧波束上侧扇弧中点连线a00a02,中间波束上侧扇弧中点a01位于连线a00a02的下方,即可确定地形类型为凹陷地形类型。
4.根据权利要求1所述的基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法,其特征在于:所述步骤S6中,将得到的大量的空间散点数据通过克里金插值法进行插值拟合处理,实现桥墩基础冲刷区地形三维重构,得到基础冲刷等高线图和三维可视化图,能在三维空间中直观查看桥墩基础冲刷区的地形地貌信息。
CN202110813425.4A 2021-07-19 2021-07-19 基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法 Active CN113376639B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202110813425.4A CN113376639B (zh) 2021-07-19 2021-07-19 基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法
US17/700,293 US11841426B2 (en) 2021-07-19 2022-03-21 Three-dimensional (3D) terrain reconstruction method for scoured area around bridge pier foundation based on mechanical scanned imaging sonar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110813425.4A CN113376639B (zh) 2021-07-19 2021-07-19 基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法

Publications (2)

Publication Number Publication Date
CN113376639A CN113376639A (zh) 2021-09-10
CN113376639B true CN113376639B (zh) 2022-08-05

Family

ID=77582490

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110813425.4A Active CN113376639B (zh) 2021-07-19 2021-07-19 基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法

Country Status (2)

Country Link
US (1) US11841426B2 (zh)
CN (1) CN113376639B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113376639B (zh) * 2021-07-19 2022-08-05 福州大学 基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法
CN116844142B (zh) * 2023-08-28 2023-11-21 四川华腾公路试验检测有限责任公司 一种桥梁基础冲刷的识别和评估方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101290237A (zh) * 2008-01-18 2008-10-22 天津大学 建筑物理数据快速三维采样方法
CN102636784A (zh) * 2012-03-31 2012-08-15 上海大学 单次声波入射快速预报声呐目标散射指向分布图的方法
CN106707265A (zh) * 2016-12-26 2017-05-24 江苏中海达海洋信息技术有限公司 一种声呐扇形图的转换更新方法
KR101802972B1 (ko) * 2016-06-29 2017-12-01 이너스페이스원정 주식회사 유인 수중건설 작업 로봇 시스템
CN107632305A (zh) * 2017-09-11 2018-01-26 河海大学 一种基于剖面声纳扫测技术的海底局部地形自主感知方法及装置
CN109765562A (zh) * 2018-12-10 2019-05-17 中国科学院声学研究所 一种三维前视声像声纳系统和方法
CN111579639A (zh) * 2020-05-19 2020-08-25 福州大学 一种成像声呐扫描桥梁水下桩柱缺陷的测点布设方法
CN112198516A (zh) * 2020-11-10 2021-01-08 广东智芯光电科技有限公司 一种水下图像采集的方法和系统
CN112285682A (zh) * 2020-10-20 2021-01-29 水利部交通运输部国家能源局南京水利科学研究院 水工工程涵洞环境的360°多波束声呐扫描装置及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113376639B (zh) * 2021-07-19 2022-08-05 福州大学 基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101290237A (zh) * 2008-01-18 2008-10-22 天津大学 建筑物理数据快速三维采样方法
CN102636784A (zh) * 2012-03-31 2012-08-15 上海大学 单次声波入射快速预报声呐目标散射指向分布图的方法
KR101802972B1 (ko) * 2016-06-29 2017-12-01 이너스페이스원정 주식회사 유인 수중건설 작업 로봇 시스템
CN106707265A (zh) * 2016-12-26 2017-05-24 江苏中海达海洋信息技术有限公司 一种声呐扇形图的转换更新方法
CN107632305A (zh) * 2017-09-11 2018-01-26 河海大学 一种基于剖面声纳扫测技术的海底局部地形自主感知方法及装置
CN109765562A (zh) * 2018-12-10 2019-05-17 中国科学院声学研究所 一种三维前视声像声纳系统和方法
CN111579639A (zh) * 2020-05-19 2020-08-25 福州大学 一种成像声呐扫描桥梁水下桩柱缺陷的测点布设方法
CN112285682A (zh) * 2020-10-20 2021-01-29 水利部交通运输部国家能源局南京水利科学研究院 水工工程涵洞环境的360°多波束声呐扫描装置及方法
CN112198516A (zh) * 2020-11-10 2021-01-08 广东智芯光电科技有限公司 一种水下图像采集的方法和系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于声呐成像的水下桩墩表观病害深度学习与智能检测;姜绍飞;《土木工程学报》;20210715;全文 *
水下无人航行器前视声呐系统研究;王玉双;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20200615;全文 *

Also Published As

Publication number Publication date
US20230014144A1 (en) 2023-01-19
US11841426B2 (en) 2023-12-12
CN113376639A (zh) 2021-09-10

Similar Documents

Publication Publication Date Title
CN113376639B (zh) 基于扫描声呐成像的桥墩基础冲刷区地形三维重构方法
CN107688180B (zh) 基于探地雷达的活动断裂浅地表层空间分布探测方法
Wilson et al. Sediment scour and deposition within harbors in California (USA), caused by the March 11, 2011 Tohoku-oki tsunami
Puente et al. Novel image analysis approach to the terrestrial LiDAR monitoring of damage in rubble mound breakwaters
CN114460554A (zh) 基于探地雷达多属性融合的岩溶洞穴三维成像方法
JP2009075017A (ja) 表層流速推定方法、装置並びにプログラム
Vetter et al. Estimating changes of riverine landscapes and riverbeds by using airborne LiDAR data and river cross-sections
Foster-Smith et al. Mapping seabed biotopes at two spatial scales in the eastern English Channel. Part 2. Comparison of two acoustic ground discrimination systems
CN108460422B (zh) 基于深度分布特征的海底地貌类型识别方法
CN114549751A (zh) 箱梁生产的模板监控系统及方法
CN112907615B (zh) 基于区域生长的海底地貌单元轮廓及细节识别方法
CN106228521A (zh) 一种基于薄板样条插值的障碍物特征提取方法
CN108896997B (zh) 复杂地形条件下侧扫声呐探测结果改正方法
CN113283437B (zh) 一种海底沙波特征识别方法
Shen et al. Rapid Three-Dimensional Reconstruction of Underwater Defective Pile Based on Two-Dimensional Images Obtained Using Mechanically Scanned Imaging Sonar
CN113009575A (zh) 一种地下水中污染物的监测方法及系统
Mirdan et al. Digital shoreline analysis system techniques for stability detection: An applied case study on Port Said, Egypt
CN113313651A (zh) 一种基于周围变化的侧扫声呐图像纹理失真区域修复方法
Goseberg et al. Highly-resolved numerical modeling of tsunami run-up and inundation scenario in the city of Padang, West Sumatra
CN113566753B (zh) 基于机械扫描成像声呐扫描桥墩基础冲刷的测点布设方法
Goseberg et al. Relevant factors on the extent of inundation based on tsunami scenarios for the city of Padang, West Sumatra
DOBROCHINSKI et al. SEDIMENT MANAGEMENT AND NAVIGATION OPTIMIZATION AT COLLIER CREEK, FLORIDA
CN117824618A (zh) 一种海底矿区可通行性地图获取方法
CN115048806B (zh) 风浪数值模拟方法
CN112487663B (zh) 能够考虑透浪和越浪的大尺度港池波要素高效模拟方法

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