CN110568417A - 基于多物理场的机载激光测深回波信号的处理方法 - Google Patents
基于多物理场的机载激光测深回波信号的处理方法 Download PDFInfo
- Publication number
- CN110568417A CN110568417A CN201910812249.5A CN201910812249A CN110568417A CN 110568417 A CN110568417 A CN 110568417A CN 201910812249 A CN201910812249 A CN 201910812249A CN 110568417 A CN110568417 A CN 110568417A
- Authority
- CN
- China
- Prior art keywords
- airborne laser
- taking
- processing
- echo signals
- formula
- 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.)
- Pending
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 35
- 238000000034 method Methods 0.000 claims abstract description 23
- 230000008878 coupling Effects 0.000 claims abstract description 10
- 238000010168 coupling process Methods 0.000 claims abstract description 10
- 238000005859 coupling reaction Methods 0.000 claims abstract description 10
- 238000013178 mathematical model Methods 0.000 claims abstract description 6
- 239000012530 fluid Substances 0.000 claims description 24
- 238000009792 diffusion process Methods 0.000 claims description 14
- 230000000704 physical effect Effects 0.000 claims description 11
- 230000035699 permeability Effects 0.000 claims description 5
- 238000005259 measurement Methods 0.000 abstract description 11
- 230000008901 benefit Effects 0.000 abstract description 3
- 238000001914 filtration Methods 0.000 description 17
- 230000000694 effects Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000013507 mapping Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 239000013535 sea water Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011160 research Methods 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/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
-
- 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
-
- 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/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/4802—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于多物理场的机载激光测深回波信号的处理方法,为解决现有技术不能精确的找到海面以及海底回波信号的位置问题,是将Darcy渗流方程与非线性扩散方程耦合,用回波信号的物理性质代替回波信号的几何性质,利用局部物理特性构造点扩散函数,接收到的回波信号作为多孔介质的渗透率。在数学模型上表现为满足初边值条件下多个方程的耦合,信号处理的结果是演化方程组在空间上的整体弱解。具有能精确地确定海面以及海底回波信号位置,从而实现测深的优点。
Description
技术领域
本发明涉及一种激光测深回波信号处理方法,特别是涉及一种基于多物理场的机载激光测深回波信号的处理方法。
背景技术
机载激光测深系统(Airborne LIDAR Bathymetry,ALB)因其灵活性好、速度快和测量精度高等特点被广泛应用于海洋测绘领域,而近岸浅海的水深、地形测量是海洋测绘最基本的任务之一。ALB主要由机上蓝绿激光发射/接收系统,以及地面数据后处理系统构成。通过激光器发射器发射大功率、窄脉冲激光来探测海底深度和水下目标。通过发射对水体穿透能力较强的蓝绿光(波长532nm)获取海底反射信号,发射对海面强反射的红外光(波长1064nm)获取高程;结合GPS、区域无线定位系统和姿态传感器所提供的空间位置和姿态信息,最终获取海底激光脚点的三维坐标。
机载激光水深测量的硬件设备已经相当成熟,我国诸多机构利用机载激光测深系统开展了一系列的研究工作,并且在很多领域都得到具体的应用。机载激光测深的最终目的是通过对接收的回波信号进行处理从而获得测量海域的水深信息。因此,机载激光测深的关键问题是准确的找到海面以及海底回波信号的位置,这就对回波信号数据处理提出了很高的要求。目前,机载激光测深系统常用的数据处理方法有回波探测、波形分解以及去卷积方法等。面对结构复杂的现实情况,各种信号滤波算法的判断条件也无法一一适应各种地貌特征且各种数据处理算法也有各自的不足之处。因此,采用较为完善的数据处理方法对机载激光测深回波信号进行处理是激光测深系统实现高精度测量的关键所在。
在机载激光测深的过程中,激光信号需经过大气、海水、海底往返等八种传输过程,因此在接收的回波信号中会引入大量的噪声,导致难以准确的确定海底、海面回波信号的位置,这也是当今海洋测绘领域的热点和难点。
发明内容
本发明目的在于克服现有技术的上述缺陷,提供一种能精确地确定海面以及海底回波信号位置的基于多物理场的机载激光测深回波信号的处理方法。也就是针对机载激光测深回波信号,提出了一种基于多物理场的处理方法,实现对机载激光回波信号的处理,从而精确的找到海面以及海底回波信号的位置,以达到实现信号处理的目的。
为实现上述目的,本发明基于多物理场的机载激光测深回波信号的处理方法将Darcy渗流方程与非线性扩散方程耦合,用回波信号的物理性质代替回波信号的几何性质,利用局部物理特性构造点扩散函数,接收到的回波信号作为多孔介质的渗透率。具有能精确地确定海面以及海底回波信号位置的优点。
作为优化,在数学模型上表现为满足初边值条件下多个方程的耦合,信号处理的结果是演化方程组在空间上的整体弱解。
作为优化,包括以下步骤:步骤一、利用探测器接收机载激光回波信号u0;步骤二、利用公式(1)中多物理场耦合模型进行处理:
u为机载激光测深回波信号,t为附加的迭代时间,为渗流场流体速度,P为渗流场流体压力,ρ为流体密度,μ为粘度流体系数,α、β、a、b为边界函数,K(u)为可选择的渗透函数,g为点扩散函数。
将上述(1)式中的(c)式代入(b)加上(d)式的边界条件,初次迭代K(u)初值取K(u0),得到:求出渗流场中流体的压力P0。
将P0代入到(c)式:
求出渗流场流体的速度场
将求出的压力P0和速度场代入(a)式:
得到第一次滤波后的信号u1。
步骤三、循环操作步骤二,依次得到u2,u3,…。
本发明提出的一种基于多物理场的机载激光测深的处理方法,具有明确的物理背景,利用回波信号的局部物理特性构造点扩散函数,用回波信号的物理性质代替回波信号的几何性质,实现回波信号的处理。可以精确地确定海面以及海底回波信号的位置,从而实现测深目的。
采用上述技术方案后,本发明基于多物理场的机载激光测深回波信号的处理方法具有能精确地确定海面以及海底回波信号位置,从而实现测深的优点。
附图说明
图1是本发明基于多物理场的机载激光测深回波信号的处理方法的数据处理流程框图;图2和3分别为本发明基于多物理场的机载激光测深回波信号的处理方法实施例一的接收信号u0图和接收信号u0处理后的效果图;图4和5分别为本发明基于多物理场的机载激光测深回波信号的处理方法实施例二的接收信号u′0图和接收信号u′0处理后的效果图。图6-8分别是本发明基于多物理场的机载激光测深回波信号的处理方法实施例二的接收信号u′0利用现有移动均值滤波处理方法、中值滤波处理方法、加权的线性最小二乘滤波(均值滤波)处理方法的处理效果图。
具体实施方式
下面结合附图和具体实施例对本发明技术方案作进一步详细描述,所描述的具体实施例仅对本发明进行解释说明,并不用以限制本发明。
本发明提出的一种基于多物理场的机载激光测深的处理方法,其设计思路是:将Darcy渗流方程与非线性扩散方程耦合,用回波信号的物理性质代替回波信号的几何性质,利用局部物理特性构造点扩散函数,接收到的回波信号作为多孔介质的渗透率,提出一种基于多物理场的信号处理方法,具有明确的物理背景。在数学模型上表现为满足初边值条件下多个方程的耦合,信号处理的结果是演化方程组在空间上的整体弱解。本发明可以很好地进行处理,精确地确定海面以及海底回波信号的位置,从而实现测深目的。
实施例一,如图1-3所示,本发明基于多物理场的机载激光测深回波信号的处理方法是将Darcy渗流方程与非线性扩散方程耦合,用回波信号的物理性质代替回波信号的几何性质,利用局部物理特性构造点扩散函数,接收到的回波信号作为多孔介质的渗透率。在数学模型上表现为满足初边值条件下多个方程的耦合,信号处理的结果是演化方程组在空间上的整体弱解。
本实施例接收信号u0图像如图2所示,数据处理效果图如图3所示,根据处理结果中海面、海底回波时刻位置就可以实现测深目的。具体实现的过程如图1所示,具体步骤如下。
步骤一、利用探测器接收机载激光回波信号u0;步骤二、利用公式(1)中多物理场耦合模型进行处理:
u为机载激光测深回波信号,t为附加的迭代时间,为渗流场流体速度,P为渗流场流体压力,ρ为流体密度,μ为粘度流体系数,α、β、a、b为边界函数,K(u)为可选择的渗透函数,g为点扩散函数。
将上述(1)式中的(c)式代入(b)加上(d)式的边界条件,初次迭代K(u)初值取K(u0),得到:
求出渗流场中流体的压力P0。
将P0代入到(c)式:
求出渗流场流体的速度场
将求出的压力P0和速度场代入(a)式:
得到第一次滤波后的信号u1。
步骤三、循环操作步骤二,依次得到u2,u3,...。
通过比较图2与图3发现,处理后的回波信号可以精确找到海面、海底回波信号的位置,从而计算出海水深度m:
本发明可以很好地进行处理,精确地确定海面以及海底回波信号的位置,从而实现测深目的。
实施例二,如图1和3-4,本发明基于多物理场的机载激光测深回波信号的处理方法是是将Darcy渗流方程与非线性扩散方程耦合,用回波信号的物理性质代替回波信号的几何性质,利用局部物理特性构造点扩散函数,接收到的回波信号作为多孔介质的渗透率。在数学模型上表现为满足初边值条件下多个方程的耦合,信号处理的结果是演化方程组在空间上的整体弱解。
本实施例接收信号u′0图像如图4所示,数据处理效果图如图5所示,根据处理结果中海面、海底回波时刻位置就可以实现测深目的。具体实现的过程如图1所示,具体步骤如下。
步骤一、利用探测器接收机载激光回波信号u′0;步骤二、利用公式(1)中多物理场耦合模型进行处理:
u为机载激光测深回波信号,t为附加的迭代时间,为渗流场流体速度,P为渗流场流体压力,ρ为流体密度,μ为粘度流体系数,α、β、a、b为边界函数,K(u)为可选择的渗透函数,g为点扩散函数。
将上述(1)式中的(c)式代入(b)加上(d)式的边界条件,初次迭代K(u)初值取K(u′0),得到:
求出渗流场中流体的压力P0。
将P0代入到(c)式:
求出渗流场流体的速度场
将求出的压力P0和速度场代入(a)式:
得到第一次滤波后的信号u′1。
步骤三、循环操作步骤二,依次得到u′2,u′3,...。
通过比较图4与图5发现,处理后的回波信号可以精确找到海面、海底回波信号的位置,从而计算出海水深度m′:
本发明可以很好地进行处理,精确地确定海面以及海底回波信号的位置,从而实现测深目的。
如图6-8所示,是基于实施例二的基于多物理场处理结果与现有的几种处理算法比较,基于多物理场的机载激光测深接收信号u′0图像如图4所示,数据处理效果图如图5所示,根据处理结果中海面、海底回波时刻位置就可以精确实现测深目的。移动均值滤波、中值滤波、加权的线性最小二乘滤波(均值滤波)处理效果如图6、图7、图8所示,具体结果见附图。
针对五组数据分别对多物理场处理方法与移动均值滤波、中值滤波、加权的线性最小二乘滤波(均值滤波)比较方差σ、标准差s以及滤波后峰值位置的确定得到如下结果。
表1五组数据分别对应的方差、标准差
表2五组数据分别对应海面、海底回波信号位置
通过观察表格数据可以得到一下结论:1)中值滤波虽然比多物理场处理方法得到的方差、标准差小,但该方法并不能准确的找到海面以及海底回波信号的位置。2)移动均值滤波、加权的线性最小二乘滤波(均值滤波)可以较好的的确定海面、海底回波信号的位置,但由上表可以观察到多物理场处理方法的方差、标准差更小。
综上分析,数值实验结果表明多物理场处理方法精度较高。与现有技术相比,本发明的有益效果是:本发明测量精度较高,本发明提出的基于多物理场的机载激光测深回波信号的处理方法,经过处理回波信号可以精确的找到海面、海底回波信号的位置,从而达到测量水深的目的。
Claims (3)
1.一种基于多物理场的机载激光测深回波信号的处理方法,其特征在于将Darcy渗流方程与非线性扩散方程耦合,用回波信号的物理性质代替回波信号的几何性质,利用局部物理特性构造点扩散函数,接收到的回波信号作为多孔介质的渗透率。
2.根据权利要求1所述基于多物理场的机载激光测深回波信号的处理方法,其特征在于在数学模型上表现为满足初边值条件下多个方程的耦合,信号处理的结果是演化方程组在空间上的整体弱解。
3.根据权利要求1-2任一所述基于多物理场的机载激光测深回波信号的处理方法,其特征在于包括以下步骤:
步骤一、利用探测器接收机载激光回波信号u0;
步骤二、利用公式(1)中多物理场耦合模型进行处理:
u为机载激光测深回波信号,t为附加的迭代时间,为渗流场流体速度,P为渗流场流体压力,ρ为流体密度,μ为粘度流体系数,α、β、a、b为边界函数,K(u)为可选择的渗透函数,g为点扩散函数;
将上述(1)式中的(c)式代入(b)加上(d)式的边界条件,初次迭代K(u)初值取K(u0),得到:
求出渗流场中流体的压力P0;
将P0代入到(c)式:
求出渗流场流体的速度场
将求出的压力P0和速度场代入(a)式:
得到第一次滤波后的信号u1;
步骤三、循环操作步骤二,依次得到u2,u3,…。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910812249.5A CN110568417A (zh) | 2019-08-30 | 2019-08-30 | 基于多物理场的机载激光测深回波信号的处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910812249.5A CN110568417A (zh) | 2019-08-30 | 2019-08-30 | 基于多物理场的机载激光测深回波信号的处理方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110568417A true CN110568417A (zh) | 2019-12-13 |
Family
ID=68776985
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910812249.5A Pending CN110568417A (zh) | 2019-08-30 | 2019-08-30 | 基于多物理场的机载激光测深回波信号的处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110568417A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116609758A (zh) * | 2023-07-17 | 2023-08-18 | 山东科技大学 | 一种机载激光测深波形旅行时提取方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105954732A (zh) * | 2016-06-14 | 2016-09-21 | 深圳大学 | 一种激光雷达测深数据的海底底质反射率提取方法及系统 |
EP3342474A1 (en) * | 2017-01-03 | 2018-07-04 | Institute of Nuclear Energy Research, Atomic | A composite material type oxygen transmission film and its preparation method |
-
2019
- 2019-08-30 CN CN201910812249.5A patent/CN110568417A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105954732A (zh) * | 2016-06-14 | 2016-09-21 | 深圳大学 | 一种激光雷达测深数据的海底底质反射率提取方法及系统 |
EP3342474A1 (en) * | 2017-01-03 | 2018-07-04 | Institute of Nuclear Energy Research, Atomic | A composite material type oxygen transmission film and its preparation method |
Non-Patent Citations (1)
Title |
---|
孙明明: "一类耦合偏微分方程问题及其在图像处理中的应用", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116609758A (zh) * | 2023-07-17 | 2023-08-18 | 山东科技大学 | 一种机载激光测深波形旅行时提取方法 |
CN116609758B (zh) * | 2023-07-17 | 2023-10-27 | 山东科技大学 | 一种机载激光测深波形旅行时提取方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Johnson‐Roberson et al. | High‐resolution underwater robotic vision‐based mapping and three‐dimensional reconstruction for archaeology | |
CN105193445B (zh) | 一种超声探头标定体模、超声探头标定系统及其标定方法 | |
Roman et al. | Application of structured light imaging for high resolution mapping of underwater archaeological sites | |
CN110207695B (zh) | 一种适用于深海auv的无速度辅助单信标定位方法 | |
CN109343022B (zh) | 估测层间土壤含水量的方法 | |
CN105954732A (zh) | 一种激光雷达测深数据的海底底质反射率提取方法及系统 | |
CN107664758B (zh) | 基于长基线或超短基线组网的深海导航定位系统及方法 | |
Wright et al. | Depth calibration and validation of the Experimental Advanced Airborne Research LIDAR, EAARL-B | |
CN107479045B (zh) | 基于全波形激光雷达点云数据剔除低矮植被的方法及系统 | |
CN110568417A (zh) | 基于多物理场的机载激光测深回波信号的处理方法 | |
CN109212598B (zh) | 基于直达波反演的三维空间二次定位方法 | |
Klopfer et al. | Assessing the potential of a low-cost 3-D sensor in shallow-water bathymetry | |
CN112907615B (zh) | 基于区域生长的海底地貌单元轮廓及细节识别方法 | |
Diamanti et al. | Geometric documentation of underwater archaeological sites | |
CN113030919A (zh) | 一种基于模型拟合的波形检测方法及系统 | |
CN117710458A (zh) | 一种基于双目视觉的舰载机着舰过程相对位置测量方法及系统 | |
CN117237430A (zh) | 一种高精度多时序水深反演方法、计算设备及存储介质 | |
CN110057990A (zh) | 一种多参数水质剖面仪的pH校正方法 | |
CN112540340B (zh) | 精度误差补偿方法和基于该方法的自校准声信标定位设备 | |
CN109816713B (zh) | 一种非自主测绘海图的水深注记精度定量评估方法 | |
CN208043067U (zh) | 水下地形及水流速度和温盐度综合自动测量装置 | |
CN113379710A (zh) | 一种水下目标声呐精准测量系统和方法 | |
Specht | Multi-sensor integration of hydroacoustic and optoelectronic data acquired from UAV and USV vehicles on the inland waterbody | |
TW201430366A (zh) | 應用微波雷達於海岸線及潮間帶地形測量之技術 | |
CN204165553U (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20191213 |