Volume 50 Issue 6
Jun.  2021
Turn off MathJax
Article Contents

E Xiangying, Dai Guangyao, Wu Songhua. ICESat-2 ATL03 data preprocessing and correction method[J]. Infrared and Laser Engineering, 2021, 50(6): 20211032. doi: 10.3788/IRLA20211032
Citation: E Xiangying, Dai Guangyao, Wu Songhua. ICESat-2 ATL03 data preprocessing and correction method[J]. Infrared and Laser Engineering, 2021, 50(6): 20211032. doi: 10.3788/IRLA20211032

ICESat-2 ATL03 data preprocessing and correction method

doi: 10.3788/IRLA20211032
  • Received Date: 2021-05-09
  • Rev Recd Date: 2021-05-17
  • Publish Date: 2021-06-30
  • On September 15, 2018, NASA successfully launched the ICESat-2 satellite, which carries the Advanced Topographic Laser Altimeter System (ATLAS). ATLAS/ICESat-2 has been widely used in many research fields, such as polar ice-sheet elevation measurement, sea-ice thickness estimation, land elevation measurement, and vegetation measurement. Unlike GLAS/ICESat, ATLAS uses a more sensitive single-photon detector in the signal detection system. The system emits 10 000 pulses per second (GLAS/ICESat system emits 40 pulses per second), which enables high-density photon data with higher repetition frequency and smaller light spots, thus enabling high-precision global observation. Due to the influence of the detector photomultiplier tube (PMT) after-pulse used by ATLAS, small amplitude pulse interference often occurs after the output of echo signal pulse, which cannot accurately reflect the original signal waveform. Therefore, the measured signal must be corrected in the practical application process. Firstly, based on ATL03 product, the photon counting echo data were preprocessed to obtain the profiles of sea ice, ocean, land and desert. The results show that there are two small pulse echoes (about 2.3 m and 4.2 m) after the main signal echo, and there is a longer delay pulse echo between land and desert (about 10-45 m). It is confirmed that different reflective surfaces are affected by after-pulse of PMT. Secondly, the impulse response functions under different beams were calculated by using the night observation data of ICESat-2 over the Sahara Desert. The overall shape of the impulse response functions obtained from the six beams was similar. After the main signal echo, there were three pulse echoes of different amplitude (about 2.3 m, 4.2 m and 6.5 m). Finally, the impulse response functions under six beams were used to correct the same ocean profile. The results show that except ATLAS Beam2, all the other beams can effectively filter out the small amplitude pulse after the main pulse, and the real underwater echo signal can be picked up.

  • 加载中
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Article Metrics

Article views(1622) PDF downloads(218) Cited by()

Related
Proportional views

ICESat-2 ATL03 data preprocessing and correction method

doi: 10.3788/IRLA20211032
  • 1. College of Marine Technology, Faculty of Information Science and Engineering, Ocean University of China, Qingdao 266100, China
  • 2. Institute for Advanced Ocean Study, Ocean University of China, Qingdao 266100, China
  • 3. Ocean Remote Sensing Institute, Ocean University of China, Qingdao 266100, China
  • 4. Laboratory for Regional Oceanography and Numerical Modeling, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China

Abstract: 

On September 15, 2018, NASA successfully launched the ICESat-2 satellite, which carries the Advanced Topographic Laser Altimeter System (ATLAS). ATLAS/ICESat-2 has been widely used in many research fields, such as polar ice-sheet elevation measurement, sea-ice thickness estimation, land elevation measurement, and vegetation measurement. Unlike GLAS/ICESat, ATLAS uses a more sensitive single-photon detector in the signal detection system. The system emits 10 000 pulses per second (GLAS/ICESat system emits 40 pulses per second), which enables high-density photon data with higher repetition frequency and smaller light spots, thus enabling high-precision global observation. Due to the influence of the detector photomultiplier tube (PMT) after-pulse used by ATLAS, small amplitude pulse interference often occurs after the output of echo signal pulse, which cannot accurately reflect the original signal waveform. Therefore, the measured signal must be corrected in the practical application process. Firstly, based on ATL03 product, the photon counting echo data were preprocessed to obtain the profiles of sea ice, ocean, land and desert. The results show that there are two small pulse echoes (about 2.3 m and 4.2 m) after the main signal echo, and there is a longer delay pulse echo between land and desert (about 10-45 m). It is confirmed that different reflective surfaces are affected by after-pulse of PMT. Secondly, the impulse response functions under different beams were calculated by using the night observation data of ICESat-2 over the Sahara Desert. The overall shape of the impulse response functions obtained from the six beams was similar. After the main signal echo, there were three pulse echoes of different amplitude (about 2.3 m, 4.2 m and 6.5 m). Finally, the impulse response functions under six beams were used to correct the same ocean profile. The results show that except ATLAS Beam2, all the other beams can effectively filter out the small amplitude pulse after the main pulse, and the real underwater echo signal can be picked up.

  • 近些年来,星载激光雷达凭借其全天时、高空间分辨率以及覆盖范围广的特点,成为多个领域研究的有效工具。2003年1月,美国NASA发射激光测高仪系统Geoscience Laser Altimeter System(GLAS),搭载于以冰、云和陆地高程卫星(Ice, Cloud and land Elevation Satellite,ICESat)平台。该系统被用于测量冰被地形,同时也监测云层和大气的特性,提供云层高度和厚度,有助于提高短期天气及气候预报能力。2003年到2009年在轨期间,收集的测高数据促使科学家对格陵兰岛和南极冰盖的总体质量变化有了更为精确的估计[1]。为了延续ICESat的后续观测任务,NASA于2018年9月15日发射了地形激光测高仪系统(Advanced Topographic Laser Altimeter System,ATLAS),搭载于以冰、云和陆地高程卫星2号(The Ice, Cloud, and Land Elevation Satellite-2,ICESat-2)平台。ATLAS/ICESat-2将帮助科学家研究地球冰冻圈发生变化的原因和幅度;测量地球温带和热带区域的高程信息,并对全球森林植被信息进行评估[2]。不同于GLAS/ICESat全波形记录,ATLAS/ICESat-2采用更为灵敏的单光子探测器,以10 kHz的重复频率发射532 nm激光脉冲,通过获取高密度光子计数回波数据,实现高精度全球观测[3]

    ATLAS/ICESat-2已经广泛应用于多个研究领域[4]。Brunt等[5]测试并验证了ATL03产品和ATL06产品获取冰盖高程的精度,利用强波束和弱波束结合可以估算冰盖的坡度变化[6];Kwok[7]等研究表明ATL07产品测量海冰高程精度约2 cm;Dandabathula等[8]利用差分全球定位系统评价了ATL08产品的地面高程精度;Horvat等[9]基于ATLAS/ICESat-2数据提出了观测海冰区域海洋表面波的新技术;Lu等[10]利用ATL03产品观测极地海洋的浮游生物垂直分布情况,有效补充了被动水色遥感所缺失的夜间数据,对于海洋浮游生物循环研究具有重要意义。ATLAS/ICESat-2已经成为诸多领域研究的有效手段,展现出了强大的应用潜力,但在数据处理过程中发现,由于光电倍增管(Photomultiplier Tube,PMT) after-pulse的影响,会使得测量信号存在误差,所以在实际应用过程中,必须对原始测量信号进行校正。

    ATLAS/ICESat-2使用PMT作为探测器来进行光子计数,经不同表面反射的单个光子将被接收器接收并被标记时间和地理位置。然而由于PMT after-pulse影响,在主信号脉冲输出后常会出现小幅度脉冲回波,这会对大幅度信号情况下精确测量低电平信号产生干扰[11]。信号强度更弱的海洋次表层信号会受到来自海面回波噪声拖尾的污染,不能准确反演水体光学参数,所以研究系统的脉冲响应特性,对测量信号的校正及水下剖面信息的正确反演具有重要意义。

    文中基于ATLAS/ICESat-2的二级产品ATL03,对观测到的光子数据进行预处理,按照垂直方向上15 cm和沿轨道方向0.001 s的精度来统计,以获取不同反射表面(海冰、海洋、陆地和沙漠)的廓线信息,结果表明:不同反射表面的廓线均会受到PMT after-pulse影响。为了消除这种影响,必须了解ATLAS/ICESat-2的脉冲响应特性。文中使用撒哈拉沙漠地区夜间观测数据来研究不同波束下的脉冲响应函数,然后对同一海洋剖面信号进行退卷积处理,结果显示:除Beam2外,其余波束的退卷积结果差异较小,均可较好的校正PMT after-pulse影响,得到更加可靠的海洋次表层剖面信息。文章分为五部分内容,第一部分主要介绍文中的研究背景及研究意义;第二部分主要介绍ATLAS/ICESat-2的基本信息以及数据产品;第三部分主要介绍了信号预处理和校正的相关原理及处理方法,包括信号的预处理、脉冲响应函数的获取以及原始信号的校正;第四部分主要介绍了信号预处理及校正的结果,包括探究不同反射表面PMT after-pulse的影响、不同波束下脉冲响应函数的获取以及海洋水体信号的校正结果;第五部分是对全文的总结以及下步工作的展望。

  • 2018年9月15日,NASA在美国加州范登堡空军基地成功发射ICESat-2卫星,ICESat-2仅搭载了地形激光测高仪系统Advanced Topographic Laser Altimeter System (ATLAS),为科学家提供全球高度测量数据,可以精确估计冰川、海冰以及陆地等地形的变化特征。其主要科学目标包括:定量评估冰盖融化对海平面上升的影响;量化冰盖和冰川变化的特征,研究引发变化的驱动机制并提高冰盖预测模型的精确性;估算海冰厚度,研究海冰/海洋/大气之间的物质能量交换;测量植被高度,研究大尺度生物量现状及变化规律[2]

    ICESat-2轨道高度约500 km,重复周期为91天,每个周期包含1387个轨道。ATLAS以10 kHz的重复频率发射532 nm激光脉冲,可获取沿轨道方向间隔0.7 m,足印尺寸约17 m的光斑。ATLAS发射的每个激光脉冲被光学元件分裂产生6束单独的脉冲,排列成三对进行探测。每对包含一个强光束和一个弱光束(强弱能量比约1∶4),在跨轨方向上相距90 m,沿轨方向强弱波束间距约2.5 km,跨轨方向波束对间距约为3.3 km[2]表1给出了ATLAS/ICESat-2的关键参数。

    ParametersValuesParametersValues
    Orbit inclination and coverage92°;88°S~88°NFootprint size/m~17
    Altitude/km500Along-track separation/m~0.7
    Track repeat period/day91Telescope diameter/m0.8
    Nominal duration of mission/year3Number of beams6 organized in 3 pairs
    Laser wavelength/nm532Distance with a pair/m90
    Transmitted pulse width/ns1.5Distance between beam pairs/km~3.3
    Pulse repetition rate/kHz10Beam energy per pulse/μJ175±17,45±5

    Table 1.  Configuration parameters of ATLAS/ICESat-2[2]

    ATLAS/ICESat-2数据分为Level0、Level1、Level2和Level3共四级,包括ATL00-ATL21共21种产品,数据产品分类详细描述见https://icesat-2.gsfc.nasa.gov/science/data-products。ATL00是0级产品,提供原始测量数据;ATL01与ATL02是1级产品,在ATL00基础上进行了格式转换和仪器校正;ATL03提供了接收光子的经纬度以及椭球高度,并被作为生成其他产品的基础数据,如陆地冰(ATL06)[12],海冰(ATL07)[13],植被冠层高度(ATL08)[14]。公开数据产品均采用HDF-5格式通过美国冰雪数据中心提供,文中使用全球定位光子数据(ATL03 Version 3),数据可由https://nsidc.org/data/atl03免费下载。

    ATLAS/ICESat-2基于光电倍增管探测光子,并将其转换为电子信号,然后计时、存储并遥测以确定地面的高程信息。然而由于PMT存在非理想的after-pulse效应,当探测器接收到海表或者地表的信号时,测量信号总会大于实际信号[15]。PMT after-pulse包含有两种类型:一种是在主信号脉冲输出之后带有非常短的延迟(几纳秒到几十纳秒),另一种则表现为更长时间的延迟(几微秒)。延时较长的后脉冲是由PMT中残留气体的电离而产生的正离子引起,这些正离子返回到光电阴极并产生许多光电子,从而形成后脉冲[11]。研究ATLAS/ICESat-2系统的瞬态响应特性,了解在某一典型信号输入下,其输出量从初始状态到稳定状态的过程,对于准确获取真实信号波形具有重要意义。

  • 文中基于ATLAS/ICESat-2的二级产品ATL03(其中接收光子高度为WGS-84椭球高度),对观测到的光子计数回波数据进行筛选、下载、预处理(按照垂直方向上15 cm和沿轨道方向0.001 s的精度统计)并通过研究其不同波束脉冲响应特性来校正海洋剖面原始观测数据,以消除PMT after-pulse影响(处理流程见图1),进而获取到更为可靠的水下剖面信号。

    Figure 1.  Flow chart of ATLAS/ICESAT-2 data preprocessing and calibration

    首先选取并下载了四种不同反射表面区域的数据个例,对四种表面类型的原始光子数据分别进行预处理,目的是统计原始的光子分布得到廓线信息。按照垂向15 cm和沿轨方向0.001 s的精度统计,统计过程中需要根据峰值bin的位置对齐。处理结果表明:不同反射表面的观测数据均会受到PMT after-pulse影响。实验为分析探测器系统脉冲响应函数影响,并考虑降低太阳背景光等其他噪声影响,文中选取了2019年1月1日至2020年11月11日卫星过境撒哈拉沙漠地区的夜间观测数据,将6个波束的数据分别进行预处理。对于每个波束来说,经预处理会得到若干条20 m以内的廓线,再按照公式(1)[16]通过每条廓线的后向散射信号来计算脉冲响应曲线$ f\left({\textit{z}}\right) $,该波束下多条$ f\left({\textit{z}}\right) $沿垂直高度方向平均得到这一波束下的脉冲响应函数$ F\left({\textit{z}}\right) $[16]

    式中:$\; \beta \left({{\textit{z}}}_{j}\right) $表示第$ j $个bin后向散射值的大小;$ j $表示每个bin对应的位置;$ p $表示地面信号所在位置即峰值信号出现的位置。

    在无风或者风速很小的情况下,海表会发生镜面反射,海面回波之后的信号被认为是海表信号缓慢衰减导致,此时水下的信号可以被忽略。文中不考虑镜面反射情况,选取了海表回波光子数在8/bin~20/bin间的数据[17],对数据进行预处理,获取归一化回波信号剖面$ S\left({\textit{z}}\right) $,使用退卷积方法[18]来校正PMT的不良响应,以获取可靠的海洋次表层信号(见公式(3))。

    ICESat-2观测到的后向散射信号为$ S\left({\textit{z}}\right) $,是正确的后向散射信号$ {S}_{c}\left({\textit{z}}\right) $与探测器的脉冲响应函数$ F\left({\textit{z}}\right) $的卷积,即$ S\left({\textit{z}}\right)=F\left({\textit{z}}\right)*{S}_{c}\left({\textit{z}}\right) $,可被表示为[18]

    所以可利用不同波束的脉冲响应函数经过退卷积处理获取正确的后向散射信号:

  • 研究选取了海冰、海洋、陆地和沙漠等不同下垫面观测数据个例,对其进行预处理。图2使用了2018年10月14日观测到的海冰回波数据,其中图2(a)给出了该数据记录的250~270 s的光子分布,图2(b)给出了由250~270 s累加的归一化光子垂直剖面分布;图3使用了2018年10月16日观测到的海洋回波数据,其中图3(a)给出了该数据记录的40~100 s的光子分布,图3(b)给出了由40~100 s累加的归一化光子垂直剖面分布;图4使用了2018年10月16日观测到的陆地回波数据,其中图4(a)给出了该数据记录的220~240 s的光子分布,图4(b)给出了由220~240 s累加的归一化光子垂直剖面分布;图5使用了2019年7月24日观测到的沙漠回波数据,其中图5(a)给出了该数据记录的280~310 s的光子分布,图5(b)给出了由280~310 s累加的归一化光子垂直剖面分布。

    Figure 2.  (a) Photon distribution from sea ice surface (Height of each received photon, relative to the WGS-84 ellipsoid); (b) The photons vertical profile distribution, x-axis is normalized photon counts per bin and y-axis is altitude in meter. The altitude of peak surface return is set to 0 meter in (b)

    Figure 3.  (a) Photon distribution from ocean surface (Height of each received photon, relative to the WGS-84 ellipsoid); (b) The photons vertical profile distribution, x-axis is normalized photon counts per bin and y-axis is altitude in meter. The altitude of peak surface return is set to 0 meter in (b)

    Figure 4.  (a) Photon distribution from land surface (Height of each received photon, relative to the WGS-84 ellipsoid); (b) The photons vertical profile distribution, x-axis is normalized photon counts per bin and y-axis is altitude in meter. The altitude of peak surface return is set to 0 meter in (b)

    Figure 5.  (a) Photon distribution from desert surface (Height of each received photon, relative to the WGS-84 ellipsoid); (b) The photons vertical profile distribution, x-axis is normalized photon counts per bin and y-axis is altitude in meter. The altitude of peak surface return is set to 0 meter in (b)

    图2~图5显示了四种反射表面(海冰、海洋、陆地和沙漠)的光子分布和归一化光子垂直剖面分布,可以看出不同反射表面得到的垂直剖面在主回波后均会存在两个小的脉冲回波(2.3 m和4.2 m处),冰面和海面之后的两个小回波数量级均在10−3~10−2间,而陆地和沙漠地表主回波后的第二个小回波(4.2 m)在10−4~10−3间。此外从图4(b)图5(b)中可以观察到,对于陆地和沙漠来说,主脉冲之后在10~45 m间存在一个更长时间延迟的脉冲回波。因此可以判断,对于不同反射表面的观测数据来说,均会受到PMT after-pulse影响。

  • 由于ATLAS/ICESat-2 PMT after-pulse影响,使得海洋次表层信号会受到来自海面回波噪声拖尾的污染。因此,要想获取更加可靠的海洋次表层信号,必须清楚系统的脉冲响应特性。由于532 nm的激光很难穿透地表,理想情况下激光雷达的后向散射信号在地面强回波之后应该迅速衰减为零,但从图4(b)图5(b)中观察到在主脉冲之后约2.3 m、4.2 m以及10~45 m间存在有三个不同幅度的脉冲回波,可以认为地面之后的后向散射信号是由接收系统本身after-pulse响应造成,所以可使用地面回波信号对水体信号进行校正。文中选取了撒哈拉沙漠地区夜间观测数据来研究不同波束下系统的脉冲响应函数。

    ATLAS/ICESat-2通过发射3对波束进行全球范围观测,每对波束中强弱能量比约为4∶1 (强波束:Beam1、Beam3以及Beam5;弱波束:Beam2、Beam4以及Beam6)。文中选取了2019年1月1日至2020年11月11日撒哈拉沙漠地区夜间观测数据,筛选出高置信度的地面回波光子数据。图6(a)~(f)给出了6个波束的脉冲响应函数$ F\left({\textit{z}}\right) $(绿色曲线),分别为该波束下多条脉冲响应曲线$ f\left({\textit{z}}\right) $沿垂直高度方向平均得到。Beam1-Beam6的脉冲响应曲线对比如图7所示,结果显示不同波束得到的脉冲响应函数整体形态相似,均在主脉冲后出现三个不同幅度的脉冲回波,三个小峰出现位置一致(约2.3 m、4.2 m和6.5 m),前两个小峰峰值大小在10−4~10−3间(其中Beam2第二个小峰略小,与最大波束相差2×10−4),第三个小峰数量级在10−5。由此可知,不同波束下脉冲响应函数出现的三个小脉冲大小相近、位置一致,PMT after-pulse影响普遍存在。

    Figure 6.  Impulse response curves of six beams, in which the green curve is the impulse response function under the beam

    Figure 7.  Impulse response of six beams

  • 由于PMT after-pulse的影响,会使得水体的测量信号大于实际信号,不能准确反映真实的剖面信息,所以在反演工作之前必须对观测数据进行校正。文中选取了海表回波光子数在8/bin-20/bin间的数据,所选数据地理位置如图8所示,图9给出了所选数据的光子分布以及归一化回波信号垂直剖面。使用3.2节中由沙漠观测数据得到的Beam1~Beam6脉冲响应函数$ F\left({\textit{z}}\right) $,对2018年10月16日观测到的同一海洋剖面$ S\left({\textit{z}}\right) $进行退卷积处理,以获取正确的后向散射信号$ {S}_{c}\left({\textit{z}}\right) $

    Figure 8.  Location of marine data used in this paper

    Figure 9.  (a) Photon distribution from ocean surface (Height of each received photon, relative to the WGS-84 ellipsoid); (b) The photons vertical profile distribution, x-axis is normalized photon counts per bin and y-axis is altitude in meter. The altitude of peak surface return is set to 0 meter in (b)

    对于同一海洋剖面选取水下15 m以内数据,利用不同波束的脉冲响应函数进行退卷积,处理结果如图10所示,图11给出了6个波束的退卷积结果对比,可以看出除Beam2外,其余波束的退卷积结果差异较小,均可有效退去主脉冲之后的回波小峰,得到更加可靠的海洋次表层信号。Beam2下脉冲响应函数无法明显退掉水下信号的回波小峰,推测由于其地面主回波后的第二个回波小峰数值稍小(与其他波束最大相差2×10−4)所致。

    Figure 10.  (a)-(f):Calibration results were obtained by using impulse response models under six beams, where the black curve represents the observed ocean profile, the blue curve represents the impulse response model, and the red curve represents the correction results

    Figure 11.  The correction results of the observed ocean profile, where the black curve represents the observed ocean profile and the other color curve represents the deconvolution results of different beams

  • (1)针对ICESat-2 ATLAS星载激光雷达高度计的海冰、海水、陆地以及沙漠的观测数据进行预处理,结果显示不同反射表面的廓线均会在主回波后存在两个小的脉冲回波(约2.3 m和4.2 m处)。此外,陆地和沙漠在10~45 m间存在一个更长时间延迟的脉冲回波。因此可以判断,对于不同反射表面的观测数据来说,PMT after-pulse影响普遍存在。

    (2)选取了ICESat-2卫星过境撒哈拉沙漠的夜间观测数据,经过预处理和模型计算,获取到6个波束下的脉冲响应函数。对比发现6条曲线整体形态相似,均在主脉冲后出现三个不同幅度的脉冲回波,且出现位置一致(约2.3 m、4.2 m和6.5 m处)、大小相近。系统脉冲响应函数的精确提取,对提高海洋光学参数反演精度具有重要意义。

    (3)利用6个波束下得到的脉冲响应函数分别对同一水体后向散射信号进行退卷积处理,结果显示除Beam2外(推测由于其地面主信号回波后的第二个回波小峰数值过小所致),其余5个波束均有效滤除了主脉冲之后的小幅度脉冲,校正了PMT after-pulse影响。校正后的水体回波信号衰减趋势较好,说明了利用沙漠地区回波数据得到的脉冲响应函数进行水体信号校正的有效性。

    (4)后续工作将进一步研究分析Beam2下脉冲响应函数在海洋水体信号校正中失效的原因,并将基于校正后的海洋水体信号反演垂直剖面的多项光学参数(漫射衰减系数$ {K}_{d} $、后向散射系数$ {b}_{b} $以及颗粒物后向散射系数$ {b}_{bp} $等),同时与生物地球化学剖面浮标(Biogeochemical Argo float, BGC-Argo)和中分辨率成像光谱仪(Moderate Resolution Imaging Spectrometer, MODIS)提供的光学参数进行交叉验证,分析ICESat-2 ATLAS反演海洋浮游植物生物量的垂向分布特性的可行性。

Reference (18)
WeChat followshare

Top

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return