留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于机载LiDAR测深水体波形的漫衰减系数提取方法

亓超 周丰年 吴敬文 宿殿鹏 王贤昆 阳凡林

亓超,周丰年,吴敬文,等. 基于机载LiDAR测深水体波形的漫衰减系数提取方法[J]. 海洋学报,2021,43(1):147–154 doi: 10.12284/hyxb2021005
引用本文: 亓超,周丰年,吴敬文,等. 基于机载LiDAR测深水体波形的漫衰减系数提取方法[J]. 海洋学报,2021,43(1):147–154 doi: 10.12284/hyxb2021005
Qi Chao,Zhou Fengnian,Wu Jingwen, et al. Extraction method for diffuse attenuation coefficient based on airborne LiDAR bathymetric water column waveform [J]. Haiyang Xuebao,2021, 43(1):147–154 doi: 10.12284/hyxb2021005
Citation: Qi Chao,Zhou Fengnian,Wu Jingwen, et al. Extraction method for diffuse attenuation coefficient based on airborne LiDAR bathymetric water column waveform [J]. Haiyang Xuebao,2021, 43(1):147–154 doi: 10.12284/hyxb2021005

基于机载LiDAR测深水体波形的漫衰减系数提取方法

doi: 10.12284/hyxb2021005
基金项目: 国家自然科学基金(41930535, 41830540, 52001189);国家重点研发计划(2017YFC1405006, 2018YFF0212203, 2018YFC1405900, 2016YFC1401210);山东科技大学科研创新团队支持计划(2019TDJH103)。
详细信息
    作者简介:

    亓超(1994-),男,山东省莱芜市人,博士生,主要从事机载激光测深及其数据处理方面的研究工作。E-mail:qichoice007@sdust.edu.cn

    通讯作者:

    阳凡林(1974-),男,湖北省荆州市人,教授,主要从事海底地形测量和海洋定位导航方面的研究工作。E-mail:flyang@sdust.edu.cn

  • 中图分类号: P229;P714.1

Extraction method for diffuse attenuation coefficient based on airborne LiDAR bathymetric water column waveform

  • 摘要: 漫衰减系数是一个重要的海洋光学参数,能够为水体环境变化、水质分析以及水产养殖等方面提供基础性数据。针对目前船载实地测量效率与分辨率低、卫星遥感反演精度与分辨率较低的局限性,本文提出一种基于机载LiDAR测深水体波形的漫衰减系数提取方法。该方法首先通过分层异构模型的机载LiDAR波形分解算法得到水体散射回波,利用激光在水体中的衰减特性,构建漫衰减系数提取模型,最终获取大面积水域漫衰减系数的空间分布。采用西沙甘泉岛与江苏连云港两个航次的实测数据对所提算法进行了验证,本算法无需每个测深点的水底底部回波强度和深度即可反演得到漫衰减系数,并且在浑浊水域也可取得较好的效果,表明在中国近海利用机载LiDAR测深系统能够有效获取高精度的漫衰减系数。
  • 漫衰减系数(Diffuse Attenuation Coefficient, Kd)是水体的固有光学属性,反映了激光在水体中受散射、吸收所引起的能量衰减变化,是分析海水水质的重要参数之一[1]。漫衰减系数的大小主要取决于水体中的浮游生物含量(叶绿素浓度)、悬浮泥沙含量(浑浊度)、营养盐含量(黄色物质、溶解有机物质、盐度指标)以及其他污染物等因素,研究漫衰减系数可对水体富营养化、水体浑浊度骤增以及水体污染等环境变化实现实时、有效监测,在海洋科学研究、海洋工程应用、海洋渔业研究和海岸带环境监测等方面发挥重要的作用[2-3]

    由于水体的流动性,导致了水中物质的非均匀分布,增加了水体中漫衰减系数的获取难度。目前,国内外针对漫衰减系数的提取主要通过实地测量、被动遥感以及主动遥感3种方式进行。实地测量利用赛克盘和透射计等工具测量水体透明度,进而反演得到漫衰减系数,这些方式需以测量船为载体,效率较低且覆盖范围小,难以快速获取大面积水域的漫衰减系数空间分布[2, 4-5];通过卫星遥感影像数据可间接反演漫衰减系数,这种被动遥感方式可对大范围海域的漫衰减系数进行求取,但是无法在云雾天气、夜晚实施测量,且由于卫星设备本身的限制,测量精度难以保证,成果分辨率较低[6-8];全波形机载LiDAR测深(Airborne LiDAR Bathymetry, ALB)技术的出现弥补了上述方法的缺陷,运用水体散射回波的强度变化可实现漫衰减系数的快速反演[4, 9-13]。Billard等[9]利用ALB技术在不同深度下的辐照度进行指数拟合,进而反演得到漫衰减系数的方法,该方法需要水底底部回波强度和相应深度信息;Ding等[10]提出了一种基于单波长ALB全波形数据的漫衰减系数计算方法,克服了传统方法的局限性,但未对水体散射回波进行拟合分解,易受噪声影响。

    因此,针对目前船载实地测量效率与分辨率低、卫星遥感反演精度与分辨率较低的局限性,本文结合ALB技术高精度、高分辨率、灵活机动、快速高效的特点,提出了一种基于机载LiDAR测深水体波形的漫衰减系数提取方法。ALB原始波形数据经小波去噪之后,利用基于分层异构模型的波形分解算法[14-15](水面−高斯函数、水体−双指数函数及水底−高斯函数)分解得到水体散射回波信号,进而获取大面积水域漫衰减系数的空间分布。

    机载LiDAR测深水体波形的漫衰减系数提取方法充分利用了ALB技术高精度、高分辨率、灵活机动、快速高效的独特优势,考虑到激光在水体中的衰减特性以及计算效率,通过分层异构模型的机载LiDAR波形分解算法分解得到水体散射回波信号,并用其与水体散射回波强度方程构建相应模型,用以计算每个激光测深点的漫衰减系数。该算法的步骤为:首先获取机载LiDAR测深原始波形数据,对其进行去噪;然后利用分层异构模型对去噪后的波形数据进行分解,分解得到水体散射回波信号;再通过水体散射回波信号分别计算得到双段指数所对应的初始漫衰减系数;最后根据双段指数的采样时间确定相应权重,对初始漫衰减系数取加权平均,得到该激光测深点最终的漫衰减系数值。本文算法具体流程如图1所示。

    图  1  本文方法流程图
    Figure  1.  Flow diagram of the new method in this study

    机载LiDAR测深系统利用具有较强透水能力的蓝绿色波段(532 nm)激光进行探测,其激光接收器能够以数字化的形式记录每个激光脉冲的全部回波波形信号(图2[15-16],包括水面反射回波、水体散射回波、水底反射回波以及噪声4部分[17],可描述为

    图  2  机载LiDAR测深波形示意图
    Figure  2.  Airborne LiDAR bathymetric waveform
    $$f\left( t \right) = {f_{\rm{s}}}(t) + {f_{\rm{c}}}(t) + {f_{\rm{b}}}(t) + {f_{\rm{N}}}(t),$$ (1)

    式中,t表示时间刻度;fs(t)表示水面反射回波;fc(t)表示水体散射回波;fb(t)表示水底反射回波;fN(t)表示原始波形的噪声;f(t)为t时刻回波的振幅。

    其中,噪声包括背景噪声和传感器内部噪声,其存在会导致波形分解函数初始参数估计不准等问题,容易获得含有粗差的波形,影响波形分解精度。因此,在波形分解之前,需对原始回波信号中的噪声进行相应的滤除,文中使用小波算法对噪声进行处理[15, 18]

    波形分解是机载LiDAR测深数据处理的重要环节,是漫衰减系数提取的基础。根据机载LiDAR测深回波信号各部分的波形特征,选择合适的函数构建分层异构模型(水面−高斯函数、水体−双指数函数及水底−高斯函数)进行分解,从而将叠加在水面反射回波和水底反射回波中的水体散射回波精确地分解出来,得到水体散射回波所对应的函数。

    考虑到ALB系统发射和接收的激光脉冲信号近似服从高斯分布[19-20],因此,根据水面和水底反射回波相应特性分析,文中采用高斯函数(公式(2))分解得到水面、水底反射回波信号,

    $${f_i}(t; {A_i}, {\mu _i}, {\sigma _i}) = {A_i}\exp \left[ { - \frac{{{{\left( {t - {\mu _i}} \right)}^2}}}{{2{\sigma _i}^2}}} \right],$$ (2)

    式中,i表示水面(水底)反射回波信号的分解结果;Aiμiσi分别为水面(水底)高斯函数的波峰值、波峰位置和半幅波宽。

    从去噪后的ALB波形数据中剔除水面、水底反射回波信号,即可得到初始水体散射回波信号。根据光在水体中辐射传输定理可知,激光在水体中传播时,其脉冲能量随水深的增加而呈现指数衰减。针对水体散射回波的分解,Ding等[21]在三角形函数[22]和四边形函数[23]分解方法的基础之上,提出了一种改进的四边形函数,将水体散射回波下降沿的一段采用指数函数描述。由于水体浑浊度在垂直剖面上非均匀分布导致激光脉冲的衰减程度不同[24],综合激光在水体中衰减特性和分解效率的问题,如图3所示,文中将水体沿垂直剖面分为两层,利用双指数函数ABCD(公式(3))分解得到水体散射回波信号

    图  3  双指数函数分解水体散射回波示意图
    Figure  3.  Decomposition of water column contribution by double-exponential function
    $${f_c}(t; {a_x}, {b_x}, {b_y}, {c_x}, {c_y}, {d_x}, {d_y}) = \left\{ {\begin{aligned} &0\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad{t < {a_x}} \\ &0\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad{t \geqslant {d_x}} \\ & {{b_y}\left( {\displaystyle \frac{{t - {a_x}}}{{{b_x} - {a_x}}}} \right)}\qquad\qquad\qquad\qquad\qquad\qquad\qquad{{a_x} \leqslant t < {b_x}} \\ & {{\rm{exp}}\left( {\displaystyle \frac{{{c_x}\ln {b_y} - {b_x}\ln {c_y}}}{{{c_x} - {b_x}}}} \right)\exp {{\left[ {\displaystyle \frac{{\left( {\ln {c_y} - \ln{b_y}} \right)}}{{{c_x} - {b_x}}}} \right]}^t}}\qquad\;{{b_x} \leqslant t < {c_x}} \\ & {{\rm{exp}}\left( {\displaystyle \frac{{{d_x}\ln {c_y} - {c_x}\ln {d_y}}}{{{d_x} - {c_x}}}} \right)\exp {{\left[ {\displaystyle \frac{{\left( {\ln {d_y} - \ln {c_y}} \right)}}{{{d_x} - {c_x}}}} \right]}^t}}\qquad\;{{c_x} \leqslant t < {d_x}} \end{aligned}} \right.,$$ (3)

    式中,axbxcxdx分别表示双指数函数4个顶点(A、B、C、D)的位置时刻;bycydy为B、C、D的振幅值。

    非线性拟合需要预设初值,以降低迭代次数,增加拟合速度。文中遵循非线性最小二乘(Non-linear Least Squares, NLS)方法(公式(4)),采用Levenberg-Marquardt算法对13个初始参数(水面–高斯函数3个、水底−高斯函数3个、水体−双指数函数7个)进行迭代优化,直到拟合偏差达到最小为止,以求得各参数的精确值。将优化后的水体反向散射回波参数代入公式(3),得到最优的水体散射回波分解结果,进而用于漫衰减系数的提取。

    $$\mathop {\min }\limits_\lambda S(\lambda ) = \mathop {\min }\limits_\lambda {\sum\limits_{j = 1}^m {\left[ {{F_j} - f({t_j},\lambda )} \right]} ^2},$$ (4)

    式中,λ表示所求参数;S(λ)表示拟合偏差的平方和;Fj表示波形分解结果在j时刻的振幅。

    机载LiDAR测深系统在深度为Di处的水体散射回波强度方程可由下式来描述[4, 17]

    $${P_c}({D_i}) = \frac{{{P_e}T_{{\rm{atm}}}^2{A_R}{\eta _e}{\eta _{\rm{R}}}F{{\left( {1 - {L_S}} \right)}^2}\beta (\phi ){{\cos }^2}\theta }}{{{{\left( {{n_w}H + {D_i}} \right)}^2}}}{\rm{exp}}\left( {\frac{{ - 2{K_d}{D_i}}}{{\cos {\theta _w}}}} \right),$$ (5)
    $${D_i} = \frac{c}{{{n_w}}}\left( {{t_i} - {t_s}} \right)\cos {\theta _w},$$ (6)

    式中,Di表示水层深度;Pc(Di)表示深度Di时的水体散射回波强度;Pe表示激光发射强度;$T_{{\rm{atm}}}^2$表示大气双程损失因子;AR表示激光接收器的视场面积;ηeηR分别表示激光发射器的光学效率和激光接收器的光学效率;F表示接收视场角损失因子;LS是通过表面传输的损失,即表面反照率;βϕ)表示体散射函数;θ表示激光传入水体时的入射角,即天底角;θw表示射入水体后的折射角;nw表示水体折射率;H表示飞机航高;Kd表示532 nm波段处的漫衰减系数;c表示激光在真空的传播速度;ti表示水层深度Di处的时刻;ts表示水面反射回波的峰值位置。

    对该方程分析可知,激光脉冲回波强度随水深深度和漫衰减系数增加呈指数衰减,因此,可以通过水体散射回波所对应的双指数函数反演水体的漫衰减系数。文中假设外业测量过程中飞机航高保持稳定,激光发射器的天底角以及接收视场角保持恒定,且所有损耗因素均得到较好的控制,对水体散射回波强度方程进行简化,并将公式(6)代入公式(5),并简化公式(5),得到简化后的水体散射回波强度方程为

    $${P_c}({t_i}) = {P_e}{{W}}\exp\left( {\frac{{2{K_d}c{t_s}}}{{{n_w}}}} \right){\rm{exp}}\left( {\frac{{ - 2{K_d}c{t_i}}}{{{n_w}}}} \right),$$ (7)

    式中,${{W}} = \displaystyle \frac{{T_{{\rm{atm}}}^2{A_R}{\eta _e}{\eta _R}F{{\left( {1 - {L_S}} \right)}^2}\beta (\phi ){{\cos }^2}\theta }}{{{{\left( {{n_w}H + {D_i}} \right)}^2}}}$

    由于近岸水深通常小于30 m,与飞机航高H(通常为400~500 m)相比,参数DiW的影响可以忽略不计,因此,W可认为是一个常数。但是W涉及未知参数太多,不能直接确定。通过分析发现,公式(3)和公式(7)均可描述为如下指数型函数的形式,

    $$y(t){\rm{ }} = m\exp \left( {nt} \right),$$ (8)

    式中,y(t)为指数型函数;mn分别为指数型函数的系数。

    考虑到传统漫衰减系数提取算法所存在的局限性[9-10],结合2.1节中水体散射回波双指数函数(公式(3))与简化后的水体散射回波强度方程(公式(7))构建模型,从而计算漫衰减系数Kd值。以双指数函数(图3)BC段为例,当bxt<cx时,fc(t)=Pc(t),由公式(8)可知,公式(3)与公式(7)中相应的指数型函数系数n是对应相等的,推导得到

    $$\frac{{{\ln {c_y} - \ln{b_y}} }}{{{c_x} - {b_x}}} = \frac{{ - 2{K_d}_1c}}{{{n_w}}}.$$ (9)

    对公式(9)做方程变换,得到BC段指数所对应的初始漫衰减系数Kd1,同理,计算得到CD段指数所对应的初始漫衰减系数Kd2,如下式所示:

    $${K_d}_1 = \frac{{{n_w}\left( {\ln {c_y} - \ln{b_y}} \right)}}{{2c\left( {{b_x} - {c_x}} \right)}},$$ (10)
    $${K_d}_2 = \frac{{{n_w}\left( {\ln {d_y} - \ln {c_y}} \right)}}{{2c\left( {{c_x} - {d_x}} \right)}}.$$ (11)

    为了保证所得漫衰减系数的准确性,以BC段、CD段的采样时间(Δt1=cxbx与Δt2=dxcx)确定相应权重,根据公式(12),对Kd1Kd2取加权平均,得到Kd1Kd2的加权平均值,将此加权平均值作为该激光测深点最终的漫衰减系数值

    $${K_d} = \frac{{\Delta {t_1}{K_d}_1 + \Delta {t_2}{K_d}_2}}{{\Delta {t_1} + \Delta {t_2}}}.$$ (12)

    为分析所提算法的性能,文中分别应用2013年1月西沙甘泉岛(Optech Aquarius)与2014年12月江苏连云港(Optech CZMIL)两个航次的ALB实验数据,两个航次的测区位置如图4所示。Aquarius和CZMIL均为全波形ALB系统,如表1所示,二者相应的技术参数是不同的。数据采集期间,由于西沙甘泉岛远离大陆,人类活动较少,相比于连云港沿岸海域,甘泉岛海域水质好,海水能见度高。

    图  4  研究区域位置
    a. 西沙甘泉岛;b. 江苏连云港沿岸
    Figure  4.  Location of the study areas
    a. Ganquan Island in the Xisha Archipelago;b. Jiangsu Lianyungang’s seacoast
    表  1  Aquarius与CZMIL的技术参数
    Table  1.  Technical parameters of Aquarius and CZMIL
    航次甘泉岛连云港
    ALB名称AquariusCZMIL
    波段Green: 532 nmIR:1064 nm Green: 532 nm
    脉冲频率70 kHz70 kHz (浅水),10 kHz (深水)
    作业航高300~500 m400 m
    最大探测深度20 mKd·DMAX = 3.75~4.0
    (白天,海底反射率>15%)
    最小探测深度 <0.15 m
    扫描天底角15°20°
      注:“−”表示“无数据”。
    下载: 导出CSV 
    | 显示表格

    通过2.1节分解算法分别对西沙甘泉岛和江苏连云港两个航次的ALB波形数据(西沙甘泉岛8 034个,连云港6 011个)进行处理,分别从两个航次的波形处理结果中随机选取一个进行展示,如图5所示。此外,利用分解结果与去噪后的波形数据对比计算得到均方根误差(Root Mean Square Error, RMSE)和确定系数(Coefficient of Determination, R2)[14-15],从而对文中分解算法性能加以验证,如表2所示。

    图  5  波形处理结果
    a. 连云港波形去噪结果;b. 连云港波形分解结果;c. 甘泉岛波形去噪结果;d. 甘泉岛波形分解结果
    Figure  5.  Waveform processing results
    a. Waveform denoising results of the Lianyungang’s seacoast;b. waveform decomposition results of the Lianyungang’s seacoast;c. waveform denoising results of the Ganquan Island;d. waveform decomposition results of the Ganquan Island
    表  2  分解算法性能指标
    Table  2.  Performance indexes of waveform fitting algorithm
    航次指标最小值最大值平均值
    连云港RMSE7.422315.578010.0197
    R20.86570.99750.9945
    甘泉岛RMSE1.844013.10717.4710
    R20.97990.99970.9947
    下载: 导出CSV 
    | 显示表格

    通过对波形去噪结果分析可知,采用小波去噪算法可将西沙甘泉岛和江苏连云港两个航次的原始回波信号中的噪声信号明显地滤除,去噪效果理想。连云港属于高浑浊海域,其ALB波形的水底反射回波较弱,且由于该测区水深较浅,水体散射回波与水面、水底反射回波叠加在一起,导致波形分解的难度变大,由波形分解结果可知,使用文中所提波形分解算法能够将水体散射回波精确地分解出来,其RMSE、R2的平均值分别为10.019 7、0.994 5。西沙甘泉岛周围海域水质较好,ALB最大探测深度相对较深,因而运用文中分解算法处理甘泉岛波形数据的精度相对更高,其RMSE、R2的平均值分别为7.471 0、0.994 7。综上所述,无论是水质清澈的西沙甘泉岛海域,还是较浑浊的江苏连云港沿岸,采用文中所提基于分层异构模型的波形分解算法均能够将叠加在水面反射回波和水底反射回波中的水体散射回波精确地提取出来,得到水体散射回波所对应的函数(双指数函数),为进一步提取漫衰减系数提供了可靠保障。

    为验证文中漫衰减系数提取算法的有效性,运用Billard B传统方法[9]与文中算法反演漫衰减系数,并作对比分析。文中以包含实验区域深水和浅水数据为原则,随机在西沙甘泉岛、江苏连云港各选取两个区块内的ALB测深数据,应用上述两种算法进行处理,得到漫衰减系数的平均值如表3

    表  3  漫衰减系数提取结果
    Table  3.  Extraction results of diffuse attenuation coefficient
    航次区块传统方法漫衰减系数/m−1文中算法漫衰减系数/m−1
    甘泉岛A0.206 50.215 6
    B0.205 20.245 3
    连云港A0.365 10.378 8
    B0.362 60.381 1
    下载: 导出CSV 
    | 显示表格

    通过表3分析可见,运用文中算法与传统方法在西沙甘泉岛和江苏连云港海域所得的漫衰减系数平均值相近,证明了利用机载LiDAR测深波形进行漫衰减系数提取算法的可行性。此外,通过表4可知[25],基于文中算法所得漫衰减系数分别反演连云港沿岸和甘泉岛海域的水质情况,显示出了较好的区分度。

    表  4  漫衰减系数与水质之间的关系[25]
    Table  4.  Relationship between diffuse attenuation coefficient and water quality[25]
    水质典型海域Kd(532)/m−1
    夏威夷、巴哈马群岛、佛罗里达东海岸<0.08
    较好佛罗里达西海岸、韩国东海岸0.08~0.2
    较浑浊日本沿海0.2~0.4
    非常浑浊密西西比海湾,东海>0.4
    下载: 导出CSV 
    | 显示表格

    为进一步验证文中算法的可行性,文中选择西沙甘泉岛西南角沿岸作为样本区域(图6a),利用其ALB数据进行处理,并生成该区域漫衰减系数空间分布图(图6c)与频率分布直方图(图7)、相应水深的空间分布图(图6b)。

    图  6  甘泉岛西南角样本区域空间分布
    a. 样本区域(红框);b. 水深值;c. 漫衰减系数
    Figure  6.  Distribution of the sample area in the Ganquan Island
    a. Sample area marked by the red box;b. depth;c. diffuse attenuation coefficient
    图  7  甘泉岛西南角样本区域Kd值频率分布直方图
    Figure  7.  Frequency distribution histogram of Kd in the Ganquan Island

    通过图6图7可知,运用文中算法所得漫衰减系数的分布规律与该海域相应水深的分布规律基本保持一致,即当甘泉岛沿岸水越来越来深时,漫衰减系数会越来越小;该区块漫衰减系数主要集中在0~0.5 m−1,其平均值为0.215 6 m−1。浅水处由于潮水的作用,海水会对海底底质产生冲刷,大量泥沙等底质在水中悬浮,因此,浅水处漫衰减系数较高;随着水深的增加,海底底质受潮水的冲刷作用减弱,漫衰减系数逐渐变小。此外,文中还绘制了连云港区块A的漫衰减系数频率分布直方图(图8),结果表明,该区域漫衰减系数主要集中在0.3~0.7 m−1,其平均值为0.378 8 m−1

    图  8  连云港Kd值频率分布直方图
    Figure  8.  Frequency distribution histogram of Kd in the Lianyungang’s seacoast

    综上所述,文中利用ALB测深水体波形提取漫衰减系数的算法具有较强可行性和适应性,可有效提取每一个激光测深点的漫衰减系数。

    结合ALB技术高精度、高分辨率、灵活机动、快速高效的特点,提出了一种基于机载LiDAR测深水体波形的漫衰减系数提取方法。通过分析ALB波形数据,针对水面、水体及水底3部分回波的相应特性,构建分层异构模型,分离得到混叠在水面回波和水底回波之间的水体后向散射回波,进而利用所得水体后向散射回波实现大面积水域漫衰减系数空间分布的精确获取,为快速、精细化分析大面积水域的水体环境变化提供了一种新的解决方案。

  • 图  1  本文方法流程图

    Fig.  1  Flow diagram of the new method in this study

    图  2  机载LiDAR测深波形示意图

    Fig.  2  Airborne LiDAR bathymetric waveform

    图  3  双指数函数分解水体散射回波示意图

    Fig.  3  Decomposition of water column contribution by double-exponential function

    图  4  研究区域位置

    a. 西沙甘泉岛;b. 江苏连云港沿岸

    Fig.  4  Location of the study areas

    a. Ganquan Island in the Xisha Archipelago;b. Jiangsu Lianyungang’s seacoast

    图  5  波形处理结果

    a. 连云港波形去噪结果;b. 连云港波形分解结果;c. 甘泉岛波形去噪结果;d. 甘泉岛波形分解结果

    Fig.  5  Waveform processing results

    a. Waveform denoising results of the Lianyungang’s seacoast;b. waveform decomposition results of the Lianyungang’s seacoast;c. waveform denoising results of the Ganquan Island;d. waveform decomposition results of the Ganquan Island

    图  6  甘泉岛西南角样本区域空间分布

    a. 样本区域(红框);b. 水深值;c. 漫衰减系数

    Fig.  6  Distribution of the sample area in the Ganquan Island

    a. Sample area marked by the red box;b. depth;c. diffuse attenuation coefficient

    图  7  甘泉岛西南角样本区域Kd值频率分布直方图

    Fig.  7  Frequency distribution histogram of Kd in the Ganquan Island

    图  8  连云港Kd值频率分布直方图

    Fig.  8  Frequency distribution histogram of Kd in the Lianyungang’s seacoast

    表  1  Aquarius与CZMIL的技术参数

    Tab.  1  Technical parameters of Aquarius and CZMIL

    航次甘泉岛连云港
    ALB名称AquariusCZMIL
    波段Green: 532 nmIR:1064 nm Green: 532 nm
    脉冲频率70 kHz70 kHz (浅水),10 kHz (深水)
    作业航高300~500 m400 m
    最大探测深度20 mKd·DMAX = 3.75~4.0
    (白天,海底反射率>15%)
    最小探测深度 <0.15 m
    扫描天底角15°20°
      注:“−”表示“无数据”。
    下载: 导出CSV

    表  2  分解算法性能指标

    Tab.  2  Performance indexes of waveform fitting algorithm

    航次指标最小值最大值平均值
    连云港RMSE7.422315.578010.0197
    R20.86570.99750.9945
    甘泉岛RMSE1.844013.10717.4710
    R20.97990.99970.9947
    下载: 导出CSV

    表  3  漫衰减系数提取结果

    Tab.  3  Extraction results of diffuse attenuation coefficient

    航次区块传统方法漫衰减系数/m−1文中算法漫衰减系数/m−1
    甘泉岛A0.206 50.215 6
    B0.205 20.245 3
    连云港A0.365 10.378 8
    B0.362 60.381 1
    下载: 导出CSV

    表  4  漫衰减系数与水质之间的关系[25]

    Tab.  4  Relationship between diffuse attenuation coefficient and water quality[25]

    水质典型海域Kd(532)/m−1
    夏威夷、巴哈马群岛、佛罗里达东海岸<0.08
    较好佛罗里达西海岸、韩国东海岸0.08~0.2
    较浑浊日本沿海0.2~0.4
    非常浑浊密西西比海湾,东海>0.4
    下载: 导出CSV
  • [1] Mobley C D. Light and Water: Radiative Transfer in Natural Waters[M]. San Diego, CA, USA: Academic Press, 1994.
    [2] Philips D M, Abbot R H, Penny M F. Remote sensing of sea water turbidity with an airborne laser system[J]. Journal of Physics D: Applied Physics, 1984, 17(8): 1749−1758.
    [3] Gomes A C C, Bernardo N, Do Carmo A C, et al. Diffuse attenuation coefficient retrieval in CDOM dominated inland water with high chlorophyll-a concentrations[J]. Remote Sensing, 2018, 10(7): 1063.
    [4] Guenther G C. Airborne laser hydrography: system design and performance factors, NOAA Professional Paper Series National Ocean Service 1[R]. Rockville, MD: National Oceanic and Atmospheric Administration, 1985.
    [5] Smart J H, Kwon K H. Comparisons between in situ and remote sensing estimates of diffuse attenuation profiles[C]//CIS Selected Papers: Laser Remote Sensing of Natural Waters: From Theory to Practice. Washington, DC, USA: International Society for Optics and Photonics, 1996, 2964: 100-109.
    [6] Gordon H R, Smith R C, Ronald J, et al. Introduction to ocean optics[C]//Proceedings of SPIE 0489, Ocean Optics VI. Monterey, CA: SPIE, 1980: 14−55.
    [7] 李凯, 童晓冲, 张永生, 等. 黄海、东海区域漫衰减系数光谱遥感反演及激光测深性能评估[J]. 遥感学报, 2015, 19(5): 761−769.

    Li Kai, Tong Xiaochong, Zhang Yongsheng, et al. Inversion of diffuse attenuation coefficient spectral in the Yellow Sea/East China Sea and evaluation of laser bathymetric performance[J]. Journal of Remote Sensing, 2015, 19(5): 761−769.
    [8] 蒋兴伟, 何贤强, 林明森, 等. 中国海洋卫星遥感应用进展[J]. 海洋学报, 2019, 41(10): 113−124.

    Jiang Xingwei, He Xianqiang, Lin Mingsen, et al. Progresses on ocean satellite remote sensing application in China[J]. Haiyang Xuebao, 2019, 41(10): 113−124.
    [9] Billard B, Abbot R H, Penny M F. Airborne estimation of sea turbidity parameters from the WRELADS laser airborne depth sounder[J]. Applied Optics, 1986, 25(13): 2080−2088.
    [10] Ding Kai, Wang Chisheng, Tao Ming, et al. A new algorithm for retrieving diffuse attenuation coefficient based on big LiDAR bathymetry data[C]//International Symposium on Cyberspace Safety and Security. Cham: Springer, 2019: 133−142.
    [11] Yang Anxiu, Wu Ziyin, Yang Fanlin, et al. Filtering of airborne LiDAR bathymetry based on bidirectional cloth simulation[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2020, 163: 49−61.
    [12] 吴自银, 阳凡林, 罗孝文, 等. 高分辨率海底地形地貌−探测处理理论与技术[M]. 北京: 科学出版社, 2017.

    Wu Ziyin, Yang Fanlin, Luo Xiaowen, et al. High-Resolution Submarine Topography: Theory and Technology for Surveying and Post-Processing[M]. Beijing: Science Press, 2017.
    [13] 吴自银, 阳凡林, 李守军, 等. 高分辨率海底地形地貌−可视计算与科学应用[M]. 北京: 科学出版社, 2017.

    Wu Ziyin, Yang Fanlin, Li Shoujun, et al. High-Resolution Submarine Topography: Visual Computation and Scientific Applications[M]. Beijing: Science Press, 2017.
    [14] Su Dianpeng, Yang Fanlin, Ma Yue, et al. Classification of coral reefs in the South China Sea by combining airborne LiDAR bathymetry bottom waveforms and bathymetric features[J]. IEEE Transactions on Geoscience and Remote Sensing, 2019, 57(2): 815−828.
    [15] 亓超, 宿殿鹏, 王贤昆, 等. 基于分层异构模型的机载激光测深波形拟合算法[J]. 红外与激光工程, 2019, 48(2): 206004.

    Qi Chao, Su Dianpeng, Wang Xiankun, et al. Fitting algorithm for airborne laser bathymetric waveforms based on layered heterogeneous model[J]. Infrared and Laser Engineering, 2019, 48(2): 206004.
    [16] Rees W G. Physical Principles of Remote Sensing[M]. Cambridge: Cambridge University Press, 2012.
    [17] Abdallah H, Baghdadi N, Bailly J S, et al. Wa-LiD: a new LiDAR simulator for waters[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(4): 744−748.
    [18] 孙磊, 张志利, 谭立龙, 等. 采用小波阈值的时变光栅莫尔信号去噪方法[J]. 红外与激光工程, 2010, 39(3): 576−580. doi: 10.3969/j.issn.1007-2276.2010.03.040

    Sun Lei, Zhang Zhili, Tan Lilong, et al. Denoising method of dynamic grating Moiré signal based on wavelet threshold[J]. Infrared and Laser Engineering, 2010, 39(3): 576−580. doi: 10.3969/j.issn.1007-2276.2010.03.040
    [19] 刘永明, 邓孺孺, 秦雁, 等. 机载激光雷达测深数据处理与应用[J]. 遥感学报, 2017, 21(6): 982−995.

    Liu Yongming, Deng Ruru, Qin Yan, et al. Data processing methods and applications of airborne LiDAR bathymetry[J]. Journal of Remote Sensing, 2017, 21(6): 982−995.
    [20] Hofton M A, Minster J B, Blair J B. Decomposition of laser altimeter waveforms[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(4): 1989−1996.
    [21] Ding Kai, Li Qingquan, Zhu Jiasong, et al. An improved quadrilateral fitting algorithm for the water column contribution in airborne bathymetric lidar waveforms[J]. Sensors, 2018, 18(2): 552.
    [22] Abdallah H, Bailly J S, Baghdadi N N, et al. Potential of space-borne LiDAR sensors for global bathymetry in coastal and inland waters[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2013, 6(1): 202−216.
    [23] Abady L, Bailly J S, Baghdadi N, et al. Assessment of quadrilateral fitting of the water column contribution in lidar waveforms on bathymetry estimates[J]. IEEE Geoscience and Remote Sensing Letters, 2014, 11(4): 813−817.
    [24] Klett J D. Stable analytical inversion solution for processing lidar returns[J]. Applied Optics, 1981, 20(2): 211−220.
    [25] Ramnath V, Feygels V, Kalluri H, et al. CZMIL (Coastal Zone Mapping and Imaging Lidar) bathymetric performance in diverse littoral zones[C]//OCEANS 2015-MTS/IEEE Washington. Washington, DC, USA: IEEE, 2015: 1−10.
  • 期刊类型引用(3)

    1. 刘超. 基于卷积神经网络的机载激光海洋测深波形分类. 自动化技术与应用. 2024(03): 48-51+56 . 百度学术
    2. 宿殿鹏,黄昱,阳凡林,赵荻能,杨安秀,刘骄阳. 顾及特征优选的机载LiDAR测深海底点云底质分类. 遥感学报. 2023(09): 2219-2228 . 百度学术
    3. 张一衡,于孝林,亓超,宿殿鹏,王志良,任国贞. 基于海底回波增强的机载LiDAR测深水体旅行时提取算法. 海洋学报. 2023(12): 145-155 . 本站查看

    其他类型引用(1)

  • 加载中
图(8) / 表(4)
计量
  • 文章访问数:  315
  • HTML全文浏览量:  119
  • PDF下载量:  39
  • 被引次数: 4
出版历程
  • 收稿日期:  2019-10-14
  • 修回日期:  2020-06-28
  • 网络出版日期:  2020-12-30
  • 刊出日期:  2021-01-25

目录

/

返回文章
返回