留言板

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

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

基于全光纤光子计数激光雷达的高分辨率三维成像

郭静菁 费晓燕 葛鹏 周安然 王磊 李正琦 盛磊

郭静菁, 费晓燕, 葛鹏, 周安然, 王磊, 李正琦, 盛磊. 基于全光纤光子计数激光雷达的高分辨率三维成像[J]. 红外与激光工程, 2021, 50(7): 20210162. doi: 10.3788/IRLA20210162
引用本文: 郭静菁, 费晓燕, 葛鹏, 周安然, 王磊, 李正琦, 盛磊. 基于全光纤光子计数激光雷达的高分辨率三维成像[J]. 红外与激光工程, 2021, 50(7): 20210162. doi: 10.3788/IRLA20210162
Guo Jingjing, Fei Xiaoyan, Ge Peng, Zhou Anran, Wang Lei, Li Zhengqi, Sheng Lei. High-resolution three-dimensional imaging based on all-fiber photon-counting Lidar system[J]. Infrared and Laser Engineering, 2021, 50(7): 20210162. doi: 10.3788/IRLA20210162
Citation: Guo Jingjing, Fei Xiaoyan, Ge Peng, Zhou Anran, Wang Lei, Li Zhengqi, Sheng Lei. High-resolution three-dimensional imaging based on all-fiber photon-counting Lidar system[J]. Infrared and Laser Engineering, 2021, 50(7): 20210162. doi: 10.3788/IRLA20210162

基于全光纤光子计数激光雷达的高分辨率三维成像

doi: 10.3788/IRLA20210162
基金项目: 军队后勤开放研究科研项目(BY119C009);安徽省自然科学基金(2008085QF315)
详细信息
    作者简介:

    郭静菁,女,工程师,博士,主要从事单光子三维成像方面的研究

    通讯作者: 葛鹏,男,高级工程师,博士,主要从事单光子三维成像方面的研究。
  • 中图分类号: TN29

High-resolution three-dimensional imaging based on all-fiber photon-counting Lidar system

  • 摘要: 为了实现全天时、高分辨率三维成像探测,文中提出了一种工作波长为1064 nm的光子计数激光雷达系统。其中,光学系统采用全光纤结构,增强了系统稳定性。雷达通过整机扫描方式对远距离目标进行成像探测,扩大扫描视场,水平方向达到360°,俯仰方向达到±30°,同时避免了由摆镜扫描引起的几何畸变;结合亚像素扫描方法,提高空间分辨率。最后通过自适应噪声阈值的多距离重建算法实现了目标的三维重建。实验结果表明,此系统在白天成功实现了3.1 km以外目标的三维重建,重建图像目标特征清晰,距离精度为0.11 m,空间分辨率为0.11 m,超过光学系统衍射极限。
  • 图  1  光子计数激光雷达系统原理图

    Figure  1.  Block diagram of photon-counting Lidar system

    图  2  亚像素扫描方法

    Figure  2.  Sub-pixel scanning method

    图  3  摆镜反射原理

    Figure  3.  Beam reflection by a mirror

    图  4  摆镜引起的几何畸变

    Figure  4.  Geometry distortion caused by a swing mirror

    图  5  整机扫描转台设计

    Figure  5.  Design of the rotation platform

    图  6  塔底幅度像成像结果。(a)可见光相机拍摄的被测目标照片;(b)噪声常数b = 2.5 Eyb;(c)噪声常数根据噪声均值分段赋值;(d)噪声抑制结果

    Figure  6.  Intensity image of the bottom of the poly. (a) Photo of the target taken by a visible-band camera; (b) Noise threshold b = 2.5 Eyb; (c) Self-adaptive noise threshold; (d) Noise censoring results

    图  7  天线塔重建成像结果。 (a)幅度像;(b)距离像;(c)三维点云重建

    Figure  7.  Image reconstruction of the poly. (a) Intensity image; (b) Range image; (c) 3D point cloud reconstruction

    图  8  可见光相机拍摄的被测场景照片

    Figure  8.  Photo of the target scene taken by a visible-band camera

    图  9  多距离目标重建结果。(a)幅度像E1;(b)幅度像E2;(c)距离像L1;(d)距离像L2

    Figure  9.  Multi-range image reconstruction results. (a) Intensity image E1; (b) Intensity image E2; (c) Depth image L1; (d) Depth Estimation L2

    图  10  三维点云重建。(a) 多距离目标重建;(b) 3.1 km处建筑物局部放大图;(c) 1.1 km处电线的局部放大图

    Figure  10.  3D point cloud reconstruction. (a) Multi-range object reconstruction; (b) Enlarged image of the building 3.1 km away; (c) Enlarged image of the cables 1.1 km away

    表  1  光学系统主要器件参数

    Table  1.   Key parameters of the optical system

    DevicesParametersValue
    LaserWavelength/nm1064
    Pulse repetition frequency/kHz40
    Peak energy/μJ6
    Pulse width/ns5
    Emitting collimatorFocus length/mm80
    Receiving telescopeDiameter/mm20.5
    Focus length/mm40
    FiberFiber diameter/μm9/125
    FilterBand width/nm1
    SPADDark count/cps100
    Dead time/ns22
    Detection efficiency2%@1064 nm
    下载: 导出CSV
  • [1] McCarthy A, Collins R J, Krichel N J, et al. Long-range time-of-flight scanning sensor based on high-speed time-correlated single-photon counting [J]. Applied Optics, 2009, 48(32): 6241-6251. doi:  10.1364/AO.48.006241
    [2] Li Z H. Multi-beam photon-counting laser imaging [D]. Shanghai: East China Normal University, 2017. (in Chinese)
    [3] Ge Peng, Guo Jingjing, Chen Cong, et al. Photon-counting 3D imaging based on Geiger-mode APD array [J]. Infrared and Laser Engineering, 2020, 49(3): 0305007. (in Chinese) doi:  0305007
    [4] Du B C, Pang C K, Wu D, et al. High-speed photon-counting laser ranging for broad range of distances [J]. Scientific Reports, 2018, 8(4198): 10.1038/s41598-018-22675-1. doi:  https://doi.org/10.1038/s41598-018-22675-1
    [5] Mccarthy A, Krichel N, Gemmell N, et al. Kilometer-range, high resolution depth imaging via 1560 nm wavelength single-photon detection [J]. Optics Express, 2013, 21(7): 8904-8915. doi:  10.1364/OE.21.008904
    [6] Shangguan M J, Xia H Y, Wang C, et al. All-fiber upconversion high spectral resolution wind [J]. Optics Express, 2016, 24(17): 19322-19336. doi:  https://doi.org/10.1364/OE.24.019322
    [7] Clifton W E, Steeleb B, Nelsonb G, et al. Medium altitude airborne Geiger-mode mapping lidar system[C]//Laser Radar Technology and Applications XX; and Atmospheric Propagation XII. International Society for Optics and Photonics, 2015: 946506.
    [8] Degnan J J. Scanning, multibeam, single photon lidars for rapid, large scale, high resolution, topographic and bathymetric mapping [J]. Remote Sensing, 2016, 8(11): 958. doi:  10.3390/rs8110958
    [9] Marino R M, Davis W R. Jigsaw: A foliage-penetrating 3D imaging laser radar system [J]. Lincoln Laboratory Journal, 2005, 15(1): 23-36.
    [10] Dumanis D. Airborne Optical Systems Testbed(AOSTB)[R]. US: MIT Lincoln Laboratory, 2016.
    [11] Albota M A, Gurjar R, Mangognia A, et al. Contributed Review: Advanced three-dimensional laser radar imaging with the airborne optical systems testbed [J]. Review of Scientific Instruments, 2018, 89(10): 101502. doi:  10.1063/1.5017855
    [12] Bahr T, Smith P. Airborne Geiger-mode lidar for large-scale, high-resolution wide-area mapping[C]//GI Forum, 2016, 1: 85-93.
    [13] Yu A W, Krainak M A, Harding D J, et al. A 16-beam non-scanning swath mapping laser altimeter instrument[C]//Proc SPIE, 2013, 8599: 85990P.
    [14] Li M. Research on technologies of photon counting Lidar based on fiber optics [D]. Shanghai: Shanghai Institute of Technical Physics, Chinese Academy of Sciences, 2017. (in Chinese)
    [15] Li Z P, Huang X, Cao Y, et al. Single-photon computational 3D imaging at 45 km [J]. arXiv, 2019: 1904.10341.
    [16] Li Z P, Huang X, Jiang P Y, et al. Super-resolution single-photon imaging at 8.2 kilometers [J]. Optics Express, 2020, 28(3): 4076-4087. doi:  10.1364/OE.383456
    [17] Qiu J. Airborne polarization lidar with 1.5 μm single-photon detectors [D]. Hefei: University of Science and Technology of China, 2020. (in Chinese)
    [18] Shangguan M J. Laser remote sensing with 1.5 μm single photon detectors [D]. Hefei: University of Science and Technology of China, 2017. (in Chinese)
    [19] Shin D, Xu F, Wong F N, et al. Computational multi-depth single-photon imaging [J]. Optics Express, 2016, 24(3): 253588.
    [20] Tobin R, Halimi A, McCarthy A, et al. Three-dimensional single-photon imaging through obscurants [J]. Optics Express, 2019, 27(4): 4590-4611. doi:  10.1364/OE.27.004590
  • [1] 冀鹏州, 穆郁, 张晨钟, 孟军合, 赵侃, 徐大维.  用于亚公里级三维成像的全光相机主物镜设计 . 红外与激光工程, 2022, 51(3): 20210493-1-20210493-7. doi: 10.3788/IRLA20210493
    [2] 徐国权, 李广英, 万建伟, 许可, 董光焰, 程光华, 王兴, 韩文杰, 马燕新.  脉冲调制激光雷达水下成像系统 . 红外与激光工程, 2022, 51(3): 20210204-1-20210204-8. doi: 10.3788/IRLA20210204
    [3] 徐璐, 刘霞霞, 杨旭, 张一嘉, 吴龙.  宽脉冲光子计数偏振激光雷达探测浅水层研究(特邀) . 红外与激光工程, 2021, 50(3): 20200452-1-20200452-10. doi: 10.3788/IRLA20200452
    [4] 张晨阳, 王春辉, 战蓝, 齐明, 蒋硕.  星载微脉冲光子计数激光雷达的探测性能仿真分析 . 红外与激光工程, 2021, 50(11): 20200502-1-20200502-9. doi: 10.3788/IRLA20200502
    [5] 朱笑笑, 王成, 习晓环, 聂胜, 杨学博, 黎东.  ICESat-2星载光子计数激光雷达数据处理与应用研究进展 . 红外与激光工程, 2020, 49(11): 20200259-1-20200259-10. doi: 10.3788/IRLA20200259
    [6] 曹杰, 郝群, 张芳华, 徐辰宇, 程阳, 张佳利, 陶禹, 周栋, 张开宇.  APD三维成像激光雷达研究进展 . 红外与激光工程, 2020, 49(9): 20190549-1-20190549-10. doi: 10.3788/IRLA20190549
    [7] 杨程, 鄢秋荣, 祝志太, 王逸凡, 王明, 戴伟辉.  基于深度学习的压缩光子计数激光雷达 . 红外与激光工程, 2020, 49(S2): 20200380-20200380. doi: 10.3788/IRLA20200380
    [8] 康岩, 薛瑞凯, 李力飞, 张同意, 张勇, 赵卫.  基于像素复用的SPAD阵列连续扫描三维成像 . 红外与激光工程, 2020, 49(S2): 20200375-20200375. doi: 10.3788/IRLA20200375
    [9] 王帅, 孙华燕, 赵延仲, 曾海瑞, 刘田间.  基于光学相控阵的提高APD阵列三维成像分辨率方法 . 红外与激光工程, 2019, 48(4): 406003-0406003(8). doi: 10.3788/IRLA201948.0406003
    [10] 王海伟, 丁宇星, 黄庚华, 侯佳, 舒嵘.  轻小型全天时远程光子计数激光雷达系统技术 . 红外与激光工程, 2019, 48(1): 106005-0106005(7). doi: 10.3788/IRLA201948.0106005
    [11] 张河辉, 丁宇星, 黄庚华.  光子计数激光测深系统 . 红外与激光工程, 2019, 48(1): 106002-0106002(5). doi: 10.3788/IRLA201948.0106002
    [12] 李潇, 石柱, 代千, 覃文治, 寇先果, 袁鎏, 刘期斌, 黄海华.  64×64 InGaAs/InP三维成像激光焦平面探测器 . 红外与激光工程, 2018, 47(8): 806004-0806004(5). doi: 10.3788/IRLA201847.0806004
    [13] 李铭, 候佳, 周成林, 舒嵘.  多波束激光雷达的高精度收发配准方法的实验验证 . 红外与激光工程, 2017, 46(7): 730001-0730001(9). doi: 10.3788/IRLA201746.0730001
    [14] 赵远, 张子静, 马昆, 徐璐, 吕华, 苏建忠.  高灵敏度的光子偏振激光雷达系统 . 红外与激光工程, 2016, 45(9): 902001-0902001(6). doi: 10.3788/IRLA201645.0902001
    [15] 张月, 张琢, 苏云, 郑国宪.  宽谱段高分辨率低温成像光谱仪制冷系统设计 . 红外与激光工程, 2016, 45(3): 323001-0323001(8). doi: 10.3788/IRLA201645.0323001
    [16] 胡杨, 朱鹤元.  1.55μm全光纤相干多普勒激光测风雷达 . 红外与激光工程, 2016, 45(S1): 71-75. doi: 10.3788/IRLA201645.S130001
    [17] 唐禹, 秦宝, 晏芸, 汪路锋, 邢孟道.  多发多收合成孔径激光雷达高分辨率宽测绘带成像 . 红外与激光工程, 2016, 45(8): 830001-0830001(8). doi: 10.3788/IRLA201645.0830001
    [18] 张振振, 杨爱玲, 赵扬, 南钢洋.  人工缺陷铝块试样的激光超声三维成像 . 红外与激光工程, 2015, 44(S1): 57-62.
    [19] 李小珍, 吴玉峰, 郭亮, 曾晓东.  合成孔径激光雷达下视三维成像构型及算法 . 红外与激光工程, 2014, 43(10): 3276-3281.
    [20] 杨芳, 张鑫, 贺岩, 陈卫标.  采用高速伪随机码调制和光子计数技术的光纤激光测距系统 . 红外与激光工程, 2013, 42(12): 3234-3238.
  • 加载中
图(11) / 表(1)
计量
  • 文章访问数:  572
  • HTML全文浏览量:  141
  • PDF下载量:  84
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-03-13
  • 修回日期:  2021-04-08
  • 刊出日期:  2021-07-23

基于全光纤光子计数激光雷达的高分辨率三维成像

doi: 10.3788/IRLA20210162
    作者简介:

    郭静菁,女,工程师,博士,主要从事单光子三维成像方面的研究

    通讯作者: 葛鹏,男,高级工程师,博士,主要从事单光子三维成像方面的研究。
基金项目:  军队后勤开放研究科研项目(BY119C009);安徽省自然科学基金(2008085QF315)
  • 中图分类号: TN29

摘要: 为了实现全天时、高分辨率三维成像探测,文中提出了一种工作波长为1064 nm的光子计数激光雷达系统。其中,光学系统采用全光纤结构,增强了系统稳定性。雷达通过整机扫描方式对远距离目标进行成像探测,扩大扫描视场,水平方向达到360°,俯仰方向达到±30°,同时避免了由摆镜扫描引起的几何畸变;结合亚像素扫描方法,提高空间分辨率。最后通过自适应噪声阈值的多距离重建算法实现了目标的三维重建。实验结果表明,此系统在白天成功实现了3.1 km以外目标的三维重建,重建图像目标特征清晰,距离精度为0.11 m,空间分辨率为0.11 m,超过光学系统衍射极限。

English Abstract

    • 光子计数激光雷达系统(Lidar)可以实现光子量级的高灵敏度探测、μrad级的高空间分辨率,近年来受到广泛关注[1]。激光能量在远距离大气传输过程中衰减严重,回波信号甚至达到光子量级,传统探测方式难以识别。光子计数激光雷达发射高重频、窄脉宽的脉冲激光,利用单光子探测方式响应光子能量级别的微弱信号,通过计算激光飞行时间,获得几百甚至几十公里外的目标距离信息,进而重建三维图像[2-3]。目前,其应用已覆盖远距离成像[4-5]、大气探测[6]、地形测绘[7-8]等领域。

      在光子计数激光雷达远距离成像领域,林肯实验室自2001年开展了大量研究工作,先后研制了多型机载激光雷达原理样机[9],在2~9 km空域完成飞行试验,验证了全天时三维成像能力[10-11]。2015年,Harris公司推出了世界上第一款商用盖革模式激光雷达系统,探测距离达到6 km[12]。这些研究均是基于盖革模式单光子雪崩二极管(GM-APD)阵列,受到阵列规模、暗计数、视场角等性能制约,难以进一步提升视场。2011年,NASA验证了基于4×4光纤阵列的机载测绘激光雷达,工作波长为1030 nm,探测距离达到10 km[13]。国内在远距离成像方面的起步较晚,中国科学院上海技术物理研究所搭建的基于51波束光纤阵列的1064 nm光子计数激光雷达,可实现3 km以内的机载三维测绘[14]。华东师范大学搭建的1550 nm收发一体式光子计数激光雷达,可以实现21 km的测距[4]。中国科学技术大学搭建的1550 nm光子计数激光雷达系统,实现了45 km远距离成像[15]。在这些已报道的光子计数激光雷达系统中,收发光学系统大多采用空间光传输方式,或空间光与光纤组合传输方式,增加了系统不稳定性。此外,扫描方式大多为摆镜扫描。为实现大范围扫描,摆镜口径通常较大,增大了系统尺寸;另一方面,摆镜旋转会导致反射图像发生几何畸变,影响成像质量。

      文中设计了一种工作波长为1064 nm的全光纤光子计数激光成像雷达系统,通过整机扫描方式和亚像素扫描方法对3 km以上目标进行成像探测,空间分辨率超越衍射极限,并通过噪声阈值自适应的图像重建算法,重建多距离目标信息。该系统适用于远距离探测、地形测绘等领域。

    • 文中设计的光子计数激光雷达系统,通过发射脉冲激光,记录脉冲飞行时间来探测目标距离,通过二维扫描来实现大视场成像。系统工作原理见图1,激光器在内部电路驱动下产生高重频、窄脉宽的脉冲激光,经准直扩束后出射。发射激光被目标反射后,进入接收望远镜。最终回波信号进入硅基单光子探测器,完成光电转换。时间相关单光子计数(TCSPC)模块记录激光发射时刻t0和回波到达时刻t1,则目标距离R为:

      图  1  光子计数激光雷达系统原理图

      Figure 1.  Block diagram of photon-counting Lidar system

      $$R = {{\left( {{t_1} - {t_0}} \right)c} / 2}{\rm{ }}$$ (1)

      式中:c为光速。

      利用TCSPC技术对单光子信号进行探测,并且记录光子响应的时间,具有极高的时间分辨率。该技术将时间轴分为离散的时间区间,时间区间的大小决定了系统的精度。当探测器探测到一个或多个光子,产生一个响应输出,记录下该次响应发生的时间,并对该时间区间内的光子计数值加1。当进行大量重复脉冲探测后,可以统计得到响应时间对应光子计数数目的统计直方图,直方图峰值表征回波幅度,峰值所在的时刻表征回波到达时间。直方图波形存在着一定的波形宽度,对应的是一定的时间跨度,即距离误差,所以需要选取一定的算法来得到最准确的距离信息。

    • (1)静态视场角

      通常,静止状态下光纤光学系统的空间分辨率由发射视场角FOV决定,可近似表达为:

      $${\rm{FOV}} \approx \frac{{{\rm{MFD}}}}{f}$$ (2)

      式中:MFD为光纤模场直径;f为望远镜焦距。

      (2)亚像素扫描成像

      在远距离成像条件下,即使是很小的发散角也会在目标位置产生较大的光斑。为提高空间分辨率,文中采用如图2所示的亚像素扫描方法[16]。利用高精度转台,控制扫描角度步进小于FOV,对静态视场覆盖区域进行nn>1)次采样。整个场景扫描完成后,通过图像重建算法实现三维成像,分辨率为静态视场角的1/n

      图  2  亚像素扫描方法

      Figure 2.  Sub-pixel scanning method

    • 文中设计的光子计数激光雷达系统主要由发射光学系统、接收光学系统、时间相关单光子计数模块、扫描系统和上位机组成。如图1所示,发射和接收光学系统为全光纤结构。发射光学系统由光纤激光器、准直镜构成。接收光学系统由滤波器、接收望远镜和单光子探测器组成。伺服系统由二维转台、伺服控制系统、角度传感器、支架等构成。上位机负责控制激光器的脉宽、重频、泵浦功率、TCSPC采集时间和采集模式,以及二维转台的转速和步进,并记录角度传感器的实时角度,完成数据采集和图像重建。

    • 光子计数激光成像雷达的典型工作波长包括532 nm、1064 nm、1030 nm和1550 nm等。波长对系统性能起到关键作用,例如太阳在532 nm的辐射量是1064 nm的三倍以上,因此大多数光子计数激光成像雷达工作于近红外波段。综合考虑太阳背景辐射、单光子探测器的光谱响应效率和技术成熟度,文中系统采用工作波长1064.4 nm的脉冲光纤激光器,重频为40 kHz。激光器采用主振荡器加功率放大器(MOPA)结构设计,经过优化设计的多级光放大,实现峰值功率6 μJ,脉宽5 ns的激光输出,线宽约为0.04 nm。激光器输出尾纤为芯径9/125 μm的SMF-28e光纤。发射准直镜焦距80 mm。根据公式(2)计算可得,系统的静态视场角FOV约为131 μrad。

      接收光学系统中,设计接收视场角为发射光束发散角的两倍。接收望远镜焦距为40 mm,孔径为20.5 mm。为了降低背景辐射噪声,接收准直镜前安装了带宽1 nm的滤光片,通过螺纹固紧于接收望远镜的前端。接收信号直接耦合进SMF-28e光纤,进入探测器,完成光电转换。

      目前,能实现光子量级响应的探测器主要包括光电倍增管(PMT)、盖革模式单光子雪崩二极管(GM-APD)和超导探测器。其中,PMT对近红外响应较低,体积较大,不适用于光子计数激光雷达。盖革模式单光子雪崩二极管(SPADs)是一种概率器件,可以一定的概率对单个光子发生雪崩。现有SPAD单光子探测器的货架产品主要采用Si基或InGaAs探测器[17]。InGaAs探测器的死时间和暗计数都较大。Si基探测器较为成熟,响应波段主要为可见光,对1064 nm波长的探测效率为2%,对1550 nm波长几乎无响应。其他报道的提高探测效率的方法包括超导探测器和上转换技术等[18],但仍处于实验室研究阶段,尚不成熟。因此,文中采用Si基单光子探测器,暗计数为100 cps,死时间为22 ns。

      光学系统主要参数见表1

      表 1  光学系统主要器件参数

      Table 1.  Key parameters of the optical system

      DevicesParametersValue
      LaserWavelength/nm1064
      Pulse repetition frequency/kHz40
      Peak energy/μJ6
      Pulse width/ns5
      Emitting collimatorFocus length/mm80
      Receiving telescopeDiameter/mm20.5
      Focus length/mm40
      FiberFiber diameter/μm9/125
      FilterBand width/nm1
      SPADDark count/cps100
      Dead time/ns22
      Detection efficiency2%@1064 nm
    • 为扩大成像视场,现有光子计数激光雷达系统基本采用摆镜扫描方法:在出射/接收光路中加入反射镜,通过旋转反射镜,实现水平和俯仰方向的光束扫描。该方法存在以下两个问题:一是反射镜口径较大,二是反射图像存在几何畸变。反射镜示意图见图3,对于收发同轴光学系统,设收发望远镜直径为w1,反射镜与激光发射光路的夹角为θ;对于收发分置系统,w1为发射和接收望远镜的最大包络距离,其他与收发同轴光学系统相同。则反射镜直径w2应大于w1/sinθ,以避免边缘能量损失。因此,望远镜口径越大,扫描角度越大,需要的反射镜口径越大。

      图  3  摆镜反射原理

      Figure 3.  Beam reflection by a mirror

      此外,摆镜引起的几何畸变如图4所示,摆镜的两个相互垂直的旋转轴为r轴和p轴,p轴在xy平面内,目标在xz平面内。如图4(a)所示,当r轴与z轴平行,p轴与x轴呈45°时,目标经摆镜反射后,反射像上下不变,左右翻转,与目标镜像对称。当r轴保持与z轴平行,摆镜围绕r轴转动时,如图4(b)所示,与z轴平行的目标的反射像在z轴上的投影长度不变;但与x轴平行的目标反射后长度发生拉伸或收缩。当p轴保持与x轴呈45°时,摆镜围绕p轴转动,如图4(c)所示,与x轴平行的目标的反射像在y轴上的投影长度不变;与z轴平行的目标经反射镜反射后,在yz平面内发生旋转。

      图  4  摆镜引起的几何畸变

      Figure 4.  Geometry distortion caused by a swing mirror

      为了避免反射像的几何畸变,系统采用整机扫描方案,即光学系统发射的激光直接照射目标,通过将光学系统整体旋转,实现水平和俯仰方向的扫描。如图5所示,将一个一维转台水平放置,以实现水平方向扫描;将另一个一维转台垂直安装于水平转台上,以实现俯仰方向扫描。光学系统中,发射准直镜、接收望远镜、滤波片及其调整、夹持结构件安装于俯仰转台转动面的转接板上,激光器和单光子探测器安装于俯仰转台底座的转接板上。转台水平扫描范围为360°,俯仰扫描范围为±30°。俯仰和水平方向的角度分辨率均为0.002°(约等于35 μrad)。因此,扫描模式的空间分辨率理论上小于凝视模式的1/3。

      图  5  整机扫描转台设计

      Figure 5.  Design of the rotation platform

      伺服控制过程如下:第一步,开启激光器和TCSPC,将激光器内部时钟作为TCSPC同步时钟;第二步,转台开始转动,当按照设定角度步进值转动一定角度后,转台停止转动;第三步,TCSPC按照设定的采集时间采集回波光子数据;第四步,判断转动角度是否达到设定转动范围,若达到,则扫描结束,关闭激光器、TCSPC和转台,若未达到,则跳转至第二步,继续扫描。

    • 为了精确提取距离信息,系统采用多距离图像重建算法[19]。对于某一像素,在连续时间域上的光子探测率函数λ(t)可以表达为:

      $$\lambda \left( t \right) = \eta \left[ {h\left( t \right) * s\left( t \right)} \right] + b$$ (3)

      式中:*表示卷积;h(t)为目标的响应函数;s(t)为探测器对发射脉冲的归一化响应函数;η为量子效率;b为噪声常数,是背景辐射噪声和暗计数的总和。

      对于真实场景来说,若在一个脉冲周期T内接收到K个不同距离的回波信号,则在光子计数直方图中将出现多个峰值[20]。该像素点接收的目标的脉冲响应函数h(t)可以表达为:

      $$\begin{array}{*{20}{c}} {h\left( t \right) = \displaystyle\sum\limits_{i = 1}^K {{r_i}\delta \left( {t - {{2{d_i}} / c}} \right),} }&{t \in \left[ {0,T} \right)} \end{array}$$ (4)

      式中:δ(·)为冲激函数;ridi分别为第i个反射体的反射率和距离;t为时间。

      当探测器接收信号之后,以一定的时间分辨率对接收到的信号离散采样,经过N个脉冲周期的积累,光子计数直方图中第k个时间单元的光子数yk服从泊松分布,即

      $$\begin{array}{*{20}{c}} {y_k} \sim {\rm{Possion}}\left( {{{\left( {{{{{Sx}}}} + b{{{{{I}}}}_n}} \right)}_k}} \right),\;\;\;\;\;k = 1,2, \cdots ,n \end{array}$$ (5)

      式中:x表示连续函数h(t)采样后的离散形式,是一个n×1向量;Sx表示s(t)离散采样后与x的卷积,S为一个n×n矩阵;b为常数;In为常数1组成的n×1向量。

    • 目标图像重建问题即从光子直方图数据y求解目标响应x的问题。为了方便计算,可以将y的概率分布函数写成关于x的负对数似然函数L(x:y,S, b)。因此,多距离图像重建问题可以转化成在x的非零元素等于K的条件下,最小化负对数似然函数的问题:

      $$ \mathop {\mathop {\min }\limits_{{x}} }\limits_{} L\left( {{{x}}:{{y}},{{S}},b} \right) + \tau {\left\| {{x}} \right\|_1} \\\begin{array}{*{20}{c}} {{{{x}}_k} \geqslant 0,}&{k = 1,2, \cdots ,n} \end{array} $$ (6)

      式中:τ > 0;L(x:y,S,b)和l1范数为关于x的凸函数,所以,最小化负对数似然函数也是一个凸优化问题。公式(6)中,各输入变量如下:

      (1) y为测量目标物获得的光子直方图数据。

      (2) S为先通过测量获得探测器对发射脉冲的响应直方图,再根据直方图拟合成归一化响应函数s(t),最终建立的S矩阵。文中拟合的归一化响应函数为:

      $$s\left( t \right) = \exp \left( { - \frac{{{x^2}}}{{{\mu ^2}}}} \right)$$ (7)

      式中:μ为发射激光脉宽。

      (3) b为根据y计算的噪声常数。对于全天时、大幅宽的单光子成像场景,背景噪声会随着扫描时长的增加而变化。定义系统盲区距离Lb内的光子直方图数据矩阵yb为:

      $${{{y}}_{{b}}} = \left\{ {{y_1},{y_2}, \cdots ,{y_m}} \right\},m = \frac{{2{L_b}}}{{\Delta c}}$$ (8)

      在评价噪声水平时,需要计算系统盲区接收光子数的平均值Eyb,用于表征噪声均值;计算系统盲区接收光子数方差σyb,用于表征噪声波动情况。实际上,由于接收噪声光子数也服从泊松分布,Eyb近似等于σyb,为了节约计算资源,可以只计算Eyb。根据噪声均值的高低,文中对噪声常数b采取分段赋值的方法:

      $$b = \left\{ {\begin{array}{*{20}{c}} {3.5{E_{{{y}}_b}}{,_{}}\begin{array}{*{20}{c}} {{E_{{{y}}_b}} \leqslant {{0.1}_{}}{{_{}}_{}}}&{} \end{array}} \\ {2.5{E_{{{y}}_b}}{,_{}}0.1 < {E_{{{y}}_b}} \leqslant 0.2} \\ {1.5{E_{{{y}}_b}}{,_{}}0.2 < {E_{{{y}}_b}} \leqslant 0.3} \end{array}} \right.$$ (9)

      (4) τ = 3b

      利用增加非零限制的迭代收缩阈值算法[19],可以求解该问题。文中对单个像素最多提取2个距离信息,即K = 2。对所有M×N个像素运行迭代收缩阈值算法,若存在目标,则将最大响应幅度存入原始幅度像矩阵E1,对应的目标距离存入原始距离像矩阵L1;若存在2个目标,则将次级响应幅度存入原始幅度像矩阵E2,对应的目标距离存入原始距离像矩阵L2

      在获得原始距离像L1L2和原始幅度像E1E2后,文中根据相邻像素的关联性进一步剔除图像噪声,算法流程如下:

      (1) 将所有E1中的非零像素按照幅值从小到大排列,定义排在前1/q的像素的最大幅值为噪声幅度阈值;将E1中所有小于噪声幅度阈值的像素赋值为0;

      (2) 以E1中任一非零像素为中心,若至少有7个相邻像素为0,则将该像素赋值为0,

      (3) 若E1中某个像素的值为0,则将E2L1L2中对应像素赋值为0;

      (4) 输出优化后的Lopt1Lopt2、Eopt1Eopt2

      根据优化后的距离像Lopt1Lopt2和转台扫描的相对角度信息,即可重建激光三维点云。

    • 实验采用逐行扫描模式,某一行水平扫描结束后,俯仰方向步进0.002°,水平方向回到初始角度,开始下一行扫描。扫描按照从左到右,从上到下的顺序进行,每个像素的采集时间为100 ms。整个过程由角度传感器记录实时角度信息。

      实验被测目标为2.1 km以外某山顶天线塔和3.1 km以外某大楼。扫描天线塔当日天气为晴天,能见度8 km;扫描大楼当日天气为晴天,能见度14 km。实验装置安装在实验室内。

      图6给出了天线塔塔底部分的幅度像成像结果。图6(a)为可见光相机拍摄的目标照片。作为对比,首先将噪声常数b设为噪声均值的固定倍数,即b = 2.5Eyb,原始幅度像重建结果见图6(b)。分析每个像素的光子直方图数据可得,从第一行像素至最后一行,噪声均值从0.28逐渐降低至0.09。这是因为随着扫描时间的增加,太阳逐渐落山,导致背景辐射强度逐渐降低。图6(b)中,在绝大部分没有目标存在的天空区域,噪声依然作为目标被提取了出来,在空间上呈连续分布,塔身的爬架几乎湮没在噪声中,图像质量较差。图6(c)为采用文中的噪声自适应方法获得的原始幅度像,天空区域的噪声大幅下降,目标和噪声的幅度对比度增大。进一步剔除图像噪声,设置幅度阈值为1,大部分噪点被剔除,但尚有少数离散的噪点存在。根据相邻像素的关联性,可以将离散分布的噪点剔除,最终获得的幅度像如图6(d)所示,可以较为清晰地辨认包括攀爬架在内的天线塔结构,噪点抑制效果明显。

      图  6  塔底幅度像成像结果。(a)可见光相机拍摄的被测目标照片;(b)噪声常数b = 2.5 Eyb;(c)噪声常数根据噪声均值分段赋值;(d)噪声抑制结果

      Figure 6.  Intensity image of the bottom of the poly. (a) Photo of the target taken by a visible-band camera; (b) Noise threshold b = 2.5 Eyb; (c) Self-adaptive noise threshold; (d) Noise censoring results

      消除图像噪点后的天线塔重建图像见图7,天线塔框架结构清晰,背景噪声抑制效果较好。该实验证明了文中提出的噪声常数自适应的图像重建算法可以有效提取目标特征,抑制背景噪声。

      图  7  天线塔重建成像结果。 (a)幅度像;(b)距离像;(c)三维点云重建

      Figure 7.  Image reconstruction of the poly. (a) Intensity image; (b) Range image; (c) 3D point cloud reconstruction

      为了进一步验证多距离图像的重建效果,笔者对存在前后遮挡关系的目标进行了成像实验。如图8所示,视场范围内1.1 km处的电线对3.1 km处的某建筑存在部分遮挡。

      图  8  可见光相机拍摄的被测场景照片

      Figure 8.  Photo of the target scene taken by a visible-band camera

      重建结果见图9E1为提取最大幅值信号重建的图像,对应的目标像为L1E2为提取次级幅值信号重建的图像,对应的距离像为L2。点云重建结果见图10。从重建结果可以看出,位于前景的电线和被电线遮挡的建筑物能够同时被提取出来,实验结果验证了文中的图像重建算法可以重建多距离目标。

      图  9  多距离目标重建结果。(a)幅度像E1;(b)幅度像E2;(c)距离像L1;(d)距离像L2

      Figure 9.  Multi-range image reconstruction results. (a) Intensity image E1; (b) Intensity image E2; (c) Depth image L1; (d) Depth Estimation L2

      进一步对建筑物进行5次重复测距,每次累计时间均为100 ms。5次测距结果为3 135.05 m、3 134.82 m、 3134.90 m、3 134.98 m和3 135.12 m。根据均方根误差公式

      $$S = \sqrt {\dfrac{{\displaystyle\sum\nolimits_{i = 1}^n {{{\left( {{R_i} - \overline R} \right)}^2}} }}{n}} $$ (11)

      式中:n为测量次数;Ri为每次距离测量值;$\overline R$为测量平均值。计算可得5次测量的均方根误差约为0.11 m。因此,文中系统的测距精度为0.11 m。

      图  10  三维点云重建。(a) 多距离目标重建;(b) 3.1 km处建筑物局部放大图;(c) 1.1 km处电线的局部放大图

      Figure 10.  3D point cloud reconstruction. (a) Multi-range object reconstruction; (b) Enlarged image of the building 3.1 km away; (c) Enlarged image of the cables 1.1 km away

    • 文中提出了一种1064 nm全光纤光子计数激光雷达系统,介绍了系统工作原理和硬件构成。然后研究了光子计数图像重建算法,给出了自适应噪声阈值的计算方法。实验结果证明:系统在白天成功重建了3.1 km以外的目标,空间分辨率为0.11 m,超越衍射极限,成功提取多距离目标信息,测距精度为0.11 m,最终重建的图像特征清晰,噪声抑制效果较好。该系统适用于地面固定式或车载式激光雷达,用于远距离目标探测。

      但是,文中搭建的系统尚存在扫描成像时间较长、视场较小、探测距离较短等问题。在接下来的工作中,笔者将继续研究基于光纤阵列的光子计数激光雷达系统,进一步缩短成像时间、扩大静态视场范围和提高探测距离,推动该项技术在机载测绘领域的发展。

参考文献 (20)

目录

    /

    返回文章
    返回