-
卫星激光测高作为星载激光雷达的一个子类，因其能获得亚米级甚至厘米级高程测量精度，已成为一种新型对地观测技术并受到广泛的关注[1-2]。2003、2018年美国国家航空航天局（National Aeronautics and Space Administration，NASA）先后发射了冰、云和陆地高程卫星ICESat （Ice、Cloud and land Elevation Satellite）和ICESat-2，在极地冰盖厚度变化、全球树高和生物量反演、全球高程控制点获取、湖泊水位监测等方面取得了很多科研应用成果[3-4]，其中ICESat-2搭载了“先进地形激光测高系统”（Advanced Topographic Laser Altimeter System，ATLAS），采用了新型光子体制，波束数量和激光点密度均比GLAS有较大提高[5-6]。NASA还在2018年12月发射了“全球生态系统动力学调查”(Global Ecosystem Dynamics Investigation Lidar，GEDI)[7]，成功将其布置在国际空间站（International Space Station，ISS）上，主要用于获取热带和温带区域地表高程、森林垂直结构信息。我国先后在2016、2019年发射的资源三号02星、高分七号卫星上配备了激光测高仪[8-9]，主要用于高程控制点的获取[10-11]，并计划在2021年发射的陆地生态系统碳监测卫星上搭载多波束的线性体制全波形激光测高仪。
不同于LOLA (Lunar Observer Laser Altimeter)、MOLA (Mars Orbiter Laser Altimeter)对月球及火星等地外星球激光测高探测[12-13]，对地观测卫星激光测高除受复杂地形影响外，大气对激光测距及最终的测高精度影响不容忽视[14-15]。目前，对于大气折射延迟分析及其改正已经开展了很多研究，相关算法基本成熟且能工程化应用[16-17]。但大气散射的影响及改正目前还处于探索阶段，离业务化应用还有较大差距。在大气传输过程中，云、气溶胶等粒子的吸收、散射作用，除导致激光束质心发生一定变形外，还会产生测距偏差[18]，而目前大气散射改正没有实现业务化，导致因大气散射而发生质量退化的激光点，要么被直接抛弃降低了总的可利用率，要么被应用但降低了精度的可靠性。因此，开展大气散射改正对于有效保证激光点的可靠性、提高总的可利用率具有重要价值。
文中对国内外激光测高载荷的发展现状及技术特点进行了描述，梳理了这些卫星所采用的大气探测算法、大气散射延迟改正方案，在此基础上进行对比分析和总结归纳，为我国后续激光测高卫星大气散射改正技术的发展提出了建议。
-
在大气传输过程中，卫星对地发射的激光脉冲受到云、气溶胶等粒子散射影响产生测距偏差。对大气散射影响的探索是一个漫长的过程，Platt等人[19]在1981年首次强调了卷云中多次散射影响，在此基础上，Eloranta[20]发展了计算激光雷达回波信号中多次散射影响的实际模型。Bissonette[21]利用辐射传输模型和实验室测量数据，模拟了多次散射的影响。Bruscaglioni等人[22-23]对多次散射进行了详细的实验论证，表明了大气散射在激光雷达回波信号中的重要作用。Mahesh等人[24-25]利用蒙特卡洛模型模拟计算了激光测高过程中云对光子多次散射的影响，结果表明，对于薄云和低层云，大气多次散射导致的路径延迟大于20 cm，但是仅涉及仿真结果的模拟及其分析。Brenner等人[26]通过对比雷达测高系统Envisat和地学激光测高系统GLAS的测量数据发现，大气的多次散射效应是影响激光测距精度的重要因素。渠丽新等人[27]模拟分析了云对激光脉冲展宽效应的影响，实验结果表明，与云层厚度、单次散射比等因素有关。宋雪平等人[28]对云雾多次散射条件下所接收的激光数据进行了评估，发现多次散射对系统稳定性有较大的干扰。李颖颖等人[29-30]通过蒙特卡洛模拟反演多次散射条件下激光雷达消光系数等大气参数，结果表明，大气参数反演必须要考虑散射的影响。陈舒杭等人[31]基于半解析蒙特卡洛模拟反演特定条件下云多次散射对激光测高数据精度的影响，结果表明，测高偏差与云层高度、云有效粒子半径、云光学厚度、浓度等多种因素有关，必须基于实时的大气观测参数和半解析蒙特卡洛法修正测高偏差，以提高测距值，但实际观测中很难同步获取实时大气相关数据。周辉等人[32-33]通过蒙特卡洛模拟激光在云中单次散射情况下的回波波形，发现波形会发生拖尾变形而影响测距精度。王婷等人[34]基于蒙特卡洛模拟分析多个红外波长激光在雾中的衰减效应，结果表明，1 064 nm激光在平流雾中具有很好的传输性，雾对激光的衰减效应与消光系数、散射系数等多个因素密切相关。么嘉棋等人[35-36]对ICESat的大气范围波形数据进行分析并反演了云光学厚度等参数，为后续开展大气散射延迟改正提供了一定参考。Duda、Palm等人[37]基于GLAS参数的蒙特卡洛模拟实验表明，云光学厚度（Cloud Optical Depth，COD）为0.2的云会引起20 cm的测高偏差。Yang等人[38-39]基于ATLAS参数的三维蒙特卡洛模拟实验表明，COD为1的云会引起4~6 cm的测高偏差。
从上述研究可以看出，大气散射延迟影响对于分米级甚至厘米级的高精度激光测量而言，是不容忽视的、必须加以考虑并修正。进入21世纪，国内外先后发射了多种型号的对地观测激光测高载荷，表1列出了相关的技术指标。国产的资源三号02星和高分七号上的激光载荷由于设计及数传等原因，没有获取大气段回波数据，因此不具备大气探测的可能性。CALIOP、CATS获取了大气传输范围内数据，由于主要探测对象并不是地表高程，并没有进行大气散射校正，但采用了各种大气参数进行数据质量控制。GLAS采用动态阈值法实现大气探测，并基于蒙特卡洛模拟结果实现了大气散射误差校正。GEDI通过多参数约束方案，将受大气散射影响的数据准确识别，实现数据质量控制。对于光子计数体制激光雷达而言，多波束、高重频是非常典型的优势，但也为大气散射改正提出新的挑战：如何在兼具技术特点、大气探测范围等因素的基础上实现误差改正[40]。对此，CATS、ATLAS采用了不同的解决方案。CATS基于阈值法实现大气探测，并基于多参数约束算法获取云相态、气溶胶类型、光学厚度等大气参数。由于高激光重频，数据本身无法实现大气探测，ATLAS采用第三方地表反射率数据集实现了大气散射探测，但是反演结果不包含大气垂直分布、层高度数据，对大气散射改正造成很大的困扰。
表 1 星载激光测高载荷主要技术指标
Table 1. Main technical parameters of space-borne laser altimeter
Payload GLAS CALIOP CATS ELA GEDI ATLAS SLA CASAL Year 2003 2006 2015 2016 2018 2018 2019 2021 Satellite ICESat-1 CALIPSO ISS ZY3-02 ISS ICESat-2 GF-7 TECSIMS Detection target Surface elevation Cloud/Aerosol Cloud/Aerosol Elevation reference point Surface elevation and Forest biomass Surface elevation Elevation reference point Height of forest and atmosphere Orbit height/km 600 700 405 505 405 500 500 506 Telescope diameter/m 1.0 1.0 0.6 0.21 0.8 0.8 0.6 — Wavelength/nm 532,1064 532,1064 355,532,1064 1064 1064 532 1064 532,1064 Laser repetition frequency/kHz 0.04 0.02 4 0.002 0.242 10 0.0030.006 VEG：0.04
ATM：0.02Single pulse energy/mJ 35-75 110 2-3 200 10 0.04-0.12 100-180 VEG：70
ATM：110Laser footprint diameter/m 70 — — 75 25 17.5 <30 25-30 Number of beams 1 1 1 1 4 6 2 VEG：5
ATM：1Data acquisition mode Linear mode Linear mode Single-photon mode Linear mode Linear mode Single-photon mode Linear mode Linear mode Resolution of atmospheric
data /m76.8 30 60 — 0.15 30 — 30 Atmospheric detection capability Yes Yes Yes No Yes Yes No Yes Optical depth inversion algorithm Threshold method +slope method Threshold method + slope method Threshold method + slope method None Threshold method + slope method Other data Laser footprint image Threshold method + slope method Correction method of atmospheric scattering error Monte Carlo simulation Monte Carlo simulation Monte Carlo simulation None Monte Carlo simulation Other data Pending research Pending research -
目前的对地观测激光载荷在大气散射探测方面，主要有基于后向散射强度的动态阈值探测、多参数约束的直接大气探测以及基于表观反射率的间接大气探测等3种方法。
-
由于后向散射效应的叠加性，云或气溶胶区域内的整体后向散射系数远远大于天气晴朗区域。基于这种特性，可以通过设置一个合理的阈值进行云或气溶胶探测。以GLAS为例，GLAS大气通道以76.8 m的垂直分辨率记录了从地面到海拔约41 km的后向散射廓线数据。基于后向散射系数经验阈值判断层结构为云或气溶胶，之后定量反演光学厚度，最后通过蒙特卡洛模拟实现大气散射误差校正。GLAS采用不同的阈值算法将数据分别保存成4种产品（40、5、1、0.25 Hz），此节主要介绍了5 Hz产品所采用的大气探测算法，算法流程如图1所示。
图 1 云/气溶胶层搜索算法流程
Figure 1. Cloud & aerosol search algorithm flow
如果3个连续的后向散射信号大于阈值，则认为在该高度存在云、气溶胶。计算激光脉冲在当前高度的后向散射系数（βm）和该海拔高度的散射常数（Ψ）之和，并乘以与当前高度有关的因子（F），通过下式可以计算出层高（Zt）：
$$ {Z}_{t}=\left({\beta }_{m}\right({\textit{z}})+{\textit{Ψ}})\cdot F $$ (3) 式中：Ψ=1.0×10−6；F=(0.70+(i_bin /2)/180.0)，i_bin指按大气探测范围每隔76.8 m所划分的节点编号。斜率法通过计算层高度范围内有效消光系数的积分，得到激光光斑范围内柱状界面的光学厚度。斜率法计算流程如下：
斜率法反演大气光学厚度步骤如下：
$$ {S_\alpha }({\textit{Z}})={\alpha _\alpha }({\textit{Z}}){\beta _\alpha }({\textit{Z}}) $$ (4) 式中：z为激光脉冲所在海拔高度；Sα(Z)为大气对激光的削弱系数（或称有效消光系数），表示在Z高度散射的全部能量）；βα(Z)为大气后向散射系数；αα(Z)为大气消光系数。
$$ \tau =\int\limits_{{{\textit{Z}}_{{\rm{bottom}}}}}^{{{\textit{Z}}_{{\rm{top}}}}} {{S_\alpha }({\textit{z}})} {\rm{d}}{\textit{z}} $$ (5) 式中：τ为光学厚度；Ztop为层顶高度；Zbottom为穿透层时的高度。至此，实现了基于回波波形数据的光学厚度参数反演，反演流程见图2。
图 2 光学厚度反演流程
Figure 2. Optical depth inversion process
-
GEDI以缓冲区的形式保存大气传输范围内的回波数据，以区别于地物回波并保存在L2产品数据集中。通过结合波形灵敏度、回波能量、土地覆盖类型等多个参数，实现大气探测。
由于大气条件从完全晴朗到多云，GEDI返回信号强度会有很大的变化，穿透遮挡的能力取决于返回信号的强度。当激光穿透云层、比较密集的树冠时，信号强度会迅速下降，甚至衰减至噪声水平。为了识别没有检测到真实地面的波形，GEDI提出了波形灵敏度（Waveform Sensitivity）指标，由公式（4）得到。
$$ {\rm{WS}}=1.0 - \frac{{min\_{\rm{detection\_energy}}}}{{rx\_{\rm{energy}}}} $$ (6) 式中：WS表示波形敏感度；min_detection_energy表示回波高斯拟合最小积分面积；rx_energy表示激光脉冲的累积能量。为了保证使用用户可以轻松的判断错误数据、较低质量数据，GEDI使用几个单独的质量评价标准组合得到总体质量控制标准，判断数据质量良好（quality_flag=1）的流程如图3所示[41]。
图 3 GEDI大气探测算法
Figure 3. Atmospheric detection algorithm of GEDI
图中，Sensitivity_a<n>表示波形敏感度。rx_processing_a<n>是指有效波形数据，其中<n>表示用于检测信号的不同参数。min_detection_energy是指通过有效波形计算的高斯分布下可检测最小积分面积，是Sensitivity_a<n>下的子类别。rx_assess是指接收激光脉冲的特征参数（例如，最小和最大幅度、能量），保存在L2A中。rx_energy是指接收激光脉冲的累积能量，是rx_assess下的子类别。在满足上述条件的基础上，根据激光点所在土地类型，进行进一步判断。如果满足上述条件，那么该数据质量判断为良好。图4是GEDI的一段真实数据，黑色为原始数据中经质量控制后筛掉的点，红色为经筛选后质量良好的点（quality_flag=1），可以看出云散射数据、背景噪声被有效过滤掉。
图 4 GEDI数据质量控制示意图
Figure 4. Data quality control diagram of GEDI
-
线性体制激光雷达可以通过接收大气后向廓线信号的变化判断是否受到云散射影响，但是对于ICESat-2/ATLAS光子计数体制而言，激光重频达到10 KHz，会出现激光脉冲混叠现象，无法采用常规的大气探测算法[42]。ATLAS采用基于表观反射率的大气探测算法[25]，当激光受到大气散射影响时，表观地表反射率数据与清洁大气时有明显的区别，通过设置合适的阈值可以实现大气探测，进而分析散射影响。
基于表观反射率ASR （Apparent Surface Reflectance，ASR）的大气探测算法流程如下：从25 Hz观测数据中分割出对应激光点每四秒数据的ASR；依据公式（5）计算大气探测阈值；依据公式（6）计算云、气溶胶现象因子，如果P大于40%则存在云、气溶胶，如果P小于等于40%则为清洁大气。
$$ T=AC{L_{{\rm{Rtru}}}}\cdot {\textit{Φ}} $$ (7) 式中：T为大气探测阈值；ACLRtru为校正之后的晴空ASR；Φ是用来校正潜在晴空ASR偏差的校正系数，如果地表覆盖类型为陆地则取值为0.9，如果覆盖类型为海洋则取值为1.0。
$$ P=\left( {1 - ASR/T} \right) \cdot 100 $$ (8) P为大气现象因子，根据该参数可实现大气探测。在算法流程中还存在其他大气探测方法，如果使用ATLAS数据使用该算法实现大气探测，则将其记录到ATL09参数Use_ATLAS中，0代表未使用该算法，1代表使用该算法实现云检测。在已知表面反射率的前提下，可计算大气柱状截面总体光学厚度，前提是返回信号强度不能被完全衰减。ICESat-2无法穿透光学厚度为3的云或者气溶胶，因此，该技术仅限于大气总光学深度之和小于3的情况[38]，超过极限阈值，返回信号将太小或完全衰减。
-
光子的平均路径延迟，与散射能量的分布息息相关，而且光子可能在某一高度发生多次散射，在多种延迟路径下返回接收视场角范围内。结合光子散射、衍射理论以及卫星技术参数，Duda等人[37]提出大气散射条件下平均延迟距离计算方案，如公式(7)~(9)所示。公式（7）表示光子散射理论框架下能量分布方程，公式（8）表示光子能量高斯分量方程，公式（9）表示光子能量平行分量方程。
$$ S={S_1}{e^{ - 2\tau }}(1 + {f_{\rm{i}}}\tau + {f_{\rm{g}}}\tau) $$ (12) $$ {f_{\rm{g}}}=1 - \exp \left( - \frac{{2{\delta _{\rm{m}}}}}{{{\textit{z}}\theta _{\rm{s}}^2}}\right) $$ (13) $$ {f_{\rm{i}}}=\frac{1}{2}(1 - \cos ({\theta _{\rm{m}}})) $$ (14) 式中：S为激光初始能量值；S1为散射之后能量大小；τ为大气光学厚度；θm在高度为Z时视场角所能允许的最大散射角度；θs向前衍射峰的宽度；δm指在当前轨道高度、接收视场角下理论上最大的延迟距离。fg、fi分别表示在垂直方向、同向的能量分量。
综合考虑垂直方向上高斯分量、同向分量，可以将平均延迟距离公式（2）改写为：
$$ \overline {{\delta _{\rm{t}}}} =\tau \frac{{\overline {{\delta _{\rm{g}}}} {f_{\rm{g}}} + \overline {{\delta _{\rm{i}}}} {f_{\rm{i}}}}}{{1 + {f_{\rm{g}}}\tau + {f_{\rm{i}}}\tau }} $$ (15) Duda等人在上述原理的基础上，经过多次蒙特卡洛模拟实验，分析得到了GLAS数据多次散射误差与多个参数之间的关系，并将OD范围在0.45~1.1的数据归纳为多个查找表。通过数据反演得到的云光学厚度值在查找表中寻找对应的改正值，实现大气散射误差校正。
-
由大气前向散射引起的激光测距误差改正一般是结合实验室标定的蒙特卡洛模拟结果进行修正，但是蒙特卡洛模拟涉及多个参数的选取，实际生产中无法同步获取这些参数，不能实现误差改正业务化。除了蒙特卡洛模拟之外，在GLAS实际案例中，统计结果表明：在低光学厚度范围内，光学厚度与大气散射引起的测距偏差统计分布规律符合指数函数。基于此，笔者提出了一种基于指数函数的大气散射改正经验模型，初步验证表明该模型能有效提高受大气散射影响的数据精度和数据的可利用率，能为激光测高数据业务化提出一种新的解决方案。
笔者选择青海湖区域的线性体制ICESat/GLAS激光测高数据，初步构建了指数函数大气散射误差改正模型并拟合出相关参数，如公式（11）所示。
$$ P={\rm{ - 0}}{\rm{.275\;8}} + {\rm{0}}{\rm{.231\;1}}{{\rm{e}}^{{\rm{1}}{{.615\;3\;{OD}}}}} $$ (16) 式中：OD为光学厚度；P代表测距改正值。
利用兴凯湖以及德国北威州区域数据进行交叉验证，三个区域的云光学厚度与测高偏差的斯皮尔曼相关系数依次是0.925、0.916、0.927，可以看出两个变量之间高度相关。如图5所示，随着光学厚度OD增加，拟合模型的精度逐渐降低。初步统计表明，对于水体区域，该拟合模型可将测高偏差修正到5 cm以内。
图 5 误差改正模型及实验选区内光学厚度与测距偏差分布散点
Figure 5. Error correction model and dispersion points of optical depth and ranging deviation distribution in experimental area
经过大气散射校正之后，如果与真实值之间的误差在10 cm (水体区域)/15 cm （陆地区域）以内，将其视为修正后可利用激光点，统计结果如表2所示，实验结果表明，对于线性体制ICESat/GLAS激光测高数据，受大气散射影响的点约占总数的20%，其中OD≤2占比约50%，提出的指数函数模型对OD≤2引起的测高偏差有较好的校正作用，数据可利用率达到90%以上，总体可利用率提高了约10%。
表 2 散射改正后激光点可用率提高情况统计表
Table 2. Statistical table of laser points availability improvement after scattering correction
Experimental area Total points Atmospheric scattering OD≤2 Availability after correction （OD≤2） Availability improvement （OD≤2） Nordrhein-Westfalen 21 440 19.32% 47.29% 90.33% 8.25% Qinghai Lake 11 240 19.29% 51.71% 96.54% 9.63% Xingkai Lake 5 921 20.26% 55.17% 93.82% 10.49% -
大气散射引起的测距偏差主要是取决于以下几个因素：云高度、云光学厚度、云粒子大小、粒子形状和接收场。无论是蒙特卡洛模拟还是基于指数函数的大气散射改正方法，均需要与激光脉冲同步获取云高、云光学厚度、气溶胶等参数。在这种情况下，仅仅基于单通道、单角度激光脉冲数据反演精度有限，因此，有必要通过辅助设备同步进行大气探测，提高反演大气参数精度，进而提高大气散射误差改正精度。除此之外，受限于没有同步观测数据，反演过程中很多关键参数无法精确确定，比如激光雷达比。比如，采用蒙特卡洛模拟改正法的GLAS将激光雷达比S1作为定值（S1=17.8 sr），尽管这一值是在地基、Mie散射雷达数据都广泛使用的值，但是激光雷达比是与激光雷达探测波长、云（或气溶胶）类型、有效粒子半径等因素紧密相关的参数，仅采用固定值对不同类型云、气溶胶反演光学厚度具有较大的误差，进而导致大气散射改正不够精确。
而对于光子计数体制的ATLAS虽然基于第三方数据实现了云检测，但是该算法仅仅判断数据是否受到云影响，而没有实现云垂直分布判断。在没有与数据同步的第三方大气垂直分布数据源的情况下，ATLAS大气散射误差校正研究陷入了瓶颈。因此，同步获取激光测高卫星工作时段的大气相关的有效粒子半径、云相态、垂直分布等参数，进而选择更合理的激光雷达比、精确反演OD值，提高大气散射改正精度显得非常重要。
-
云、气溶胶等粒子是变化频繁的大气成分，变化特征观测主要围绕着粒子大小、形态、吸收和散射特性等因素进行，偏振测量是现有比较可靠的探测技术，李正强等人[43]对星载大气探测设备现状与发展进行了梳理，如表3所示。结合前述表1中CALIPSO和CATS的激光偏振观测通道，对于激光测高卫星而言，配置偏振通道既可以提高大气散射改正精度，同时能拓展激光数据在大气方面的应用。
表 3 星载大气探测载荷主要技术参数
Table 3. Main technical parameters of space-borne atmospheric observation payload
Payload Satellite Time Organization/Country Orbit height/km Technical parameters (Polarization is referred to as P) POLDER-1 ADEOS I 1996 CNES/France 797 Wavelength/nm：443 （P）, 490, 565, 670 （P）,763, 765,
865 （P）, 910APS Glory Mission 2011 NASA/USA 705 Wavelength/nm：410(P), 443(P), 555(P), 670(P), 865(P), 910(P), 1370(P), 1610(P), 2200(P) CALIOP CALIPSO 2006 NASA/USA 700 Wavelength/nm：532(P), 1064 MAI TG-2 2016 China ~400 Wavelength/nm：565(P), 670(P), 763, 765, 865(P), 910 DPC GF-5 2018 China 705 Wavelength/nm：443, 490(P), 565, 670(P),763, 765,
865(P), 910HARP CubeSat 2018 UMBC/USA ~400 Wavelength/nm：440(P), 550(P), 670(P), 870(P) HARP2 PACE 2022 UMBC/USA 675 Wavelength/nm：440(P), 550(P), 670(P), 870(P) POSP HJ-2 2019 China 644 Wavelength/nm：410(P), 443(P), 555(P), 670(P), 865(P), 910(P), 1380(P), 1610(P), 2250(P) PCF GF-5(02) 2020 China 705 Wavelength of DPC/nm：443, 490(P), 565, 670(P), 763, 76, 865(P), 910
Wavelength of POSP /nm：380(P), 410(P), 443(P), 490(P), 670(P), 865(P), 1 380(P), 1 610(P), 2 250(P)MAIA OTB-2 2022 NASA/USA Lower orbit Wavelength/nm：365, 391, 415, 444(P), 550, 646(P), 750, 763, 866, 943, 1 044(P), 1 610, 1 886, 2 126 -
近几年来，对地观测卫星激光测高技术飞速发展，由线性计数体制转变为光子体制，由单波束测量发展为多波束探测，由低重频发展为高重频观测。对于对地观测卫星激光测高而言，大气散射引起的测高误差是不可忽视的，分米级甚至厘米级的精度需求对高精度的大气散射误差改正提出了更高要求。梳理激光测高卫星大气散射延迟改正技术的现状及面临的难点问题，基本可以得出如下结论：
（1）目前大气散射改正算法还处于研究阶段，并没有得到业务化应用，无论是GLAS、GEDI还是ALTAS基本采用对大气散射影响的激光点进行质量控制来确保测量精度，导致点的可利用率受到一定影响。
（2）基于指数函数的大气散射改正算法，在有效反演光学厚度OD值的前提下，可实现规模化应用，是未来业务化应用的一个重要选择。
（3）为了实现激光测高大气散射改正的业务化，同步探测大气参数非常必要且很关键。
大气散射误差改正的前提是具有较高精度的大气层结构数据，GLAS以532 nm单通道获取了75 m分辨率的层高度信息，CALIPSO以532 nm偏振双通道获取了30、60 m分辨率的层高度信息，CATS以532 nm偏振双通道获取了60 m分辨率的层高度信息。而偏振通道能提高退偏比、激光雷达比等大气参数反演精度，这些大气参数与大气散射改正息息相关。因此，后续发展国产激光测高卫星时，建议在常规激光探测通道外，可考虑补充激光大气探测通道或偏振探测模式，如同步配备不超过75 m垂直分辨率的大气探测激光通道和偏振探测通道，提高大气参数反演精度，既提高大气散射改正精度、也进一步拓展激光雷达大气探测应用。
实现大范围三维地形、森林树高、极地冰盖、湖泊水位、浅海水下地形等信息的有效获取，是未来激光测高卫星发展的必然趋势。在以陆海激光测高卫星为代表的国产激光测高卫星观测体系中，大气散射改正是卫星立项论证与后续处理应用时必须考虑的问题，希望文中的有关研究能引起更多学者的关注，最终有效提高国产激光测高卫星的设计和应用水平。
Progress and prospect of atmospheric scattering correction for laser altimetry satellite
-
摘要: 激光测高卫星可以在大范围内获得亚米甚至厘米级的地表高程信息，但不可避免受云、气溶胶等粒子引起的散射影响，其中前向散射引起的激光测距和最终测高误差不容忽视。文中系统梳理激光测高卫星大气散射误差改正技术，介绍了国内外卫星激光测高系统参数、大气探测及散射改正算法，有别于已有的蒙特卡洛模拟改正方法，提出了一种基于指数函数模型的大气散射改正算法，选择青海湖等区域的ICESat/GLAS开展试验，结果表明，光学厚度小于2时能有效提升受大气散射影响的数据精度，同时提升数据可利用率约9%，该算法更易于实现业务化应用。最后针对大气参数同步探测的必要性，结合星载大气参数探测设备对后续国产激光测高卫星大气散射改正提出了若干建议。Abstract: Laser altimetry satellite can obtain the surface elevation information of sub-meter or even centimeter-level in a wide range, but it is inevitably affected by the scattering caused by particles such as clouds and aerosols. The laser ranging and final height measurement errors caused by forward scattering of the cloud or fog can not be ignored. In this paper, the atmospheric scattering error correction technology of laser altimetry satellite was systematically reviewed, and the satellite laser altimeter system parameters, atmospheric detection and scattering correction algorithm at home and abroad were introduced. Different from the Monte Carlo simulation correction method theory, an atmospheric scattering correction algorithm based on exponential function model was proposed. The data of GLAS (Geo-science Laser Altimeter System) on the ICESat (Ice, Cloud and land Elevation Satellite) in Qinghai Lake and other regions was selected for the experiment, and the experimental results show that the algorithm can effectively improve the accuracy of the altimetry data affected by atmospheric scattering when the optical thickness is less than 2, and the data availability rate can be improved by about 9%. The algorithm is easier to realize operation application. Finally, according to the necessity of synchronous detection of atmospheric parameters, some suggestions for atmospheric scattering correction of domestic laser altimetry satellites were put forward in combination with onboard atmospheric parameter detection equipment.
-
图 5 误差改正模型及实验选区内光学厚度与测距偏差分布散点
Figure 5. Error correction model and dispersion points of optical depth and ranging deviation distribution in experimental area
表 1 星载激光测高载荷主要技术指标
Table 1. Main technical parameters of space-borne laser altimeter
Payload GLAS CALIOP CATS ELA GEDI ATLAS SLA CASAL Year 2003 2006 2015 2016 2018 2018 2019 2021 Satellite ICESat-1 CALIPSO ISS ZY3-02 ISS ICESat-2 GF-7 TECSIMS Detection target Surface elevation Cloud/Aerosol Cloud/Aerosol Elevation reference point Surface elevation and Forest biomass Surface elevation Elevation reference point Height of forest and atmosphere Orbit height/km 600 700 405 505 405 500 500 506 Telescope diameter/m 1.0 1.0 0.6 0.21 0.8 0.8 0.6 — Wavelength/nm 532,1064 532,1064 355,532,1064 1064 1064 532 1064 532,1064 Laser repetition frequency/kHz 0.04 0.02 4 0.002 0.242 10 0.0030.006 VEG：0.04
ATM：0.02Single pulse energy/mJ 35-75 110 2-3 200 10 0.04-0.12 100-180 VEG：70
ATM：110Laser footprint diameter/m 70 — — 75 25 17.5 <30 25-30 Number of beams 1 1 1 1 4 6 2 VEG：5
ATM：1Data acquisition mode Linear mode Linear mode Single-photon mode Linear mode Linear mode Single-photon mode Linear mode Linear mode Resolution of atmospheric
data /m76.8 30 60 — 0.15 30 — 30 Atmospheric detection capability Yes Yes Yes No Yes Yes No Yes Optical depth inversion algorithm Threshold method +slope method Threshold method + slope method Threshold method + slope method None Threshold method + slope method Other data Laser footprint image Threshold method + slope method Correction method of atmospheric scattering error Monte Carlo simulation Monte Carlo simulation Monte Carlo simulation None Monte Carlo simulation Other data Pending research Pending research 
下载: 导出CSV
表 2 散射改正后激光点可用率提高情况统计表
Table 2. Statistical table of laser points availability improvement after scattering correction
Experimental area Total points Atmospheric scattering OD≤2 Availability after correction （OD≤2） Availability improvement （OD≤2） Nordrhein-Westfalen 21 440 19.32% 47.29% 90.33% 8.25% Qinghai Lake 11 240 19.29% 51.71% 96.54% 9.63% Xingkai Lake 5 921 20.26% 55.17% 93.82% 10.49%
下载: 导出CSV
表 3 星载大气探测载荷主要技术参数
Table 3. Main technical parameters of space-borne atmospheric observation payload
Payload Satellite Time Organization/Country Orbit height/km Technical parameters (Polarization is referred to as P) POLDER-1 ADEOS I 1996 CNES/France 797 Wavelength/nm：443 （P）, 490, 565, 670 （P）,763, 765,
865 （P）, 910APS Glory Mission 2011 NASA/USA 705 Wavelength/nm：410(P), 443(P), 555(P), 670(P), 865(P), 910(P), 1370(P), 1610(P), 2200(P) CALIOP CALIPSO 2006 NASA/USA 700 Wavelength/nm：532(P), 1064 MAI TG-2 2016 China ~400 Wavelength/nm：565(P), 670(P), 763, 765, 865(P), 910 DPC GF-5 2018 China 705 Wavelength/nm：443, 490(P), 565, 670(P),763, 765,
865(P), 910HARP CubeSat 2018 UMBC/USA ~400 Wavelength/nm：440(P), 550(P), 670(P), 870(P) HARP2 PACE 2022 UMBC/USA 675 Wavelength/nm：440(P), 550(P), 670(P), 870(P) POSP HJ-2 2019 China 644 Wavelength/nm：410(P), 443(P), 555(P), 670(P), 865(P), 910(P), 1380(P), 1610(P), 2250(P) PCF GF-5(02) 2020 China 705 Wavelength of DPC/nm：443, 490(P), 565, 670(P), 763, 76, 865(P), 910
Wavelength of POSP /nm：380(P), 410(P), 443(P), 490(P), 670(P), 865(P), 1 380(P), 1 610(P), 2 250(P)MAIA OTB-2 2022 NASA/USA Lower orbit Wavelength/nm：365, 391, 415, 444(P), 550, 646(P), 750, 763, 866, 943, 1 044(P), 1 610, 1 886, 2 126
下载: 导出CSV
-
[1] 李国元, 唐新明, 张重阳, 等. 多准则约束的ICESat/GLAS高程控制点筛选[J]. 遥感学报, 2017, 21(1): 96-104. doi: 10.11834/jrs.20175269 Li Guoyuan, Tang Xinming, Zhang Chongyang, et al. Multi-criteria constraint algorithm for selecting ICESat/GLAS data as elevation control points [J]. Journal of Remote Sensing, 2017, 21(1): 96-104. (in Chinese doi: 10.11834/jrs.20175269 [2] 唐新明, 李国元. 激光测高卫星的发展与展望[J]. 国际太空, 2017(11): 13-18. doi: 10.3969/j.issn.1009-2366.2017.11.004 Tang Xinming, Li Guoyuan. Development and prospect of laser altimetry satellite [J]. Space International, 2017(11): 13-18. (in Chinese doi: 10.3969/j.issn.1009-2366.2017.11.004 [3] Zwally H Jay, Schutz B, Abdalati W, et al. ICESat's laser measurements of polar ice, atmosphere, ocean, and land [J]. Journal of Geodynamics, 2002, 34(3): 405-445. [4] Wang Xianwei, Cheng Xiao, Gong Peng, et al. Earth science applications of ICESat/GLAS [J]. International Journal of Remote Sensing, 2011, 32(23): 8837-8864. doi: 10.1080/01431161.2010.547533 [5] 杨帆, 温家洪. ICESat与ICESat-2应用进展与展望[J]. 极地研究, 2011, 23(2): 138-148. Yang Fan, Wen Jiahong. Application progress and prospect of ICESat and ICESat-2 [J]. Chinese Journal of Polar Research, 2011, 23(2): 138-148. (in Chinese [6] Yu W Anthoy, Stephen M et al. Space Laser Transmitter Development for ICESat-2 Mission[C]//Proceedings of SPIE - The International Society for Optical Engineering, 2010, 7578(2): 757809. [7] 谢栋平, 李国元, 唐新明, 等. 美国GEDI天基激光测高系统及其应用[J]. 国际太空, 2018(12): 39-42. doi: 10.3969/j.issn.1009-2366.2018.12.008 Xie Dongping, Li Guoyuan, Tang Xinming, et al. GEDI space-based laser altimetry system and its application in the United States [J]. Space International, 2018(12): 39-42. (in Chinese doi: 10.3969/j.issn.1009-2366.2018.12.008 [8] 李国元, 唐新明. 资源三号02星激光测高精度分析与验证[J]. 测绘学报, 2017, 46(12): 1939-1949. doi: 10.11947/j.AGCS.2017.20170174 Li Guoyuan, Tang Xinming. Analysis and validation of ZY-3 02 satellite laser altimetry data [J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(12): 1939-1949. (in Chinese doi: 10.11947/j.AGCS.2017.20170174 [9] 李国元, 高小明, 陈继溢, 等. 资源三号02星激光测高数据质量分析[J]. 遥感学报, 2019, 23(6): 1159-1166. Li Guoyuan, Gao Xiaoming, Chen Jiyi, et al. Data quality analysis of ZY-3 02 satellite laser altimeter [J]. Journal of Remote Sensing, 2019, 23(6): 1159-1166. (in Chinese [10] Li Guoyuan, Tang Xinming, Gao Xiaoming, et al. ZY-3 block adjustment supported by GLAS laser altimetry data [J]. Photogrammetric Record, 2016, 31(153): 88-107. doi: 10.1111/phor.12138 [11] Li Guoyuan, Tang Xinming, Gao Xiaoming, et al. Integration of ZY3-02 satellite laser altimetry data and stereo images for high-accuracy mapping [J]. Photogrammetric Engineering & Remote Sensing, 2018, 84(9): 569-578. [12] Maria T Zuber, Gregory A Neumann, Frank G Lemoine, et al. Initial observations from the lunar orbiter laser altim-eter (LOLA) [J]. Geophysical Research Letters, 2010, 37(18): 1-6. [13] James W Head, Mikhail K, Harald H, et al. Oceans in the past history of Mars: Tests for their presence using Mars Orbiter Laser Altimeter (MOLA) data [J]. Geophysical Research Letters, 2013, 25(24): 4401-4404. [14] 李国元. 对地观测卫星激光测高数据处理方法与工程实践[D]. 武汉: 武汉大学, 2017. Li Guoyuan. Earth observing satellite laser altimeter data processing method and engineer practice[D]. Wuhan: Wuhan University, 2017. (in Chinese) [15] 张玉发, 孙晓泉. 激光主动探测影响因素分析[J]. 半导体光电, 2016, 37(6): 882-885. Zhang Yufa, Sun Xiaoquan. Analysis on the influence factors of laser active detection technology [J]. Semiconductor Optoelectronics, 2016, 37(6): 882-885. (in Chinese [16] 李松, 马跃, 周辉, 等. 星载激光测高系统大气折射延迟的影响分析[J]. 应用光学, 2014, 35(3): 520-524. Li Song, Ma Yue, Zhou Hui, et al. Analysis on atmospheric refraction delay for satellite laser altimeter system [J]. Journal of Applied Optics, 2014, 35(3): 520-524. (in Chinese [17] 李国元, 唐新明, 樊文锋, 等. 基于地面红外探测器的星载激光测高仪在轨几何定标[J]. 红外与激光工程, 2017, 46(11): 1117004. Li Guoyuan, Tang Xinming, Fan Wenfeng, et al. On-orbit geometric calibration of satellite laser altimeter using ground-based IR detectors [J]. Infrared and Laser Engineering, 2017, 46(11): 1117004. (in Chinese [18] 杨瑞科, 马春林, 韩香娥, 等. 激光在大气中传输衰减特性研究[J]. 红外与激光工程, 2007, 36(z2): 415-418. doi: 10.3969/j.issn.1007-2276.2007.z2.105 Yang Ruike, Ma Chunlin, Han Xiang’e, et al. Study of the attenuation characteristics of laser propagation in the atmosphere [J]. Infrared and Laser Engineering, 2007, 36(z2): 415-418. (in Chinese doi: 10.3969/j.issn.1007-2276.2007.z2.105 [19] Platt M R C. Remote sounding of high clouds. III: Monte Carlo calculations of multiple-scattered lidar returns [J]. Journal of the Atmospheric Sciences, 1981, 38(1): 156-167. doi: 10.1175/1520-0469(1981)038<0156:RSOHCI>2.0.CO;2 [20] Edwin W Eloranta. Practical model for the calculation of multiply scattered lidar returns [J]. Applied Optics, 1998, 37(12): 2464-2472. doi: 10.1364/AO.37.002464 [21] Bissonnette L R, Bruscaglioni P, Ismaelli A, et al. LIDAR multiple scattering from clouds [J]. Applied Physics, 1995, B60(4): 355-362. [22] Bruscaglioni P, Flesia C, Ismaelli A, et al. Multiple scattering and lidar returns [J]. Pure Applied Optics, 1999, 7(6): 1273. [23] 刘建斌, 吴健. 群体雾粒子对激光束散射光强的角谱分析[J]. 光学技术, 2005, 31(5): 659-661. doi: 10.3321/j.issn:1002-1582.2005.05.045 Liu Jianbin, Wu Jian. Analysis of angular spectrum of scattering intensity with particles for laser beam [J]. Optical Technique, 2005, 31(5): 659-661. (in Chinese doi: 10.3321/j.issn:1002-1582.2005.05.045 [24] Mahesh A, Spinhirne D James, Duda P David, et al. Atmospheric multiple scattering effects on GLAS altimetry. II. Analysis of expected errors in Antarctic altitude measurements [J]. Geoscience & Remote Sensing IEEE Transactions, 2001, 40(11): 2353-2362. [25] 么嘉棋,唐新明,李国元,等. 激光测高卫星ICESat-2云检测及其相关算法研究[J]. 激光与光电子学进展, 2020, 57(13): 240-248. Yao Jiaqi, Tang Xinming, Li Guoyuan, et al. Cloud detection of laser altimetry satellite icesat-2 and the related alogorithm [J]. Laser & Optoelectronics Progress, 2020, 57(13): 240-248. (in Chinese [26] Brenner A C, DiMarzio J P, Zwally H J. Precision and accuracy of satellite radar and laser altimeter data over the continental ice sheets [J]. IEEE Transactions on Geoscience & Remote Sensing, 2007, 45(2): 321-331. [27] 渠丽新, 李晓峰. 大气信道对激光通信信号散射脉冲展宽的研究[C]//中国宇航学会深空探测技术专业委员会学术会议, 2005. Qu Lixin, Li Xiaofeng. Study on the broadening of scattering pulse of laser communication signal in atmospheric channel[C]//Academic Meeting of Deep Space Exploration Technical Committee of China Aerospace Society, 2005. (in Chinese) [28] 宋雪平, 刘锋, 覃一凡. 云雾多次散射对激光散射信号生成的影响[C]//先进激光技术发展与应用研讨会, 2007. Song Xueping, Liu Feng, Qin Yifan. The influence of multiple scattering of cloud and fog on the generation of laser scattering signal[C]//Seminar on the development and application of advanced laser technology, 2007. (in Chinese) [29] 李颖颖, 孙东松, 沈法华, 等. 多次散射的激光雷达消光系数反演方法研究[C]//2007年先进激光技术发展与应用研讨会论文集, 2007. Li Yingying, Sun Dongsong, Shen Fahua, et al. Research on lidar extinction coefficient inversion method based on multiple scattering[C]//Seminar on The Development and Application of Advanced Laser Technology, 2007. (in Chinese) [30] 李颖颖, 孙东松, 王珍珠, 等. 云的多次散射对激光雷达测量结果影响的研究[J]. 激光技术, 2008, 32(6): 611-613, 638. Li Yingying, Sun Dongsong, Wang Zhenzhu, et al. Study of clouds multiple-scattering influence on lidar measurement [J]. Laser Technology, 2008, 32(6): 611-613, 638. (in Chinese [31] 陈舒杭, 李子乐, 陈梦竹, 等. 大气多次散射效应对星载激光测高仪测距偏差值的影响[J]. 红外与激光工程, 2012, 41(9): 2522-2526. doi: 10.3969/j.issn.1007-2276.2012.09.050 Chen Shuhang, Li Zile, Chen Mengzhu, et al. Influence of atmospheric multiple scattering effects on the range bias for satellite laser altimeter [J]. Infrared and Laser Engineering, 2012, 41(9): 2522-2526. (in Chinese doi: 10.3969/j.issn.1007-2276.2012.09.050 [32] 周辉, 李松, 王良训, 等. 单次大气散射效应对星载激光测高仪接收脉冲回波的影响[J]. 红外与激光工程, 2016, 45(1): 0106002. doi: 10.3788/irla201645.0106002 Zhou Hui, Li Song, Wang Liangxun, et al. Influence of single atmospheric scattering effect on received pulse waveform of satellite laser altimeter [J]. Infrared and Laser Engineering, 2016, 45(1): 0106002. (in Chinese doi: 10.3788/irla201645.0106002 [33] 赵欣, 张毅, 张黎明, 等. 激光测高仪高斯回波分解算法[J]. 红外与激光工程, 2012, 41(3): 101-106. Zhao Xin, Zhang Yi, Zhang Liming, et al. Algorithm of laser altimeter reflecting waveform based on gaussian decomposition [J]. Infrared and Laser Engineering, 2012, 41(3): 101-106. (in Chinese [34] 王婷, 徐军, 赵选科, 等. 基于Monte Carlo方法的雾红外传输仿真及分析[J]. 光子学报, 2018, 47(12): 122-131. Wang Ting, Xu Jun, Zhao Xuanke, et al. Simulation and analysis of infrared radiation transmission in fog based on Monte Carlo method [J]. Acta Photonica Sinica, 2018, 47(12): 122-131. (in Chinese [35] 么嘉棋, 高小明, 李国元, 等. ICESat/GLAS回波能量数据的云光学厚度反演[J]. 红外与激光工程, 2019, 48(S1): S117009. Yao Jiaqi, Gao Xiaoming, Li Guoyuan, et al. Cloud optical depth inversion based on ICESat/GLAS echo energy data [J]. Infrared and Laser Engineering, 2019, 48(S1): S117009. (in Chinese [36] Tang Xinming, Yao Jiaqi, Li Guoyuan. Influence of cloud scattering on satellite laser altimetry data and its correction [J]. Applied Optics, 2020, 59(13): 4064-4075. doi: 10.1364/AO.385808 [37] Duda P David, Spinhirne D James, Edwin W Eloranta. Atmospheric multiple scattering effects on GLAS altimetry-Part I: Calculations of single pulse bias [J]. IEEE Transactions on Geoscience and Remote Sensing, 1999, 39(1): 92-101. [38] Yang Yuekui, Marshak A, Palm P Stephen, et al. Assessment of cloud screening with apparent surface reflectance in support of the ICESat-2 Mission [J]. IEEE Transactions on Geoscience & Remote Sensing, 2013, 51(2): 1037-1045. [39] Michelle H, Story S, Yi D. Algorithm theoretical basis document for GEDI transmit and receive waveform processing for L1 and L2 products[EB/OL]. 2019.https://gedi.umd.edu/. [40] Palm P Stephen, Yang Yuekui, Herzfeld U, et al. ATLAS/ICESat-2 L3A calibrated backscatter profiles and atmospheric layer characteristics[EB/OL]. 2019.https://doi.org/10.5067/ATLAS/ATL09.001. [41] Palm P Stephen, Yang Yuekui, Herzfeld U, et al. ATLAS/ICESat-2 L2A normalized relative backscatter profiles[EB/OL].2019.https://doi.org/10.5067/ATLAS/ATL04.001. [42] Yang Yuekui, Marshak A, Palm P Sstephen, et al. Cloud impact on surface altimetry from a spaceborne 532-nm micropulse photon-counting lidar: system modeling for cloudy and clear atmospheres [J]. IEEE Transactions on Geoscience & Remote Sensing, 2011, 49(12): 4910-4919. [43] Dubovik O, Li Zhengqiang, Mishchenko I Mechael, et al. Polarimetric remote sensing of atmospheric aerosols: instruments, methodologies, results, and perspectives [J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 2018, 224: 474-511. -
点击查看大图
图(5) / 表(3)
计量
- 文章访问数: 568
- HTML全文浏览量: 183
- PDF下载量: 79
- 被引次数: 0
