留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

ICESat-2/ATLAS数据反演林下地形精度验证

黄佳鹏 邢艳秋 秦磊 夏婷婷

黄佳鹏, 邢艳秋, 秦磊, 夏婷婷. ICESat-2/ATLAS数据反演林下地形精度验证[J]. 红外与激光工程, 2020, 49(11): 20200237. doi: 10.3788/IRLA20200237
引用本文: 黄佳鹏, 邢艳秋, 秦磊, 夏婷婷. ICESat-2/ATLAS数据反演林下地形精度验证[J]. 红外与激光工程, 2020, 49(11): 20200237. doi: 10.3788/IRLA20200237
Huang Jiapeng, Xing Yanqiu, Qin Lei, Xia Tingting. Accuracy verification of terrain under forest estimated from ICESat-2/ATLAS data[J]. Infrared and Laser Engineering, 2020, 49(11): 20200237. doi: 10.3788/IRLA20200237
Citation: Huang Jiapeng, Xing Yanqiu, Qin Lei, Xia Tingting. Accuracy verification of terrain under forest estimated from ICESat-2/ATLAS data[J]. Infrared and Laser Engineering, 2020, 49(11): 20200237. doi: 10.3788/IRLA20200237

ICESat-2/ATLAS数据反演林下地形精度验证

doi: 10.3788/IRLA20200237
基金项目: 国家重点研发计划(2017YFD060090402);中央高校基本科研业务费专项资金项目(2572019AB18);卫星测绘技术与应用国家测绘地理信息局重点实验室项目(KLSMTA-201706)
详细信息
    作者简介:

    黄佳鹏(1993-),男,博士生,主要从事星载激光雷达理论和应用方面的研究。Email:huangjp@nefu.edu.cn

    邢艳秋(1970-),女,教授,博士,主要从事林业遥感应用方面的研究。Email:yanqiuxing@nefu.edu.cn

  • 中图分类号: S75

Accuracy verification of terrain under forest estimated from ICESat-2/ATLAS data

  • 摘要: 针对星载激光雷达数据反演林下数字地面模型(Digital Terrain Model,DTM)存在困难的问题,研究了冰云陆地高程卫星-2 (Ice, Cloud, and land Elevation Satellite,ICESat-2)/先进地形激光高度计系统(Advanced Opographic Laser Altimeter System,ATLAS)的强弱光束数据反演林下地形的精度,并探究了冠层高度及植被覆盖率对于ICESat-2/ATLAS反演林下DTM精度的影响。研究结果表明:强波束反演林下DTM的精度为R2=1,RMSE=0.74 m,弱波束反演林下DTM的精度为R2=1,RMSE=0.76 m,强波束相对弱波束表现出更优的反演精度,但是,强光束与弱光束的光子云数据均可为反演林下DTM提供科学数据。从研究区植被的整体情况来看,随冠层高度及植被覆盖率的增加,不同激光类型数据均出现误差逐步增加的情况。
  • 图  1  研究区位置示意图

    Figure  1.  Location of study area

    图  2  ATLAS数据示意图[17]

    Figure  2.  Schematic of ATLAS data

    图  3  研究流程图

    Figure  3.  Research flow chart

    图  4  ATLAS数据与G-LiHT数据重叠示意图

    Figure  4.  Schematic diagram of overlap between ATLAS data and G-LiHT data

    图  5  ATLAS与G-LiHT数据反演林下地形示意图

    Figure  5.  Schematic diagram of ground elevation under forest by using ATLAS and G-LiHT data

    图  6  不同激光类型条件下ATLAS数据与G-LiHT数据散点图。(a)弱波束;(b)强波束

    Figure  6.  Scatter diagram of ATLAS data and G-LiHT data under different laser types. (a) Weak beam; (b) Strong beam

    图  7  不同冠层高度条件下ATLAS数据反演精度。(a)强波束;(b)弱波束

    Figure  7.  Estimation accuracy of ATLAS data under different canopy height conditions. (a) Strong beam; (b) Weak beam

    表  1  研究使用的ATLAS参数

    Table  1.   ATLAS parameters used in the paper

    ParametersDescription of parametersSource of parameters
    lat_phLatitude of each received photonATL03/gtx/heights/lat_ph
    lon_phLongitude of each received photonATL03/gtx/heights/lon_ph
    h_phHeight of each received photonATL03/gtx/heights/h_ph
    geoidGeoid height above WGS-84 reference ellipsoidATL03/gtx/geophys_corr/
    geoid
    segment_idA 7 digit number identifiying the along-track geolocation segment numberATL03/gtx/ geolocation/segment_id
    ph_index_begIndex (1-based) within the photon-rate data of the first photon within this segmentATL03/gtx/geolocation/
    ph_index_beg
    classed_pc_flagLand vegetation ATBD classification flag for each photon as either noise, ground, canopy,
    and top of canopy. 0=noise, 1=ground, 2=canopy, or 3=top of canopy
    ATL08/gtx/ signal_photons/
    classed_pc_flag
    classed_pc_indxIndex (1-based) of the ATL08 classified signal photon from the start of the ATL03 geolocation
    segment specified on the ATL08 product at the photon rate in the
    corresponding parameter, ph_segment_id
    ATL08/gtx/signal_photons/
    classed_pc_indx
    ph_segment_idSegment ID of photons tracing back to specific 20 m segment_id on ATL03ATL08/gtx/signal_photons/
    ph_segment_id
    下载: 导出CSV

    表  2  不同激光类型条件下ATLAS数据反演精度统计

    Table  2.   Estimation accuracy statistics of ATLAS data under different laser types

    Different laser intensities typesRMSE/mMAE/mR2ME/m
    Weak beam0.760.541.000.33
    Strong beam0.740.511.000.27
    Mean of different laser intensities types0.750.531.000.30
    下载: 导出CSV

    表  3  不同冠层高度条件下强、弱波束数据反演精度统计

    Table  3.   Estimation accuracy statistics of strong beam and weak beam data under different canopy heights

    Canopy heights/mRMSE/mMAE/mR2ME/mNpearsonspearmankendallVegetation coverage
    0-50.49/0.900.41/0.571.00/0.990.33/0.3221.39%/20.48%1.00/1.001.00/0.990.95/0.9327%(28.5%)/25%(35.7%)
    6-100.47/0.500.41/0.381.00/1.000.29/0.2522.37%/24.34%1.00/1.001.00/0.990.97/0.9427%(33.0%)/66%(36.4%)
    11-150.52/0.820.42/0.681.00/0.990.33/0.5627.20%/33.01%1.00/1.001.00/0.990.96/0.9436%(30.4%)/66%(61.8%)
    16-200.66/0.970.47/0.761.00/0.990.35/0.5013.28%/11.57%1.00/1.000.99/0.940.93/0.8436%(50.0%)/66%(51.0%)
    21-1.25/0.660.99/0.580.98/0.980.77/0.4115.76%/10.60%0.99/1.000.84/0.680.67/0.5236%(42.0%)/40%(30.9%)
    注:表格中数据分别为强波束数据/弱波束数据的精度统计数据
    下载: 导出CSV
  • [1] Long Tengfei, Zhang Zhaoming, He Guojin, et al. 30 m resolution global annual burned area mapping based on Landsat images and Google earth engine [J]. Remote Sensing, 2019, 11: 489. doi:  10.3390/rs11050489
    [2] Puliti S, Hauglin M, Breidenbach J, et al. Modelling above-ground biomass stock over Norway using national forest inventory data with ArcticDEM and Sentinel-2 data [J]. Remote Sensing of Environment, 2020, 236: 111501. doi:  10.1016/j.rse.2019.111501
    [3] Sarma C P, Dey A, Krishna A M. Influence of digital elevation models on the simulation of rainfall-induced landslides in the hillslopes of Guwahati, India [J]. Engineering Geology, 2020, 268: 105523. doi:  10.1016/j.enggeo.2020.105523
    [4] Daniela A, Ranjith G, Mikko K, et al. A method for vertical adjustment of digital aerial photogrammetry data by using a high-quality digital terrain model [J]. International Journal of Applied Earth Observation and Geoinformation, 2020, 84: 101954. doi:  10.1016/j.jag.2019.101954
    [5] Ai H, Kentaro T, Ram A, et al. Synthesis of L-band SAR and forest heights derived from TanDEM-X DEM and 3 digital terrain models for biomass mapping [J]. Remote Sensing, 2020, 12: 349. doi:  10.3390/rs12030349
    [6] 马跃, 阳凡林, 王明伟, 等. 利用GLAS激光测高仪计算格陵兰冰盖高程变化[J]. 红外与激光工程, 2015, 44(12): 65-69.

    Ma Yue, Yang Fanlin, Wang Mingwei, et al. Calculation of elevation changing of Greenland’s ice sheet using GLAS laser altimeter [J]. Infrared and Laser Engineering, 2015, 44(12): 65-69. (in Chinese
    [7] 李国元, 唐新明, 樊文锋, 等. 基于地面红外探测器的星载激光测高仪在轨几何定标[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
    [8] Wang Xianwei, David M Holland, Gudmundsson G Hilmar, et al. Accurate coastal DEM generation by merging ASTER GDEM and ICESat/GLAS data over Mertz Glacier, Antarctica [J]. Remote Sensing of Environment, 2018, 206: 218-230. doi:  10.1016/j.rse.2017.12.041
    [9] Nie Sheng, Wang Cheng, Dong Pinliang, et al. A novel model for terrain slope estimation using ICESat/GLAS waveform data [J]. IEEE Transactions on Geoence and Remote Sensing, 2018, 50, 1: 217-227.
    [10] Lamsters K, Karuss J, Krievans M, et al. High-resolution orthophoto map and digital surface models of the largest Argentine Islands (the Antarctic) from unmanned aerial vehicle photogrammetry [J]. Journal of Maps, 2020, 16(2): 335-347. doi:  10.1080/17445647.2020.1748130
    [11] Li Hui, Zhao Jiayang. Evaluation of the newly released worldwide AW3D30 DEM over typical landforms of China using two global DEMs and ICESat/GLAS data [J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2018, 11: 4430-4440. doi:  10.1109/JSTARS.2018.2874361
    [12] Yang Xuebo, Wang Cheng, Nie Sheng, et al. application and validation of a model for terrain slope estimation using space-borne LiDAR waveform data [J]. Remote Sensing, 2018, 10: 1691. doi:  10.3390/rs10111691
    [13] Irfan A, Jadunandan D, Saleem U, et al. A novel approach to estimate canopy height using ICESat/GLAS data: A case study in the New Forest National Park, UK [J]. International Journal of Applied Earth Observation and Geoinformation, 2013, 23: 109-118. doi:  10.1016/j.jag.2012.12.009
    [14] Yue Linwei, Shen Huanfeng, Zhang Liangpei, et al. High-quality seamless DEM generation blending SRTM-1, ASTER GDEMv2 and ICESat/GLAS observations [J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2017, 123: 20-34. doi:  10.1016/j.isprsjprs.2016.11.002
    [15] Neuenschwander A L, Magruder L A. Canopy and terrain height retrievals with ICESat-2: a first look [J]. Remote Sensing, 2019, 11: 1721. doi:  10.3390/rs11141721
    [16] Wang Cheng, Zhu Xiaoxiao, Nie Sheng, el al. Ground elevation accuracy verification of ICESat-2 data: a case study in Alaska, USA [J]. Optics Express, 2019, 27: 38168-38179. doi:  10.1364/OE.27.038168
    [17] Huang Jiapeng, Xing Yanqiu, You Haotian, et al. Particle swarm optimization-based noise filtering algorithm for photon cloud data in forest area [J]. Remote Sensing, 2019, 11(8): 980. doi:  10.3390/rs11080980
    [18] 黄佳鹏, 邢艳秋, 秦磊, 等. 弱光束条件下森林区域光子云去噪算法精度研究[J]. 农业机械学报, 2020, 51(4): 164-172. doi:  10.6041/j.issn.1000-1298.2020.04.019

    Huang Jiapeng, Xing Yanqiu, Qin Lei, et al. Accuracy of photon cloud noise filtering algorithm in forest area under weak beam conditions [J]. Transactions of the Chinese Society for Agricultural Machinery, 2020, 51(4): 164-172. (in Chinese doi:  10.6041/j.issn.1000-1298.2020.04.019
    [19] Chen Jun, Cao Xin, Peng Shu, et al. Analysis and applications of GlobeLand30: A review [J]. ISPRS International Journal of Geo-Information, 2017, 6(8): 230. doi:  10.3390/ijgi6080230
    [20] Hansen M, Song X P. Vegetation Continuous Fields (VCF) yearly global 0.05 deg [DB/OL]. [2020-08-10]https://doi.org/10.5067/MEaSUREs/VCF/VCF5KYR.001.
    [21] Neuenschwander Amy, Pitts Katherine, Jelley Benjamin, et al. ICE, CLOUD, and Land Elevation Satellite (ICESat-2) algorithm theoretical basis document (ATBD) for land-vegetation along-track products [EB/OL]. [2020-01-15] https://nsidc.org/data/atl08.
  • [1] 李亚茹, 周亮, 刘朝晖, 折文集.  天基数字成像链路数理模型构建方法 . 红外与激光工程, 2023, 52(12): 20230351-1-20230351-12. doi: 10.3788/IRLA20230351
    [2] 刘佳琪, 程用志, 陈浮, 罗辉, 李享成.  基于几何相位超表面的高效独立双频点圆偏振太赫兹波束调控 . 红外与激光工程, 2023, 52(2): 20220377-1-20220377-11. doi: 10.3788/IRLA20220377
    [3] 周智标, 周辉, 马跃, 宋越, 李松.  ICESat-2激光雷达海面信号提取和海浪要素计算 . 红外与激光工程, 2023, 52(2): 20220366-1-20220366-9. doi: 10.3788/IRLA20220366
    [4] 隋钧铖, 任登凤, 韩玉阁.  地面目标红外建模缺失参数反演与模型验证方法 . 红外与激光工程, 2022, 51(10): 20220033-1-20220033-12. doi: 10.3788/IRLA20220033
    [5] 邢晖, 张岩岫, 雷萍, 马娜, 王娟锋, 康大勇.  对光学侦察卫星的地面探测模型与试验研究 . 红外与激光工程, 2022, 51(9): 20210978-1-20210978-7. doi: 10.3788/IRLA20210978
    [6] 俄相颖, 戴光耀, 吴松华.  ICESat-2 ATL03数据预处理及校正方法 . 红外与激光工程, 2021, 50(6): 20211032-1-20211032-11. doi: 10.3788/IRLA20211032
    [7] 汪自军, 张扬, 刘东, 王晓波, 袁金如, 潘超, 赵一鸣, 韩晓爽, 周雨迪, 刘群, 王成.  新型多波束陆-海激光雷达探测卫星技术发展研究 . 红外与激光工程, 2021, 50(7): 20211041-1-20211041-11. doi: 10.3788/IRLA20211041
    [8] 邵光灏, 刘昂, 翟计全, 张国强.  光波束形成中色散延时的非线性修正 . 红外与激光工程, 2021, 50(7): 20210235-1-20210235-6. doi: 10.3788/IRLA20210235
    [9] 李泓成, 周辉, 李松, 马跃, 王玥.  面向ICESat-2单光子激光雷达的陆地目标测距误差分析 . 红外与激光工程, 2020, 49(11): 20200247-1-20200247-8. doi: 10.3788/IRLA20200247
    [10] 李彬彬, 谢欢, 童小华, 叶丹, 孙凯鹏, 李铭.  基于随机森林的ICESat-2卫星数据地表覆盖分类 . 红外与激光工程, 2020, 49(11): 20200292-1-20200292-7. doi: 10.3788/IRLA20200292
    [11] 张岩, 李春, 卞博锐, 张文, 蒋玲.  新型太赫兹波束分离器的设计 . 红外与激光工程, 2020, 49(5): 20190290-20190290-7. doi: 10.3788/IRLA20190290
    [12] 朱笑笑, 王成, 习晓环, 聂胜, 杨学博, 黎东.  ICESat-2星载光子计数激光雷达数据处理与应用研究进展 . 红外与激光工程, 2020, 49(11): 20200259-1-20200259-10. doi: 10.3788/IRLA20200259
    [13] 杨书娟, 张珂殊, 邵永社.  地面激光扫描仪的灰度优先距离检校模型 . 红外与激光工程, 2019, 48(1): 106007-0106007(6). doi: 10.3788/IRLA201948.0106007
    [14] 谭敏, 谢晨波, 王邦新, 吴德成, 马晖, 刘东, 王英俭.  北京2014年冬季边界层高度与颗粒物浓度的相关性研究 . 红外与激光工程, 2018, 47(7): 717007-0717007(8). doi: 10.3788/IRLA201847.0717007
    [15] 侯重远, 李恒年, 杨元, 路毅.  地面强激光操控空间碎片避碰的控制策略 . 红外与激光工程, 2017, 46(3): 329003-0329003(6). doi: 10.3788/IRLA201746.0329003
    [16] 李铭, 候佳, 周成林, 舒嵘.  多波束激光雷达的高精度收发配准方法的实验验证 . 红外与激光工程, 2017, 46(7): 730001-0730001(9). doi: 10.3788/IRLA201746.0730001
    [17] 廖天河, 刘伟, 高穹.  不同形状激光波束在大气中传输的湍流效应 . 红外与激光工程, 2015, 44(S1): 41-45.
    [18] 卢发铭, 张盛, 夏元钦, 陈德应.  双色飞秒强激光作用下的CO2分子高次谐波 . 红外与激光工程, 2014, 43(1): 77-80.
    [19] 刘永川, 宋恩民, 金人超, 许向阳, 刘宏, 马光志.  一种基于强散射介质的断层重建模型 . 红外与激光工程, 2014, 43(9): 3094-3098.
    [20] 屈檀, 吴振森, 牟媛, 李正军.  在轴零阶贝塞尔波束对手征介质球的散射特性分析 . 红外与激光工程, 2014, 43(9): 2867-2872.
  • 加载中
图(7) / 表(3)
计量
  • 文章访问数:  772
  • HTML全文浏览量:  590
  • PDF下载量:  122
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-06-09
  • 修回日期:  2020-07-10
  • 刊出日期:  2020-11-25

ICESat-2/ATLAS数据反演林下地形精度验证

doi: 10.3788/IRLA20200237
    作者简介:

    黄佳鹏(1993-),男,博士生,主要从事星载激光雷达理论和应用方面的研究。Email:huangjp@nefu.edu.cn

    邢艳秋(1970-),女,教授,博士,主要从事林业遥感应用方面的研究。Email:yanqiuxing@nefu.edu.cn

基金项目:  国家重点研发计划(2017YFD060090402);中央高校基本科研业务费专项资金项目(2572019AB18);卫星测绘技术与应用国家测绘地理信息局重点实验室项目(KLSMTA-201706)
  • 中图分类号: S75

摘要: 针对星载激光雷达数据反演林下数字地面模型(Digital Terrain Model,DTM)存在困难的问题,研究了冰云陆地高程卫星-2 (Ice, Cloud, and land Elevation Satellite,ICESat-2)/先进地形激光高度计系统(Advanced Opographic Laser Altimeter System,ATLAS)的强弱光束数据反演林下地形的精度,并探究了冠层高度及植被覆盖率对于ICESat-2/ATLAS反演林下DTM精度的影响。研究结果表明:强波束反演林下DTM的精度为R2=1,RMSE=0.74 m,弱波束反演林下DTM的精度为R2=1,RMSE=0.76 m,强波束相对弱波束表现出更优的反演精度,但是,强光束与弱光束的光子云数据均可为反演林下DTM提供科学数据。从研究区植被的整体情况来看,随冠层高度及植被覆盖率的增加,不同激光类型数据均出现误差逐步增加的情况。

English Abstract

    • 截止2018年年底,全球森林总面积为38.15亿hm2,约占全球陆地总面积的25.60% (按全球陆地面积149亿hm2计算)[1],故林下数字地面模型(Digital Terrain Model,DTM)是构建全球DTM模型的重要部分[2]。为获取全球范围内的大尺度、多时相数据,通常选择卫星观测数据作为主要研究对象[3-4]。不同于光学卫星遥感技术及微波卫星遥感技术反演林下地形过程中易受云层及森林信号饱和的影响,卫星激光测高技术可以获取分米级甚至厘米级的林下DTM数据[5-7]

      目前为止,全球共有4个已经成功发射的以激光测高技术为主载荷的对地观测卫星,分别为:冰云陆地高程卫星(Ice, Cloud, and land Elevation Satellite,ICESat)/ 地学激光高度计系统(Geoscience Laser Altimeter System,GLAS),ICESat-2/先进地形激光高度计系统(Advanced Topographic Laser Altimeter System,ATLAS),高分七号(GF-7)和全球生态系统动力学研究卫星(Global Ecosystem Dynamics Investigation,GEDI)。其中,ICESat,GF-7,GEDI均以线性体制的波形方式记录激光测高数据,而ICESat-2是以光子计数方式记录激光测高数据。基于ICESat/GLAS数据已经分别完成不同类型DTM数据的反演研究。其中,Wang、Nie、Lamsters等将ICESat/GLAS数据应用于冰川的DTM数据反演,并得出ICESat/GLAS数据与实测值具有较高相关性[8-10]。Li等[11]讨论了ICESat/GLAS数据与全球先进陆地观测卫星(Advance Land Observing Satellite,ALOS)的ALOS world 3D–30 m(AW3D30)、航天飞机雷达地形任务(Shuttle Radar Topography Mission,SRTM)的SRTM1数据及Global DEM version 2(GDEM2)的相关性,通过试验证明,3种数据存在明显相关性,其中,AW3D30相关性最高,RMSE为4.81 m。Yang等[12]利用ICESat/GLAS波形数据量化了地形坡度,试验结果表明:ICESat/GLAS数据反演地面坡度的精度分别为R2=0.83,RMSE=3.6°。Irfan等[13]利用ICESat/GLAS数据完成森林冠层高度信息反演研究,通过试验结果证明ICESat/GLAS数据可解释研究区内的79%冠层高度数据,且ICESat/GLAS数据可用于反演林下地形。但是,由于森林区域波形数据的特性,使得如何利用星载波形数据正确获取林下DTM数据存在困难。故利用ICESat/GLAS数据反演的DTM研究更多集中于非林下的DTM探究,因此如何利用星载激光雷达数据反演林下DTM将是星载激光波形数据应用的主要难点[9, 14]。由于ICESat-2/ATLAS采用不同于ICESat/GLAS的光子计数形式记录激光数据,所以探究光子数据反演林下DTM数据精度对于星载激光雷达数据应用具有重要意义。Neuenschwander等[15]通过试验证明,芬兰林区机载激光雷达数据反演的林下DTM与对应的ICESat-2/ATLAS数据之间存在良好的相关性(R2=0.99,RMSE=0.85 m)。Wang等[16]通过试验证明在阿拉斯加林区,ICESat-2/ATLAS数据与机载激光雷达获取的DTM之间存在良好相关性(R2=0.99,RMSE=1.96 m)。与ICESat/GLAS采用的单波束激光发射方式不同,ICESat-2/ATLAS采用3对6波束(每对激光束由一强一弱激光束组成,强弱光束能量比为4∶1)的激光发射方式,而以上针对ICESat-2/ATLAS数据的研究均以其强光束类型数据作为主要研究对象,并未分析ICESat-2/ATLAS弱光束数据反演林下DTM的精度[17-18]

      为探究ATLAS/ICESat-2强、弱光束数据反演林下DTM的精度,研究拟以ATL03、ATL08数据作为研究对象,探究ICESat-2/ATLAS的强、弱光束数据反演林下DTM的精度,并讨论不同冠层高度及植被覆盖率对于ICESat-2/ATLAS数据反演林下DTM的精度影响。

    • 研究选择的研究区位于美国的南卡罗莱纳州Aiken郡,其经纬度范围为(33.56°N~33.78°N,81.69°W~81.73°W),如图1所示。Aiken郡位于南卡罗莱纳州西南部,属于亚热带季风湿润气候。研究区内植被类型包括:耕地(0.04%)、森林(88.52%)、灌木(0.51%)、湿地(6.97%)、人工地表(3.95%)[19]。研究区林下地形缓坡居多,海拔范围为91~164 m。植被覆盖率范围为25%~66%,大部分足印覆盖区的植被覆盖率为36% (29.5%),66% (26.2%)[20]。冠层高度范围为0.1~34.8 m。

      图  1  研究区位置示意图

      Figure 1.  Location of study area

    • ICESat-2采用的激光发射频率为10 kHz,结合ICESat-2的卫星轨道高度及卫星运动速度可知,ICESat-2的足印直径约为17 m,沿轨足印间距0.7 m。每对波束垂轨间距约为3 km,每对波束内强弱光束垂轨方向间距约为90 m,沿轨方向间距约为2.5 km。ATLAS数据足印分布形式如图2所示[17]

      图  2  ATLAS数据示意图[17]

      Figure 2.  Schematic of ATLAS data

      根据美国国家冰雪数据中心(National Snow and Ice Data Center,NSIDC)公布的ICESat-2数据产品形式,ICESat-2有3种级别共21种产品,分别为ATL01~ATL21,文件均采用hdf5文件格式存储数据。其中与林下DTM相关的数据产品包括:ATL03、ATL08。其中ATL03为全球定位光子数据产品(Global Geolocated Photon Data),ATL03产品包括每个光子事件的时间、纬度、经度和高程信息。在ATL03文件中,每个光子采用2种索引方法,一种是按照光子的传输时间顺序进行编号,称为光子数据传输时间序号,光子数据的经纬度、高程等信息按照该序号被记录在ATL03/gtx/heights组内。另一种索引方式为按沿轨距离每20 m划分为一个区段,并按区段进行数据处理生成ATL08级数据。对各组区段依次进行编号,称为区段号,参数名称为segment_id,每组区段的起始光子序号参数名称为ph_index_beg。

      ATL08为陆地植被数据产品(Land and Vegetation Height),ATL08产品为利用差分、回归和高斯自适应最近邻(Differential, Regressive, and Gaussian Adaptive Nearest Neighbor,DRAGANN)对ATL03数据进行去噪,并利用NASA官方分类算法完成光子云的分类,将光子云划分为冠层顶光子、冠层光子、地面光子和噪声光子[15],光子分类参数名称为:classed_pc_flag,ATL08数据的光子分类参数为按ATL03文件的20 m区段进行储存,光子分类参数区段号的参数名称为ph_segment_id,区段内相对序号的参数名称为classed_pc_indx。研究使用的ATL03、ATL08数据均可在网站https://search.earthdata.nasa.gov下载。

    • 为确定研究区范围,利用机载激光雷达(Goddard's LiDAR, Hyperspectral Thermal Imager, G-LiHT)公布的机载轨迹标记语言 (Keyhole Markup Language, KML)文件联合ATL03文件的卫星运行轨迹,确定两数据重合位置,作为研究区。为验证ATLAS数据反演林下DTM的精度,研究选择G-LiHT数据发布的数字地面模型(DTM)、冠层高度模型(Canopy Height Model,CHM)产品作为验证数据。该类数据的分辨率为1 m,数据格式为TIFF。研究使用的G-LiHT相关数据来自于https://glihtdata.gsfc.nasa.gov

    • 为探究ATLAS强、弱光束数据反演林下DTM精度,研究选择对ATL03、ATL08数据进行关联,着重探究ATL08标定的ATL03中的林下地面光子与G-LiHT的DTM数据的关系。研究方法流程如图3所示,研究主要包括以下步骤:

      图  3  研究流程图

      Figure 3.  Research flow chart

      (1)根据ATLAS轨道范围结合G-LiHT的KML文件,选择ATLAS轨道与G-LiHT轨迹重合区域作为研究区。

      (2)提取ATL03数据产品中光子云数据的经度(lat_ph)、纬度(lon_ph)、高程(h_ph)和大地水准面改正信息(geoid)。由于光子云数据为基于WGS84坐标系统,需将光子云高程信息减去大地水准面改正信息获得带有大地水准面改正信息的光子高程数据。研究使用的ATL03文件与ATL08文件中的参数汇总如表1所示。

      表 1  研究使用的ATLAS参数

      Table 1.  ATLAS parameters used in the paper

      ParametersDescription of parametersSource of parameters
      lat_phLatitude of each received photonATL03/gtx/heights/lat_ph
      lon_phLongitude of each received photonATL03/gtx/heights/lon_ph
      h_phHeight of each received photonATL03/gtx/heights/h_ph
      geoidGeoid height above WGS-84 reference ellipsoidATL03/gtx/geophys_corr/
      geoid
      segment_idA 7 digit number identifiying the along-track geolocation segment numberATL03/gtx/ geolocation/segment_id
      ph_index_begIndex (1-based) within the photon-rate data of the first photon within this segmentATL03/gtx/geolocation/
      ph_index_beg
      classed_pc_flagLand vegetation ATBD classification flag for each photon as either noise, ground, canopy,
      and top of canopy. 0=noise, 1=ground, 2=canopy, or 3=top of canopy
      ATL08/gtx/ signal_photons/
      classed_pc_flag
      classed_pc_indxIndex (1-based) of the ATL08 classified signal photon from the start of the ATL03 geolocation
      segment specified on the ATL08 product at the photon rate in the
      corresponding parameter, ph_segment_id
      ATL08/gtx/signal_photons/
      classed_pc_indx
      ph_segment_idSegment ID of photons tracing back to specific 20 m segment_id on ATL03ATL08/gtx/signal_photons/
      ph_segment_id

      (3)通过ATL03数据产品中光子云数据的经纬度信息获取对应G-LiHT的DTM数据和CHM数据作为验证数据,如图4所示。图4(a)、4(b)分别为ATL03数据与对应研究区G-LiHT的DTM数据和CHM数据重合示意图。图中红色轨迹为gt2r轨迹,绿色轨迹为gt2l轨迹。图中ATLAS文件名为:ATL03_20181226053112_13530106_001_01.h5,G-LiHT的DTM、CHM文件名分别为:AMIGACarb_Augusta_FIA_Sep2011_ l47s557_DTM.jpg、AMIGACarb_ Augusta _FIA_Sep2011_l47s557_CHM.jpg。

      图  4  ATLAS数据与G-LiHT数据重叠示意图

      Figure 4.  Schematic diagram of overlap between ATLAS data and G-LiHT data

      (4)提取ATLAS数据产品中的光子分类参数(classed_pc_flag),关联参数(segment_id、ph_index_beg、classed_pc_indx、ph_segment_id)。关联方式为:遍历ATL03、ATL08的光子云数据,通过ATL03中的segment_id与ATL08中的ph_segment_id匹配两数据文件相同区段,获得该组的起始光子序号(ph_index_beg),光子的相对序号(classed_pc_indx)与所在组的起始光子序号相加,可得到其在ATL03文件中光子数据传输时间序号,即光子数据传输时间序号[21]。由此,可将ATL08中的光子分类参数与ATL03中的光子数据经纬度及高程进行关联,按照ATL08的分类参数确定ATL03光子云数据中的地面光子。

      (5) classed_pc_flag为ATL08产品官方提供的NASA官方的光子云数据的分类参数,classed_pc_flag的值可能为0、1、2、3,分别对应噪声光子,地面光子,冠层光子,冠层顶光子。研究选择利用classed_pc_flag=1即ATL08分类为地面的光子作为研究使用的林下地面光子。ATL03与ATL08数据反演林下地形结果如图5所示。图5(a)、5(b)分别为ATLAS的弱光束与强光束对应的反演林下地形示意图,其中红色点为ATL03的数据,蓝线曲线为带有ATL08分类参数的ATL03光子数据所生成的DTM,绿线曲线为G-LiHT数据所生成的DTM。

      图  5  ATLAS与G-LiHT数据反演林下地形示意图

      Figure 5.  Schematic diagram of ground elevation under forest by using ATLAS and G-LiHT data

      (6)比较分析ATL03产品中的地面光子与机载G-LiHT数据对应DTM数值的对应关系。若地面光子的高度与机载G-LiHT数据相应高度之间的绝对差大于20 m,则认为该光子被错误分类并剔除该光子。

      (7)为探究不同冠层高度及植被覆盖率对于ATLAS强、弱激光束反演林下DTM的影响程度,研究选择以G-LiHT数据的CHM数据作为冠层高度验证数据,并以5 m为间隔,将研究区足印覆盖的冠层高度分为:0~5 m,6~10 m,11~15 m,16~20 m,大于21 m,统计不同冠层区段内,两种数据反演林下DTM相关性评价数据。

    • 为评价ATLAS产品的地面光子高程与机载G-LiHT产品中的DTM数据的关系。研究选择统计了两种数据的均方根误差(RMSE),绝对平均误差(MAE),决定系数(R2),平均误差(ME)定量评价两种数据之间的关系。为探究不同冠层高度对于ATLAS数据反演林下DTM的影响程度,研究拟计算了ATLAS地面光子高程减去机载G-LiTH的对应高程数据,得到对应的高程误差,并统计不同冠层高度内,两组数据的误差统计指标(RMSE、MAE、R2、ME、不同冠层高度数量占比(N)、pearson相关系数、spearman相关系数、kendall相关系数),不同冠层高度内对应的高程误差及植被覆盖率。

    • 为探究强弱激光类型的ATLAS数据反演林下DTM精度,统计了强弱激光类型的ATLAS数据反演林下DTM精度如表2所示,散点图如图6所示。

      表 2  不同激光类型条件下ATLAS数据反演精度统计

      Table 2.  Estimation accuracy statistics of ATLAS data under different laser types

      Different laser intensities typesRMSE/mMAE/mR2ME/m
      Weak beam0.760.541.000.33
      Strong beam0.740.511.000.27
      Mean of different laser intensities types0.750.531.000.30

      图  6  不同激光类型条件下ATLAS数据与G-LiHT数据散点图。(a)弱波束;(b)强波束

      Figure 6.  Scatter diagram of ATLAS data and G-LiHT data under different laser types. (a) Weak beam; (b) Strong beam

      试验结果的定量结果如图6表2所示。其中,R2值均为1,说明ATL03数据反演的林下DTM数据与机载G-LiHT的林下DTM高度一致。ATLAS数据的RMSE平均值和MAE平均值分别为0.75 m和0.53 m。强、弱光束条件下,RMSE分别为0.74 m、0.76 m,MAE分别为0.51 m、0.54 m。研究结果显示,该研究的强弱光束反演结果相对参考文献[15-16]表现出更高的R2,与更低的RMSEMAE。试验结果显示,该研究中强光束类型相对弱光束类型表现出更优的反演精度。试验结果表明,激光强度对于ATLAS数据反演林下DTM存在影响,由于强、弱波束类型的标定激光发射强度比为4∶1[18],当测量同一类植被类型时,强波束类型可以获得更多的光子数据,用以描述被测地物,但是此次研究结果表明两种激光波束都可以为反演林下DTM提供科学数据。

    • 为探究冠层高度对于ATLAS强、弱激光束反演林下DTM的影响程度,统计不同冠层高度条件下强、弱波束数据反演精度评价标准及植被覆盖情况如表3所示,其中表3中的植被覆盖率为该区段冠层高度情况下覆盖占比最大的植被覆盖率[20]。不同冠层高度条件下两组数据误差的箱形图如图7所示。

      表 3  不同冠层高度条件下强、弱波束数据反演精度统计

      Table 3.  Estimation accuracy statistics of strong beam and weak beam data under different canopy heights

      Canopy heights/mRMSE/mMAE/mR2ME/mNpearsonspearmankendallVegetation coverage
      0-50.49/0.900.41/0.571.00/0.990.33/0.3221.39%/20.48%1.00/1.001.00/0.990.95/0.9327%(28.5%)/25%(35.7%)
      6-100.47/0.500.41/0.381.00/1.000.29/0.2522.37%/24.34%1.00/1.001.00/0.990.97/0.9427%(33.0%)/66%(36.4%)
      11-150.52/0.820.42/0.681.00/0.990.33/0.5627.20%/33.01%1.00/1.001.00/0.990.96/0.9436%(30.4%)/66%(61.8%)
      16-200.66/0.970.47/0.761.00/0.990.35/0.5013.28%/11.57%1.00/1.000.99/0.940.93/0.8436%(50.0%)/66%(51.0%)
      21-1.25/0.660.99/0.580.98/0.980.77/0.4115.76%/10.60%0.99/1.000.84/0.680.67/0.5236%(42.0%)/40%(30.9%)
      注:表格中数据分别为强波束数据/弱波束数据的精度统计数据

      图  7  不同冠层高度条件下ATLAS数据反演精度。(a)强波束;(b)弱波束

      Figure 7.  Estimation accuracy of ATLAS data under different canopy height conditions. (a) Strong beam; (b) Weak beam

      由于研究区处于亚热带季风湿润气候,森林冠层高度范围为0~34.8 m,故研究选择将植被冠层高度值分为5组,每组间隔5 m,如表3N所示每组冠层高度数据分布较为平均。

      表3中统计了各冠层高度区段内占比最高的植被覆盖率,由表3植被覆盖率情况可知,从研究区冠层高度整体变化趋势来看,随着冠层高度的增加,植被覆盖率也随之增加。但是,在弱光束冠层高度大于21 m的情况下,出现了植被覆盖率降低的情况。

      表3可知,从不同统计指标的整体变化趋势来看,在研究区冠层高度内,随冠层高度及植被覆盖率的增加,强波束激光束的MSE、MAE、ME也逐步增加。不同类型相关系数随植被冠层高度及植被覆盖率的增加,出现冠层高度与高程误差的相关性逐步降低的情况,该情况与参考文献[16]表现类似。其原因可能是由于随植被冠层高度及植被覆盖率的增加,ICESat-2/ATLAS发射的光子通过森林植被进而撞击地表和经过地表反射回ICESat-2接收器的概率降低,进而造成较大的高程误差。

      表3还可知,在研究区冠层高度范围为6~20 m的情况下,虽然主要植被覆盖率保持在66%,但是随冠层高度的增加,弱波束激光束的MSE、MAE、ME逐步增加。不同类型相关系数随植被冠层高度的增加,出现冠层高度与高程误差的相关性逐步降低的情况。其原因可能是由于随冠层高度的增加,地表反射回ICESat-2接收器的概率降低。这种情况表明,即使在植被覆盖率一定的情况下,随冠层高度的增加,误差也会随之增大。在冠层高度范围为0~5 m的条件下,出现RMSE与MAE不符合强波束激光类型的变化情况,其原因可能是由于弱激光束能量相对较低,在冠层高度为0~5 m的情况下,ICESat-2/ATLAS的光子较少,地面信号光子与冠层信号光子距离较近,难以通过ATL08的分类算法区分冠层信号光子与地面光子,故出现较大误差。在研究区冠层高度大于21 m的情况下,出现RMSE与MAE不符合误差随冠层高度增加而增加的情况,综合考虑冠层高度与植被覆盖率可知,虽然冠层高度增加,但是植被覆盖率出现降低情况,这种情况表明,冠层植被覆盖率也是影响光子云数据反演林下地形的一个重要因素。

    • 为探究ICESat-2/ATLAS不同激光类型数据反演林下DTM的精度,以美国Aiken郡的林区作为研究区,以G-LiHT数据作为验证数据,分析了ICESat-2/ATLAS的不同激光类型数据的反演精度,并着重讨论了冠层高度及植被覆盖率对于反演精度的影响得出以下结论:

      (1)研究结果显示,强波束反演精度为R2=1,RMSE=0.74 m,弱波束反演精度为R2=1,RMSE=0.76 m。虽然强波束反演林下DTM的精度表现优于弱波束,但是,强光束与弱光束的数据均可为反演林下DTM提供科学数据。

      (2)从研究区植被整体情况来看,随冠层高度及植被覆盖率的增加,ICESat-2/ATLAS的不同激光类型数据均出现误差逐步增加的情况。

      总体而言,研究针对ATLAS反演林下地面高程进行了初步探索,并探讨了冠层高度对地面高程误差的影响,这将为ATL03数据反演林下地面高程提供基础。但是,该研究主要分析了冠层高度对于ATLAS数据反演林下地面高程的影响,下一步的研究重点为植被覆盖率、森林类型等因素对于ATLAS数据反演林下地形的影响。

参考文献 (21)

目录

    /

    返回文章
    返回