相干调频连续波激光雷达畸变补偿技术研究 下载: 936次
1 引言
无人驾驶是未来汽车领域的主要发展方向和趋势之一。激光雷达具有主动全天时、可三维(3D)成像、抗干扰能力强等优点,在无人驾驶汽车众多传感器中的地位越来越突出。调频连续波(FMCW)激光雷达是无线电线性调频探测与激光探测结合的产物,在探测虚警率、探测精度等方面具有明显的优势。车载激光雷达包括面成像激光雷达和点扫描激光雷达两种类型,目前市场上的主流产品均为点扫描激光雷达[1-5]。随着激光雷达发射和探测技术的提升,测距精度越来越高,一般均达到厘米级[6-11]。
采用点扫描成像方法,目标上不同的点探测的时间不同,整个时间差最长为帧频的倒数。对于静态或慢速运动场景,运动产生的畸变并不突出。但对于车载应用领域,汽车在行驶中存在非常大的速度,常用的方法是在激光雷达上安装加速度计测量激光雷达自身的运动速度,通过算法进行运动畸变的补偿[12-14]。1996年,Harbor等[15]将惯性测量单元与激光雷达结合,利用集成在惯性测量单元内的加速度计进行速度测量。2015年,Marlton等[16]使用加速度计对激光雷达回波数据进行校准。但是,通过积分求解速度,加速度计方法会带来累积误差。另外,加速度计方法只能解决激光雷达运动、目标静止带来的运动畸变,而无法解决目标运动产生的畸变。
三维相干调频连续波激光雷达是无人驾驶领域的新技术,可以解决远距离、人眼安全、高分辨率的应用需求。由于采用相干体制,可以实现目标距离和速度的同时测量,实现更多维目标信息的获取。
本文利用三维相干调频连续波激光雷达的特点,提出了一种不依赖于外部传感器进行测量的运动畸变逐点补偿方法,采用逐点补偿可以大幅提升激光雷达对运动速度范围的适应性。该方法不仅可以解决激光雷达运动产生的畸变,还可以解决目标运动产生的畸变,对于车载激光雷达的实际应用具有重要意义。
2 相干调频连续波激光雷达原理
相干调频连续波激光雷达的基本原理如
调频连续波激光器发射的光信号经1×2光纤分束器,一路由光纤环形器1端进入,2端出射,经光纤准直器发射出去,照射到目标表面的后向散射光进入准直器,进入光纤环形器2端,并由3端出射,进入2×2光纤耦合器,与1×2光纤分束器的本振光信号相干后由平衡探测器进行光电转换,再由数据采集卡进行信息采集,采集的信息由上位机进行数据处理。
调频连续波激光雷达的发射波形如
式中:Δfdopp为回波信号与发射信号间的多普勒频移;ΔfFMu和ΔfFMd分别为激光雷达上升沿和下降沿解线调频频率。具体公式表示为[17]
式中:c为光速。
图 2. 相干调频连续波激光雷达发射、本振、回波信号瞬时频率随时间变化关系图
Fig. 2. Relationship of instantaneous frequency of coherent frequency modulation continuous wave lidar emission, local oscillator, and echo signal with time
假设
式中:f为调频连续波激光器的发射激光的平均频率。
从(6)、(7)式可以看出,相干调频连续波激光雷达在测量目标与激光雷达相对距离的同时可以获取目标与激光雷达的相对速度,二者是同时获得的,利用这一点可以对成像过程中产生的运动畸变进行逐点补偿。
3 逐点运动畸变补偿理论
实际车载激光雷达工作模式如
按照
对于每个时刻,激光雷达成像系统都将获取该时刻的激光雷达与目标点距离L,出射光束的方位角θ,以及出射光束的俯仰角φ。在激光雷达运动坐标系x'ly'lz'lo'l下,对于t时刻扫描目标上某一位置点P在激光雷达运动坐标系x'ly'lz'lo'l上的坐标(x'P,y'P,z'P),有
则在0时刻,该点在激光雷达初始静止坐标系下的坐标为
式中:
所以有
式中:vlx、vly、vlz分别为vl在三个坐标轴方向的分量;vox、voy、voz分别为vo在三个坐标轴方向的分量。由(10)式可知,本文所述基于相干调频连续波激光雷达速度测量对畸变进行逐点补偿方法无需使用加速度计,其应用场景为激光雷达与目标仅在z轴方向存在相对运动。而对于激光雷达与目标的相对运动在x轴或y轴方向存在速度分量的情况,则需要借助其他算法进行补偿。
当激光雷达静止时,(10)式可简化为
此时,传统加速度计方法无法获取目标速度,进而也就无法进行畸变补偿。
当目标静止时,(11)式可简化为
从(12)式可以看出,此时激光雷达测量的速度即为车载平台的运动速度,这种情况下的补偿效果与采用传统加速度计测量方法相同。
对于某个时间段τ,获取到的空间点集可写为
而点集∑P就组成了最终的点云图像。
4 数值模拟仿真分析
本文仿真了激光雷达和目标只沿z轴方向运动的成像情形。
1)不进行速度补偿
假设激光雷达以速度vl=60 km/h沿z轴正方向匀速运动,靠近目标物体,目标物体静止,不进行速度补偿时,激光雷达点云成像如
由图可见,由于存在相对运动,不同时刻在激光雷达坐标系下采集到的点坐标与静态坐标系下的实际位置坐标存在差异,不进行速度补偿则会造成成像畸变,无法正常成像。
2)使用传统加速度计方法进行速度补偿
假设激光雷达以速度vl=60 km/h沿z轴正方向匀速运动,靠近目标物体,目标物体静止,激光雷达点云成像如
图 7. 激光雷达运动且目标物体静止点云图(加速度计方法)
Fig. 7. Point cloud image of lidar moving and target stationary (accelerometer method)
由
而在激光雷达以速度vl=60 km/h沿z轴负方向匀速运动时,目标物体以速度20 km/h沿z轴负方向与激光雷达相向运动时,激光雷达点云成像如
图 8. 激光雷达运动且目标物体运动点云图(加速度计方法)
Fig. 8. Point cloud image of lidar moving and target moving (accelerometer method)
由
3)使用本文所述方法进行速度补偿
假设激光雷达以速度vl=60 km/h沿z轴负方向匀速运动,同时,目标物体以速度为20 km/h沿z轴负方向与激光雷达相向运动。使用本文方法进行激光雷达点云成像的结果如
由
图 9. 激光雷达运动且目标物体运动点云图(本文方法)
Fig. 9. Point cloud image of lidar moving and target moving (proposed method)
使用补偿误差δ作为成像畸变补偿效果的评价指标,
式中:Δv为补偿后激光雷达与目标间的剩余速度差;f为扫描成像帧频。在扫描成像帧频不变的情况下,补偿后激光雷达与目标间的剩余速度差越小,补偿误差越小,精度越高。而使用加速度计方法补偿时,激光雷达与目标的相对速度仍有20 km/h,所以补偿后仍然出现一定程度的成像畸变;使用本文方法补偿后激光雷达与目标理论上不存在相对速度,所以精度较高。
5 实验验证
设计实验来验证激光雷达和目标只沿z轴方向运动的成像情况。基于本文方法进行速度补偿的三维相干调频连续波激光雷达成像技术可以有效地为激光雷达运动推扫成像提供准确的三维点云数据。如
表 1. 基于速度补偿的三维相干调频连续波激光雷达推扫成像实验参数
Table 1. Experimental parameters of 3D coherent FMCW lidar push-scan imaging based on velocity compensation
|
经本文方法速度补偿后的三维激光雷达点云图如
图 13. 使用本文方法进行速度补偿的三维激光雷达点云
Fig. 13. 3D lidar point cloud with velocity compensation by proposed method
忽略墙体表面起伏,将墙面近似为平面,则认为墙体成像产生的厚度即为速度补偿成像距离误差。测量成像墙体三维点云厚度峰峰值,计算得距离误差δ=±4.76 cm。激光雷达成像推扫时间t=50.95 s,速度误差δvelocity=0.19 cm/s,约占平均速度的0.52%。经分析,此误差的主要来源为对频率信号进行数字采样所带来的量化误差,所以提高精度的主要方法在于提高信号接收系统的采样率和量化位数。
6 结论
本文提出了一种基于相干调频连续波激光雷达速度测量对于运动产生畸变进行逐点补偿的方法,同时给出了车载三维相干调频连续波激光雷达成像点云坐标方程,以及该方程在不同情形下的讨论。分不同情形对激光雷达扫描成像进行仿真,最后进行了实际成像的实验室验证,获取了清晰的三维成像结果。仿真与实验结果表明,相干调频连续波激光雷达速度测量产生的相对误差为0.52 %,进行逐点补偿后,车载激光雷达点云图像的距离误差为±4.76 cm,与传统加速度计方法相比精度更高,普适性更强,并且不依赖于加速度计等外部测量,非常适合应用在车载或动目标检测等场合。在理论分析基础上,在实验室进行了实验验证,通过圆形扫描方式得到了三维目标图像,下一步将继续提升测距、测角精度,提升成像质量。
[1] ManandharD, ShibasakiR. Vehicle-borne laser mapping system (VLMS) for 3D GIS[C]∥IEEE International Geoscience & Remote Sensing Symposium. IEEE, 2001.
[2] 张立朔, 程效军. 基于激光点云的隧道形变分析方法[J]. 中国激光, 2018, 45(4): 0404004.
[3] Uchino O, Tabata I. Mobile lidar for simultaneous measurements of ozone, aerosols, and temperature in the stratosphere[J]. Applied Optics, 1991, 30(15): 2005-2012.
[4] Weibring P, Edner H, Svanberg S. Versatilemobile lidar system for environmental monitoring[J]. Applied Optics, 2003, 42(18): 3583-3594.
[5] Carmer D C, Peterson L M. Laser radar in robotics[J]. Proceedings of the IEEE, 1996, 84(2): 299-320.
[6] Hall D S. High definition lidar system: US7969558B2. [P].2008-01-17.
[7] PacalaA, Yu, TY, Eldada L. Cost-effective lidar sensor for multi-signal detection, weak signaldetection and signaldisambiguation and method of using same: US20140211194[P].2014-07-31.
[8] Matsubara H, Soga M, Niclass C, et al. Development of next generation LIDAR[J]. R&D Review of Toyota CRDL, 2012, 43(1): 7-12.
[9] 谢庚承, 叶一东, 李建民, 等. 脉冲激光测距回波特性及测距误差研究[J]. 中国激光, 2018, 45(6): 0610001.
[10] 冯长中, 吴松华, 刘秉义. 相干多普勒激光雷达风场反演方法研究与实验印证[J]. 中国激光, 2018, 45(4): 0410001.
[11] Mao XS, InoueD, MatsubaraH, et al., 2012, E95. B( 8): 2631- 2637.
[12] Moreira J R. A new method of aircraft motion error extraction from radar raw data for real time motion compensation[J]. IEEE Transactions on Geoscience and Remote Sensing, 1990, 28(4): 620-626.
[13] IlasC. Electronic sensing technologies for autonomous ground vehicles: a review[C]∥2013 8th International Symposium on Advanced Topics in Electrical Engineering (Atee). 23-25 May 2013, Bucharest, Romania. New York: IEEE Press, 2013: 1- 6.
[14] TeichmanA, LevinsonJ, ThrunS. Towards 3D object recognition via classification of arbitrary object tracks[C]∥2011 IEEE International Conference on Robotics and Automation. 9-13 May 2011, Shanghai, China.New York: IEEE Press, 2011: 4034- 4041.
[15] Harbor L K. Space-based interceptor hardware-in-the-loop simulation featuring IR scene projection and dual-sensor IR and lidar signal injection[J]. Proceedings of SPIE, 1996, 2741: 439-446.
[16] Marlton G J, Harrison R G, Nicoll K A, et al. Note: a balloon-borne accelerometer technique for measuring atmospheric turbulence[J]. The Review of Scientific Instruments, 2015, 86(1): 016109.
[17] 孙建锋, 周煜, 卢智勇, 等. 基于调频连续波的扫描激光雷达技术[ C]∥2017年光学技术研讨会暨交叉学科论坛论文集, 2017: 16- 21.
Sun JF, ZhouY, Lu ZY, et al. Scanning lidar technology based on FMCW[ C]∥2017 Optical Technology Symposium and Interdisciplinary Forum, 2017: 16- 21.
Article Outline
蔡新雨, 孙建锋, 卢智勇, 李跃新, 从海胜, 韩荣磊. 相干调频连续波激光雷达畸变补偿技术研究[J]. 中国激光, 2020, 47(9): 0910003. Cai Xinyu, Sun Jianfeng, Lu Zhiyong, Li Yuexin, Cong Haisheng, Han Ronglei. Distortion Compensation Technology of Coherent Frequency Modulation Continuous Wave Lidar[J]. Chinese Journal of Lasers, 2020, 47(9): 0910003.