CN114674292A - 一种基于机载多波长激光雷达探测海洋光学剖面的系统及反演方法 - Google Patents

一种基于机载多波长激光雷达探测海洋光学剖面的系统及反演方法 Download PDF

Info

Publication number
CN114674292A
CN114674292A CN202111591030.0A CN202111591030A CN114674292A CN 114674292 A CN114674292 A CN 114674292A CN 202111591030 A CN202111591030 A CN 202111591030A CN 114674292 A CN114674292 A CN 114674292A
Authority
CN
China
Prior art keywords
light beam
wavelength
laser radar
band
light
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.)
Granted
Application number
CN202111591030.0A
Other languages
English (en)
Other versions
CN114674292B (zh
Inventor
毛志华
袁大鹏
陈鹏
贺岩
张镇华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiaotong University
Second Institute of Oceanography MNR
Original Assignee
Shanghai Jiaotong University
Second Institute of Oceanography MNR
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Jiaotong University, Second Institute of Oceanography MNR filed Critical Shanghai Jiaotong University
Priority to CN202111591030.0A priority Critical patent/CN114674292B/zh
Publication of CN114674292A publication Critical patent/CN114674292A/zh
Application granted granted Critical
Publication of CN114674292B publication Critical patent/CN114674292B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection
    • G01N21/49—Scattering, i.e. diffuse reflection within a body or fluid
    • G01N21/53—Scattering, i.e. diffuse reflection within a body or fluid within a flowing fluid, e.g. smoke
    • 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
    • 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/483—Details of pulse systems
    • G01S7/484—Transmitters
    • 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
    • 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/4868—Controlling received signal intensity or exposure of sensor
    • 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
    • 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/497—Means for monitoring or calibrating
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Health & Medical Sciences (AREA)
  • Hydrology & Water Resources (AREA)
  • Optical Radar Systems And Details Thereof (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开一种基于机载多波长激光雷达探测海洋光学剖面的系统，该系统包括多波长激光发射系统、望远镜接收系统、分视场镜系统、四个分光系统、七个光电探测系统、偏振分光系统、数据采集与处理系统。反演方法采用高重复频率的多波长激光器作为发射光源，基于斜率‑扰动法反演海洋水体光学剖面参数。针对深海和浅海区域水体的光学特性分别采用大视场和小视场，重点探测海洋水体光学剖面。本发明的系统及方法具有探测深海和浅海区域的可靠性能，适应一类水体和二类水体探测的最优波长，精确反演海洋光学剖面以及浮游植物垂直分层结构。

Description

一种基于机载多波长激光雷达探测海洋光学剖面的系统及反 演方法
技术领域
本发明涉及海洋光学探测领域，具体涉及一种基于机载多波长激光雷达探测海洋光学剖面的系统及反演方法。
背景技术
随着激光技术的发展，主动式激光雷达具有高时空分辨率、可探测水体光学剖面和昼夜持续观测等优点。不仅在大气云层和吸收性气溶胶方面得到广泛的应用，而且在海洋浮游植物层探测、鱼群探测、水深反演等方面有较大的应用前景。海洋激光雷达探测原理是利用激光在海水中产生的后向散射信号反演海洋环境参数，在激光雷达系统硬件方面，目前仍存在着许多未解决的技术问题：如激光雷达同时探测浅海和深海区域时，浅海区域回波信号过饱和；针对不同类型水体，输出单一波长的激光雷达无法达到理想的探测性能。在反演算法方面，存在基于激光雷达回波信号反演海水光学参数精度较低的问题。因此在海洋激光雷达探测领域，亟需发展新型探测技术和较高精度的海洋光学参数反演方法。
发明内容
针对现有技术的不足，本发明公开一种基于机载多波长激光雷达探测海洋光学剖面的系统及反演方法，该系统采用高重复频率的多波长激光器作为发射光源，通过分视场镜把接收光路分为大视场和小视场。反演方法基于斜率-扰动法反演海洋水体光学剖面参数，针对深海和浅海区域水体的光学特性分别采用大视场和小视场，重点探测海洋水体光学剖面。
本发明的目的通过如下的技术方案来实现：
一种基于机载多波长激光雷达探测海洋光学剖面的系统，该系统包括：
多波长激光发射系统、望远镜接收系统、分视场镜系统、分光系统一、分光系统二、分光系统三、分光系统四、光电探测系统一、光电探测系统二、光电探测系统三、光电探测系统四、光电探测系统五、光电探测系统六、光电探测系统七、偏振分光系统、数据采集与处理系统；
所述多波长激光发射系统包括多波长激光器和激光扩束器，多波长激光器用于发射四种不同波长的激光，激光扩束器用于出射激光扩束并压缩其发散角；
所述望远镜接收系统用于接收激光在海水中产生的后向散射信号；
所述分视场镜系统包括聚焦透镜一和分视场镜，所述聚焦透镜一将所述后向散射信号聚焦到分视场镜，分视场镜将光束分成两束，一束进入分光系统一，另一束进入分光系统二；
分光系统一包括准直透镜一、分光镜一和反射镜，准直透镜一将光束准直后，分光镜一将准直后的光束分为两束，一束进入光电探测系统一，另一束经反射镜入射到光电探测系统二；
光电探测系统一包括带通滤光片一、聚焦透镜二和光电探测器一，所述带通滤光片一将光束滤光后，聚焦透镜二将滤光后的光束聚焦，光电探测器一接收聚焦后的光束信号；
光电探测系统二包括带通滤光片二、聚焦透镜三和光电探测器二，所述带通滤光片二将光束滤光后，聚焦透镜三将滤光后的光束聚焦，光电探测器二接收聚焦后的光束信号；
分光系统二包括准直透镜二和分光镜二，准直透镜二将光束准直后，分光镜二将光束分为两束，一束进入光电探测系统三，另一束进入分光系统三；
光电探测系统三包括带通滤光片三、聚焦透镜四和光电探测器三，带通滤光片三将光束进行滤光后，聚焦透镜四将滤光后的光束聚焦，光电探测器三接收聚焦后的光束信号；
分光系统三为分光镜三，其将光束分为两束，一束进入光电探测系统四，另一束进入分光系统四；
光电探测系统四包括带通滤光片四、聚焦透镜五和光电探测器四，带通滤光片四将光束滤光后，聚焦透镜五对滤光后的光束聚焦，光电探测器四接收聚焦后的光束信号；
分光系统四包括分光镜四和带通滤光片五，分光镜四将光束分为两束，一束进入光电探测系统五，另一束经过带通滤光片五滤光后进入偏振分光系统；
光电探测系统五包括带通滤光片六、聚焦透镜六和光电探测器五，带通滤光片六将光束滤光后，聚焦透镜六对滤光后的光束聚焦，光电探测器五接收聚焦后的光束信号；
偏振分光系统为偏振分光镜，其将光束分成两束，一束进入光电探测系统六，另一束进入光电探测系统七；
光电探测系统六包括聚焦透镜七和光电探测器六，聚焦透镜七将光束聚焦后，光电探测器六接收聚焦后的光束信号；
光电探测系统七包括聚焦透镜八和光电探测器七，聚焦透镜八将光束聚焦后，光电探测器七接收聚焦后的光束信号；
数据采集与处理系统包括多通道数据采集卡和计算机，两者均与所有的光电探测器电连接，多通道数据采集卡用于机载多波长海洋激光雷达数据的采集，计算机用于激光雷达数据的处理分析。
进一步地，所述望远镜接收系统采用200mm直径的卡式望远镜，焦距为700mm。
进一步地，所述带通滤光片一的中心波长为500nm,透过率大于90％，短波介质范围： 200-472nm，长波截止范围：527-1200nm；所述带通滤光片二的中心波长为532nm，透过率大于90％，短波截止范围：200-512nm，长波截止范围：552-1200nm；所述带通滤光片三的中心波长为1064nm，透过率大于90％，短波截止范围：200-1039nm，长波截止范围： 1089-1200nm；所述带通滤光片四的中心波长为355nm，透过率大于90％，短波截止范围： 200-345nm，长波截止范围：365-1200nm；所述带通滤光片五的中心波长为500nm,透过率大于90％，短波介质范围：200-472nm，长波截止范围：527-1200nm；所述带通滤光片六的中心波长为532nm，透过率大于90％，短波截止范围：200-512nm，长波截止范围：552-1200nm。
一种基于上述的机载多波长激光雷达系统探测海洋光学剖面的反演方法，该方法包括如下步骤：
S1：所述多波长激光器发射1064nm、532nm、486nm和355nm的激光，分视场镜将光束分为大视场角的光束和小视场角的光束，大视场角的光束进入分光系统一，小视场角的光束进入分光系统二；
S2：计算机对七个光电探测器采集到的机载多波长激光雷达数据进行预处理；
(1)基于飞机POS数据中的飞行时间、经纬度和激光雷达数据采集时间对七个光电探测器采集到的机载多波长激光雷达数据进行分离，保留多波长海洋激光雷达数据Sraw；
(2)将多波长海洋激光雷达数据Sraw中的脉冲数据进行平均处理，得到Smean；
(3)对Smean中脉冲数据的背景噪声点进行平均处理，得到Snoise，然后基于公式Sdenoise＝Smean-Snoise得到去噪后的脉冲数据Sdenoise；
(4)利用1064nm通道采集的激光雷达数据计算机载激光雷达系统与海表之间的高度，并对所有通道的激光雷达数据进行距离校正，得到S(λ，z)；其中，λ为激光波长，z为海水深度；
S3:计算海洋水体激光雷达衰减系数Klidar(λ，z)；
假设海洋水体后向散射系数与激光雷达衰减系数满足β(λ，z)＝AKlidar(λ，z)R，取激光雷达探测的某边界值深度为zB，则激光雷达衰减系数表示为：
Figure BDA0003429913730000031
其中，A为相关参数，R为激光雷达后向散射消光对数比；
S4：计算海洋水体激光雷达后向散射系数β(λ，z)；
(1)假设激光雷达后向散射系数是不随深度变化的均匀水体光学参数和随深度变化的非均匀水体光学参数之和，则激光雷达后向散射系数β(λ，z)表述为
Figure BDA0003429913730000041
其中，
S(λ，z)＝Sh(λ，z)+Sinh(λ，z)
Sh(λ，z)＝k·βh(λ，z)·exp[-2Klidar(λ，z)]
Sinh(λ，z)＝k·βinh(λ，z)·exp[-2Klidar(λ，z)]
lnSh(λ，z)＝lnk·βh(λ，0)-2Klidar(λ，0)
其中，Sh(λ，z)为均匀水体部分产生的激光雷达回波信号，Sinh(λ，z)为非均匀水体产生的激光雷达回波信号，k为激光雷达系统常数，βh(λ，z)为均匀水体产生的后向散射系数，βinh(λ，z)为非均匀水体产生的后向散射系数；
S5：分别计算大视场和小视场各个波段的海洋水体光学剖面参数Klidar(532nm，z)，Klidar(486nm，z)，Klidar(355nm，z)，β(532nm，z)，β(486nm，z)和β(355nm，z)。
进一步地，zB选取激光雷达的最大有效探测深度。
进一步地，所述大视场角的范围为6-25mrad，小视场角的范围为0-6mrad。
本发明的有益效果如下：
(1)本发明的基于机载多波长激光雷达探测海洋光学剖面的系统，具有探测海洋中一类水体和二类水体的能力，通过接收到的多波长激光雷达信号，可用于探测海洋水体光学剖面、次表层信息、浮游植物垂直分层结构、以及海洋混合层深度的时空变化。
(2)通过分视场镜系统把光路分为小视场和大视场，能更有效的探测浅海区域和深海区域。
(3)本发明的反演方法首先通过1064nm通道对所有通道的激光雷达数据进行距离校正，然后选取激光雷达最大有效探测深度作为光学参数反演的边界条件，最后利用扰动法假设激光雷达后向散射系数是不随深度变化的均匀水体光学参数和随深度变化的非均匀水体光学参数之和。在距离校正方面，提升了不同深度处的光学参数的反演精度；在最大有效探测深度范围内，降低了因为信号突变而造成的反演误差；扰动法提升了激光雷达在水体光学特性变化不剧烈时光学参数的反演精度。通过该反演方法获得的海洋光学剖面参数，能够用于评估浮游植物粒径分布，并进一步评估海洋浮游植物群落分布。
附图说明
图1为本发明的基于机载多波长激光雷达探测海洋光学剖面的系统的示意图；
图2为本发明的反演方法的流程图；
图3为反演的海洋水体激光雷达衰减系数结果图；
图4为反演的海洋水体激光雷达后向散射系数结果图。
图中，多波长激光器1、激光扩束器2、望远镜接收系统3、聚焦透镜一4、分视场镜5、准直透镜一6、分光镜一7、带通滤光片一8、聚焦透镜二9、光电探测器一10、反射镜11、带通滤光片二12、聚焦透镜三13、光电探测器二14、准直透镜二15、分光镜二16、带通滤光片三17、聚焦透镜四18、光电探测器三19、分光镜三20、带通滤光片四21、聚焦透镜五 22、光电探测器四23、分光镜四24、带通滤光片六25、聚焦透镜六26、光电探测器五27、带通滤光片五28、偏振分光镜29、聚焦透镜七30、光电探测器六31、聚焦透镜八32、光电探测器七33、多通道数据采集卡34、计算机35。
具体实施方式
下面根据附图和优选实施例详细描述本发明，本发明的目的和效果将变得更加明白，应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
如图1所示，本发明的基于机载多波长激光雷达探测海洋光学剖面的系统包括多波长激光器1、激光扩束器2、望远镜接收系统3、聚焦透镜一4、分视场镜5、准直透镜一6、分光镜一7、带通滤光片一8、聚焦透镜二9、光电探测器一10、反射镜11、带通滤光片二12、聚焦透镜三13、光电探测器二14、准直透镜二15、分光镜二16、带通滤光片三17、聚焦透镜四18、光电探测器三19、分光镜三20、带通滤光片四21、聚焦透镜五22、光电探测器四 23、分光镜四24、带通滤光片六25、聚焦透镜六26、光电探测器五27、带通滤光片五28、偏振分光镜29、聚焦透镜七30、光电探测器六31、聚焦透镜八32、光电探测器七33、多通道数据采集卡34、计算机35。
多波长激光器1和激光扩束器2组成多波长激光发射系统，多波长激光器1用于发射四种不同波长的激光，作为优选，发射1064nm、532nm、486nm和355nm的激光，1064nm激光通道用来校正机载激光雷达系统与海面之间的高度，532nm激光通道优先探测二类水体光学特性，486nm激光通道优先探测一类水体光学特性，355nm激光通道优先探测海洋水体的荧光特性。多波长激光器1波长355nm激光输出的平均功率为3.8W，486nm激光输出的平均功率为0.3W，532nm波长输出的平均功率为1.7W，1064nm激光输出的平均功率为4.2W。激光扩束器2用于出射激光扩束并压缩其发散角。
望远镜接收系统3用于接收激光在海水中产生的后向散射信号，为了更好的提升机载多波长激光雷达的探测性能，其采用的是200mm直径的卡式望远镜，焦距为700mm。
聚焦透镜一4和分视场镜5组成分视场镜系统，聚焦透镜一4将望远镜接收系统3接收到的海水中产生的后向散射信号聚焦到分视场镜5，分视场镜将光束分成两束，一束进入分光系统一，另一束进入分光系统二。
准直透镜一6、分光镜一7和反射镜11组成分光系统一，准直透镜一6将光束准直后，分光镜一7将准直后的光束分为两束，一束进入光电探测系统一，另一束经反射镜11入射到光电探测系统二；
带通滤光片一8、聚焦透镜二9和光电探测器一10组成光电探测系统一，带通滤光片一 8将光束滤光后，聚焦透镜二9将滤光后的光束聚焦，光电探测器一10接收聚焦后的光束信号。作为优选，带通滤光片一8的中心波长为500nm，透过率大于90％，短波截止范围：200-472nm，长波截止范围：527-1200nm；光束经带通滤光片一8滤光后，后向散射光中只有486nm波长的散射光通过，486nm波长的散射光经聚焦透镜二9后由光电探测器一10接收。
带通滤光片二12、聚焦透镜三13和光电探测器二14组成光电探测系统二，带通滤光片二12将光束滤光后，聚焦透镜三13将滤光后的光束聚焦，光电探测器二14接收聚焦后的光束信号。作为优选，带通滤光片二12的中心波长为532nm，透过率大于90％，短波截止范围： 200-512nm，长波截止范围：552-1200nm；光束经过带通滤光片二12，后向散射光中只有532nm 波长的散射光通过，532nm波长的散射光经聚焦透镜三13后由光电探测器二14接收。
准直透镜二15和分光镜二16组成分光系统二，准直透镜二15将光束准直后，分光镜二 16将光束分为两束，一束进入光电探测系统三，另一束进入分光系统三。
带通滤光片三17、聚焦透镜四18和光电探测器三19组成光电探测系统三，带通滤光片三17将光束进行滤光后，聚焦透镜四18将滤光后的光束聚焦，光电探测器三19接收聚焦后的光束信号。作为优选，带通滤光片三17的中心波长为1064nm，透过率大于90％，短波截止范围：200-1039nm，长波截止范围：1089-1200nm。光束经过带通滤光片三17后，后向散射光中只有1064nm波长的散射光通过，1064nm波长的散射光经聚焦透镜四18后由光电探测器三19接收。
分光镜三20将光束分为两束，一束进入光电探测系统四，另一束进入分光系统四。
带通滤光片四21、聚焦透镜五22和光电探测器四23组成光电探测系统四，带通滤光片四21将光束滤光后，聚焦透镜五22对滤光后的光束聚焦，光电探测器四23接收聚焦后的光束信号。作为优选，带通滤光片四21的中心波长为355nm，透过率大于90％，短波截止范围： 200-345nm，长波截止范围：365-1200nm。光束经过带通滤光片四21后，后向散射光中只有 355nm波长的散射光通过，355nm波长的散射光经聚焦透镜五22后由光电探测器四23接收。
分光镜四24和带通滤光片五28组成分光系统四，分光镜四将光束分为两束，一束进入光电探测系统五，另一束经过带通滤光片五28滤光后进入偏振分光系统。作为优选，带通滤光片五28的中心波长为500nm，透过率大于90％，短波介质范围：200-472nm，长波截止范围：527-1200nm。光束经带通滤光片五28滤光后，后向散射中只有486nm波长的散射光通过，486nm波长的散射光进入偏振分光系统。
带通滤光片六25、聚焦透镜六26和光电探测器五27组成光电探测系统五，带通滤光片六将光束滤光后，聚焦透镜六对滤光后的光束聚焦，光电探测器五接收聚焦后的光束信号。作为优选，带通滤光片六25的中心波长为532nm，透过率大于90％，短波截止范围：200-512nm，长波截止范围：552-1200nm。光束经过带通滤光片六25后，后向散射光中只有532nm波长的散射光通过，532nm波长的散射光经聚焦透镜六26后由光电探测器五27接收。
偏振分光镜29将光束分成两束，一束进入光电探测系统六，另一束进入光电探测系统七。
聚焦透镜七30和光电探测器六31组成光电探测系统六，聚焦透镜七30将垂直偏振光聚焦后，光电探测器六31接收聚焦后的光束信号。
聚焦透镜八32和光电探测器七33组成光电探测系统七，聚焦透镜八32将光束聚焦后，光电探测器七33接收聚焦后的光束信号。
数据采集与处理系统包括多通道数据采集卡34和计算机35，两者均与所有的光电探测器电连接，多通道数据采集卡用于机载多波长海洋激光雷达数据的采集，计算机用于激光雷达数据的处理分析。
如图2所示，本发明的基于机载多波长激光雷达系统探测海洋光学剖面的反演方法，该方法包括如下步骤：
步骤一：多波长激光器发射1064nm、532nm、486nm和355nm的激光，为了更好的探测浅海区域和深海区域，分视场镜将光束分为大视场角的光束和小视场角的光束，大视场角的光束进入分光系统一，小视场角的光束进入分光系统二。作为优选，小视场光路的接收视场角范围为：0-6mrad，大视场光路的接收视场角范围为：6-25mrad。
步骤二：计算机对七个光电探测器采集到的机载多波长激光雷达数据进行预处理；
(1)基于飞机POS数据中的飞行时间、经纬度和激光雷达数据采集时间对七个光电探测器采集到的机载多波长激光雷达数据进行分离，保留多波长海洋激光雷达数据Sraw；
(2)将多波长海洋激光雷达数据Sraw中的脉冲数据进行平均处理，得到Smean；例如，将每50个脉冲数据进行平均处理，从而降低异常信号对反演结果的影响。
(3)对Smean中脉冲数据的背景噪声点进行平均处理，得到Snoise，然后基于公式Sdenoise＝Smean-Snoise得到去噪后的脉冲数据Sdenoise。例如，将后200个背景噪声点进行平均处理得到Snoise。
(4)利用1064nm通道采集的激光雷达数据计算机载激光雷达系统与海表之间的高度，并对所有通道的激光雷达数据进行距离校正，得到S(λ，z)；其中，λ为激光波长，z为海水深度；
步骤三：计算海洋水体激光雷达衰减系数Klidar(λ，z)；
(1)假设海洋水体后向散射系数与激光雷达衰减系数满足
β(λ，z)＝AKlidar(λ，z)R
其中A为相关参数，R为激光雷达后向散射消光对数比；
(2)对激光雷达信号S(λ，z)求导可以得到：
Figure BDA0003429913730000081
(3)把公式(1)带入公式(2)可得：
Figure BDA0003429913730000082
(4)取激光雷达探测的某边界值深度为zB，则激光雷达衰减系数表示为：
Figure BDA0003429913730000083
步骤四：计算海洋水体激光雷达后向散射系数β(λ，z)；
(1)假设激光雷达后向散射系数是不随深度变化的均匀水体光学参数和随深度变化的非均匀水体光学参数之和，则激光雷达回波信号和后向散射系数可以表示为：
S(λ，z)＝Sh(λ，z)+Sinh(λ，z), (5)
β(λ，z)＝βh(λzz)+βinh(λ，z), (6)
Sh(λ，z)＝K·βh(λ，z)·exp[-2Klidar(λ，z)] (7)
Sinh(λ，z)＝K·βinh(λ，z)·exp[-2Klidar(λ，z)] (8)
其中Sh(λ，z)为均匀水体部分产生的激光雷达回波信号，Sinh(λ，z)为非均匀水体产生的激光雷达回波信号，k为激光雷达系统常数，βh(λ，z)为均匀水体产生的后向散射系数，βinh(λ，z)为非均匀水体产生的后向散射系数。
均匀水体中激光雷达回波信号和衰减系数、后向散射系数之间的关系表示为：
lnSh(λ，z)＝lnK·βh(λ，0)-2Klidar(λ，0) (9)
因此，结合公式(5)-(9)，激光雷达后向散射系数表述为：
Figure BDA0003429913730000091
步骤五：分别计算大视场和小视场各个波段的海洋水体光学剖面参数Klidar(532nm，z)， Klidar(486nm，z)，Klidar(355nm，z)，β(532nm，z)，β(486nm，z)和β(355nm，z)。
根据上述计算得到的海洋水体光学剖面参数，可用于评估浮游植物粒径分布，并进一步评估海洋浮游植物群落分布。
图3是采用本发明的基于机载多波长激光雷达探测海洋光学剖面的系统及方法所反演得到的海洋水体激光雷达衰减系数结果，从图中可以清晰的看出海洋水体激光雷达衰减系数剖面分布。
图4是采用本发明提出的基于机载多波长激光雷达探测海洋光学剖面的系统及方法所反演得到的海洋水体激光雷达后向散射系数结果，从图中可以清晰的看出海洋水体激光雷达后向散射系数剖面分布。
本领域普通技术人员可以理解，以上所述仅为发明的优选实例而已，并不用于限制发明，尽管参照前述实例对发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实例记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在发明的精神和原则之内，所做的修改、等同替换等均应包含在发明的保护范围之内。

Claims (6)

1.一种基于机载多波长激光雷达探测海洋光学剖面的系统，其特征在于，该系统包括：
多波长激光发射系统、望远镜接收系统、分视场镜系统、分光系统一、分光系统二、分光系统三、分光系统四、光电探测系统一、光电探测系统二、光电探测系统三、光电探测系统四、光电探测系统五、光电探测系统六、光电探测系统七、偏振分光系统、数据采集与处理系统；
所述多波长激光发射系统包括多波长激光器和激光扩束器，多波长激光器用于发射四种不同波长的激光，激光扩束器用于出射激光扩束并压缩其发散角；
所述望远镜接收系统用于接收激光在海水中产生的后向散射信号；
所述分视场镜系统包括聚焦透镜一和分视场镜，所述聚焦透镜一将所述后向散射信号聚焦到分视场镜，分视场镜将光束分成两束，一束进入分光系统一，另一束进入分光系统二；
分光系统一包括准直透镜一、分光镜一和反射镜，准直透镜一将光束准直后，分光镜一将准直后的光束分为两束，一束进入光电探测系统一，另一束经反射镜入射到光电探测系统二；
光电探测系统一包括带通滤光片一、聚焦透镜二和光电探测器一，所述带通滤光片一将光束滤光后，聚焦透镜二将滤光后的光束聚焦，光电探测器一接收聚焦后的光束信号；
光电探测系统二包括带通滤光片二、聚焦透镜三和光电探测器二，所述带通滤光片二将光束滤光后，聚焦透镜三将滤光后的光束聚焦，光电探测器二接收聚焦后的光束信号；
分光系统二包括准直透镜二和分光镜二，准直透镜二将光束准直后，分光镜二将光束分为两束，一束进入光电探测系统三，另一束进入分光系统三；
光电探测系统三包括带通滤光片三、聚焦透镜四和光电探测器三，带通滤光片三将光束进行滤光后，聚焦透镜四将滤光后的光束聚焦，光电探测器三接收聚焦后的光束信号；
分光系统三为分光镜三，其将光束分为两束，一束进入光电探测系统四，另一束进入分光系统四；
光电探测系统四包括带通滤光片四、聚焦透镜五和光电探测器四，带通滤光片四将光束滤光后，聚焦透镜五对滤光后的光束聚焦，光电探测器四接收聚焦后的光束信号；
分光系统四包括分光镜四和带通滤光片五，分光镜四将光束分为两束，一束进入光电探测系统五，另一束经过带通滤光片五滤光后进入偏振分光系统；
光电探测系统五包括带通滤光片六、聚焦透镜六和光电探测器五，带通滤光片六将光束滤光后，聚焦透镜六对滤光后的光束聚焦，光电探测器五接收聚焦后的光束信号；
偏振分光系统为偏振分光镜，其将光束分成两束，一束进入光电探测系统六，另一束进入光电探测系统七；
光电探测系统六包括聚焦透镜七和光电探测器六，聚焦透镜七将光束聚焦后，光电探测器六接收聚焦后的光束信号；
光电探测系统七包括聚焦透镜八和光电探测器七，聚焦透镜八将光束聚焦后，光电探测器七接收聚焦后的光束信号；
数据采集与处理系统包括多通道数据采集卡和计算机，两者均与所有的光电探测器电连接，多通道数据采集卡用于机载多波长海洋激光雷达数据的采集，计算机用于激光雷达数据的处理分析。
2.根据权利要求1所述的基于机载多波长激光雷达探测海洋光学剖面的系统，其特征在于，所述望远镜接收系统采用200mm直径的卡式望远镜，焦距为700mm。
3.根据权利要求1所述的基于机载多波长激光雷达探测海洋光学剖面的系统，其特征在于，所述带通滤光片一的中心波长为500nm，透过率大于90％，短波介质范围：200-472nm，长波截止范围：527-1200nm；所述带通滤光片二的中心波长为532nm，透过率大于90％，短波截止范围：200-512nm，长波截止范围：552-1200nm；所述带通滤光片三的中心波长为1064nm，透过率大于90％，短波截止范围：200-1039nm，长波截止范围：1089-1200nm；所述带通滤光片四的中心波长为355nm，透过率大于90％，短波截止范围：200-345nm，长波截止范围：365-1200nm；所述带通滤光片五的中心波长为500nm，透过率大于90％，短波介质范围：200-472nm，长波截止范围：527-1200nm；所述带通滤光片六的中心波长为532nm，透过率大于90％，短波截止范围：200-512nm，长波截止范围：552-1200nm。
4.一种基于权利要求1所述的机载多波长激光雷达系统探测海洋光学剖面的反演方法，其特征在于，该方法包括如下步骤：
S1：所述多波长激光器发射1064nm、532nm、486nm和355nm的激光，分视场镜将光束分为大视场角的光束和小视场角的光束，大视场角的光束进入分光系统一，小视场角的光束进入分光系统二；
S2：计算机对七个光电探测器采集到的机载多波长激光雷达数据进行预处理；
(1)基于飞机POS数据中的飞行时间、经纬度和激光雷达数据采集时间对七个光电探测器采集到的机载多波长激光雷达数据进行分离，保留多波长海洋激光雷达数据Sraw；
(2)将多波长海洋激光雷达数据Sraw中的脉冲数据进行平均处理，得到Smean；
(3)对Smean中脉冲数据的背景噪声点进行平均处理，得到Snoise，然后基于公式Sdenoise＝Smean-Snoise得到去噪后的脉冲数据Sdenoise；
(4)利用1064nm通道采集的激光雷达数据计算机载激光雷达系统与海表之间的高度，并对所有通道的激光雷达数据进行距离校正，得到S(λ，z)；其中，λ为激光波长，z为海水深度；
S3：计算海洋水体激光雷达衰减系数Klidar(λ，z)；
假设海洋水体后向散射系数与激光雷达衰减系数满足β(λ，z)＝AKlidar(λ，z)R，取激光雷达探测的某边界值深度为zB，则激光雷达衰减系数表示为：
Figure FDA0003429913720000031
其中，A为相关参数，R为激光雷达后向散射消光对数比；
S4：计算海洋水体激光雷达后向散射系数β(λ，z)；
(1)假设激光雷达后向散射系数是不随深度变化的均匀水体光学参数和随深度变化的非均匀水体光学参数之和，则激光雷达后向散射系数β(λ，z)表述为
Figure FDA0003429913720000032
其中，
S(λ，z)＝Sh(λ，z)+Sinh(λ，z)
Sh(λ，z)＝k·βh(λ，z)·exp[-2Klidar(λ，z)]
Sinh(λ，z)＝k·βinh(λ，z)·exp[-2Klidar(λ，z)]
lnSh(λ，z)＝Ink·βh(λ，0)-2Klidar(λ，0)
其中，Sh(λ，z)为均匀水体部分产生的激光雷达回波信号，Sinh(λ，z)为非均匀水体产生的激光雷达回波信号，k为激光雷达系统常数，βh(λ，z)为均匀水体产生的后向散射系数，βinh(λ，z)为非均匀水体产生的后向散射系数；
S5：分别计算大视场和小视场各个波段的海洋水体光学剖面参数Klidar(532nm，z)，Klidar(486nm，z)，Klidar(355nm，z)，β(532nm，z)，β(486nm，z)和β(355nm，z)。
5.根据权利要求4所述的基于权利要求1所述的机载多波长激光雷达系统探测海洋光学剖面的反演方法，其特征在于，zB选取激光雷达的最大有效探测深度。
6.根据权利要求4所述的基于权利要求1所述的机载多波长激光雷达系统探测海洋光学剖面的反演方法，其特征在于，所述大视场角的范围为6-25mrad，小视场角的范围为0-6mrad。
CN202111591030.0A 2021-12-23 2021-12-23 基于机载多波长激光雷达探测海洋光学剖面的系统及方法 Active CN114674292B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111591030.0A CN114674292B (zh) 2021-12-23 2021-12-23 基于机载多波长激光雷达探测海洋光学剖面的系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111591030.0A CN114674292B (zh) 2021-12-23 2021-12-23 基于机载多波长激光雷达探测海洋光学剖面的系统及方法

Publications (2)

Publication Number Publication Date
CN114674292A true CN114674292A (zh) 2022-06-28
CN114674292B CN114674292B (zh) 2024-04-26

Family

ID=82070373

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111591030.0A Active CN114674292B (zh) 2021-12-23 2021-12-23 基于机载多波长激光雷达探测海洋光学剖面的系统及方法

Country Status (1)

Country Link
CN (1) CN114674292B (zh)

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105403877A (zh) * 2015-11-12 2016-03-16 中国科学院上海光学精密机械研究所 大动态范围光学分视场探测激光雷达
CN105486664A (zh) * 2015-12-31 2016-04-13 浙江大学 探测海洋浮游植物生物量和poc的激光雷达装置及方法
US20160274026A1 (en) * 2015-03-16 2016-09-22 Fredrick S. Solheim Characterizing tropospheric boundary layer thermodynamic and refractivity profiles utilizing multiband infrared observations
CN106199562A (zh) * 2016-07-06 2016-12-07 山东省科学院海洋仪器仪表研究所 基于机载激光雷达测量海底地形的海面误差校正方法
US20170030830A1 (en) * 2013-01-28 2017-02-02 Woods Hole Oceanographic Institution Submersible n-wavelength interrogation system and method for multiple wavelength interferometers
CN109187365A (zh) * 2018-08-22 2019-01-11 中国科学院上海光学精密机械研究所 海洋剖面双波长偏振光学参数观测仪
US20190025423A1 (en) * 2017-07-19 2019-01-24 Rezatec Limited Water network monitoring system
CN110673157A (zh) * 2019-11-07 2020-01-10 自然资源部第二海洋研究所 一种探测海洋光学参数的高光谱分辨率激光雷达系统
CN110673108A (zh) * 2019-09-25 2020-01-10 自然资源部第二海洋研究所 一种基于迭代Klett的机载海洋激光雷达信号处理方法
US20200209020A1 (en) * 2017-07-26 2020-07-02 Terra15 Pty Ltd Distributed optical sensing systems and methods
CN111595365A (zh) * 2020-07-06 2020-08-28 山东省科学院激光研究所 一种用于海洋温度及压力同步监测的多波长激光器
AU2020101889A4 (en) * 2020-08-19 2020-09-24 Guilin University Of Technology A three-channel optical receiver for airborne dual-frequency lidar
CN112034480A (zh) * 2020-07-16 2020-12-04 浙江大学 一种双波长海洋激光雷达探测的波长择优方法
CN112526482A (zh) * 2020-11-30 2021-03-19 北京遥测技术研究所 一种星载激光近海岸地形探测激光雷达及探测方法
CN112731341A (zh) * 2020-11-20 2021-04-30 自然资源部第二海洋研究所 基于双级虚像相位阵列的双波长高光谱分辨率激光雷达系统

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170030830A1 (en) * 2013-01-28 2017-02-02 Woods Hole Oceanographic Institution Submersible n-wavelength interrogation system and method for multiple wavelength interferometers
US20160274026A1 (en) * 2015-03-16 2016-09-22 Fredrick S. Solheim Characterizing tropospheric boundary layer thermodynamic and refractivity profiles utilizing multiband infrared observations
CN105403877A (zh) * 2015-11-12 2016-03-16 中国科学院上海光学精密机械研究所 大动态范围光学分视场探测激光雷达
CN105486664A (zh) * 2015-12-31 2016-04-13 浙江大学 探测海洋浮游植物生物量和poc的激光雷达装置及方法
CN106199562A (zh) * 2016-07-06 2016-12-07 山东省科学院海洋仪器仪表研究所 基于机载激光雷达测量海底地形的海面误差校正方法
US20190025423A1 (en) * 2017-07-19 2019-01-24 Rezatec Limited Water network monitoring system
US20200209020A1 (en) * 2017-07-26 2020-07-02 Terra15 Pty Ltd Distributed optical sensing systems and methods
CN109187365A (zh) * 2018-08-22 2019-01-11 中国科学院上海光学精密机械研究所 海洋剖面双波长偏振光学参数观测仪
CN110673108A (zh) * 2019-09-25 2020-01-10 自然资源部第二海洋研究所 一种基于迭代Klett的机载海洋激光雷达信号处理方法
CN110673157A (zh) * 2019-11-07 2020-01-10 自然资源部第二海洋研究所 一种探测海洋光学参数的高光谱分辨率激光雷达系统
CN111595365A (zh) * 2020-07-06 2020-08-28 山东省科学院激光研究所 一种用于海洋温度及压力同步监测的多波长激光器
CN112034480A (zh) * 2020-07-16 2020-12-04 浙江大学 一种双波长海洋激光雷达探测的波长择优方法
AU2020101889A4 (en) * 2020-08-19 2020-09-24 Guilin University Of Technology A three-channel optical receiver for airborne dual-frequency lidar
CN112731341A (zh) * 2020-11-20 2021-04-30 自然资源部第二海洋研究所 基于双级虚像相位阵列的双波长高光谱分辨率激光雷达系统
CN112526482A (zh) * 2020-11-30 2021-03-19 北京遥测技术研究所 一种星载激光近海岸地形探测激光雷达及探测方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
JAMES H.CHURNSIDE 等: "Inversion of oceanographic profiling lidars by a perturbation to a linear regression", 《APPLIED OPTICS》, vol. 56, no. 18, pages 5228 - 5233 *
LIU Q 等: "Desigh of the ship-borne multi-wavelength polarization ocean lidar system and measurement of seawater optical properties", 《29TH INTERNATIONAL LASER RADAR CONFERENCE》, vol. 237, 31 December 2020 (2020-12-31), pages 07007 - 1 *
刘金涛;张凯临;卢渊;郭金家;郑荣儿;: "用于浮游植物探测的三维激光诱导荧光光谱系统", 激光生物学报, no. 06, 15 December 2014 (2014-12-15), pages 70 - 74 *
唐军武 等: "海洋三维遥感与海洋剖面激光雷达", 《遥感学报》, vol. 25, no. 01, 25 January 2021 (2021-01-25), pages 460 - 500 *
李凯鹏 等: "双波长海洋激光雷达探测近岸到大洋水体的叶绿素剖面", 《中国激光》, vol. 48, no. 20, pages 2010002 - 1 *

Also Published As

Publication number Publication date
CN114674292B (zh) 2024-04-26

Similar Documents

Publication Publication Date Title
CN110095477B (zh) 测定微尘的分析仪
US7656526B1 (en) Lidar system for remote determination of calibrated, absolute aerosol backscatter coefficients
US4537507A (en) Dual beam maximum intensity laser sizing system
WO2020063073A1 (zh) 多波长测量大气二氧化碳浓度及气溶胶垂直廓线的激光雷达系统装置
Guo et al. Development of a single-wavelength airborne bathymetric LiDAR: System design and data processing
CN106772438A (zh) 一种全天时准确测量大气温度和气溶胶参数的激光雷达系统
CN109946710B (zh) 一种双波长多偏振激光成像装置
CN1987520A (zh) 气象与大气环境观测拉曼散射激光雷达系统
CN110673157B (zh) 一种探测海洋光学参数的高光谱分辨率激光雷达系统
CN112731341B (zh) 基于双级虚像相位阵列的双波长高光谱分辨率激光雷达系统
IL260362B2 (en) Device and method for increasing the dynamic range of a particle sensor
CN107976686B (zh) 一种多视场角海洋激光雷达及其视场角择优方法
CN110940992B (zh) 可提高激光雷达探测距离和精度的信号检测方法和系统
CN113281774B (zh) 高效紧凑的高光谱分辨率激光雷达系统、气溶胶后向散射系数和消光系数的获取方法
Li et al. Co-path full-waveform LiDAR for detection of multiple along-path objects
CN113219480B (zh) 基于扫描激光雷达的气溶胶消光系数三维空间分布反演方法
CN114167437B (zh) 测水激光雷达多通道设计方法
Wang et al. Development of ZJU high-spectral-resolution lidar for aerosol and cloud: Feature detection and classification
CN115856933A (zh) 高通量光子计数激光雷达成像装置及方法
CN114674292B (zh) 基于机载多波长激光雷达探测海洋光学剖面的系统及方法
CA2228499C (en) Optical measurement of marine conditions
Cester Atmospheric aerosol monitoring at the Pierre Auger Observatory
CN115773864A (zh) 一种基于光腔衰荡技术的高反射光学元件总积分散射的测量方法
Nicolae et al. Laser remote sensing of tropospheric aerosol
Somekawa et al. Depolarization light detection and ranging using a white light LIDAR system

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