留言板

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

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

面向ICESat-2单光子激光雷达的陆地目标测距误差分析

李泓成 周辉 李松 马跃 王玥

李泓成, 周辉, 李松, 马跃, 王玥. 面向ICESat-2单光子激光雷达的陆地目标测距误差分析[J]. 红外与激光工程, 2020, 49(11): 20200247. doi: 10.3788/IRLA20200247
引用本文: 李泓成, 周辉, 李松, 马跃, 王玥. 面向ICESat-2单光子激光雷达的陆地目标测距误差分析[J]. 红外与激光工程, 2020, 49(11): 20200247. doi: 10.3788/IRLA20200247
Li Hongcheng, Zhou Hui, Li Song, Ma Yue, Wang Yue. Analysis on ranging error of terrain targets for ICESat-2 single-photon lidar[J]. Infrared and Laser Engineering, 2020, 49(11): 20200247. doi: 10.3788/IRLA20200247
Citation: Li Hongcheng, Zhou Hui, Li Song, Ma Yue, Wang Yue. Analysis on ranging error of terrain targets for ICESat-2 single-photon lidar[J]. Infrared and Laser Engineering, 2020, 49(11): 20200247. doi: 10.3788/IRLA20200247

面向ICESat-2单光子激光雷达的陆地目标测距误差分析

doi: 10.3788/IRLA20200247
基金项目: 国家自然科学基金(41971302);装备预研航天科技联合基金(6141B06110102);国防科工局十三五民用航天技术预先研究项目(250000887)
详细信息
    作者简介:

    李泓成(1995-),男,硕士生,主要从事单光子激光雷达数据处理等方面的研究。Email:18883285813@163.com

    周辉(1979-),男,副教授,博士,主要从事激光遥感与空间探测等方面的研究。Email:zhouhui@whu.edu.cn

  • 中图分类号: TN958.98

Analysis on ranging error of terrain targets for ICESat-2 single-photon lidar

  • 摘要: 单光子激光雷达所获取的光子事件存在随机分布的特点,使得其激光测距值出现不确定性,从而降低了单光子激光雷达的测距精度。在不减小采样分辨率的情况下,采用累积邻近多光斑的光子事件来构建光子累积直方图,并基于反算的目标响应函数的时间重心来确定激光测距值。针对ICESat-2单光子激光雷达,以测距均方根误差(RMSE)和平均绝对误差(MAE)为性能指标,构建出一种综合考虑陆地地形、累积光斑数目和回波光子数等多因素对其测距误差影响的评估方法。同时,选取ICESat-2过境美国犹他州西瓦利城的某观测条带的随机地形数据进行验证分析。结果表明,所提出的激光测距值解算方法能够使该条带的测距RMSE值由114.25 cm降低到63.84 cm,MAE值由70.97 cm降低到48.52 cm,均优于ICESat-2数据产品提供的137.96 cm RMSE值和97.24 cm MAE值,这对提升单光子激光雷达在陆地区域的测距精度具有一定的借鉴作用。
  • 图  1  单光子激光雷达测距误差分析流程

    Figure  1.  Analysis process of ranging error for single-photon lidar

    图  2  三种坡度区域累积多光斑前后的测距误差图

    Figure  2.  Range errors before and after accumulating multiple spots at three slope areas

    图  3  三种坡度区域MAE和RMSE值之和分布图

    Figure  3.  Distribution of the sum values of MAE and RMSE at three slope areas

    图  4  三种坡度区域的测距误差起伏分布

    Figure  4.  Fluctuation distributions of range error at three slope areas

    图  5  最优累积光斑数目与平均回波光子数的关系图

    Figure  5.  Relationship between the optimal number of accumulative spots and the average echo photons

    图  6  平均回波光子数与测距误差的关系图

    Figure  6.  Relationship between the average echo photons and range error

    图  7  研究区域与对应的实际高程和坡度分布图

    Figure  7.  Study area with corresponding true elevation and slope distribution

    图  8  累积光斑前后的高程对比图和ATL08产品的条带高程

    Figure  8.  Comparison of elevation before and after the accumulative spots and the stripe elevation label of ATL08 products

    表  1  ICESat-2单光子激光雷达的主要技术参数

    Table  1.   Main technological parameters of ICESat-2 single-photon lidar

    Parameters of ICESat-2Value
    Detector channels 16 (strong)/4 (weak)
    Pulse width/ns 0.64
    Beam numbers 6 (3 strong/3 weak)
    Energy ratio 4∶1 (strong: weak)
    RMS radius of spot/m 4.375
    Along-track distance of spot/m 0.7
    Dead time/ns 3.2
    Time resolution/ps 200
    下载: 导出CSV

    表  2  机载激光雷达点云数据的主要参数

    Table  2.   Main parameters of the lidar point clouds

    TerrainSlope/(°)Latitude and longitudePoint density/pts·m−2
    Low0-5(172°38′17″E, 40°42′37″S)9.54
    Moderate5-15(107°13′45″W, 37°59′33″N)15.83
    High20-35(115°68′50″W, 44°20′60″N)13.71
    下载: 导出CSV

    表  3  累积光斑前后所对应的不同坡度区域的MAERMSE

    Table  3.   MAE and RMSE at the test areas with different slopes before and after accumulating spots

    TypeTerrainMAE/cmRMSE/cm
    No accumulationLow12.0717.21
    Moderate30.5142.25
    High99.80127.69
    AccumulationLow4.906.66
    Moderate10.9513.64
    High35.4343.78
    下载: 导出CSV
  • [1] Zhou H, Chen Y, Hyyppä J, et al. An overview of the laser ranging method of space laser altimeter [J]. Infrared Physics & Technology, 2017, 86: 147-158.
    [2] 卜禹铭, 杜小平, 曾朝阳, 等. 无扫描激光三维成像雷达研究进展及趋势分析[J]. 中国光学, 2018, 11(5): 711-727. doi:  10.3788/co.20181105.0711

    Bu Yuming, Du Xiaoping, Zeng Zhaoyang, et al. Research progress and trend analysis of non-scanning laser 3D imaging radar [J]. Chinese Optics, 2018, 11(5): 711-727. (in Chinese doi:  10.3788/co.20181105.0711
    [3] 赵泉华, 陈为多, 王玉, 等. 偏正态全波激光雷达数据的可变分量波形分解[J]. 光学 精密工程, 2018, 26(1): 161-171. doi:  10.3788/OPE.20182601.0161

    Zhao Quanhua, Chen Weiduo, Wang Yu, et al. Full-waveform LiDAR data decomposition based on skew-normal distribution with unknown number of components [J]. Optics and Precision Engineering, 2018, 26(1): 161-171. (in Chinese doi:  10.3788/OPE.20182601.0161
    [4] Wang Fei, Zhao Yuan, Zhang Yu, et al. Range accuracy limitation of pulse ranging systems based on Geiger mode single-photon detectors [J]. Applied Optics, 2010, 49(29): 5561-5566. doi:  10.1364/AO.49.005561
    [5] Zhou Xin, Sun Jianfeng, Jiang Peng, et al. Influence investigation on ranging performance for range-gated Geiger-mode avalanche photodiode ladar [J]. Applied Optics, 2018, 57(10): 2667-2674. doi:  10.1364/AO.57.002667
    [6] Oh M S, Kong H J, Kim T H, et al. Reduction of range walk error in direct detection laser radar using a Geiger mode avalanche photodiode [J]. Optics Communications, 2010, 283(2): 304-308. doi:  10.1016/j.optcom.2009.10.009
    [7] 黄科, 李松, 马跃, 等. 单光子激光测距的漂移误差理论模型及补偿方法[J]. 物理学报, 2018, 67(6): 064205.

    Huang Ke, Li Song, Ma Yue, et al. Theoretical model and correction method of range walk error for single-photon laser ranging [J]. Acta Physica Sinica, 2018, 67(6): 064205. (in Chinese
    [8] Popescu S C, Zhou T, Nelson R, et al. Photon counting LiDAR: An adaptive ground and canopy height retrieval algorithm for ICESat-2 data [J]. Remote Sensing of Environment, 2018, 208: 154-170. doi:  10.1016/j.rse.2018.02.019
    [9] Neuenschwander A, Pitts K. The ATL08 land and vegetation product for the ICESat-2 Mission [J]. Remote Sensing of Environment, 2019, 221: 247-259. doi:  10.1016/j.rse.2018.11.005
    [10] Kwok R, Markus T, Kurtz N T, et al. Surface height and sea ice freeboard of the Arctic Ocean from ICESat-2: Characteristics and early results [J]. Journal of Geophysical Research: Oceans, 2019, 124(10): 6942-6959. doi:  10.1029/2019JC015486
    [11] 杨子健, 陈锋, 李抄, 等. 微脉冲激光雷达中的光子计数死区时间瞬态效应[J]. 光学 精密工程, 2015, 23(2): 408-414. doi:  10.3788/OPE.20152302.0408

    Yang Zijian, Chen Feng, Li Chao, et al. Transient effect of dead time of photon-counting in micro-pulse lidar [J]. Optics and Precision Engineering, 2015, 23(2): 408-414. (in Chinese doi:  10.3788/OPE.20152302.0408
    [12] 陈露, 高志山, 袁群, 等. 星载激光测高仪距离参数地面标定方法[J]. 中国光学, 2019, 12(4): 896-904. doi:  10.3788/co.20191204.0896

    Chen Lu, Gao Zhishan, Yuan Qun, et al. Ground calibration method for distance parameters of satellite-borne laser altimeter [J]. Chinese Optics, 2019, 12(4): 896-904. (in Chinese doi:  10.3788/co.20191204.0896
    [13] 沈姗姗, 陈钱, 何伟基, 等. 单光子测距系统性能优化研究和实现[J]. 红外与激光工程, 2016, 45(2): 0217001.

    Shen Shanshan, Chen Qian, He Weiji, et al. Research and realization on performance of single photon counting ranging system optimizing[J]. Infrared and Laser Engineering, 2016, 45(2): 0217001.
    [14] 许艺腾, 李国元, 邱春霞, 等. 基于地形相关和最小二乘曲线拟合的单光子激光数据处理技术[J]. 红外与激光工程, 2019, 48(12): 1205004. doi:  10.3788/IRLA201948.1205004

    Xu Yiteng, Li Guoyuan, Qiu Chunxia, et al. Single photon laser data processing technology based on Terrain Correlation and least square curve fitting [J]. Infrared and Laser Engineering, 2019, 48(12): 1205004. (in Chinese doi:  10.3788/IRLA201948.1205004
    [15] Markus T, Neumann T, Martino A, et al. The Ice, Cloud, and land Elevation Satellite-2 (ICESat-2): science requirements, concept, and implementation [J]. Remote Sensing of Environment, 2017, 190: 260-273. doi:  10.1016/j.rse.2016.12.029
    [16] Harding D J, Bufton J L, Frawley J J. Satellite laser altimetry of terrestrial topography: vertical accuracy as a function of surface slope, roughness, and cloud cover [J]. IEEE Transactions on Geoscience and Remote Sensing, 1994, 32(2): 329-339. doi:  10.1109/36.295048
    [17] 张文豪, 李松, 张智宇, 等. 利用波形匹配实现卫星激光测高脚点精确定位的方法[J]. 红外与激光工程, 2018, 47(11): 1117007. doi:  10.3788/IRLA201847.1117007

    Zhang Wenhao, Li Song, Zhang Zhiyu, et al. Using waveform matching to precisely locate footprints of a satellite laser altimeter [J]. Infrared and Laser Engineering, 2018, 47(11): 1117007. (in Chinese doi:  10.3788/IRLA201847.1117007
  • [1] 刘博, 蒋贇, 王瑞, 陈臻, 赵彬, 黄风云, 杨玉强.  全天时单光子激光雷达技术进展与系统评价 . 红外与激光工程, 2023, 52(1): 20220748-1-20220748-15. doi: 10.3788/IRLA20220748
    [2] 蒋筱朵, 赵晓琛, 冒添逸, 何伟基, 陈钱.  采用传感器融合网络的单光子激光雷达成像方法 . 红外与激光工程, 2022, 51(2): 20210871-1-20210871-7. doi: 10.3788/IRLA20210871
    [3] 李雨佳, 周晓青, 李国元, 郭金权, 马跃, 谌一夫.  星载单光子激光雷达浅水测深技术研究进展和展望 . 红外与激光工程, 2022, 51(10): 20220003-1-20220003-10. doi: 10.3788/IRLA20220003
    [4] 邬淼, 陆俣, 冒添逸, 何伟基, 陈钱.  单光子激光雷达的时间相关多深度估计 . 红外与激光工程, 2022, 51(2): 20210885-1-20210885-9. doi: 10.3788/IRLA20210885
    [5] 李致廷, 刘长明, 王与烨, 常继英, 陈锴, 李吉宁, 钟凯, 徐德刚, 姚建铨.  基于单光子探测的目标光学散射特性研究 . 红外与激光工程, 2022, 51(9): 20210825-1-20210825-8. doi: 10.3788/IRLA20210825
    [6] 岳广, 孙文邦, 张星铭, 李铜哨.  航空面阵图像拼接的累积误差消除方法 . 红外与激光工程, 2021, 50(9): 20200529-1-20200529-9. doi: 10.3788/IRLA20200529
    [7] 刘心溥, 元志安, 王玲, 许可, 万建伟.  副载波调制水下激光雷达测距性能仿真 . 红外与激光工程, 2020, 49(S2): 20200193-20200193. doi: 10.3788/IRLA20200193
    [8] 侯阿慧, 胡以华, 赵楠翔, 董骁, 曾祥熙.  光子探测距离漂移误差与大气湍流效应建模分析 . 红外与激光工程, 2020, 49(S2): 20200192-20200192. doi: 10.3788/IRLA20200192
    [9] 胡恒, 杜小平, 张朋, 渠晋.  调制函数对偏振调制激光测距精度影响分析 . 红外与激光工程, 2018, 47(10): 1006004-1006004(7). doi: 10.3788/IRLA201847.1006004
    [10] 黄民双, 刘晓晨, 马鹏.  脉冲飞行时间激光测距系统中周期误差补偿 . 红外与激光工程, 2018, 47(3): 317004-0317004(5). doi: 10.3788/IRLA201847.0317004
    [11] 姚宗辰, 张合, 张祥金, 黄正祥, 岳连永.  探测距离对激光定距引信探测能力影响的分析与评价 . 红外与激光工程, 2018, 47(3): 303005-0303005(9). doi: 10.3788/IRLA201847.0303005
    [12] 许晓飞, 张丽, 王森.  光子计数探测器的响应函数标定及校正方法 . 红外与激光工程, 2017, 46(12): 1217005-1217005(6). doi: 10.3788/IRLA201746.1217005
    [13] 熊显名, 李三龙, 张文涛, 张良, 李鹏飞.  激光雷达测距系统光源驱动设计 . 红外与激光工程, 2016, 45(2): 206004-0206004(5). doi: 10.3788/IRLA201645.0206004
    [14] 都琳, 孙华燕, 张廷华, 王帅.  基于PCA的相机响应函数模型标定算法 . 红外与激光工程, 2016, 45(10): 1026001-1026001(9). doi: 10.3788/IRLA201645.1026001
    [15] 杨驰, 胡文怡, 罗敏, 孙莹莹, 周辉.  椭圆高斯足印对星载激光测高仪测距值及其误差的影响 . 红外与激光工程, 2016, 45(7): 717003-0717003(7). doi: 10.3788/IRLA201645.0717003
    [16] 马跃, 阳凡林, 卢秀山, 冯成凯, 李松.  对地观测星载激光测高系统高程误差分析 . 红外与激光工程, 2015, 44(3): 1042-1047.
    [17] 周辉, 李松, 王良训, 郑国兴.  噪声对星载激光测高仪测距误差的影响 . 红外与激光工程, 2015, 44(8): 2256-2261.
    [18] 王良训, 周辉, 李子乐, 刘国根, 王虹, 王雅培.  面向星载激光测高仪的陆地目标响应函数仿真 . 红外与激光工程, 2015, 44(11): 3424-3430.
    [19] 冯晓宇, 宗肖颖.  一种去除拼接干涉图中累积误差的简单方法 . 红外与激光工程, 2014, 43(3): 997-1001.
    [20] 王竹萍, 刘汝卿, 郭仙菊, 张珂殊.  应用于便携式激光雷达测距仪的单芯片全集成前置放大电路系统 . 红外与激光工程, 2013, 42(8): 2050-2054.
  • 加载中
图(8) / 表(3)
计量
  • 文章访问数:  620
  • HTML全文浏览量:  154
  • PDF下载量:  84
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-06-05
  • 修回日期:  2020-10-01
  • 刊出日期:  2020-11-25

面向ICESat-2单光子激光雷达的陆地目标测距误差分析

doi: 10.3788/IRLA20200247
    作者简介:

    李泓成(1995-),男,硕士生,主要从事单光子激光雷达数据处理等方面的研究。Email:18883285813@163.com

    周辉(1979-),男,副教授,博士,主要从事激光遥感与空间探测等方面的研究。Email:zhouhui@whu.edu.cn

基金项目:  国家自然科学基金(41971302);装备预研航天科技联合基金(6141B06110102);国防科工局十三五民用航天技术预先研究项目(250000887)
  • 中图分类号: TN958.98

摘要: 单光子激光雷达所获取的光子事件存在随机分布的特点,使得其激光测距值出现不确定性,从而降低了单光子激光雷达的测距精度。在不减小采样分辨率的情况下,采用累积邻近多光斑的光子事件来构建光子累积直方图,并基于反算的目标响应函数的时间重心来确定激光测距值。针对ICESat-2单光子激光雷达,以测距均方根误差(RMSE)和平均绝对误差(MAE)为性能指标,构建出一种综合考虑陆地地形、累积光斑数目和回波光子数等多因素对其测距误差影响的评估方法。同时,选取ICESat-2过境美国犹他州西瓦利城的某观测条带的随机地形数据进行验证分析。结果表明,所提出的激光测距值解算方法能够使该条带的测距RMSE值由114.25 cm降低到63.84 cm,MAE值由70.97 cm降低到48.52 cm,均优于ICESat-2数据产品提供的137.96 cm RMSE值和97.24 cm MAE值,这对提升单光子激光雷达在陆地区域的测距精度具有一定的借鉴作用。

English Abstract

    • 星载单光子激光雷达是近几年发展起来的新一代对地观测载荷,它采用低能量、高重频和多波束激光脉冲[1],通过采集地面目标的单光子事件来获取高精度和高分辨率目标信息,从而使星载单光子激光雷达的使用性能得到大幅提升,它也成为未来对地观测的一个重要发展方向。

      单光子激光雷达只能探测到信号的有无,不能探测到信号的大小,无法像传统线性激光雷达一样能够采集到完整的回波波形后通过求时间重心来确定激光测距值[2-3],因此,其测距精度会受到光子事件随机性的影响。Wang和Zhou通过减小发射脉冲宽度,增加发射激光能量以及使用窄带滤光片等硬件设计措施来提升激光测距精度[4-5]。Oh和黄科针对静态目标,通过累积几千次的光子事件产生光子直方图,并根据光子直方图的时间重心来估算激光测距值,从而修正光子事件随机性引入的误差[6-7]。然而,星载单光子激光雷达始终处于高速运动状态,每次测量的目标都是在动态变化的,各种噪声会影响激光测距值的解算。因而,目前提升星载单光子激光雷达测距精度的主要手段还是集中在其点云数据的精细化处理方式上[8]。但是,经过单光子数据处理后得到的有效信号的点云事件仍然存在一定的随机性。因此,ICESat-2研究团队在牺牲空间分辨率的情况下,针对陆地和植被目标,选择100 m范围内目点云高程的中值作为目标高程值[9],而海冰高程产品则是基于累积150个左右的光子事件所产生的光子直方图来获取初始海冰面高程值[10],不过没有提供相关的选取依据。

      针对光子事件随机性对测距精度影响的问题,在不牺牲测量空间分辨率的情况下,通过累积邻近有限数量的光子事件产生光子直方图,并通过理论推导反算出目标响应函数获取激光测距结果。同时,以测距均方根误差和平均绝对误差为性能指标,构建一种综合考虑陆地地形、累积光斑数目和回波光子数对测距误差影响的分析方法。采用ICESat-2单光子激光雷达的实测数据和产品结果,选择随机起伏的陆地条带,验证激光测距值解算及其误差分析方法,这为单光子激光雷达测距值的获取提供技术支撑。

    • 星载激光雷达所接收的脉冲回波信号可以近似为发射脉冲信号与目标响应函数的卷积。对于常见的线性陆地目标而言,其脉冲回波信号的理论表达式为:

      $$ \begin{split} s\left(t\right)=& g\left(t\right)\cdot m\left(t\right) =g\left(t\right)\cdot \\ &\dfrac{\gamma }{\sqrt{2\pi }\delta }exp\left[-\dfrac{{\left(t-{t}_{h}\right)}^{2}}{2{\delta }^{2}}\right] \end{split} $$ (1)

      式中:$ g\left(t\right) $$ m\left(t\right) $分别表示发射激光脉冲和目标响应函数;$ \gamma $$ {t}_{h} $$ \delta $则是目标响应函数的面积、时间重心和均方根脉宽,其中$ {t}_{h} $直接反映激光测距值的大小。

      单光子激光雷达仅能捕获到离散的光子事件,无法获取脉冲回波信号。然而,单光子激光雷达的重复频率很高,相邻光斑的间距很小,使得一段沿轨距离范围内的目标具有一定的相关性,因而,可以将当前光斑作为中心,按照一定的时间(高程)分辨率,通过累积前后相邻的多个光斑的光子事件来形成光子计数直方图。

      将光子直方图进行归一化运算,归一化的光子直方图直接反映出光子事件的概率分布,基于单光子探测模型[11],该概率近似满足以下关系式:

      $${P_r}\left( i \right) = \left[ {1 - \mathop \sum \limits_{j = i - d}^{i - 1} {P_r}\left( j \right)} \right]\left\{ {1 - exp\left[ { - K\left( i \right)} \right]} \right\} $$ (2)

      式中:$ {P}_{r}\left(i\right) $$ K\left(i\right) $表示在第$ i $个时间片的光子事件概率和光子数;$ d $表示死区时间所对应的时间片数目。因此,由光子直方图反算可以得到脉冲回波分布为:

      $$ K\left(i\right)=-{\rm{ln}}\left[1-\dfrac{{P}_{r}\left(i\right)}{1-\displaystyle\sum _{j=i-d}^{i-1}{P}_{r}\left(j\right)}\right] $$ (3)

      考虑到累积光斑数是有限的,则反算得到的脉冲回波分布与理论的脉冲回波之间会存在一定的区别。

      基于公式(1),可以反算得到目标响应函数为:

      $$ M\left(i\right)={f}^{-1}\left\{\dfrac{f\left[K\left(i\right)\right]}{f\left[g\left(i\right)\right]}\right\} $$ (4)

      式中:$ f $$ {f}^{-1} $分别表示傅里叶变换和傅里叶逆变换。

      由于常见的线性目标所对应的目标响应函数在理论上满足高斯分布,则采用高斯模型对反算的目标响应函数$ M $进行非线性最小二乘拟合得到高斯函数$ {L}_{n} $。通过解算拟合高斯函数的时间重心[12],即可得到激光测距值的表达形式:

      $${{{R}}_{{n}}} = \frac{{{{c}}\tau }}{2}\frac{{\displaystyle\sum\nolimits_{{{i}} = 1}^{{m}} {{i}} {{{L}}_{{n}}}\left( {{i}} \right)}}{{\displaystyle\sum\nolimits_{{{i}} = 1}^{{m}} {{{{L}}_{{n}}}} \left( {{i}} \right)}} $$ (5)

      式中:$ n $表示累积光斑的个数;$ c $表示光速;$ \tau $表示单光子激光雷达的时间片长度;$ {L}_{n}\left(i\right) $表示拟合高斯函数在第$ i $个时间片对应的数值;$ m $表示时间片的总数。

      当然,激光光斑内的实际陆地地形可能会偏离线性分布,导致目标响应函数无法采用高斯函数来表示,因此,在解算激光测距值时,可以将公式(5)中的$ {L}_{n} $替换为公式(4)中的$ M $,这样得到光斑内的陆地目标的平均距离。

      星载单光子激光雷达的地形起伏、累积光斑数和平均回波光子数是影响陆地目标测距误差的三个重要因素。考虑到单光子激光雷达的测距误差无法采用解析形式来描述,则采用如图1所示的测距误差评估流程。

      图  1  单光子激光雷达测距误差分析流程

      Figure 1.  Analysis process of ranging error for single-photon lidar

      单光子激光雷达测距误差的评估分析主要包括四个过程:单光子点云仿真[13]、激光测距值解算、参考距离提取和测距误差分析。其中,单光子点云仿真是对机载激光雷达点云数据所生成的伪回波进行随机抽样后,再经过死区时间剔除所产生的。测距值是由反算的目标响应函数的时间重心解算得到的,参考距离是光斑内目标的加权平均高程。测距误差即为解算测距值与参考距离之间的差异。

      参考距离是在考虑激光光斑的高斯空间分布和机载激光点云的强度情况下得到的。以目标的相对高程来表征参考距离,则其表达式为:

      $$ H=\dfrac{\displaystyle\sum _{j=1}^{l}{\rm{exp}}\left[-\frac{{d}^{2}\left(j\right)}{2{r}^{2}}\right]e\left(j\right)h\left(j\right)}{\displaystyle\sum _{j=1}^{l}{\rm{exp}}\left[-\dfrac{{d}^{2}\left(j\right)}{2{r}^{2}}\right]e\left(j\right)} $$ (6)

      式中:$ e\left(j\right) $$ h\left(j\right) $$ d\left(j\right) $和分别为第$ j $个机载激光点的强度、相对高程及其与光斑中心的水平距离;$ l $表示当前光斑中机载激光点云的总个数,$ r $表示当前光斑的均方根半径。

      采用均方根误差(RMSE)和平均绝对误差(MAE)形式来表征某个机载条带目标的测距误差[14],其中,RMSE和MAE分别代表条带目标测距值的随机误差和系统误差。根据累积多光斑事件得到的测距值以及参考距离的表达式,可以得到测距值的RMSE和MAE的形式为:

      $$ {{RMSE}}\left(n\right)=\sqrt{\frac{1}{p}\sum\limits _{k=1}^{p}{\left|H\left(k\right)-{R}_{n}\left(k\right)\right|}^{2}} $$ (7)
      $$ {{MAE}}\left(n\right)=\frac{1}{p}\sum\limits _{k=1}^{p}\left|H\left(k\right)-{R}_{n}\left(k\right)\right| $$ (8)

      式中:$ p $表示沿轨光斑点的总个数。

    • ICESat-2星载单光子激光雷达是目前国际上唯一一个对地观测的光子计数激光雷达,其主要目的在于利用单光子激光雷达的高精度测距结果监测极区冰层的消融情况。单光子激光雷达发射系统和接收系统的主要技术参数如表1所示[15]

      表 1  ICESat-2单光子激光雷达的主要技术参数

      Table 1.  Main technological parameters of ICESat-2 single-photon lidar

      Parameters of ICESat-2Value
      Detector channels 16 (strong)/4 (weak)
      Pulse width/ns 0.64
      Beam numbers 6 (3 strong/3 weak)
      Energy ratio 4∶1 (strong: weak)
      RMS radius of spot/m 4.375
      Along-track distance of spot/m 0.7
      Dead time/ns 3.2
      Time resolution/ps 200

      不同地形起伏会影响激光测距结果,因而选取平坡、缓坡和陡坡三种目标类型[16]作为分析对象。三种目标所对应的机载激光点云数据分别位于新西兰金海湾区域,美国科罗拉达州西南部滑坡区域和美国爱达荷州岩基泥石流区域,其选取的区域面积大小均为150 m×750 m,三种区域的点密度已经足够描述伪回波波形[17],相关数据的主要参数如表2所示。

      表 2  机载激光雷达点云数据的主要参数

      Table 2.  Main parameters of the lidar point clouds

      TerrainSlope/(°)Latitude and longitudePoint density/pts·m−2
      Low0-5(172°38′17″E, 40°42′37″S)9.54
      Moderate5-15(107°13′45″W, 37°59′33″N)15.83
      High20-35(115°68′50″W, 44°20′60″N)13.71
    • 以ICESat-2单光子激光雷达的强激光作为测量光束,预估陆地目标的平均回波光子数为3,累积光斑数目选取为11,仿真出三种地形的测距误差分布图。

      图2中可以看出,累积光斑前后不同坡度的地形都会对测距误差造成影响。在累积多光斑前,三种不同坡度区域的测距误差值波动性很大,平坡、缓坡和陡坡的沿轨光斑点的最大测距误差值分别达到了100 cm、200 cm和400 cm;在累积多光斑后,误差值的波动趋近于平稳,测距误差值得到明显减小,其沿轨光斑点的最大测距误差值也分别减小到50 cm、100 cm和200 cm。通过统计分析得到,三种坡度区域所对应的均方根误差和平均绝对误差分别见表3所示。

      图  2  三种坡度区域累积多光斑前后的测距误差图

      Figure 2.  Range errors before and after accumulating multiple spots at three slope areas

      表 3  累积光斑前后所对应的不同坡度区域的MAERMSE

      Table 3.  MAE and RMSE at the test areas with different slopes before and after accumulating spots

      TypeTerrainMAE/cmRMSE/cm
      No accumulationLow12.0717.21
      Moderate30.5142.25
      High99.80127.69
      AccumulationLow4.906.66
      Moderate10.9513.64
      High35.4343.78

      表3中的结果可以得出,累积光斑前后,随着坡度的增加,条带的测距误差值均会相应增加。且无论平坡、缓坡还是陡坡区域,累积光斑后的RMSEMAE均能得到改善,三种地形的RMSE值分别减小10.55 cm、28.61 cm和83.91 cm,MAE值分别减小7.17 cm、19.56 cm和64.37 cm。

    • 理想情况下,对于短死区的单光子激光雷达系统而言,当光斑大小固定不变且累积光斑数目趋近于无穷大时,目标响应的时间重心对应于实际的激光测距值。然而,对于有限的累积数,累积直方图分布与回波光子数分布之间出现偏差,其分布形式与累积数有关,这使得不同累积数条件下的激光测距值也会发生改变。

      $\Delta {\rm{\rho }}\left( {{n}} \right) = {{RMSE}}\left( {{n}} \right) + {{MAE}}\left( {{n}} \right)$来表征激光测距误差的大小,模拟出三种目标区域所对应的测距误差与累积光斑数目之间的关系曲线,如图3所示。

      图  3  三种坡度区域MAE和RMSE值之和分布图

      Figure 3.  Distribution of the sum values of MAE and RMSE at three slope areas

      图3的仿真结果中可以看出,随着累积光斑数目的增加,三种目标区域的测距误差值都会逐渐减小,但是当累积光斑数目超过某个值时,其绝对误差值又会逐渐增大。三种目标区域所对应的最优累积光斑数目取值不同,对于平坡、缓坡和陡坡,最优的累积光斑数分别为19、23和21。考虑到星载单光子激光雷达测量过程中,无法预知被测陆地区域目标的坡度先验知识,因此,采用以下原则来确定最终的累积光斑数目优化值:

      $$ {{N}} = \left\{ {n|\mathop \sum \nolimits_{\mu = 1}^3 \Delta {{\rm{\rho }}_\mu }\left( n \right) = min} \right\} $$ (9)

      式中:$ \Delta {{\rm{\rho }}_\mu }\left( {\mu = 1,2,3} \right)$分别表示平坡、缓坡和陡坡目标的测距误差。基于该原则得到,对于陆地目标,当累积光斑数为21个时,ICESat-2单光子激光雷达的测距误差达到最小。

    • 考虑到单光子事件仿真随机性会引起激光测距误差的起伏,在固定光斑累积数为21,回波光子数为3的情况下,针对平坡、缓坡和陡坡三种地形,重复100次仿真模拟,得到其测距误差的起伏分布,如图4所示。

      图  4  三种坡度区域的测距误差起伏分布

      Figure 4.  Fluctuation distributions of range error at three slope areas

      图4可以看出,随着仿真次数的变化,三种地形所对应的测距误差会发生微小起伏,其平均值分别为10.86 cm、16.75 cm和66.40 cm,其标准差分别为0.70 cm、1.25 cm和4.22 cm,对应的变异系数(标准差/平均值)分别为6.45%,7.46%和6.36%。上述结果表明,测距误差的起伏远小于其均值水平,即累积光斑数目后的测距误差具备收敛性。

    • 若陆地目标的回波光子数发生改变,则最终的优化累积光斑数目和对应的测距误差也会相应地出现变化,图5仿真模拟了不同回波光电子数条件下的最优累积光斑数目的分布。

      图5的结果可以看出,随着平均回波光子数的增加,最优累积光斑数会逐渐减小,当平均回波光子数超过三个时,累加光斑数趋于稳定。当平均回波光子数为10时,对应的累积光斑数目为15,这与ICESat-2提供的ATL07产品中海冰累积150个光子是相吻合的,这也证明了文中累积光斑数优化方法的正确性。

      图  5  最优累积光斑数目与平均回波光子数的关系图

      Figure 5.  Relationship between the optimal number of accumulative spots and the average echo photons

      在最优累积光斑数目的情况下,不同回波光子数对应的最小测距误差也是不同的,图6给出了三种区域平均回波光子数与测距误差的关系图。

      图  6  平均回波光子数与测距误差的关系图

      Figure 6.  Relationship between the average echo photons and range error

      图6的结果表明,随着回波光子数的增加,相应的测距误差值会逐渐减小,当回波光子数超过四个时,三种地形的测距误差会趋于稳定,因此,适当地增加回波光子数量有利于提高测距精度。

    • 以ICESat-2卫星过境美国犹他州西瓦利城附近陆地区域(40°42′ N, 112°7′ W)的一个条带数据作为分析对象,其条带的总长度为420 m,对应的光斑数为600个。其中,滤噪后的光子事件数据及其对应的目标高程数据分别存储在ATL03和ATL08产品中。根据ALT03产品中光子事件的经纬度信息,选取相应的机载激光雷达点云数据作为参考源,其点密度为11.93 pts/m2,如图7所示。

      图  7  研究区域与对应的实际高程和坡度分布图

      Figure 7.  Study area with corresponding true elevation and slope distribution

      通过对图7沿轨光斑的坡度统计分析,光斑内目标坡度为0°~5°、5°~20°和20°~35°的光斑数分别占总光斑数的3%、70%和27%,其地形存在不规则的随机起伏特点。

      为减少噪声对测距误差的影响,仅选取ALT03产品中置信度标记为4的光子事件进行处理分析,该条带上的平均回波光子数为2.2个,对应的累积光斑数为21。采用累积21个光斑的光子事件所得到的光子直方图,解算得到沿轨条带上的目标高程。将该高程与累加前的目标高程以及ICESat-2提供的ATL08产品中的目标高程进行比较,如图8所示。

      图  8  累积光斑前后的高程对比图和ATL08产品的条带高程

      Figure 8.  Comparison of elevation before and after the accumulative spots and the stripe elevation label of ATL08 products

      通过图8显示的结果计算分析,累积光斑前后RMSE值由114.25 cm减小到了63.84 cm,MAE值由70.97 cm减小到了48.52 cm,两者之和减小的幅度达到72.86 cm。将ATL08产品提供的间隔100 m范围的A、B、C、D和E这五个高程值点与目标参考高程值比较,计算得到其RMSE值和MAE值分别是137.96 cm和97.24 cm,两者的值均大于累积21个光斑后该条带测出的结果值。

      相比于仿真的误差分析结果,实际地形的测距误差结果略偏大,这主要是由于ICESat-2过境的这片条带地形的相关性不足所导致的,使得某些位置处的目标回波偏离理想的单高斯回波信号,存在多个模式回波的混叠。

    • 文中从星载单光子激光雷达测距模型出发,提出了一种通过累积直方图反算出目标响应函数的时间重心来确定激光测距值的方法,分析了陆地地形、累积光斑数目和回波光子数对测距误差的影响。从数值模拟和实际数据验证的角度,论证了在不降低采样分辨率的条件下,该方法的正确性与可行性。

      通过仿真分析表明,测距精度的提升与目标起伏程度有关,目标越平坦,测距精度提升效果就越好,反之亦然;累积光斑数的优化结果与目标平均回波光子数存在反比关系,当平均回波光子数超过三个时,累积光斑数的优化值趋近于平稳。随着平均回波光子数的增加,测距误差会逐渐减小。同时,多次重复仿真表明,所提出的激光测距解算方法具有较好的收敛性。然而,激光测距精度提升效果受地形起伏的制约。若地形无法采用单一的平面形式来描述,其测距精度的提升效果就会下降。

      所提出的激光测距解算方法有利于提高单光子激光雷达在陆地目标的测距精度,这对我国未来单光子激光雷达数据处理具有一定的参考意义。

参考文献 (17)

目录

    /

    返回文章
    返回