CN113759387B - 一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法 - Google Patents
一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法 Download PDFInfo
- Publication number
- CN113759387B CN113759387B CN202111019743.XA CN202111019743A CN113759387B CN 113759387 B CN113759387 B CN 113759387B CN 202111019743 A CN202111019743 A CN 202111019743A CN 113759387 B CN113759387 B CN 113759387B
- Authority
- CN
- China
- Prior art keywords
- wave
- resistant building
- embankment
- laser radar
- amount
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000001914 filtration Methods 0.000 claims abstract description 15
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 18
- 238000012545 processing Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000003746 surface roughness Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 4
- 235000014653 Carica parviflora Nutrition 0.000 description 2
- 241000243321 Cnidaria Species 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法,包括如下步骤:采用激光雷达获取堤顶波面的三维坐标点云;对三维坐标点云进行滤波处理;计算防浪建筑物的堤顶压强和堤顶波峰处的流速分布;基于获取的堤顶压强和堤顶波峰处的流速分布,计算防浪建筑物堤顶流量;根据防浪建筑物堤顶流量,计算获得防浪建筑物堤顶总越浪量。本发明可以实时获取防浪建筑物堤顶的流量和总越浪量,不但测量精度高,而且可以反映整个越浪过程,为实验室和现场的越浪量分析提供了切实可行的方法。
Description
技术领域
本发明属于海洋观测技术领域,具体涉及一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法。
背景技术
越浪量是防浪建筑物设计的重要参数,对防浪建筑物堤身稳定和堤顶安全具有重要影响。流量和总越浪量是反映越浪量的重要参数,流量是指单位时间通过堤顶的水量,总越浪量是一段时间内通过海堤的总水量。
目前,越浪量的测量主要有两种方法。第一种方法是通过收纳越浪水体计算越浪量,另一种方法是通过视频和图像处理技术。第一种方法测量较为准确,但是这种方法不能反映越浪过程。第二种方法虽然能反映越浪过程,但是图像和视频的处理难度较大,并且越浪量测量误差较大。
发明内容
发明目的:为了克服现有技术中存在的不足,提供一种基于三维激光雷达的珊瑚礁海岸防浪建筑物越浪量测量方法,可以实时获取防浪建筑物堤顶的流量和总越浪量,不但测量精度高,而且可以反映整个越浪过程,为实验室和现场的越浪量分析提供了切实可行的方法。
技术方案:为实现上述目的,本发明提供一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法,包括如下步骤:
S1:采用激光雷达获取堤顶波面的三维坐标点云;
S2:对三维坐标点云进行滤波处理;
S3:计算防浪建筑物的堤顶压强和堤顶波峰处的流速分布;
S4:基于获取的堤顶压强和堤顶波峰处的流速分布,计算防浪建筑物堤顶流量;
S5:根据防浪建筑物堤顶流量,计算获得防浪建筑物堤顶总越浪量。
进一步地,所述步骤S1中激光雷达布置于堤顶上方,所述激光雷达内置有GPS和惯性导航系统。能够获取世界标准时间和激光雷达偏转角。三维激光雷达可以实时获取堤顶波面的三维坐标点云、时间戳和反射强度信息。
所述防浪建筑物的堤顶压强主要是用于分析建筑物受力,能够更好地进行越浪量分析。
进一步地,所述步骤S2中滤波处理的具体方法为:
A1:采用直通滤波器提取堤顶范围内波面的三维坐标点云;
A2:采用半径滤波器去除波面三维坐标点云中的离群点;
A3:对每一帧的波面数据进行插值;
A4:采用移动平均滤波器对波形做平滑处理。
进一步地,所述步骤S2中直通滤波器在水平方向上设置有坐标阈值,水平方向的坐标阈值范围等于防浪建筑物所在位置的水平坐标范围。
进一步地,所述步骤S2中所述半径滤波器设定有滤波半径和半径范围内点数的阈值,其中,滤波半径的取值和半径范围内点数阈值的取值与防浪建筑物堤顶三维坐标点云的采样数有关。其主要目的是去除离群点,滤波半径取值和点数阈值通过经验调试,选取滤波效果较好时的参数。
进一步地,所述步骤S3中防浪建筑物的堤顶压强的计算方法为:P=ρgz,其中:ρ是水体密度;z是波面高度(波面高程与堤顶高程之差);P是压强,用于防浪建筑物的受力分析。
进一步地,所述步骤S3中防浪建筑物的堤顶波峰处的流速分布的计算方法为:
通过单位时间内波峰移动的水平距离得到波峰处的速度,由流速沿水深的对数分布式获取到堤顶波峰处的流速分布。
进一步地,所述步骤S3中流速沿水深的对数分布式为:
式中:κ=0.4是卡门常数;z是计算点位置距防浪建筑物顶部的高度;是速度零点高度,△是防浪建筑物的表面粗糙度;u是任意水深处的流速;u*是摩阻流速,摩阻流速可基于流速沿水深的对数分布式和波面流速计算得到。
进一步地,所述步骤S4中防浪建筑物堤顶流量的计算方法为:
其中,q是通过防浪建筑物堤顶的流量;n是这一时刻波峰线上坐标点的数量;zj是沿岸方向上第j个波峰处的波面高度(波面高程与堤顶高程之差);uj是沿岸方向上第j个波峰处的流速;Δyj是第j个波峰处所对应的波面宽度。
进一步地,所述步骤S5中防浪建筑物堤顶总越浪量的计算方法为:
其中,Q是堤顶总越浪量;qi是i时刻的堤顶流量;Δti是i时刻与前一时刻的时间差。
有益效果:本发明与现有技术相比,具备如下优点:
1、通过三维激光雷达直接获取防浪建筑物堤顶的波形数据,数据集具有高精度和高分辨率的特性,保证了测量精度。
2、对三维坐标点云进行了滤波处理,能够去除数据集中的噪声点,提升测量精度。
3、能够实时获取堤顶的流量和总越浪量,可以反映整个越浪过程。
4、本发明方法的整个测量过程实现了无接触测量,可以保护测量仪器的安全,同时保证测量结果不受波浪作用力的影响。
附图说明
图1是本发明方法的流程图;
图2是本发明具体实施例的模型实验布置图;
图3是堤顶流量计算方法简图。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本发明提供一种基于三维激光雷达的珊瑚礁海岸防浪建筑物越浪量测量方法,如图1所示,其包括如下步骤:
S1:采用激光雷达获取堤顶波面的三维坐标点云:
如图2所示,激光雷达布置在堤顶上方,内置有GPS和惯性导航系统,获取世界标准时间和激光雷达偏转角。三维激光雷达可以实时获取堤顶波面的三维坐标点云、时间戳和反射强度信息。
S2:对三维坐标点云进行滤波处理:
A1:采用直通滤波器提取堤顶所在范围内波面的三维坐标点云,其中,直通滤波器需要在水平方向上设置阈值范围,水平方向的阈值范围等于海堤位置的水平坐标范围。
A2:采用半径滤波器去除波面三维坐标点云中的离群点,半径滤波器需要设定滤波半径和半径范围内其他点个数的阈值,其中,滤波半径的取值和半径范围内其他点个数阈值的取值与堤顶三维坐标点云的采样数有关。
A3:对每一帧的波面数据进行插值;
A4:采用移动平均滤波器对波形做平滑处理。
S3:计算防浪建筑物的堤顶压强和堤顶波峰处的流速分布:
防浪建筑物的堤顶压强的计算方法为:P=ρgz,其中:ρ是水体密度;z是波面高度(波面高程与堤顶高程之差);P是压强,用于防浪建筑物的受力分析。
防浪建筑物的堤顶波峰处的流速分布的计算方法为:
通过单位时间内波峰移动的水平距离得到波峰处的速度,由流速沿水深的对数分布式获取到堤顶波峰处的流速分布。
流速沿水深的对数分布式为:
式中:κ=0.4是卡门常数;z是计算点位置距防浪建筑物顶部的高度;是速度零点高度,△是防浪建筑物的表面粗糙度;u是任意水深处的流速;u*是摩阻流速,摩阻流速可基于流速沿水深的对数分布式和波面流速计算得到。
S4:基于获取的堤顶压强和堤顶波峰处的流速分布,计算防浪建筑物堤顶流量:
如图3所示,防浪建筑物堤顶流量的计算方法为:
其中,q是通过防浪建筑物堤顶的流量;n是这一时刻波峰线上坐标点的数量;zj是沿岸方向上第j个波峰处的波面高度(波面高程与堤顶高程之差);uj是沿岸方向上第j个波峰处的流速;Δyj是第j个波峰处所对应的波面宽度。
S5:根据防浪建筑物堤顶流量,计算获得防浪建筑物堤顶总越浪量:
防浪建筑物堤顶总越浪量的计算方法为:
其中,Q是堤顶总越浪量;qi是i时刻的堤顶流量;Δti是i时刻与前一时刻的时间差。
Claims (7)
1.一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法,其特征在于,包括如下步骤:
S1:采用激光雷达获取堤顶波面的三维坐标点云;
S2:对三维坐标点云进行滤波处理;
S3:计算防浪建筑物的堤顶压强和堤顶波峰处的流速分布;
通过单位时间内波峰移动的水平距离得到波峰处的速度,由流速沿水深的对数分布式获取到堤顶波峰处的流速分布;所述流速沿水深的对数分布式为:
式中:κ是卡门常数;z是计算点位置距防浪建筑物顶部的高度;是速度零点高度,△是防浪建筑物的表面粗糙度;u是任意水深处的流速;u*是摩阻流速,摩阻流速可基于流速沿水深的对数分布式和波面流速计算得到;
S4:基于获取的堤顶波峰处的流速分布,计算防浪建筑物堤顶流量;
所述防浪建筑物堤顶流量的计算方法为:
其中,q是通过防浪建筑物堤顶的流量;n是这一时刻波峰线上坐标点的数量;zj是沿岸方向上第j个波峰处的波面高度;uj是沿岸方向上第j个波峰处的流速;Δyj是第j个波峰处所对应的波面宽度;
S5:根据防浪建筑物堤顶流量,计算获得防浪建筑物堤顶总越浪量。
2.根据权利要求1所述的一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法,其特征在于,所述步骤S1中激光雷达布置于堤顶上方,所述激光雷达内置有GPS和惯性导航系统。
3.根据权利要求1所述的一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法,其特征在于,所述步骤S2中滤波处理的具体方法为:
A1:采用直通滤波器提取堤顶范围内波面的三维坐标点云;
A2:采用半径滤波器去除波面三维坐标点云中的离群点;
A3:对每一帧的波面数据进行插值;
A4:采用移动平均滤波器对波形做平滑处理。
4.根据权利要求3所述的一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法,其特征在于,所述步骤S2中直通滤波器在水平方向上设置有坐标阈值,水平方向的坐标阈值范围等于防浪建筑物所在位置的水平坐标范围。
5.根据权利要求3所述的一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法,其特征在于,所述步骤S2中所述半径滤波器设定有滤波半径和半径范围内点数的阈值,其中,滤波半径的取值和半径范围内点数阈值的取值与防浪建筑物堤顶三维坐标点云的采样数有关。
6.根据权利要求1所述的一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法,其特征在于,所述步骤S3中防浪建筑物的堤顶压强的计算方法为:P=ρgz,其中:ρ是水体密度;z是波面高度;P是压强,用于防浪建筑物的受力分析。
7.根据权利要求1所述的一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法,其特征在于,所述步骤S5中防浪建筑物堤顶总越浪量的计算方法为:
其中,Q是堤顶总越浪量;qi是i时刻的堤顶流量;Δti是i时刻与前一时刻的时间差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111019743.XA CN113759387B (zh) | 2021-09-01 | 2021-09-01 | 一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111019743.XA CN113759387B (zh) | 2021-09-01 | 2021-09-01 | 一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113759387A CN113759387A (zh) | 2021-12-07 |
CN113759387B true CN113759387B (zh) | 2023-12-22 |
Family
ID=78792304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111019743.XA Active CN113759387B (zh) | 2021-09-01 | 2021-09-01 | 一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113759387B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115840975B (zh) * | 2022-11-15 | 2023-08-01 | 广东省水利水电科学研究院 | 一种风暴潮增水漫堤预警方法、系统、装置和存储介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012100763A1 (de) * | 2011-01-24 | 2012-08-02 | Technische Universität Dresden | Verfahren zur bestimmung von geschwindigkeiten in strömungen und phasen-frequenz-geschwindigkeits-feldsensor |
CN106595777A (zh) * | 2016-12-01 | 2017-04-26 | 广西师范大学 | 一种非接触式探测河流断面流量的计算方法 |
CN107576361A (zh) * | 2017-09-01 | 2018-01-12 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种基于激光多普勒测速的大口径管路流量测量装置 |
CN108917860A (zh) * | 2018-05-22 | 2018-11-30 | 赵彦增 | 基于航拍激光点云的河道流量测验方法 |
CN109781196A (zh) * | 2019-02-14 | 2019-05-21 | 苏州万闻信息工程有限公司 | 一种非接触式的测量带式输送机上物料流量的方法 |
CN111649791A (zh) * | 2020-06-18 | 2020-09-11 | 浙江省水利河口研究院 | 一种水面上测量多流态水体流量的装置及方法 |
CN112950610A (zh) * | 2021-03-18 | 2021-06-11 | 河海大学 | 一种裂流监测和预警方法及系统 |
CN113219481A (zh) * | 2021-03-29 | 2021-08-06 | 河海大学 | 一种基于三维激光雷达的破波带波浪水动力监测方法及其系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8903134B2 (en) * | 2010-07-21 | 2014-12-02 | Ron Abileah | Methods for mapping depth and surface current |
-
2021
- 2021-09-01 CN CN202111019743.XA patent/CN113759387B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012100763A1 (de) * | 2011-01-24 | 2012-08-02 | Technische Universität Dresden | Verfahren zur bestimmung von geschwindigkeiten in strömungen und phasen-frequenz-geschwindigkeits-feldsensor |
CN106595777A (zh) * | 2016-12-01 | 2017-04-26 | 广西师范大学 | 一种非接触式探测河流断面流量的计算方法 |
CN107576361A (zh) * | 2017-09-01 | 2018-01-12 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种基于激光多普勒测速的大口径管路流量测量装置 |
CN108917860A (zh) * | 2018-05-22 | 2018-11-30 | 赵彦增 | 基于航拍激光点云的河道流量测验方法 |
CN109781196A (zh) * | 2019-02-14 | 2019-05-21 | 苏州万闻信息工程有限公司 | 一种非接触式的测量带式输送机上物料流量的方法 |
CN111649791A (zh) * | 2020-06-18 | 2020-09-11 | 浙江省水利河口研究院 | 一种水面上测量多流态水体流量的装置及方法 |
CN112950610A (zh) * | 2021-03-18 | 2021-06-11 | 河海大学 | 一种裂流监测和预警方法及系统 |
CN113219481A (zh) * | 2021-03-29 | 2021-08-06 | 河海大学 | 一种基于三维激光雷达的破波带波浪水动力监测方法及其系统 |
Non-Patent Citations (4)
Title |
---|
Nonintrusive Measurement of Ocean Waves: Lidar Wave Gauge;Jennifer L. Irish et al.;《Journal of Atmospheric and Oceanic Technology》;第23卷(第11期);1559-1572 * |
三维激光雷达扫描技术在河流模拟中的应用;张力元 等;《水电能源科学》;第38卷(第8期);121-124 * |
浅化波浪层流边界层流速分布特性的数值分析;李诚;张弛;隋倜倜;;海洋学报(第05期);143-151 * |
海流测量技术发展及应用;周庆伟;白杨;封哲;汪小勇;王延伟;;海洋测绘(第03期);76-80 * |
Also Published As
Publication number | Publication date |
---|---|
CN113759387A (zh) | 2021-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Vousdoukas et al. | The role of combined laser scanning and video techniques in monitoring wave-by-wave swash zone processes | |
CN104156629B (zh) | 一种基于相对辐射校正的导航雷达图像反演海面风向方法 | |
CN113012195B (zh) | 基于光流计算的河流表面流速估计的方法和可读存储介质 | |
CN113759387B (zh) | 一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法 | |
van Scheltinga et al. | Observations of dune interactions from DEMs using through-water Structure from Motion | |
CN211927018U (zh) | 一种新型液位测量装置 | |
Muste et al. | Validation and extension of image velocimetry capabilities for flow diagnostics in hydraulic modeling | |
Porter et al. | Comparison of three techniques for scour depth measurement: photogrammetry, echosounder profiling and a calibrated pile | |
CN103852813A (zh) | 雨滴三维尺度检测装置及利用该装置计算雨滴体积的方法 | |
JP2008216010A (ja) | 河川流量算出装置、河川流量算出方法およびコンピュータプログラム | |
CN116222676B (zh) | 精准定位的毫米波水流量监测方法及系统 | |
Negri et al. | Self-synchronized phase averaging of PIV measurements in the base region of a rectangular cylinder | |
CN104297753A (zh) | 一种基于自适应缩减算子的导航雷达图像反演海面风向方法 | |
CN203745670U (zh) | 雨滴三维尺度检测装置 | |
CN107554719B (zh) | 一种基于声呐系统的船舶载重测量方法 | |
Maas et al. | Improvements in lidar bathymetry data analysis | |
CN110412312A (zh) | 涌潮流速垂向分布的一种计算方法 | |
Li et al. | Spatial-temporal measurement of waves in laboratory based on binocular stereo vision and image processing | |
CN113358191B (zh) | 基于条纹投影结构光的全域洪水水位实时监测方法 | |
Johnson et al. | Measurements of water surface profile and velocity field at a circular pier | |
Tuyen | Flow over oblique weirs | |
CN116203277B (zh) | 基于ptv与piv技术的海表面小尺度流场测量方法 | |
De Vries et al. | Using stereo photogrammetry to measure coastal waves | |
Chen et al. | Application of large-scale particle image velocimetry at the Hydraulics Laboratory of Colorado State University | |
CN116125459B (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 |