摘要
星载单光子激光雷达以超高灵敏度和超高重频的优势,在海洋探测领域展现了广泛的应用前景。雷达系统中的单光子探测器件具有极高的灵敏度,可以探测光子量级的回波信号,同时也极易受太阳背景光噪声影响。由于背景噪声直接影响激光雷达的工作性能,还会对星上原始的数据量产生影响,在卫星系统设计阶段,对噪声强度的准确估计至关重要。本文综合考虑大气后向散射、水面反射及水体后向散射的贡献,建立了一个星载单光子激光雷达海洋噪声的估计模型。以全球首个对地观测星载激光雷达ATLAS为例,在输入系统参数和环境参数后,模型估计的噪声与ATLAS实测噪声误差在15%以内,证实了该噪声模型的正确性。
星载激光雷达是一种通过主动发射激光脉冲并接收回波信号,得到地表或海表精确高程信息的卫星遥感测量设备。随着国外第一代激光测高卫星ICESat(The Ice, Cloud, and land Elevation Satellite)以及国内的高分七号等卫星的发
受益于单光子探测器件极高的灵敏度,星载单光子激光雷达可以发射远低于线性体制激光雷达单脉冲能量的激光脉冲,从而在总能量受限的约束下极大地提高了激光器重复频率,由此带来了极高的沿轨测量分辨率。2018年底美国国家航空航天局(NASA,National Aeronautics and Space Administration)发射了搭载新一代单光子激光雷达ATLAS(Advanced Topographic Laser Altimeter System)系统的卫星ICESat-2(Ice, Cloud, and land Elevation Satellite-2),并在海面波浪信息提取、海洋水体参数获取、近岸测深等诸多应用领域取得了巨大的成
目前,星载激光雷达的噪声估计主要采用经典的线性体制噪声模
ATLAS系统主动发射激光脉冲,并接收由地球表面反射回的激光信号,通过准确地记录发射脉冲时间和接收时间,可以准确地获取卫星与地球表面的距离。然而,在信号开窗的距离门内,太阳背景光将持续不断地进入激光雷达接收系统,并被单光子探测器响应,大大增加了信号识别和提取的难度。如
, | (1) |

图 1 星载激光雷达海洋区域噪声组成示意图
Fig 1 Schematic diagram of noise composition in oceanic spaceborne lidar
ICESat-2的单个探测器暗计数噪声率约为400 H
对于粒子大小远小于光波波长的气体分子,一般用瑞利散射近似描述其散射特征,可以用单次散射近似描述瑞利散射引入的噪声强
, | (2) |
这里F是星载激光雷达系统标定系数,η是系统整体接收效率,包含光学系统接收效率与探测器量子效率的影响,Nλ是太阳在大气层外的平均辐照度,θr是接收系统半视场角,Ar是望远镜口径大小,hv是单个光子能量,Pr(θ)为瑞利散射相位函数,θ为瑞利散射角,角度θs和ϕs分别是太阳所处位置的天顶角和方位角,θv和ϕv分别是星载激光雷达系统的天顶角和方位角。对于接近正天顶发射激光的ATLAS系统,有θ±=±θs,r(θ)是入射角为θ时水气界面的菲涅耳反射系数。ICESat-2 ATL03产品提供了卫星及太阳的天顶角和方位角。这里Pr(θ+)和Pr(θ-)分别代表散射角θ取θ+和θ-时的散射相位函数。
这里τr为瑞利光学厚度,主要取决于激光波长λ和海表大气压P。瑞利光学厚度与大气压P成正比关系,在标准大气压和标准温度下(1013.25 hPa, 288.15 K),瑞利光学厚度τr0可以表示
, | (3) |
任意大气压下的瑞利光学厚度τr可以从标准大气压瑞利光学厚度τr0转化
, | (4) |
这里P是海平面大气压,P0是标准大气压。需要强调的是,对于海洋区域,海平面的大气压P很接近标准大气压P0,因此海洋区域瑞利光学厚度τr可以近似为标准大气压下的光学厚度τr0。
对于较大的散射粒子,如尘埃、烟雾、霾、云等,一般用气溶胶散射(或米氏散射)描述其散射特征。同样地,可以用单次散射近似描述气溶胶粒子散射引入的噪声强
, | (5) |
这里气溶胶散射相位函数Pa(θ)可以用TTHG(Two-Term Henyey-Greenstein)函数表示,g1和g2分别代表前向HG函数和后向HG函数的非对称因子。而f(θ,g)代表HG散射相位函数,g大于0时为前向散射偏强,g小于0时为后向散射偏强,公式中的f(θ,g1)和f(θ,g2)分别代表g取g1和g2的HG函数。wa是气溶胶单次散射比率, τa是气溶胶光学厚度,其中单次散射比率与气溶胶类型和湿度相关,可以表示
, | (6) |
这里AM是气溶胶类型,RH为大气相对湿度。气溶胶散射与瑞利散射不同,不同地方的云、雾的厚度会极大程度上影响气溶胶散射噪声大小,一般可通过气象数据获得当地的各项气溶胶参数。

图 2 瑞利散射与气溶胶散射相位函数
Fig 2 The phase function of Rayleigh scattering and aerosol scattering
水面反射噪声主要包括:水面白沫反射产生的背景光噪声以及水面镜面反射产生的背景光噪声。由于海面存在风浪,水面上部分区域会被泡沫覆盖,一般来说,水面泡沫部分反射特性可以视为漫反射,当太阳天顶角为θs,星载激光雷达收到的白沫反射的噪声光子率可以表示为
, | (7) |
这里W为白沫覆盖比例,海洋表面白沫覆盖比例可以与风速之间建立关联[
在非白沫覆盖区域,水表面反射主要由表面随机分布的各个小面元的镜面反射组成,只有海面点刚好能将太阳入射光反射到卫星方向,才能被接收系统接收。海洋反射截面与水面粗糙度密切相关。对于ICESat-2系统,激光入射天底角(也即接收系统天底角)仅为0.33°,海表面反射的噪声可以近似表示
, | (8) |
这里ρs是水面菲涅尔反射率,Ta代表正天顶方向的大气直射透过率。
, | (9) |
在532 nm的光波波段,太阳光可以穿过水面进入水体。由于水体的后向散射,部分太阳光被再次反射出水面而被激光雷达探测器接收到。在被动遥感领域,来自水面以下的水体太阳光反射被称为离水辐亮度,一般用遥感反射率Rrs表示水面以下(包括水体和水底)反射太阳光的能力。遥感反射率Rrs要受水体光学参数及水深影响,在大洋区域,主要取决于水体后向散射系数和水体吸收系数,本研究中不作详细讨论。遥感反射率Rrs可以用来表示S水体的反射能力,该参数可以直接从被动卫星遥感数据源获取,如MODIS(Moderate Resolution Imaging Spectroradiometer)。水面以下的太阳光反射项可以表示为
, | (10) |
ICESat-2卫星发射532 nm波长的绿光,在垂轨方向分成3个强波束和3个弱波束。在500 km的轨道高度上,以10 kHz的重频产生约0.7 m的沿轨光斑间隔。如

图 3 研究区域与所选的四个ICESat-2轨道概览:(a)位于澳大利亚北部附近;(b)位于澳大利亚北部附近;(c)位于英国东部附近;(d)位于印度西部附近。红色线段代表ICESat-2轨道
Fig 3 Overview of the study areas and the corresponding four ICESat-2 tracks: (a) Near north Australia; (b) Near the north of Australia; (c) Near the eastern of the United Kingdom; (d) Near the Western of India. The red lines represent the location of the ICESat-2 tracks
在输入激光雷达系统参数,和对应的实测环境参数后,模型估计了ICESat-2在运行时沿轨分布的噪声率。
, | (11) |
这里ftheo为理论总噪声,fmea为实测总噪声,n为每一个轨道的噪声样本数量。从

图 4 海洋区域的噪声模型验证:(a-d)四个图的噪声率分别对应图3中(a-d)四条轨道的沿轨噪声率。图中蓝色实现代表ICESat-2实测噪声,红色实现代表理论总噪声,不同颜色的虚线分别代表理论不同噪声项
Fig 4 Validation of the noise model in ocean areas: (a-d) represent the result of the noise rates for four tracks in Fig 3 (a-d). The blue solid line represent the ICESat-2 measured noise rate, the red solid line represent the theoretical noise rate , and the dashed line represent the different contribution of the noise rate
从
本文综合考虑了大气瑞利散射、气溶胶散射、海面镜面反射、海面白沫反射、水体后向散射等诸多太阳光反射项,较为准确地建立了星载单光子激光雷达海洋噪声模型。该模型综合考虑了太阳天底角、气溶胶厚度、海面风速、气压、湿度等环境因素的影响,通过输入系统参数和环境参数即可完成星载单光子激光雷达的噪声估计。
以ICESat-2卫星单光子激光雷达系统的实际参数以及NCEP和MODIS提供的环境参数为输入,对ICESat-2系统4个轨道的海面沿轨噪声进行了估计。在清澈的大洋区域,大气后向散射噪声占据噪声中的主体地位,瑞利散射引入的噪声较为稳定,气溶胶散射沿轨变化较大。在近岸区域,水体变得更加浑浊,水体后向散射噪声变强,对噪声水平产生了较大影响。将模型估计的总噪声与ICESat-2卫星单光子激光雷达系统的海面实测总噪声进行了对比,实测结果与理论模型估计结果吻合较好,整体误差在10%左右,最大误差不超过15%,验证了噪声估计模型的正确性。
本文提出的星载单光子激光雷达海洋噪声模型可以计算不同系统参数和环境条件下星载单光子激光雷达的海面噪声,进而可以分析不同环境下的信噪比,评估激光载荷的海洋探测能力。另外,该模型也可以用于估计星载单光子激光雷达的数据量,协助激光载荷参数的闭合论证及优化设计。
References
TANG Xin-Yue, XIE Jun-Feng, MO Fan, et al. GF-7 dual-beam laser altimeter on-orbit geometric calibration and test verification [J]. Acta Geodaetica et Cartographica Sinica.唐新明, 谢俊峰, 莫凡,等.高分七号卫星双波束激光测高仪在轨几何检校与试验验证 [J].测绘学报, 2021, 50(03): 384-395. [百度学术]
TANG Xin-Yue, CHEN Ji-Yi, LI Guo-Yuan, et al. Error analysis and preliminary pointing angle calibration of laser altimeter on Ziyuan-3 02 satellite[J]. Geomatics and Information Science of Wuhan University.唐新明, 陈继溢, 李国元, 等.资源三号02星激光测高误差分析与指向角粗标定 [J]. 武汉大学学报: 信息科学版, 2018, 43(11): 9. [百度学术]
Ye L, Gu G, He W, et al. A Real-time restraint method for range walk error in 3-D imaging lidar via dual detection [J]. IEEE Photonics Journal, 2018, 10(2): 1-1. [百度学术]
Lancaster R S, Spinhirne J D, Palm S P. Laser pulse reflectance of the ocean surface from the GLAS satellite lidar [J]. Geophysical Research Letters, 2005, 32(22): L22S10. [百度学术]
Winker D M, Hunt W H, Mcgill M J. Initial performance assessment of CALIOP [J]. Geophysical Research Letters, 2007, 34(19): 228-262. [百度学术]
Ma Y, Zhang W, Li S, et al. A new wind speed retrieval method for an ocean surface using the waveform width of a laser altimeter[J]. Canadian Journal of Remote Sensing, 2017, 43(4): 309-317. [百度学术]
Neumann, Thomas A , Martino, et al. The ice, cloud, and land elevation satellite-2 mission: A global geolocated photon product derived from the advanced topographic laser altimeter system[J]. Remote Sensing of Environment, 2019, 233: 111325. [百度学术]
Parrish C E, Magruder L A, Neuenschwander A L, et al. Validation of ICESat-2 ATLAS bathymetry and analysis of ATLAS's bathymetric mapping performance[J]. Remote Sensing, 2019, 11(14): 1634. [百度学术]
Neuenschwander A, Guenther E, White J C, et al. Validation of ICESat-2 terrain and canopy heights in boreal forests [J]. Remote Sensing of Environment, 2020, 251: 112110. [百度学术]
Lu X, Hu Y, Yang Y, et al. Antarctic spring ice-edge blooms observed from space by ICESat-2 [J]. Remote Sensing of Environment, 2020, 245: 111827. [百度学术]
Ma Y, Xu N, Liu Z, et al. Satellite-derived bathymetry using the ICESat-2 lidar and Sentinel-2 imagery datasets [J]. Remote Sensing of Environment, 2020, 250: 112047. [百度学术]
Degnan J. Photon-counting multikilohertz microlaser altimeters for airborne and spaceborne topographic measurements [J]. Journal of Geodynamics, 2002, 34(3-4): 503-549. [百度学术]
Martino A J, Neumann T A, Kurtz N T, et al. ICESat-2 mission overview and early performance [C]. Sensors, Systems, and Next-Generation Satellites XXIII. SPIE, 2019. [百度学术]
Gardner C S. Target signatures for laser altimeters: An analysis [J]. Applied Optics, 1982, 21(3): 448-453. [百度学术]
Yang G, Martino A J, Lu W, et al. IceSat-2 ATLAS photon-counting receiver: Initial on-orbit performance[C]. Advanced Photon Counting Techniques XIII. 2019. [百度学术]
Gordon H R, Brown J W, Evans R H. Exact Rayleigh scattering calculations for use with the Nimbus-7 Coastal Zone Color Scanner [J]. Applied Optics, 1988, 27(5): 862. [百度学术]
Bodhaine B A, Wood N B, Dutton E G, et al. On Rayleigh optical depth calculations[J]. Journal of Atmospheric and Oceanic Technology, 1999, 16(11): 1854. [百度学术]
Gordon H R, Wang M. Retrieval of water-leaving radiance and aerosol optical thickness over the oceans with SeaWiFS: a preliminary algorithm [J]. Applied Optics, 1994, 33(3): 443. [百度学术]
Gregg W W, Carder K L. A simple spectral solar irradiance model for cloudless maritime atmospheres [J]. Limnology and Oceanography, 1990, 35(8): 1657-1675. [百度学术]
Monahan E C, Muircheartaigh, Iognáidó. Optimal power-law description of oceanic whitecap coverage dependence on wind speed[J]. Journal of Physical Oceanography, 1980, 10(12): 2094-2099. [百度学术]
Cox C, Munk W. Measurement of the roughness of the sea surface from photographs of the sun's glitter [J]. Journal of the Optical Society of America (1917-1983), 1954, 44(11): 838. [百度学术]
And C R J, Alpers W. The role of the critical angle in brightness reversals on sunglint images of the sea surface [J]. Journal of Geophysical Research: Oceans, 2010, 115(C9): C09019. [百度学术]