CN107831485A - 船载多视场激光雷达探测多个水体光学特征参数的方法 - Google Patents
船载多视场激光雷达探测多个水体光学特征参数的方法 Download PDFInfo
- Publication number
- CN107831485A CN107831485A CN201710977199.7A CN201710977199A CN107831485A CN 107831485 A CN107831485 A CN 107831485A CN 201710977199 A CN201710977199 A CN 201710977199A CN 107831485 A CN107831485 A CN 107831485A
- Authority
- CN
- China
- Prior art keywords
- water body
- scattering
- laser radar
- water
- coefficient
- 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
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
- 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)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明涉及船载多视场激光雷达探测多个水体光学特征参数的方法,包括:利用窄视场探测回波反演激光雷达衰减系数估计水体多次散射特征参数;基于多次前向‑单次后向散射雷达模型,建立不同视场接收时回波比值方程组,结合参数m值估计水体前向散射系数;基于Walker‑Mclean模型采用最小二乘法计算水体的吸收和散射系数;通过极限深度假设反演不同水层的向下漫射衰减系数。本发明利用视场变化影响激光雷达回波的多次散射增量,结合多次前向‑单次后向和Walker‑Mclean模型,引入多次散射因子,基于船载探测时不受参数m值影响,通过单水层极限深度假设,实现水体多个光学特征参数反演,提高了方法的水体光学特性探测能力。
Description
技术领域
本发明属于海洋光学技术领域,涉及一种水体光学特性的探测方法,具体地说是船载多视场激光雷达探测多个水体光学特征参数的方法。
背景技术
水体光学特性,是针对水域初级生产力调查、水质探测和生物资源开发等方面的重要研究对象。在诸多鉴别手段中,激光雷达可进行水域的大尺度、快速探测,能够应用于船载平台进行走航式测量,具备耗时短、抗实验环境干扰、测量效率高等优点。
目前,用于水体光学特性探测的多是基于弹性散射机制的单一视场激光雷达,仅能探测激光波长处单一的海水光学特征参数,即水体的吸收系数、衰减系数、或向下漫射衰减系数。美国NAVAIR支持的LOCO计划下完成的SeaTROLL系统,是目前成熟的激光雷达水体光学特性探测系统,主要用于测量海洋水体的光学特性结构特点和浮游植物的分层结构,但只能获得两个水体光学特征参数,即水体的衰减系数和向下漫射衰减系数。因此,利用激光雷达探测多个水体光学特征参数的方法有待于发展,以实现快速、高效地进行水体光学特性结构的探测,解决激光雷达探测水体光学特性的技术中特征参数分析不全面的问题。
发明内容
本发明针对现有激光雷达测量技术中水体光学特征参数探测不全面的问题,提供了一种船载多视场激光雷达探测多个水体光学特征参数的方法,利用接收视场变化对激光雷达回波中多次散射增量的影响,并结合多次前向-单次后向雷达模型和Walker-Mclean模型,实现水体的多个水体光学特征参数反演,提高了方法的水体光学特性探测能力。
本发明为实现上述目的所采用的技术方案是:船载多视场激光雷达探测多个水体光学特征参数的方法,包括以下步骤:
步骤1:利用船载激光雷达窄视场的探测回波反演激光雷达衰减系数,以估计水体的多次散射特征参数m的估计值;
步骤2:利用多次散射特征参数m的估计值,基于多次前向-单次后向散射雷达模型,建立不同视场接收时激光雷达回波比值的方程组,以计算水体前向散射系数bf;
步骤3:基于Walker-Mclean模型,获得不同视场的激光雷达回波比,并采用最小二乘法模型计算水体的吸收系数a和多次散射项g(z,τ);
步骤4:在Walker-Mclean模型中,进行极限深度假设,反演不同水层的向下漫射衰减系数Kd。
所述步骤1中的多次散射特征参数m的估计值取6~8之间的整数值。
所述步骤2包括:
步骤2.1:将多次散射特征参数m的估计值代入多次前向-单次后向散射雷达模型,得到多次散射项:
其中,m为多次散射特征参数,2θrcvr为接收视场角,n水体折射率,H为激光雷达距离水面的高度,bf为水体前向散射系数,z为探测水深,2θlas为激光光束发散角,x为光子散射后位置向量的空间频率,J1为1阶贝塞尔函数,rlas和Rrcvr分别为探测激光束和望远镜接收通光面的半径;
步骤2.2:基于多次前向-单次后向散射雷达模型,计算船载激光雷达的回波强度:
其中,Q为激光能量,Arcvr为接收视场入瞳面积,T2为激光波长在大气中和气海界面的双程衰减,v为光速,τd为脉冲宽度,n为水体折射率,z为探测深度,H为激光雷达与水面间距离,β(π)为水体后向180°散射系数;a为水体吸收系数,bb为后向散射系数;F(z)为多次散射项;
步骤2.3:利用相同视场入瞳面积、不同的有效接收视场i,计算船载激光雷达的回波强度:
Pi(z)=Pwater(ψi,z)=M(z)F(ψi,z),
获得回波比值方程组:
其中, 对于同一水体,该项与视场无关;
求解回波比值方程组,计算积分项
中的未知数:水体前向散射系数bf和多次散射特征参数m的准确值。
所述步骤3包括:
步骤3.1:先将水体前向散射系数bf近似为散射系数b,代入Walker-Mclean模型,则有:
其中,b为散射系数,τ为多次散射时间,g(z,τ)为对多次散射时间τ统计的Gamma分布概率密度函数,简称多次散射项;
其中,
<θ2>为散射角均方值,
θ=arccos[(0.142m2-1)/0.132m2],m为多次散射特征参数m的准确值;
步骤3.2:利用不同视场角θ1、θ2下的激光雷达回波比值,由于散射系数b取近似值,并去除水体后向180°散射系数β(π)的计算影响:
求解得到水体吸收系数a;
步骤3.3:采用最小二乘法模型,
结合实验测量回波比值和Walker-Mclean模型所计算的回波比值多次迭代计算水体的多次散射项g(z,τ):
其中,θi为不同于θ3的视场角。
所述步骤4包括:
步骤4.1:首先将Walker-Mclean模型的深度z趋近于无穷大,得到基于极限深度假设的Walker-Mclean模型:
其中,~为近似符号,
步骤4.2:将计算得到的水体吸收系数a和多次散射项g(z,τ)值代入基于极限深度假设的Walker-Mclean模型,利用limz→∞exp(-2KLidarz)=exp(-2Kdz)计算不同水层处向下漫射衰减系数Kd的趋近值。
本发明具有以下有益效果及优点:
1、通过船载多视场激光雷达探测方式,利用接收视场变化对激光雷达回波中多次散射增量的影响,可实现多个水体光学特征参数的反演,与现有激光雷达水体光学探测技术相比,获得了更为全面的水体光学特征参数。
2、通过引入多次散射因子,利用船载激光雷达探测不受多次散射特征参数m值影响的特点,可实现水体前向散射系数的反演,并结合多次前向-单次后向雷达模型和Walker-Mclean模型,通过单水层极限深度假设,可实现水体的多个水体光学特征参数反演,提高了船载激光雷达探测水体光学特性的效率。
3、本发明所述船载多视场激光雷达探测多个水体光学特征参数的方法具有探测参数全面、快速、有效的特点,尤其适用于水域光学特性的走航式测量应用。
附图说明
图1为本发明船载多视场激光雷达探测多个水体光学特征参数的方法实施流程图;
图2为不同接收视场探测情况下,清洁水体中激光雷达后向散射回波的反演值。
具体实施方式
下面结合附图和具体实施方式对本发明的技术方案作进一步详细的说明。
一种船载多视场激光雷达探测多个水体光学特征参数的方法,包括多次散射特征参数m值估计,水体前向散射系数反演,吸收和散射系数计算,以及利用极限深度假设反演向下漫射衰减系数。其中,水体前向散射系数的反演,是利用参数m的估值,基于多次前向-单次后向散射雷达模型,利用相同视场入瞳面积、不同的有效接收视场,计算船载激光雷达的回波强度Pi(z)=Pwater(ψi,z)=M(z)F(ψi,z),建立回波比值方程组,来计算未知数bf和m,求解获得其准确值。其中, 对于同一水体,该项与视场无关;bf为前向散射系数,m为多次散射特征参数。本发明结合多次前向-单次后向雷达模型和Walker-Mclean模型,利用船载激光雷达探测不受多次散射特征参数m值影响的特点,提高了船载激光雷达探测水体光学特性的效率、以及参数探测的全面性。
如上所述的船载多视场激光雷达探测方法,如果获得窄视场的激光雷达回波信号,从而估计参数m值,则进一步采用所述多次前向-单次后向雷达模型进行水体前向散射系数bf的计算。
如上所述的船载多视场激光雷达探测方法,水体前向散射系数bf作为散射系数b的近似值,代入Walker-Mclean模型中,进行水体的吸收系数和多次散射项精确值的计算。
如上所述的船载多视场激光雷达探测方法,基于Walker-Mclean模型和无穷大水深假设,计算不同水层处向下漫射衰减系数的趋近值。
如上所述的船载多视场激光雷达探测方法,为了保证参数计算的准确,最大限度地去除太阳背景及其他噪声信号的影响,用于计算的激光雷达回波为多次测量的原始回波平均、并进行探测能量归一化处理后,得到激光后向散射回波的衰减信号。
本实施例具体步骤如下:
首先,简要说明本发明所述船载多视场激光雷达探测多个水体光学特征参数方法的立足依据:常规激光雷达探测水体光学特性时仅能获得单一的光学特征参数,尤其对于船载激光雷达,只能利用激光后向散射回波反演值与水体的衰减系数或吸收系数做相关性比较。因此,基于上述问题,利用多视场的探测结构,结合不同的激光雷达模型反演多个水体光学特征参数,可以解决船载激光雷达应用中水体光学特征参数探测不全面的缺陷。
图1示出了本发明船载多视场激光雷达探测多个水体光学特征参数的方法实施流程图。
如图1所示,该实施例的船载多视场激光雷达探测多个水体光学特征参数方法,包括多次散射特征参数m值估计,水体前向散射系数反演,吸收和衰减系数计算,以及利用极限深度假设反演向下漫射衰减系数。
第一步,利用船载激光雷达窄视场的探测回波反演激光雷达衰减系数,从而判断水体的浑浊程度,以对参数m值进行估计。根据窄视场的探测回波NWwater,激光雷达衰减系数KLidar=Δlog[NWwater×(nH+z)2]/Δz。m为小角度散射情况下的特性参数:
而其中<θ2>为散射角均方值。该近似关系满足前向性强的散射物特性(水体特性表现为对于船载平台探测水体情况来讲,m取整数值在6~8之间,如:当KLidar小于0.1m-1时,m可取估值6;当KLidar大于2m-1时,m取估值8。
第二步,利用估计的参数m值,基于多次前向-单次后向散射雷达模型,建立不同视场接收时激光雷达回波比值的方程组,以计算水体前向散射系数bf。基于多次前向-单次后向散射雷达模型,得到多次散射项
其中,bf为前向散射系数,m为多次散射特征参数,m为多次散射特征参数,2θrcvr为接收视场角,n水体折射率,H为激光雷达距离水面的高度,bf为水体前向散射系数,z为探测水深,2θlas为激光光束发散角,x为光子散射后位置向量的空间频率,J1为1阶贝塞尔函数,rlas和Rrcvr分别为探测激光束和望远镜接收通光面的半径。
则,船载激光雷达的回波强度
其中,bb为后向散射系数,Q为激光能量,Arcvr为接收视场入瞳面积,T2为激光波长在大气中和气海界面的双程衰减,v为光速,τd为脉冲宽度,n为水体折射率,z为探测深度,H为激光雷达与水面间距离,β(π)为水体后向180°散射系数。
设定相同视场入瞳面积、不同的有效接收视场i,计算船载激光雷达的回波强度Pi(z)=Pwater(ψi,z)=M(z)F(ψi,z),获得回波比值方程组
其中,积分项中包含未知数bf和m,求解获得其准确值。
第三步是基于Walker-Mclean模型,构建激光雷达回波比,并采用最小二乘法模型进行水体的吸收和散射系数计算。先将水体前向散射系数bf近似为散射系数b,代入Walker-Mclean模型,
再利用激光雷达回波比值,去除水体后向散射系数β(π)的计算影响,得到水体吸收系数a,采用最小二乘法模型
结合实验测量回波比值和Walker-Mclean模型所计算的回波比值,多次迭代计算水体的多次散射项g(z,τ),其中,为实验测量回波比值,为Walker-Mclean模型所计算的回波比值。
最后,在Walker-Mclean模型中,进行极限深度假设,反演不同水层的向下漫射衰减系数。将Walker-Mclean模型的深度z趋近于无穷大,激光雷达视场既能接收到单次散射,也能接收到多次散射的光子,因此,
将计算得到的水体吸收系数a和多次散射项g(z,τ)值,代入基于极限深度假设的Walker-Mclean模型,利用limz→∞exp(-2KLidarz)=exp(-2Kdz)计算得到不同水层处向下漫射衰减系数Kd的趋近值。
该实例中,在不同水体类型下,利用532nm激光探测实验水体,将船载平台上多个接收视场情况下采集到的激光后向散射回波进行去除背景噪声、归一化处理后,获得不同水体类型的船载多视场激光雷达探测信号,共包括5000组激光雷达回波数据。实验水体分别为清洁水体和浑浊水体,实验过程中采用原位探测设备同步测量水体的多个水体光学特征参数。
利用窄视场8.7mrad情况下的激光后向散射回波衰减系数KLidar=Δlog[NWwater×(nH+z)2]/Δz,估计清洁水体和浑浊水体的参数m值,分别为7和8。
然后,基于多次前向-单次后向散射雷达模型,利用参数m估值分别建立不同视场接收时激光雷达回波比值的方程组
得到清洁水体和浑浊水体的前向散射系数bf分别为0.378和2.045。
再基于Walker-Mclean模型,将水体前向散射系数bf近似为散射系数b后代入,构建激光雷达回波比:
其中,
θ=arccos[(0.142m2-1)/0.132m2]
采用最小二乘法模型进行
结合实验测量回波比值和Walker-Mclean模型所计算的回波比值计算清洁水体和浑浊水体的吸收和散射系数,分别为(a=0.066m-1,b=0.369m-1)和(a=0.162m-1,b=2.109m-1)。
不同接收视场探测情况下,清洁水体中激光雷达后向散射回波的反演值如图2所示,反演值随接收视场变化而明显不同。
最后,将Walker-Mclean模型的深度z趋近于无穷大,反演水体的向下漫射衰减系数得到清洁水体和浑浊水体的漫射衰减系数0.073m-1和0.268m-1。将水体的吸收、散射、衰减和漫射衰减系系数与原位设备测量值比较,误差均小于5%。
在该实施例中,基于船载多视场激光雷达,结合多次前向-单次后向雷达模型和Walker-Mclean模型探测的多个水体光学特征参数,与原位设备观测值比较的误差均小于为5%,从而印证了本发明针对船载激光雷达探测多个水体光学特征参数方法的有效性。该发明应用过程中,若激光雷达的接收视场或探测水体进行更改,均按照方法的实施流程进行分析,依次获得多次散射特征参数m值、水体前向散射系数、水体的吸收和散射系数、向下漫射衰减系数,即可获得探测水体更为全面的参数结果,提高了船载激光雷达探测水体光学特性的效率。
以上所述的实施例仅用以说明本发明技术方案,而非对其进行方法限制;尽管上述实施例对本发明进行了详述,对于本领域的技术人员来说,依然可以对该实施例所述的技术方案进行修改,或对其中的部分技术特征进行等同性地替换;而这些修改或替换,并不使相应技术方案的本质脱离本发明所要求保护的技术方案的精神和范围。
Claims (5)
1.船载多视场激光雷达探测多个水体光学特征参数的方法,其特征在于,包括以下步骤:
步骤1:利用船载激光雷达窄视场的探测回波反演激光雷达衰减系数,以估计水体的多次散射特征参数m的估计值;
步骤2:利用多次散射特征参数m的估计值,基于多次前向-单次后向散射雷达模型,建立不同视场接收时激光雷达回波比值的方程组,以计算水体前向散射系数bf;
步骤3:基于Walker-Mclean模型,获得不同视场的激光雷达回波比,并采用最小二乘法模型计算水体的吸收系数a和多次散射项g(z,τ);
步骤4:在Walker-Mclean模型中,进行极限深度假设,反演不同水层的向下漫射衰减系数Kd。
2.按照权利要求1所述的船载多视场激光雷达探测多个水体光学特征参数的方法,其特征在于,所述步骤1中的多次散射特征参数m的估计值取6~8之间的整数值。
3.按照权利要求2所述的船载多视场激光雷达探测多个水体光学特征参数的方法,其特征在于,所述步骤2包括:
步骤2.1:将多次散射特征参数m的估计值代入多次前向-单次后向散射雷达模型,得到多次散射项:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>m</mi>
<mi>&psi;</mi>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mi>f</mi>
</msub>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>&infin;</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>b</mi>
<mi>f</mi>
</msub>
<mi>z</mi>
</mrow>
<mi>x</mi>
</mfrac>
</msup>
<mo>&times;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<msup>
<mi>m</mi>
<mn>2</mn>
</msup>
</mrow>
<mn>4</mn>
</mfrac>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>r</mi>
<mrow>
<mi>l</mi>
<mi>a</mi>
<mi>s</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>R</mi>
<mrow>
<mi>r</mi>
<mi>c</mi>
<mi>v</mi>
<mi>r</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Theta;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
<mo>&rsqb;</mo>
<msub>
<mi>J</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>x</mi>
<mi>&psi;</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>x</mi>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,m为多次散射特征参数,2θrcvr为接收视场角,n水体折射率,H为激光雷达距离水面的高度,bf为水体前向散射系数,z为探测水深,2θlas为激光光束发散角,x为光子散射后位置向量的空间频率,J1为1阶贝塞尔函数,rlas和Rrcvr分别为探测激光束和望远镜接收通光面的半径;
步骤2.2:基于多次前向-单次后向散射雷达模型,计算船载激光雷达的回波强度:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>QA</mi>
<mrow>
<mi>r</mi>
<mi>c</mi>
<mi>v</mi>
<mi>r</mi>
</mrow>
</msub>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>v&tau;</mi>
<mi>d</mi>
</msub>
<mo>/</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mi>H</mi>
<mo>+</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mi>&beta;</mi>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>b</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>z</mi>
<mo>&rsqb;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,Q为激光能量,Arcvr为接收视场入瞳面积,T2为激光波长在大气中和气海界面的双程衰减,v为光速,τd为脉冲宽度,n为水体折射率,z为探测深度,H为激光雷达与水面间距离,β(π)为水体后向180°散射系数;a为水体吸收系数,bb为后向散射系数;F(z)为多次散射项;
步骤2.3:利用相同视场入瞳面积、不同的有效接收视场i,计算船载激光雷达的回波强度:
Pi(z)=Pwater(ψi,z)=M(z)F(ψi,z),
获得回波比值方程组:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>/</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>=</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
<mo>/</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>P</mi>
<mn>3</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>/</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>=</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mn>3</mn>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
<mo>/</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
其中, 对于同一水体,该项与视场无关;
求解回波比值方程组,计算积分项中的未知数:水体前向散射系数bf和多次散射特征参数m的准确值。
4.按照权利要求3所述的船载多视场激光雷达探测多个水体光学特征参数的方法,所述步骤3包括:
步骤3.1:先将水体前向散射系数bf近似为散射系数b,代入Walker-Mclean模型,则有:
其中,b为散射系数,τ为多次散射时间,g(z,τ)为对多次散射时间τ统计的Gamma分布概率密度函数,简称多次散射项;
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>,</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>&mu;</mi>
<mrow>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>&mu;</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&mu;</mi>
<mi>&tau;</mi>
</mrow>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<msup>
<mi>&mu;</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&mu;</mi>
<mi>&tau;</mi>
</mrow>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中,
<θ2>为散射角均方值,
θ=arccos[(0.142m2-1)/0.132m2],m为多次散射特征参数m的准确值;
步骤3.2:利用不同视场角θ1、θ2下的激光雷达回波比值,由于散射系数b取近似值,并去除水体后向180°散射系数β(π)的计算影响:
<mrow>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</msqrt>
<mo>=</mo>
<mfrac>
<mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>b</mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>b</mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>&times;</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>&infin;</mi>
</msubsup>
<mi>d</mi>
<mi>&tau;</mi>
<mfrac>
<mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>a</mi>
<mi>&tau;</mi>
<mi>v</mi>
<mo>/</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>,</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mfrac>
<mn>4</mn>
<mn>3</mn>
</mfrac>
<mi>z</mi>
<mi>&tau;</mi>
<mi>v</mi>
<mo>/</mo>
<mi>n</mi>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>tan&theta;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mi>H</mi>
<mo>+</mo>
<mi>z</mi>
<mo>/</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>b</mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>b</mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>&times;</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>&infin;</mi>
</msubsup>
<mi>d</mi>
<mi>&tau;</mi>
<mfrac>
<mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>a</mi>
<mi>&tau;</mi>
<mi>v</mi>
<mo>/</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>,</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mfrac>
<mn>4</mn>
<mn>3</mn>
</mfrac>
<mi>z</mi>
<mi>&tau;</mi>
<mi>v</mi>
<mo>/</mo>
<mi>n</mi>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>tan&theta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mi>H</mi>
<mo>+</mo>
<mi>z</mi>
<mo>/</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</mfrac>
</mrow>
</mfrac>
</mrow>
求解得到水体吸收系数a;
步骤3.3:采用最小二乘法模型,
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>F</mi>
<mi>u</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>g</mi>
<mo>(</mo>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>&tau;</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>&Integral;</mo>
<mo>|</mo>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>3</mn>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
</mrow>
<mi>exp</mi>
</msub>
<mo>-</mo>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>3</mn>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
</mrow>
<mi>mod</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>z</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
结合实验测量回波比值和Walker-Mclean模型所计算的回波比值多次迭代计算水体的多次散射项g(z,τ):
其中,θi为不同于θ3的视场角。
5.按照权利要求4所述的船载多视场激光雷达探测多个水体光学特征参数的方法,其特征在于,所述步骤4包括:
步骤4.1:首先将Walker-Mclean模型的深度z趋近于无穷大,得到基于极限深度假设的Walker-Mclean模型:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<munder>
<mi>lim</mi>
<mrow>
<mi>z</mi>
<mo>&RightArrow;</mo>
<mi>&infin;</mi>
</mrow>
</munder>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>M</mi>
<mn>0</mn>
</msub>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>K</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>d</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>~</mo>
<msub>
<mi>M</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>{</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mi>z</mi>
<mo>&rsqb;</mo>
<mo>+</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>a</mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>b</mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>&times;</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>&infin;</mi>
</msubsup>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>a</mi>
<mi>&tau;</mi>
<mi>v</mi>
<mo>/</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>,</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>&tau;</mi>
<mo>}</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mi>M</mi>
<mn>0</mn>
</msub>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,~为近似符号,
步骤4.2:将计算得到的水体吸收系数a和多次散射项g(z,τ)值代入基于极限深度假设的Walker-Mclean模型,利用limz→∞exp(-2KLidarz)=exp(-2Kdz)计算不同水层处向下漫射衰减系数Kd的趋近值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710977199.7A CN107831485A (zh) | 2017-10-19 | 2017-10-19 | 船载多视场激光雷达探测多个水体光学特征参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710977199.7A CN107831485A (zh) | 2017-10-19 | 2017-10-19 | 船载多视场激光雷达探测多个水体光学特征参数的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107831485A true CN107831485A (zh) | 2018-03-23 |
Family
ID=61648392
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710977199.7A Pending CN107831485A (zh) | 2017-10-19 | 2017-10-19 | 船载多视场激光雷达探测多个水体光学特征参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107831485A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110095784A (zh) * | 2019-05-09 | 2019-08-06 | 北京航空航天大学 | 一种复杂环境影响下的海洋-低层大气激光传输建模方法 |
CN110865389A (zh) * | 2019-10-29 | 2020-03-06 | 浙江大学 | 一种海洋激光雷达系统响应优化处理方法 |
CN111239713A (zh) * | 2020-03-18 | 2020-06-05 | 武汉大学 | 一种星载单光子激光雷达的最大测量深度评估方法 |
CN111965608A (zh) * | 2020-07-16 | 2020-11-20 | 自然资源部第二海洋研究所 | 一种基于水体叶绿素浓度的星载海洋激光雷达探测能力评估方法 |
CN113219496A (zh) * | 2021-03-18 | 2021-08-06 | 杭州电子科技大学 | 一种星载海洋大气参数激光雷达探测系统 |
CN113655495A (zh) * | 2021-08-31 | 2021-11-16 | 自然资源部第二海洋研究所 | 一种基于水体漫衰减系数光谱依赖性的星载海洋激光雷达最优波段评估方法 |
CN113740873A (zh) * | 2021-08-31 | 2021-12-03 | 自然资源部第二海洋研究所 | 一种基于高斯卷积的海洋激光雷达快速仿真方法 |
CN114167437A (zh) * | 2021-11-22 | 2022-03-11 | 桂林理工大学 | 测水激光雷达多通道设计方法 |
CN115615958A (zh) * | 2022-09-29 | 2023-01-17 | 哈尔滨工业大学(威海) | 闪光成像激光雷达测量水体后向衰减系数的方法及应用 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101915912A (zh) * | 2010-07-02 | 2010-12-15 | 武汉大学 | 一种全面的激光测高回波仿真方法 |
CN102012419A (zh) * | 2010-11-03 | 2011-04-13 | 浙江工业大学 | 基于视觉感知鱼类行为的生物水质监测系统 |
CN104614333A (zh) * | 2015-02-10 | 2015-05-13 | 中国人民解放军63655部队 | 一种大气水汽分子吸收系数廓线测量方法 |
US20150205895A1 (en) * | 2012-07-24 | 2015-07-23 | Navico Holding As | Reporting for an aquatic geographic information system |
CN105143914A (zh) * | 2013-04-17 | 2015-12-09 | 梅查理斯系统股份有限公司 | 用于在近程区域和远程区域中光学间距测量的方法 |
CN204855797U (zh) * | 2015-08-06 | 2015-12-09 | 中国科学院海洋研究所 | 一种用于水体光学特性测量的视场可调节光学系统 |
CN105372641A (zh) * | 2015-12-23 | 2016-03-02 | 国家海洋局第一海洋研究所 | 激光雷达水体回波信号仿真方法 |
CN105954732A (zh) * | 2016-06-14 | 2016-09-21 | 深圳大学 | 一种激光雷达测深数据的海底底质反射率提取方法及系统 |
CN106526611A (zh) * | 2016-10-31 | 2017-03-22 | 鲁东大学 | 一种利用激光雷达系统测量海洋固有光学参数的方法 |
CN106802289A (zh) * | 2017-01-20 | 2017-06-06 | 深圳大学 | 基于测深激光全波形数据的漫衰减系数提取方法及系统 |
-
2017
- 2017-10-19 CN CN201710977199.7A patent/CN107831485A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101915912A (zh) * | 2010-07-02 | 2010-12-15 | 武汉大学 | 一种全面的激光测高回波仿真方法 |
CN102012419A (zh) * | 2010-11-03 | 2011-04-13 | 浙江工业大学 | 基于视觉感知鱼类行为的生物水质监测系统 |
US20150205895A1 (en) * | 2012-07-24 | 2015-07-23 | Navico Holding As | Reporting for an aquatic geographic information system |
CN105143914A (zh) * | 2013-04-17 | 2015-12-09 | 梅查理斯系统股份有限公司 | 用于在近程区域和远程区域中光学间距测量的方法 |
CN104614333A (zh) * | 2015-02-10 | 2015-05-13 | 中国人民解放军63655部队 | 一种大气水汽分子吸收系数廓线测量方法 |
CN204855797U (zh) * | 2015-08-06 | 2015-12-09 | 中国科学院海洋研究所 | 一种用于水体光学特性测量的视场可调节光学系统 |
CN105372641A (zh) * | 2015-12-23 | 2016-03-02 | 国家海洋局第一海洋研究所 | 激光雷达水体回波信号仿真方法 |
CN105954732A (zh) * | 2016-06-14 | 2016-09-21 | 深圳大学 | 一种激光雷达测深数据的海底底质反射率提取方法及系统 |
CN106526611A (zh) * | 2016-10-31 | 2017-03-22 | 鲁东大学 | 一种利用激光雷达系统测量海洋固有光学参数的方法 |
CN106802289A (zh) * | 2017-01-20 | 2017-06-06 | 深圳大学 | 基于测深激光全波形数据的漫衰减系数提取方法及系统 |
Non-Patent Citations (5)
Title |
---|
ELEONORA ZEGE: "ELASTIC AND RAMAN LIDAR SOUNDING OF COASTAL WATERS.THEORY, COMPUTER SIMULATION, INVERSION POSSIBILITIES", 《EARSEL EPROCEEDINGS》 * |
FEYGELS: "Airborne lidar system with variable-field-of-view receiver for water optical properties measurement", 《PROCEEDINGS OF SPIE》 * |
李志刚: "一种激光雷达测量海洋光学参数的新方法", 《应用光学》 * |
李晓龙: "海洋激光雷达水体光学特性反演模型对比及分析", 《光学学报》 * |
李晓龙: "视场可调节海洋激光雷达实验系统与ICCD激光荧光实验研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110095784B (zh) * | 2019-05-09 | 2020-11-24 | 北京航空航天大学 | 一种复杂环境影响下的海洋-低层大气激光传输建模方法 |
CN110095784A (zh) * | 2019-05-09 | 2019-08-06 | 北京航空航天大学 | 一种复杂环境影响下的海洋-低层大气激光传输建模方法 |
CN110865389A (zh) * | 2019-10-29 | 2020-03-06 | 浙江大学 | 一种海洋激光雷达系统响应优化处理方法 |
CN111239713B (zh) * | 2020-03-18 | 2022-03-04 | 武汉大学 | 一种星载单光子激光雷达的最大测量深度评估方法 |
CN111239713A (zh) * | 2020-03-18 | 2020-06-05 | 武汉大学 | 一种星载单光子激光雷达的最大测量深度评估方法 |
CN111965608A (zh) * | 2020-07-16 | 2020-11-20 | 自然资源部第二海洋研究所 | 一种基于水体叶绿素浓度的星载海洋激光雷达探测能力评估方法 |
CN111965608B (zh) * | 2020-07-16 | 2024-01-12 | 自然资源部第二海洋研究所 | 一种基于水体叶绿素浓度的星载海洋激光雷达探测能力评估方法 |
CN113219496A (zh) * | 2021-03-18 | 2021-08-06 | 杭州电子科技大学 | 一种星载海洋大气参数激光雷达探测系统 |
CN113655495A (zh) * | 2021-08-31 | 2021-11-16 | 自然资源部第二海洋研究所 | 一种基于水体漫衰减系数光谱依赖性的星载海洋激光雷达最优波段评估方法 |
CN113740873A (zh) * | 2021-08-31 | 2021-12-03 | 自然资源部第二海洋研究所 | 一种基于高斯卷积的海洋激光雷达快速仿真方法 |
CN113740873B (zh) * | 2021-08-31 | 2024-05-28 | 自然资源部第二海洋研究所 | 一种基于高斯卷积的海洋激光雷达快速仿真方法 |
CN113655495B (zh) * | 2021-08-31 | 2024-06-11 | 自然资源部第二海洋研究所 | 一种基于水体漫衰减系数光谱依赖性的星载海洋激光雷达最优波段评估方法 |
CN114167437A (zh) * | 2021-11-22 | 2022-03-11 | 桂林理工大学 | 测水激光雷达多通道设计方法 |
CN114167437B (zh) * | 2021-11-22 | 2023-03-14 | 桂林理工大学 | 测水激光雷达多通道设计方法 |
CN115615958A (zh) * | 2022-09-29 | 2023-01-17 | 哈尔滨工业大学(威海) | 闪光成像激光雷达测量水体后向衰减系数的方法及应用 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107831485A (zh) | 船载多视场激光雷达探测多个水体光学特征参数的方法 | |
Arnaudon et al. | Riemannian medians and means with applications to radar signal processing | |
CN110095784B (zh) | 一种复杂环境影响下的海洋-低层大气激光传输建模方法 | |
CN105954732B (zh) | 一种激光雷达测深数据的海底底质反射率提取方法及系统 | |
CN107632308B (zh) | 一种基于递归叠加算法的车辆前方障碍物轮廓检测方法 | |
US8559753B2 (en) | Reflection removal system | |
CN113219496A (zh) | 一种星载海洋大气参数激光雷达探测系统 | |
CN115797760B (zh) | 主被动融合的水质立体遥感反演方法、系统及存储介质 | |
CN110673108B (zh) | 一种基于迭代Klett的机载海洋激光雷达信号处理方法 | |
CN110134976A (zh) | 一种机载激光测深信号提取方法及系统 | |
CN114295585B (zh) | 一种基于解析模型的多视场海洋激光雷达数据正则化反演方法 | |
Park et al. | New approach to sea surface wind retrieval from GNSS-R measurements | |
CN110376572B (zh) | 一种机载激光雷达水下探测深度模拟装置 | |
CN107219525B (zh) | Sar海面仿真方法及装置 | |
CN113310851B (zh) | 一种实时浑浊介质目标识别系统和方法 | |
CN110133680B (zh) | 一种机载激光测深接收波形有效信号初值确定方法及系统 | |
CN116609758B (zh) | 一种机载激光测深波形旅行时提取方法 | |
RU2586078C2 (ru) | Однопозиционный пассивный радиоэлектронный комплекс для определения горизонтальных координат, элементов движения цели и коэффициента километрического затухания электромагнитного излучения цели | |
CN111650128A (zh) | 一种基于地表反射率库的高分辨率大气气溶胶反演方法 | |
CN108919236B (zh) | 一种空间碎片激光测距效果仿真分析方法和装置 | |
Cheng et al. | Visibility improvement in turbid water by the fusion technology of Mueller matrix images | |
Guenther et al. | Prediction and correction of propagation-induced depth measurement biases plus signal attenuation and beam spreading for airborne laser hydrography | |
CN103593877A (zh) | 合成孔径声纳图像的仿真方法及系统 | |
CN106501790A (zh) | 一种浅海水下目标高光谱参数化特征提取方法 | |
CN116430353B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180323 |