激光雷达体回波信号仿真方法

申请号 CN201510975332.6 申请日 2015-12-23 公开(公告)号 CN105372641A 公开(公告)日 2016-03-02
申请人 国家海洋局第一海洋研究所; 发明人 马毅; 崔廷伟; 张震; 张靖宇;
摘要 激光雷达 水 体 回波 信号 仿真方法,根据接收到的激光雷达水体回波信号及其参数,分别计算出水表面回波信号、水体回波信号、水底底部回波信号、背景噪声和为接收器内部噪声;将上述参数输入激光雷达水体回波信号仿真模型,输出为激光雷达接收的总回波信号。本 发明 完善了包含纯 海水 、浮游 植物 、悬浮物和黄色物质的水体组分吸收和散射特性的参数,以构造物理意义明确的激光雷达回波信号仿真模型,用于解决 现有技术 中仅对激光雷达回波 波形 处理的问题,便于直观清晰描述激光雷达回波信号,为激光雷达测深系统的设计和应用提供技术 支撑 。
权利要求

1.激光雷达体回波信号仿真方法,其特征在于,根据接收到的激光雷达水体回波信号及其参数,分别计算出水表面回波信号、水体回波信号、水底底部回波信号、背景噪声和接收器内部噪声;
①、所述水表面回波信号Ps(t)为水表面瞬时回波能量与发射脉冲的卷积;包括首先根据水表面的双向反射率分布函数确定水表面损失因素LS,
其中,ks镜面反射系数,kd为漫反射系数,O为水面BRDF几何衰减系数,Fr为菲涅尔光学反射函数, m=cos(θ),g2=m2+n2-1,n为海水折射
率;
其次,将水表面损失因素LS输入到水表面瞬时回波能量模型中,输出水表面瞬时回波能量;
其中,Ps为发射能量, 为双程大气损失,AR传感器接受面积,ηe为脉冲发射光学效率,ηR为接收光学效率,θ为脉冲入射,H为传感器的高度,LS为水表面损失因素;
然后,再将水表面瞬时回波能量与发射脉冲进行卷积,
其中,ts=2H/ccos(θ),ts为脉冲从发射到达水表面的双程时间;
②、所述水体回波信号Pc(t)为不同水深段瞬时水体回波能量与发射脉冲的卷积;首先根据输入的激光雷达水体回波信号参数,经计算输出包括黄色物质、浮游植物、非藻类颗粒物等贡献总和的水体浑浊度参数(a(λ),b(λ)),
a(λ)=aw(λ)+ay(λ)+aph(λ)+as(λ) (4)
其中,aw(λ)为纯水的吸收系数,aph(λ)为浮游植物的吸收系数;其中黄色物质的吸收系数ay(λ)=ay(λ0)exp(-Sg(λ-λ0)),非藻类颗粒物的吸收系数as(λ)=as(λ0)exp(-Sd(λ-λ0)),散射系数 其
中,Sg为黄色物质的光谱吸收斜率,Sd为非藻类颗粒物吸收系数曲线斜率,bw(λ)纯水的散射系数,bp(λ)颗粒物的散射系数,C为叶绿素浓度;
其次,将水体浑浊度参数(a(λ),b(λ))输入漫衰减系数k模型中,经计算输出漫衰减系数k:
其中,ω0=b(λ)/c(λ),c(λ)=a(λ)+b(λ),ω0为单次散射反照率,a(λ)、b(λ)、c(λ)分别为λ波段处水体的光束吸收系数、散射系数、总衰减系数;
最后,输入的漫衰减系数k到水体回波信号模型中,输出对应水深时的水体回波信号,其中,β(φ)为散射相函数,nw为水的折射率,tc=ts+(2Z/cwcosθw),cw为光在水中的速度,θw为脉冲在水中的入射角,漫衰减系数k;
③、所述底部回波信号Pb(t)是水底部瞬时回波能量与发射脉冲的卷积,其中,Rb为底部反照率,Z为底部水深,tc=ts+(2Z/cwcosθw);
④、所述背影噪声Pbg(t)是均值为0、标准差为1的高斯白噪声与瞬时回波信号Pbg的卷积,所述瞬时回波信号 其中,Is为水体反射的太阳辐射
γγ为接收器昏暗比,△ λ为接收器最优滤波带宽;
⑤、所述接收器内部噪声PN(t)是均值为0,标准差为σN(t)的正态分布参数,所述标准差 其中,Pext(t)=Ps(t)+Pc(t)+Pb(t)+Pbg(t),e为电荷
(1.6×10-19C),B为探测器的电子带宽,G为增益噪声,Id为暗电流,Rλ为响应度。
2.根据权利要求1所述的激光雷达水体回波信号仿真方法,其特征在于,将水表面回波信号、水体回波信号、水底底部回波信号、背景噪声和接收器内部噪声数据输入到激光雷达水体回波信号仿真模型,输出激光雷达接收的总回波信号PT(t),
PT(t)=Ps(t)+Pc(t)+Pb(t)+Pbg(t)+PN(t) (8)
其中,t为时间,PT(t)是激光雷达接收的总回波信号,Ps(t)为水表面回波信号,Pc(t)为水体回波信号,Pb(t)为水底底部回波信号,Pbg(t)为背景噪声,PN(t)为接收器内部噪声。

说明书全文

激光雷达体回波信号仿真方法

技术领域

[0001] 本发明涉及一种遥感水深探测方法,尤其涉及通过辐射传输模型进行水体回波信号定量描述的激光雷达水体回波信号仿真方法。

背景技术

[0002] 激光雷达测深系统集成激光、微光探测、高精度定位及高速数据处理等技术于一体,是当今国际技术发展的前沿。激光雷达测深是一种主动式遥感测深手段,由于蓝绿光在水中的衰减系数最小、穿透能最强的特性,激光雷达测深系统中一般会选用532nm蓝绿激光作为发射波段,其测深的基本原理,是发射器以一定的度向水面发射激光,激光经过水表面反射作用一部分返回到测深系统的接收器,另一部分经过水表面透射到水底或水中再返回到接收器中,利用两者的时间差计算水深。
[0003] 分析和研究激光雷达从发射脉冲到接收回波信号的物理过程,可以对激光雷达水体回波信号的物理机制实现仿真模拟,为激光雷达测深系统的设计提供技术支撑。目前在已公开的Wa-Lid模型中,但该模型水体三要素的吸收、散射等公式未声明,缺乏对波形物理参数认识,进而无法准确描述激光雷达水体回波信号。

发明内容

[0004] 本发明提供一种激光雷达水体回波信号仿真方法,完善了包含纯海水、浮游植物、悬浮物和黄色物质的水体组分吸收和散射特性的参数;进而,可以构造出明确的激光雷达回波信号仿真模型,用于解决现有技术中仅对激光雷达回波波形处理的问题,便于直观清晰描述激光雷达回波信号。
[0005] 激光雷达水体回波信号仿真方法,根据接收到的激光雷达水体回波信号及其参数,分别计算出水表面回波信号、水体回波信号、水底底部回波信号、背景噪声和接收器内部噪声;
[0006] ①、所述水表面回波信号Ps(t)为水表面瞬时回波能量与发射脉冲的卷积;包括首先根据水表面的双向反射率分布函数确定水表面损失因素LS,
[0007]
[0008] 其中,ks镜面反射系数,kd为漫反射系数,O为水面BRDF几何衰减系数,Fr为菲涅尔光学反射函数, g2=m2+n2-1,n为海水折射率;
[0009] 其次,将水表面损失因素LS输入到水表面瞬时回波能量模型中,输出水表面瞬时回波能量;
[0010]
[0011] 其中,Ps为发射能量, 为双程大气损失,AR传感器接受面积,ηe为脉冲发射光学效率,ηR为接收光学效率,θ为脉冲入射角,H为传感器的高度,LS为水表面损失因素;
[0012] 然后,再将水表面瞬时回波能量与发射脉冲进行卷积,
[0013]
[0014] 其中,ts=2H/ccos(θ),ts为脉冲从发射到达水表面的双程时间;
[0015] ②、所述水体回波信号Pc(t)为不同水深段瞬时水体回波能量与发射脉冲的卷积;首先根据输入的激光雷达水体回波信号参数,经计算输出包括黄色物质、浮游植物、非藻类颗粒物等贡献总和的水体浑浊度参数(a(λ),b(λ)),
[0016] a(λ)=aw(λ)+ay(λ)+aph(λ)+as(λ) (4)
[0017] 其中,aw(λ)为纯水的吸收系数,aph(λ)为浮游植物的吸收系数;其中黄色物质的吸收系数ay(λ)=ay(λ0)exp(-Sg(λ-λ0)),非藻类颗粒物的吸收系数as(λ)=as(λ0)exp(-Sd(λ-λ0)),散射系数 其中,Sg为黄色物质的光谱吸收斜率,Sd为非藻类颗粒物吸收系数曲线斜率,bw(λ)纯水的散射系数,bp(λ)颗粒物的散射系数,C为叶绿素浓度;
[0018] 其次,将水体浑浊度参数(a(λ),b(λ))输入漫衰减系数k模型中,经计算输出漫衰减系数k:
[0019]
[0020] 其中,ω0=b(λ)/c(λ),c(λ)=a(λ)+b(λ),ω0为单次散射反照率,a(λ)、b(λ)、c(λ)分别为λ波段处水体的光束吸收系数、散射系数、总衰减系数;
[0021] 最后,输入的漫衰减系数k到水体回波信号模型中,输出对应水深时的水体回波信号,
[0022]
[0023] 其中,β(φ)为散射相函数,nw为水的折射率,tc=ts+(2Z/cwcosθw),cw为光在水中的速度,θw为脉冲在水中的入射角,漫衰减系数k;
[0024] ③、所述底部回波信号Pb(t)是水底部瞬时回波能量与发射脉冲的卷积,[0025]
[0026] 其中,Rb为底部反照率,Z为底部水深,tc=ts+(2Z/cwcosθw);
[0027] ④、所述背影噪声Pbg(t)是均值为0、标准差为1的高斯白噪声与瞬时回波信号Pbg的卷积,所述瞬时回波信号 其中,Is为水体反射的太阳辐射,γγ为接收器昏暗比,△ λ为接收器最优滤波带宽;
[0028] ⑤、所述接收器内部噪声PN(t)是均值为0,标准差为σN(t)的正态分布参数,所述标准差 其中,Pext(t)=Ps(t)+Pc(t)+Pb(t)+Pbg(t),e为电荷(1.6×10-19C),B为探测器的电子带宽,G为增益噪声,Id为暗电流,Rλ为响应度。
[0029] 最后,将水表面回波信号、水体回波信号、水底底部回波信号、背景噪声和接收器内部噪声数据输入到激光雷达水体回波信号仿真模型,输出激光雷达接收的总回波信号PT(t),
[0030] PT(t)=Ps(t)+Pc(t)+Pb(t)+Pbg(t)+PN(t) (8)
[0031] 其中,t为时间,PT(t)是激光雷达接收的总回波信号,Ps(t)为水表面回波信号,Pc(t)为水体回波信号,Pb(t)为水底底部回波信号,Pbg(t)为背景噪声,PN(t)为接收器内部噪声。
[0032] 激光雷达水体回波信号模型能仿真水表面、水体、海底等反射激光雷达波束所形成的回波波形,而且充分考虑了太阳噪声和传感器内部噪声的影响,因此,LiDAR测深系统接收的回波信号可认为是关于时间函数的五大部分回波信号的叠加
[0033] 本发明的有益效果:
[0034] 本发明对现有技术中的Wa-Lid模型,主要针对其中的水色三要素,包括黄色物质、沉积物、浮游植物等对激光波束的吸收和散射系数进行补充完善;进而,可以构造出明确的激光雷达回波信号辐射传输仿真模型,直观了解和掌握波形物理参数,以便更加清晰地描述激光雷达水体回波信号以及准确地处理激光雷达回波信号,为激光雷达测深系统的设计提供技术支撑。附图说明
[0035] 图1是本发明实施例中以532nm激光形成的模拟LiDAR波形图;
[0036] 具体,1a对应于水表面回波信号,1b对应于水体回波信号,1c对应于水底底部回波信号,1d对应于背景噪声,1e对应于接收器内部噪声,1f对应于激光激光雷达接收的总回波信号。

具体实施方式

[0037] 为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
[0038] 结合附图1对本发明具体实施方式进一步详细描述,激光雷达水体回波信号仿真方法,根据接收到的激光雷达水体回波信号及其参数,分别计算出水表面回波信号、水体回波信号、水底底部回波信号、背景噪声和接收器内部噪声;
[0039] ①、所述水表面回波信号Ps(t)为水表面瞬时回波能量与发射脉冲的卷积;包括首先根据水表面的双向反射率分布函数确定水表面损失因素LS,
[0040]
[0041] 其中,ks镜面反射系数,kd为漫反射系数,O为水面BRDF几何衰减系数,Fr为菲涅2 2 2
尔光学反射函数, g=m +n-1,n为海水
折射率;
[0042] 其次,将水表面损失因素LS输入到水表面瞬时回波能量模型中,输出水表面瞬时回波能量;
[0043]
[0044] 其中,Ps为发射能量, 为双程大气损失,AR传感器接受面积,ηe为脉冲发射光学效率,ηR为接收光学效率,θ为脉冲入射角,H为传感器的高度,LS为水表面损失因素;
[0045] 然后,再将水表面瞬时回波能量与发射脉冲进行卷积,
[0046]
[0047] 其中,ts=2H/ccos(θ),ts为脉冲从发射到达水表面的双程时间;
[0048] ②、所述水体回波信号Pc(t)为不同水深段瞬时水体回波能量与发射脉冲的卷积;首先根据输入的激光雷达水体回波信号参数,经计算输出包括黄色物质、浮游植物、非藻类颗粒物等贡献总和的水体浑浊度参数(a(λ),b(λ)),
[0049] a(λ)=aw(λ)+ay(λ)+aph(λ)+as(λ) (4)
[0050] 其中,aw(λ)为纯水的吸收系数,aph(λ)为浮游植物的吸收系数;其中黄色物质的吸收系数ay(λ)=ay(λ0)exp(-Sg(λ-λ0)),非藻类颗粒物的吸收系数as(λ)=as(λ0)exp(-Sd(λ-λ0)),散射系数 其中,Sg为黄色物质的光谱吸收斜率,Sd为非藻类颗粒物吸收系数曲线斜率,bw(λ)纯水的散射系数,bp(λ)颗粒物的散射系数,C为叶绿素浓度;
[0051] 其次,将水体浑浊度参数(a(λ),b(λ))输入漫衰减系数k模型中,经计算输出漫衰减系数k:
[0052]
[0053] 其中,ω0=b(λ)/c(λ),c(λ)=a(λ)+b(λ),ω0为单次散射反照率,a(λ)、b(λ)、c(λ)分别为λ波段处水体的光束吸收系数、散射系数、总衰减系数;
[0054] 最后,输入的漫衰减系数k到水体回波信号模型中,输出对应水深时的水体回波信号,
[0055]
[0056] 其中,β(φ)为散射相函数,nw为水的折射率,tc=ts+(2Z/cwcosθw),cw为光在水中的速度,θw为脉冲在水中的入射角,漫衰减系数k;
[0057] ③、所述底部回波信号Pb(t)是水底部瞬时回波能量与发射脉冲的卷积,[0058]
[0059] 其中,Rb为底部反照率,Z为底部水深,tc=ts+(2Z/cwcosθw);
[0060] ④、所述背影噪声Pbg(t)是均值为0、标准差为1的高斯白噪声与瞬时回波信号Pbg的卷积,所述瞬时回波信号 其中,Is为水体反射的太阳辐射,γγ为接收器昏暗比,△ λ为接收器最优滤波带宽;
[0061] ⑤、所述接收器内部噪声PN(t)是均值为0,标准差为σN(t)的正态分布参数,所述标准差 其中,Pext(t)=Ps(t)+Pc(t)+Pb(t)+Pbg(t),e为电荷(1.6×10-19C),B为探测器的电子带宽,G为增益噪声,Id为暗电流,Rλ为响应度。
[0062] 本实施例中的激光雷达水体回波信号模型能仿真水表面、水体、海底等反射激光雷达波束所形成的回波波形,而且考虑了太阳噪声和传感器内部噪声的影响。LiDAR测深系统接收的回波信号可以确定为和时间函数有关的上述五部分回波信号的叠加。因此,将水表面回波信号、水体回波信号、水底底部回波信号、背景噪声和接收器内部噪声数据输入到激光雷达水体回波信号仿真模型,输出激光雷达接收的总回波信号PT(t),[0063] PT(t)=Ps(t)+Pc(t)+Pb(t)+Pbg(t)+PN(t) (8)
[0064] 其中,t为时间,PT(t)是激光雷达接收的总回波信号,Ps(t)为水表面回波信号,Pc(t)为水体回波信号,Pb(t)为水底底部回波信号,Pbg(t)为背景噪声,PN(t)为接收器内部噪声。
[0065] 具体而言,本实施例选用HawkEye机载激光雷达系统,具体的测深系统参数和环境参数见表1:
[0066] 表1 HawkEye测深系统参数和环境参数
[0067]
[0068]
[0069] 选取沉积物浓度S=2.8mg/l,C=4mg/m3,r=0.2,通过设置上述几个参数以及表1的参数,对激光雷达水体信号进行仿真,当水深为2m时形成的回波波形,如图1所示。
[0070] 本发明构造了物理意义明确的激光雷达回波信号辐射传输仿真模型,直观了解和掌握波形物理参数,激光雷达水体回波信号进行了清晰描述,为激光雷达测深系统的设计提供技术支撑。
[0071] 本发明未详尽描述的技术内容均为公知技术。
QQ群二维码
意见反馈