-
星载激光雷达所接收的脉冲回波信号可以近似为发射脉冲信号与目标响应函数的卷积。对于常见的线性陆地目标而言,其脉冲回波信号的理论表达式为:
$$ \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所示的测距误差评估流程。
单光子激光雷达测距误差的评估分析主要包括四个过程:单光子点云仿真[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-2 Value 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
Terrain Slope/(°) Latitude and longitude Point density/pts·m−2 Low 0-5 (172°38′17″E, 40°42′37″S) 9.54 Moderate 5-15 (107°13′45″W, 37°59′33″N) 15.83 High 20-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 累积光斑前后所对应的不同坡度区域的MAE与RMSE
Table 3. MAE and RMSE at the test areas with different slopes before and after accumulating spots
Type Terrain MAE/cm RMSE/cm No accumulation Low 12.07 17.21 Moderate 30.51 42.25 High 99.80 127.69 Accumulation Low 4.90 6.66 Moderate 10.95 13.64 High 35.43 43.78 从表3中的结果可以得出,累积光斑前后,随着坡度的增加,条带的测距误差值均会相应增加。且无论平坡、缓坡还是陡坡区域,累积光斑后的RMSE和MAE均能得到改善,三种地形的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可以看出,随着仿真次数的变化,三种地形所对应的测距误差会发生微小起伏,其平均值分别为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的结果表明,随着回波光子数的增加,相应的测距误差值会逐渐减小,当回波光子数超过四个时,三种地形的测距误差会趋于稳定,因此,适当地增加回波光子数量有利于提高测距精度。
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值,这对提升单光子激光雷达在陆地区域的测距精度具有一定的借鉴作用。Abstract: The photon events acquired by the single-photon lidar have the characteristics of random distribution, which makes the laser ranging value be uncertain and reduces the ranging accuracy of the single-photon lidar. Under the circumstance that the sampling resolution was not decreased, an accumulative histogram of photon was constructed by accumulating photon events from the multiple contiguous spots. The range was determined on the basis of the time centroid of the derived target response function. As for the ICESat-2 single-photon lidar, a comprehensive evaluation approach of the range error was put forward by considering the effects of the terrestrial relief, the accumulative spot numbers and returned photons on the range error according to the root mean square error (RMSE) and mean absolute error (MAE). The random terrain on a strip of ICESat-2 in the West Valley City of American Utah was selected to validate the proposed method. The results demonstrate that the RMSE value is reduced from 114.25 cm to 63.84 cm and the MAE value is reduced from 70.97 cm to 48.52 cm. They outweigh the RMSE of 137.96 cm and MAE of 97.24 cm from the ICESat-2 products. The conclusions might provide some guidance for improving the ranging precision on the terrain.
-
Key words:
- single-photon lidar /
- cumulative histogram /
- target response function /
- ranging error
-
表 1 ICESat-2单光子激光雷达的主要技术参数
Table 1. Main technological parameters of ICESat-2 single-photon lidar
Parameters of ICESat-2 Value 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 表 2 机载激光雷达点云数据的主要参数
Table 2. Main parameters of the lidar point clouds
Terrain Slope/(°) Latitude and longitude Point density/pts·m−2 Low 0-5 (172°38′17″E, 40°42′37″S) 9.54 Moderate 5-15 (107°13′45″W, 37°59′33″N) 15.83 High 20-35 (115°68′50″W, 44°20′60″N) 13.71 表 3 累积光斑前后所对应的不同坡度区域的MAE与RMSE
Table 3. MAE and RMSE at the test areas with different slopes before and after accumulating spots
Type Terrain MAE/cm RMSE/cm No accumulation Low 12.07 17.21 Moderate 30.51 42.25 High 99.80 127.69 Accumulation Low 4.90 6.66 Moderate 10.95 13.64 High 35.43 43.78 -
[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