CN110865389B - 一种海洋激光雷达***响应优化处理方法 - Google Patents
一种海洋激光雷达***响应优化处理方法 Download PDFInfo
- Publication number
- CN110865389B CN110865389B CN201911038431.6A CN201911038431A CN110865389B CN 110865389 B CN110865389 B CN 110865389B CN 201911038431 A CN201911038431 A CN 201911038431A CN 110865389 B CN110865389 B CN 110865389B
- Authority
- CN
- China
- Prior art keywords
- pulse width
- system response
- response
- laser radar
- laser
- 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
- 230000004044 response Effects 0.000 title claims abstract description 118
- 238000005457 optimization Methods 0.000 title claims abstract description 18
- 238000003672 processing method Methods 0.000 title claims abstract description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 43
- 238000000034 method Methods 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 10
- 238000012937 correction Methods 0.000 claims abstract description 9
- 230000003287 optical effect Effects 0.000 claims abstract description 8
- 230000008569 process Effects 0.000 claims abstract description 5
- 238000005070 sampling Methods 0.000 claims description 27
- 238000009792 diffusion process Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 5
- 238000010187 selection method Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000001514 detection method Methods 0.000 description 8
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 4
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical group O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 4
- 229910052799 carbon Inorganic materials 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000000630 rising effect Effects 0.000 description 4
- 239000002028 Biomass Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 231100000716 Acceptable daily intake Toxicity 0.000 description 2
- 229910002092 carbon dioxide Inorganic materials 0.000 description 2
- 239000001569 carbon dioxide Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005316 response function Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000012851 eutrophication Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000243 photosynthetic effect Effects 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Images
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
-
- 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
- G01S17/10—Systems determining position data of a target for measuring distance only using transmission of interrupted, pulse-modulated waves
-
- 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
-
- 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/483—Details of pulse systems
- G01S7/486—Receivers
- G01S7/4861—Circuits for detection, sampling, integration or read-out
-
- 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/483—Details of pulse systems
- G01S7/486—Receivers
- G01S7/487—Extracting wanted echo signals, e.g. pulse detection
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)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种海洋激光雷达***响应优化处理方法,首先,优化了海洋激光雷达***响应的器件选择和计算方法:提出了一种***响应近似模型,可快速获取激光雷达***响应脉宽;提出了激光器、探测器和采集卡三个器件之间的匹配选型方法,可发挥不同器件的最佳性能;其次,对因海洋激光雷达***响应导致的误差进行了抑制处理:通过选择适合不同水体的***响应脉宽,以减小因水面反射和***展宽导致的水面失真区;通过修正方程对***响应导致的浮游植物层的测量误差进行修正。本发明通过优化处理海洋激光雷达***响应,能够便利于研制海洋激光雷达时的器件选取,降低***响应对水体光学特性反演的影响。
Description
技术领域
本发明属于激光雷达技术领域,尤其是涉及一种海洋激光雷达***响应优化处理方法。
背景技术
海洋占据了地球表面71%的面积,而人类目前已探索的海洋面积仅占总面积的5%。海洋生态***是海洋中由许多生物群落及其环境构成的复杂集合,其中,浮游植物构成大多数海洋生态***的基础,他们的净光合作用固碳量约等于所有陆生植物固碳量的总和。食物链最基础的生产力驱动着大气-海洋的二氧化碳交换和深海的燃料碳固定,因此,浮游植物生产力在地球的海气耦合***中发挥了重要作用。而且,全球海洋浮游植物的存量每周都在消耗和再生长,这种快速的转换支撑了海洋浮游植物网以及渔业资源和全球食物供应。然而,随着科技的日益发展,经济的发展需求愈演愈烈,工业地区排放污染物、局部海域富营养化等都会对海洋环境造成恶劣的影响;二氧化碳排放、全球变暖等气候问题也会造成海平面上升、海洋生态***失衡等后果。因此,发展高效的、多样化的海洋探测方法迫在眉睫,用以帮助人类更好地认识海洋、开发海洋,同时更好地保护海洋环境。
海水光学参数的测量对研究全球气候变化和物质循环具有重要的意义,能有效用于遥感探测海洋内部特性的仪器具有巨大的优势。海洋激光雷达作为一种主动式光学遥感设备,可用于探测水体的光学特性,能够对上层海洋垂直分布结构进行遥感测量,具备高时空分辨率、昼夜连续观测、全球尺度测量等优势。目前,海洋激光雷达已能够成功应用于生物探测、浮游植物层分布研究、浅海地形测绘、海洋内波探测等相关研究。
公开号为CN105486664A的中国专利文献,公开了一种探测海洋浮游植物生物量和颗粒有机碳的装置及方法,基于高光谱分辨率激光雷达,结合海洋浮游植物生物量和POC反演算法,利用激光雷达方程中散射角为π时的粒子体积散射方程和有效衰减系数同步反演浮游植物生物量和POC。
公开号为CN107976686A的中国专利文献公开了一种多视场角海洋激光雷达及其视场角择优方法,其采用工作在蓝绿光波段的多视场角海洋激光雷达,建立参数灵敏度评价方法,择优选取多视场角激光雷达***的视场角,在水体等强散射体的多次散射效应下,仍能够高精度地遥测水体参数的垂直分布特性。
海洋激光雷达信号可以通过激光雷达方程来描述,因此,从信号中反演有效的信息通常也是基于激光雷达方程,如扰动法、Fernald法和斜率法。然而,实际接收到的激光雷达回波信号由包含目标特性的激光雷达方程以及激光雷达***时间***响应卷积组成,激光器、探测器和采集卡是影响***响应的主要器件。***响应会导致原始激光雷达方程发生变化,降低反演精度。
发明内容
本发明提供了一种海洋激光雷达***响应优化处理方法,可以提高海洋激光雷达反演精度,降低***响应的影响,精确得到近水面和浮游植物层中的信息。
本发明的技术方案如下:
一种海洋激光雷达***响应优化处理方法,其中,海洋激光雷达***包括激光发射***、接收***、数据采集及处理***,数据发射***包括激光器和扩束镜,接收***包括望远镜、光阑、窄带滤光片、会聚透镜和探测器,数据采集及处理***包括采集卡和计算机,具体的优化处理方法包括:
(1)根据激光器、探测器和采集卡对***响应的贡献,建立快速获取***响应脉宽的近似模型,所述的近似模型采用高斯模型,表达式为:
其中,T′T为近似的***响应脉宽,TL为激光器脉宽,TD为探测器脉宽;
(2)根据近似模型计算得到近似的***响应脉宽,并根据近似的***响应脉宽估算所需的采集卡采样率SR;
(3)采集卡根据得到的采样率SR进行采样,***通过计算得到实测浮游植物层厚度Thk;
(4)通过修正方程对***响应导致的实测浮游植物层厚度的误差进行修正,最后得到修正后的浮游植物层厚度。
本发明基于工作在蓝绿光波段的海洋激光雷达,通过采用海洋激光雷达***响应优化处理方法,可精确得到近水面和浮游植物层中的信息。
步骤(1)中,本发明提出了一种快速获取***响应脉宽的近似模型,通过分析三种不同激光波形对***响应的展宽影响,由于高斯模型符合大多数激光波形的情况,因此采用高斯模型。采用该模型,可快速计算近似的***响应脉宽T′T,大大减小了***响应脉宽TT的计算时间。
步骤(2)中,在根据近似模型计算得到近似的***响应脉宽之后,还包括:基于水体浑浊度判断近似的***响应脉宽是否符合测量条件,若满足,则继续进行采样率估算;若不满足,则重新设定激光器和探测器的脉宽。所述的测量条件为:
其中,T′T为近似的***响应脉宽,Kd为水体漫射衰减系数。
定义在水面以下的一段距离内,反演得到的漫射衰减系数相对误差大于10%的范围为信号失真区(SDZ),在信号失真区内,回波信号会出现一定程度的畸变。信号失真区与水体类型以及所使用的激光雷达***参数有关。本发明为使信号失真区小于一定值,通过分析***响应与水体漫射衰减系数的关系,得知不同水体类型所对应的激光雷达***响应脉宽具有一定的范围。随着水体的漫射衰减系数不断变大,***响应脉宽范围逐渐变小。当***响应脉宽与水体漫射衰减系数满足上述公式时,可使信号失真区范围与最大可探测的深度的比值小于20%,能够获得较准确的激光雷达信号,以反演较准确的水体光学参数。
步骤(2)中,估算所需的采集卡采样率SR的过程如下:
通过对数据进行线性拟合得到***响应脉宽TT与采样率SR之间的反比例线性关系,表达式为:
SR=0.8234/TT+3.3×106
再将***响应脉宽TT用近似的***响应脉宽T′T代替,利用上述表达式,得到估算的采集卡采样率。得到的采样率SR需大于或等于两倍的***带宽B。
根据奈奎斯特采样定理,以奈奎斯特速率或高于奈奎斯特速率采样的信号,能保证真实的水体信号可以被准确无误地反演出来。在激光雷达***中,如果采集卡的采样率过低,导致采样时间大于总***脉宽,则激光雷达***对回波信号进行欠采样,这将导致激光雷达***无法进行探测。由所述的公式所知,通过***响应脉宽能够大致估算所需要的采集卡采样率。对于估算***采样率,真实值TT和近似模型T′T所对应的结果近乎一致。
步骤(4)中,浮游植物层厚度校正时,当***响应脉宽较大时,浮游植物层的厚度误差较大,步骤(4)中,所述的修正方程为:
其中,Thk′为修正后的浮游植物层厚度;Thk为浮游植物层的实测厚度;LT为近似的***响应脉冲长度,其由近似的***响应脉宽T′T通过光速匹配得到。
激光雷达***及其参数会对浮游植物层的探测有一定的展宽作用。展宽效应主要表现为回波信号中浮游植物层的厚度变宽,探测得到的浮游植物层的厚度会随着激光脉宽和探测器响应的增加而增加。***响应脉宽越小,浮游植物层的探测精度则越大。可以根据总时间响应来估算浮游植物层达到的精度,或者通过需要的浮游植物层精度来匹配所需要的。
本发明中,所述激光器的时域波形为高斯型,所述探测器的时域波形为重尾型。
与现有技术相比,本发明具有以下有益效果;
1、本发明优化了海洋激光雷达***响应的器件选择和计算方法:提出了一种***响应近似模型,可快速获取激光雷达***响应脉宽;提出了激光器、探测器和采集卡三个器件之间的匹配选型方法,可发挥不同器件的最佳性能;
2、本发明对因海洋激光雷达***响应导致的误差进行了抑制处理:通过选择适合不同水体的***响应脉宽,以减小因水面反射和***展宽导致的水面失真区;通过修正方程对***响应导致的浮游植物层的测量误差进行修正。
3、本发明通过优化处理海洋激光雷达***响应,能够便利于研制海洋激光雷达时的器件选取,降低***响应对水体光学特性反演的影响。
附图说明
图1为本发明一种海洋激光雷达***响应优化处理方法的流程示意图;
图2为本发明实施例中近似模型相对误差与激光器、探测器的关系;
图3为本发明实施例中激光雷达采样率与***时间响应的关系,其中,(a)提供了采集卡与激光器以及探测器的匹配情况,(b)表示TT与SR之间的关系;
图4为本发明实施例中不同水体类型所对应的激光雷达***响应脉宽;
图5为本发明实施例中浮游植物层厚度的相对误差与真实厚度以及***响应脉宽的关系,其中,(a)为实测浮游植物层厚度的相对误差与真实厚度以及LT的关系,(b)为修正后的浮游植物层厚度与真实值的相对误差分布情况。
具体实施方式
下面结合附图和实施例对本发明做进一步详细描述,需要指出的是,以下所述实施例旨在便于对本发明的理解,而对其不起任何限定作用。
如图1所示,一种海洋激光雷达***响应优化处理方法,基于工作在蓝绿光波段的海洋激光雷达,通过采用海洋激光雷达***响应优化处理方法,精确得到近水面和浮游植物层中的信息。
首先,优化了海洋激光雷达***响应的器件选择和计算方法:根据激光器、探测器和采集卡等三个关键器件对***响应的贡献,提出了一种***响应近似模型,快速获取激光雷达***响应脉宽;提出了激光器、探测器和采集卡三个器件之间的匹配选型方法,发挥不同器件的最佳性能;其次,对因海洋激光雷达***响应导致的误差进行了抑制处理:根据准单次散射近似的激光雷达方程,通过选择适合不同水体的***响应脉宽,减小因水面反射和***展宽导致的水面失真区;通过修正方程对***响应导致的浮游植物层的测量误差进行修正。
本发明涉及到的出射激光脉冲形状模型通常为高斯模型、重尾模型或截段抛物线模型。其中本实施例使用的高斯模型可以表示为式(1):
光电倍增管的下降沿时间是上升沿时间Tr的2~3倍,脉宽TD是上升沿的2.5倍左右。PMT的时间响应函数可以用重尾模型来建模,假设TD=3Tr,PMT响应函数可表示为式(2):
激光脉宽TL选取范围为2~18ns,PMT上升沿时间设置为0.6~2.2ns内,对应的探测器脉宽TD为1.8~6.6ns,采用控制变量法进行分析,并对激光雷达信号进行了归一化处理,将激光雷达信号的动态范围设置为常见的3个数量级。在仿真模型中,设置表面反射率为0.02,漫射衰减系数Kd大于等于0.1m-1,大气损耗为1,水面反射的太阳辐射为0.14mw/m2srnm。具体仿真中用到的激光雷达***参数如表1所示。
表1
本发明使用的激光雷达***响应近似模型,通过分析三种不同激光波形对***响应的展宽影响,综合考虑,由于高斯模型符合大多数激光波形的情况,采用高斯模型。激光脉宽用TL表示,探测器脉宽用TD表示,根据高斯函数卷积的特性,提出***响应脉宽的近似模型可以表示为式(3):
其中TD为探测器脉宽,TL为激光脉宽。***响应脉宽TT是激光器和探测器共同作用的结果。近似模型T′T与真实值TT之间的误差可写作式(4):
图2展示了不同激光器、探测器脉宽情况下计算得到的***响应近似模型T′T与真值TT之间的相对误差在3%以内,近似模型的相对误差分布取决于激光器脉宽与探测器脉宽的绝对差值,当绝对差值比较大时,误差较小。
本发明所述的***器件带宽匹配方法,其特征在于采集卡的采样率SR必须大于或等于两倍***带宽B(频域上的半高全宽值)才能正常探测信号。对数据进行线性拟合得到TT与SR之间的关系可表示为式(5):
SR=0.8234/TT+3.3×106 (5)
图3示出了激光雷达***采样率和***响应各部分的关系。
图3中,(a)为***响应与探测器脉宽以及激光器脉宽之间的关系,采集卡的采样率SR必须大于或等于两倍***带宽B(频域上的半高全宽值)才能正常探测信号。图中提供了采集卡与激光器以及探测器的匹配情况,虚线表示***采样率为120MHz时所对应的激光器脉宽和探测器脉宽,可以看出120MHz的采集卡满足绝大多数情况下的***配置。
图3中,(b)示出了TT与SR之间的关系,可以看出TT与SR具有极好的反比例线性关系。其中,实线对应的是近似模型计算的结果,虚线对应真实值计算的结果,可以看出对于估算***采样率,TT、T′T所对应的结果近乎一致。
本发明基于水体浑浊度的***响应脉宽选择,不同水体类型需选择合适的激光雷达***响应脉宽。水体漫射衰减系数Kd与允许的最大***响应脉宽TT满足关系式(6):
如图4所示,进一步讨论分析***响应与水体探测能力的关系,假设信号失真区与最大可探测的深度的比值为20%,探究不同水体类型所对应的激光雷达***响应脉宽TT。图中圈线表示选用任意的激光脉宽与探测器脉宽的情况,图中实线表示激光器固定为2ns的情况,图中的虚线表示探测器固定为1.8ns的情况。
分析图4可知,不同水体类型所对应的激光雷达***响应脉宽具有一定的波动范围,且随着水体的漫射衰减系数不断变大,***响应脉宽范围越来越小。波动范围的下边界为激光器脉宽固定为最小值的情况,上边界为探测器固定为最小值的情况。由于激光器和探测器的时域响应波形不同(激光器的时域波形为高斯型,探测器的时域波形为重尾型),以探测器响应为主时,探测不同水体所允许的***响应脉宽小于以激光器为主时的情况。
本发明浮游植物层厚度校正,在***响应脉宽较大时,浮游植物层的厚度误差较大,厚度可以通过实测的厚度Thk和***响应脉冲长度LT进行修正,可表示为式(7):
***响应对浮游植物层探测的影响,主要指浮游植物层的厚度,***响应脉宽TT和探测到的浮游植物层厚度相关,TT越小,浮游植物层的探测精度则越大。根据总时间响应估算浮游植物层达到的精度,或者通过需要的浮游植物层精度来匹配所需要的TT。本发明所述的***响应脉宽的近似模型T′T同样适用于评估浮游植物层的探测。
通过图5,进一步展示并分析了浮游植物层厚度误差与***响应脉宽之间的影响,其中,***响应脉宽TT设置范围为2.7~18.6ns,对应于水体中的脉冲长度LT范围为0.3~2.1m,浮游植物层的厚度范围为0.2~4.2m。
图5中,(a)展示了实测浮游植物层厚度的相对误差与真实厚度以及LT的关系,从图中可以看出,当真实厚度和LT之间相差较大时,实测值的相对误差也会越来越小。图中线性拟合出了相对误差为10%、20%、40%的情况,当相对误差为定值时,浮游植物层厚度与总***脉冲宽度TT或LT大致满足线性关系。在***响应脉宽较大时,浮游植物层的厚度误差较大。
图5中,(b)为修正后的浮游植物层厚度与真实值的相对误差分布情况,修正后的相对误差平均值由原来的55%提升至6%。采用线性回归的方法拟合出了相对误差为5%、10%、20%的情况,当Thk>0.5LT+0.5时,修正误差不超过5%,当Thk=0.4LT+0.3附近时,相对误差稍微偏大,但其分布仍在10%左右。当Thk<0.3LT+0.1时,此时修正后的厚度误差仍然大于20%。综上,本发明所述的浮游植物层修正公式在绝大多数情况下能够校准实测浮游植物层的厚度。
本发明通过优化处理海洋激光雷达***响应,能够便利于研制海洋激光雷达时的器件选取,降低***响应对水体光学特性反演的影响。
以上所述的实施例对本发明的技术方案和有益效果进行了详细说明,应理解的是以上所述仅为本发明的具体实施例,并不用于限制本发明,凡在本发明的原则范围内所做的任何修改、补充和等同替换,均应包含在本发明的保护范围之内。
Claims (6)
1.一种海洋激光雷达***响应优化处理方法,其中,海洋激光雷达***包括激光发射***、接收***、数据采集及处理***,激光发射***包括激光器和扩束镜,接收***包括望远镜、光阑、窄带滤光片、会聚透镜和探测器,数据采集及处理***包括采集卡和计算机,其特征在于,具体的优化处理方法包括:
(1)根据激光器、探测器和采集卡对***响应的贡献,建立快速获取***响应脉宽的近似模型,所述的近似模型采用高斯模型,表达式为:
其中,T′T为近似的***响应脉宽,TL为激光器脉宽,TD为探测器脉宽;
(2)根据近似模型计算得到近似的***响应脉宽,并根据近似的***响应脉宽估算所需的采集卡采样率SR;估算所需的采集卡采样率SR的过程如下:
通过对数据进行线性拟合得到***响应脉宽TT与采样率SR之间的反比例线性关系,表达式为:
SR=0.8234/TT+3.3×106
再将***响应脉宽TT用近似的***响应脉宽TT′代替,利用上述表达式,得到估算的采集卡采样率;
(3)采集卡根据得到的采样率SR进行采样,***通过计算得到实测浮游植物层厚度Thk;
(4)通过修正方程对***响应导致的实测浮游植物层厚度的误差进行修正,最后得到修正后的浮游植物层厚度;所述的修正方程为:
其中,Thk′为修正后的浮游植物层厚度;Thk为浮游植物层的实测厚度;LT为近似的***响应脉冲长度,其由近似的***响应脉宽TT′通过光速匹配得到。
2.根据权利要求1所述的海洋激光雷达***响应优化处理方法,其特征在于,步骤(2)中,在根据近似模型计算得到近似的***响应脉宽之后,还包括:基于水体浑浊度判断近似的***响应脉宽是否符合测量条件,若满足,则继续进行采样率估算;若不满足,则重新设定激光器和探测器的脉宽。
4.根据权利要求2所述的海洋激光雷达***响应优化处理方法,其特征在于,不同水体浑浊度所对应的***响应脉宽具有一定的波动范围,且随着水体漫射衰减系数不断变大,***响应脉宽的波动范围逐渐变小;所述波动范围的下边界为激光器脉宽固定为最小值的情况,上边界为探测器脉宽固定为最小值的情况。
5.根据权利要求1或4所述的海洋激光雷达***响应优化处理方法,其特征在于,步骤(2)中,采样率SR需大于或等于两倍的***带宽B。
6.根据权利要求1所述的海洋激光雷达***响应优化处理方法,其特征在于,所述激光器的时域波形为高斯型,所述探测器的时域波形为重尾型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911038431.6A CN110865389B (zh) | 2019-10-29 | 2019-10-29 | 一种海洋激光雷达***响应优化处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911038431.6A CN110865389B (zh) | 2019-10-29 | 2019-10-29 | 一种海洋激光雷达***响应优化处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110865389A CN110865389A (zh) | 2020-03-06 |
CN110865389B true CN110865389B (zh) | 2021-08-10 |
Family
ID=69654168
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911038431.6A Active CN110865389B (zh) | 2019-10-29 | 2019-10-29 | 一种海洋激光雷达***响应优化处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110865389B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113655495B (zh) * | 2021-08-31 | 2024-06-11 | 自然资源部第二海洋研究所 | 一种基于水体漫衰减系数光谱依赖性的星载海洋激光雷达最优波段评估方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5091778A (en) * | 1990-12-21 | 1992-02-25 | Kaman Aerospace Corporation | Imaging lidar systems and K-meters employing tunable and fixed frequency laser transmitters |
WO2006076474A1 (en) * | 2005-01-13 | 2006-07-20 | Arete Associates | Optical system with wavefront sensor |
CN106526611A (zh) * | 2016-10-31 | 2017-03-22 | 鲁东大学 | 一种利用激光雷达***测量海洋固有光学参数的方法 |
CN106970392A (zh) * | 2017-05-31 | 2017-07-21 | 南京先进激光技术研究院 | 高灵敏度气体污染探测激光雷达*** |
CN110268283A (zh) * | 2016-11-16 | 2019-09-20 | 创新科技有限公司 | 激光雷达***和方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110210237A1 (en) * | 2010-01-21 | 2011-09-01 | Shmuel Sternklar | Methods and systems for measuring the frequency response and impulse response of objects and media |
WO2014066645A1 (en) * | 2012-10-24 | 2014-05-01 | Naiad Maritime Group, Inc. | Predictive sea state mapping for ship motion control |
US9274221B2 (en) * | 2013-02-01 | 2016-03-01 | GM Global Technology Operations LLC | Method and apparatus for remote object sensing employing compressive sensing |
CN103308926B (zh) * | 2013-06-18 | 2015-07-22 | 浙江大学 | 一种高光谱分辨率激光雷达装置 |
US9720415B2 (en) * | 2015-11-04 | 2017-08-01 | Zoox, Inc. | Sensor-based object-detection optimization for autonomous vehicles |
CN105486664B (zh) * | 2015-12-31 | 2018-04-20 | 浙江大学 | 探测海洋浮游植物生物量和poc的激光雷达装置及方法 |
US10837773B2 (en) * | 2016-12-30 | 2020-11-17 | DeepMap Inc. | Detection of vertical structures based on LiDAR scanner data for high-definition maps for autonomous vehicles |
CN107831485A (zh) * | 2017-10-19 | 2018-03-23 | 中国科学院海洋研究所 | 船载多视场激光雷达探测多个水体光学特征参数的方法 |
CN107976686B (zh) * | 2017-10-26 | 2020-12-08 | 浙江大学 | 一种多视场角海洋激光雷达及其视场角择优方法 |
CN108459313B (zh) * | 2018-02-07 | 2023-10-20 | 苏州镭图光电科技有限公司 | 基于光学微机电***的激光雷达回波处理方法 |
CN208109385U (zh) * | 2018-05-16 | 2018-11-16 | 北京华宇德信光电技术有限公司 | 一种激光响应时间测试*** |
-
2019
- 2019-10-29 CN CN201911038431.6A patent/CN110865389B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5091778A (en) * | 1990-12-21 | 1992-02-25 | Kaman Aerospace Corporation | Imaging lidar systems and K-meters employing tunable and fixed frequency laser transmitters |
WO2006076474A1 (en) * | 2005-01-13 | 2006-07-20 | Arete Associates | Optical system with wavefront sensor |
CN106526611A (zh) * | 2016-10-31 | 2017-03-22 | 鲁东大学 | 一种利用激光雷达***测量海洋固有光学参数的方法 |
CN110268283A (zh) * | 2016-11-16 | 2019-09-20 | 创新科技有限公司 | 激光雷达***和方法 |
CN106970392A (zh) * | 2017-05-31 | 2017-07-21 | 南京先进激光技术研究院 | 高灵敏度气体污染探测激光雷达*** |
Also Published As
Publication number | Publication date |
---|---|
CN110865389A (zh) | 2020-03-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gommenginger et al. | Retracking altimeter waveforms near the coasts: A review of retracking methods and some applications to coastal waveforms | |
Abdallah et al. | Wa-LiD: A new LiDAR simulator for waters | |
Saylam et al. | Airborne lidar bathymetry: Assessing quality assurance and quality control methods with Leica Chiroptera examples | |
CN111965608B (zh) | 一种基于水体叶绿素浓度的星载海洋激光雷达探测能力评估方法 | |
CN102759731B (zh) | 基于星载激光测高仪回波的海洋表面风、浪特性反演方法 | |
Fore et al. | SMAP salinity and wind speed data user’s guide | |
Schlembach et al. | Round robin assessment of radar altimeter Low Resolution Mode and delay-Doppler retracking algorithms for significant wave height | |
Li et al. | Dual antenna space-based GNSS-R ocean surface mapping: Oil slick and tropical cyclone sensing | |
CN110673108B (zh) | 一种基于迭代Klett的机载海洋激光雷达信号处理方法 | |
CN105445751A (zh) | 一种浅水区域水深比值遥感反演方法 | |
Astudillo et al. | Surface winds off Peru-Chile: observing closer to the coast from radar altimetry | |
Guo et al. | Development of a single-wavelength airborne bathymetric LiDAR: System design and data processing | |
CN111460681B (zh) | 监测近海水体真光层深度的卫星遥感方法 | |
CN115293198A (zh) | 一种基于多隐层神经网络提高gnss-r测高反演精度的方法 | |
CN110865389B (zh) | 一种海洋激光雷达***响应优化处理方法 | |
CN111561916B (zh) | 一种基于四波段多光谱遥感图像的浅海水深无控提取方法 | |
CN114295585B (zh) | 一种基于解析模型的多视场海洋激光雷达数据正则化反演方法 | |
CN110376572B (zh) | 一种机载激光雷达水下探测深度模拟装置 | |
Durand et al. | How will radar layover impact SWOT measurements of water surface elevation and slope, and estimates of river discharge? | |
CN111426643A (zh) | 一种大气二氧化氮柱浓度反演条带噪声自适应窗口校正方法 | |
Zhang et al. | A multiband model with successive projections algorithm for bathymetry estimation based on remotely sensed hyperspectral data in Qinghai Lake | |
Desmars | Real-time reconstruction and prediction of ocean wave fields from remote optical measurements | |
Rikka et al. | Sea state in the Baltic Sea from space-borne high-resolution synthetic aperture radar imagery | |
Hiroji | Extracting sonar relative beam patterns for multi-sector multibeam sonar | |
CN117075149A (zh) | 基于ddm的星载gnss-r台风位置估计方法及*** |
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 |