CN103064076B - 光子计数三维成像激光雷达距离行走误差校正系统及方法 - Google Patents

光子计数三维成像激光雷达距离行走误差校正系统及方法 Download PDF

Info

Publication number
CN103064076B
CN103064076B CN201210574661.6A CN201210574661A CN103064076B CN 103064076 B CN103064076 B CN 103064076B CN 201210574661 A CN201210574661 A CN 201210574661A CN 103064076 B CN103064076 B CN 103064076B
Authority
CN
China
Prior art keywords
photon
laser
pulse
error
distance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201210574661.6A
Other languages
English (en)
Other versions
CN103064076A (zh
Inventor
何伟基
陈钱
顾国华
司马博羽
陈云飞
张闻文
钱惟贤
隋修宝
于雪莲
路东明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201210574661.6A priority Critical patent/CN103064076B/zh
Publication of CN103064076A publication Critical patent/CN103064076A/zh
Application granted granted Critical
Publication of CN103064076B publication Critical patent/CN103064076B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Optical Radar Systems And Details Thereof (AREA)

Abstract

本发明公开了一种光子计数三维成像激光雷达距离行走误差校正系统及方法,首先对光子计数三维成像激光雷达的距离行走误差进行多点标定,得到不同激光脉冲响应率下的距离行走误差;然后采用数学拟合的方法,获取距离行走误差关于激光脉冲响应率的函数,记为距离行走误差函数;最后,在光子计数三维成像激光雷达的工作过程中,统计原始三维距离图像的激光脉冲响应率分布,采用距离行走误差函数对距离行走误差进行预测,对原始三维距离图像应用补偿进行校正,得到校正后的三维距离图像。本发明能够消除距离行走误差。

Description

光子计数三维成像激光雷达距离行走误差校正系统及方法
技术领域
本发明属于光学、激光雷达和图像处理技术,特别是一种光子计数三维成像激光雷达距离行走误差校正系统及方法。
背景技术
距离图采用三维影像模式描述感兴趣区域,数据立方中包含了目标的几何不变特性,可以避免二维图像中的扭曲和混淆,广泛应用于工业模具设计、军事目标自动识别等领域。目前已知有若干技术途径可以用于获取距离图像。其中,光子计数三维激光成像雷达采用具有单光子灵敏度、工作于盖革模式下的雪崩二极管作为激光回波信号光子探测器,以光子飞行模式实现高时间分辨率测量。对目标进行多脉冲重复测量,再利用统计原理、光子计数技术得到目标单一像素的距离信息,通过扫描整个视场,得到目标的完整三维距离信息。
光子计数三维激光成像雷达的测量精度是最为重要的参数。精度定义为测量结果与实际距离间的差异,是探测器响应特性的随工作时间、工作状态和外部输入条件等因素改变而产生测量漂移,也称之为距离行走误差,以光子飞行时间测量误差或测量长度误差单位表示。为了克服距离行走误差,已有两种技术途径。(1)假定造成距离行走误差的原因是GmAPD输出电流脉冲上升速率因激光脉冲回波光子数多寡存在差异,提出一种采用多级阈值比较电路的校正方法[G.Kirchner,F.Koidl,et al,Proc.SPIE.3218,106-112(1997).]。(2)认为造成距离行走误差是激光脉冲回波光子数差异造成的GmAPD探测概率变化,进而提出了基于GmAPD探测概率模型的查表式校正方法[Min Seok Oh,Hong Jin Kong,Tae HoonKim,Keun Ho Hong,Byung Wook Kim,Opt.Commun.283,304-308(2010).]。
以上两种距离行走误差校正方法,能够在一定程度上校正误差,但也存在着很多不足之处。多级阈值比较电路过于复杂,难以集成到单片式读出电路中,因而应用价值有限。回波光子数差异改变探测器探测概率,难以经过实验验证,且查表式校正方法局限于单一实验系统,不具有普遍性。
在我们的研究中发现,造成距离行走误差的因素是激光脉冲回波光子数差异导致的原电子产生平均时间和GmAPD电流上升速率的变化。有必要采取更完善的校正方法来进行测量结果的校正。
发明内容
本发明的目的在于提供一种光子计数三维成像激光雷达距离行走误差校正系统及方法,能够消除距离行走误差,避免了测量误差对距离信息的干扰,提高成像质量。
实现本发明目的的技术解决方案为:一种光子计数三维成像激光雷达距离行走误差校正系统,包括光学系统、单光子探测器、脉冲激光光源、时间相关单光子计数模块和数据处理系统控制模块,光学系统由望远镜镜头、透镜、Y轴扫描转镜、X轴扫描转镜、准直透镜组、四分之一波片、第二半波片、分光棱镜、第一半波片、滤波片、光纤耦合器组成,脉冲激光光源同时发出同步起始信号和脉冲周期激光,同步起始信号输入数据处理系统控制模块,脉冲周期激光依次经过第一半波片、分光棱镜部分反射后在经过第二半波片、四分之一波片、Y轴扫描转镜、准直透镜组、X轴扫描转镜、透镜和望远镜镜头,出射到目标,经过目标漫反射被望远镜镜头接收,依次经过透镜、X轴扫描转镜、准直透镜组、Y轴扫描转镜、四分之一波片、第二半波片,在分光棱镜处部分透射,经过滤波片,达光纤耦合器,形成回波信号,被单光子探测器所探测接收,经过时间相关单光子计数模块产生截止信号输入数据处理系统控制模块;
脉冲激光光源发出周期脉冲激光,同时给予数据处理系统控制模块同步起始信号,经过光学系统后,激光反射信号到达单光子探测器,激励探测器形成响应脉冲,被时间相关单光子计数模块和数据处理系统控制模块接收,结合同步起始信号,即截止信号与同步起始信号的时间差,该时间差乘以光速,形成目标距离信息;通过固定X、Y轴扫描转镜,在数据处理系统控制模块实现单像素的先验标定;通过对X、Y轴扫描转镜的控制,控制出射激光到达目标的不同位置,从而达到对整个视场,即目标的整个平面进行二维扫描,得到每个像素的距离信息及误差信息,合成三维距离图像,对该三维距离图像在数据处理系统控制模块进行实时校正。
本发明与现有技术相比,其显著优点:(1)同时考虑GmAPD上升电流和探测率的影响,更符合实际情况,更准确,能够消除距离行走误差;(2)测量的同时进行误差校正,实时性好;(3)实施容易,无需额外增加硬件系统复杂性;(4)避免了测量误差对距离信息的干扰,提高成像质量。
下面结合附图对本发明作进一步详细描述。
附图说明
图1是本发明光子计数三维成像激光雷达距离行走误差校正系统的示意图。
图2是本发明采用的光子计数三维激光成像雷达的先验标定实验测量结果及数值拟合结果。
图3是本发明适用于光子计数三维成像激光雷达的先验标定距离行走误差校正方法的实时校正方法流程示意图。
图4是本发明实施效果验证中所采用的目标及实验场景实物图。
图5是本发明实施效果验证中光子计数三维成像激光雷达输出的原始距离图像。
图6是本发明实施效果验证中光子计数三维成像激光雷达输出的原始距离图像对应的响应率图像。
图7是本发明实施效果验证中光子计数三维成像激光雷达输出的原始距离图像对应的距离行走误差补偿图像。
图8是本发明实施效果验证中光子计数三维成像激光雷达输出的原始距离图像对应的实时校正图像。
具体实施方式
本发明适用于光子计数三维成像激光雷达的先验标定距离行走误差校正方法适用装置为光子计数三维成像激光雷达。首先通过多组实验,实验中改变出射激光的能量,从而得到不同的回波光子数,定义回波光子数与出射脉冲总数的比值为激光脉冲响应率。每次实验可测得相应脉冲响应率下的距离行走误差。然后采用数学拟合的方法,获取距离行走误差关于激光脉冲响应率的函数,记为距离行走误差函数。通过此次标定,就能得到该实验系统的距离行走误差校正函数,即校正方法。实验测量中,得到含有距离行走误差的原始三维距离图像,同时可以得到三维距离图像中每个像素点的激光脉冲响应率,根据此脉冲响应率,结合误差行走函数,能够得到该像素点的距离行走误差,然后对其进行补偿校正,就能够得到更接近真实目标的三维距离图像。
结合图1,光子计数三维成像激光雷达距离行走误差校正(实验)系统,包括光学系统、单光子探测器、脉冲激光光源、时间相关单光子计数模块(市购产品,如PicoHarp 300)和数据处理系统控制模块,光学系统由望远镜镜头1、透镜2,Y轴扫描转镜3、X轴扫描转镜4、准直透镜组5、四分之一波片6、第二半波片7、分光棱镜8、第一半波片9、滤波片10、光纤耦合器11组成,脉冲激光光源同时发出同步起始信号和脉冲周期激光,同步起始信号输入数据处理系统控制模块,脉冲周期激光依次经过第一半波片9、分光棱镜8部分反射后再经过第二半波片7、四分之一波片6、Y轴扫描转镜3、准直透镜组5、X轴扫描转镜4、透镜2和望远镜镜头1,出射到目标,经过目标漫反射被望远镜镜头1接收,依次经过透镜2、X轴扫描转镜4、准直透镜组5、Y轴扫描转镜3、四分之一波片6、第二半波片7,在分光棱镜8处部分透射,经过滤波片10,达光纤耦合器11,形成回波信号,被单光子探测器所探测接收,经过时间相关单光子计数模块产生截止信号输入数据处理系统控制模块;
脉冲激光光源发出周期脉冲激光,同时给予数据处理系统控制模块同步起始信号,经过光学系统后,激光反射信号到达单光子探测器,激励探测器形成响应脉冲,被时间相关单光子计数模块和数据处理系统控制模块接收,结合同步起始信号,即截止信号与同步起始信号的时间差,该时间差乘以光速,形成目标距离信息;通过固定X、Y轴扫描转镜,在数据处理系统控制模块实现单像素的先验标定;通过对X、Y轴扫描转镜的控制,控制出射激光到达目标的不同位置,从而达到对整个视场,即目标的整个平面进行二维扫描,得到每个像素的距离信息及误差信息,合成三维距离图像,对该三维距离图像在数据处理系统控制模块进行实时校正。
上述单像素数据处理系统控制模块的先验标定过程和实时校正过程同方法中步骤。
结合图3,本发明利用上述的光子计数三维成像激光雷达距离行走误差校正系统实现校正方法,包括先验标定和实时校正方法,其中先验标定方法包括以下步骤:
(1)事先准备高反射率平板一块作为目标,放置于光子计数三维成像激光雷达的望远镜镜头1正前方固定距离处;
(2)使得光子计数三维成像激光雷达的X轴扫描转镜4、Y轴扫描转镜3处于静止状态,设定脉冲激光光源输出脉冲能量为最大,激光脉冲能量记为E0
(3)使得脉冲激光光源处于周期工作状态,连续输出激光脉冲,激光脉冲重复次数均设定为Ntotal;同时,采用时间相关单光子计数模块分别记录每个激光脉冲中光子飞行时间,记为ttof(0,j)(1≤j≤Ntotal);
(4)计算激光脉冲能量E0对应的激光脉冲响应率为:R(0)=Npusle(0)/Ntotal,其中,表示Npusle(0)在测量持续时间内,探测器响应的脉冲总数;Ntotal表示在在测量持续时间内,激光脉冲的总数;
(5)计算激光脉冲能量E0对应的光子飞行时间测量均值:
t me ‾ ( 0 ) = 1 N total Σ j = 1 N total t tof ( 0 , j ) ;
(6)计算距离行走误差其中,ttof为恒定不变量,是光子在光子计数三维激光成像雷达与目标间光子的飞行时间;
(7)逐次减小脉冲激光光源的激光脉冲输出能量,重复步骤(2)至(6),得到不同激光脉冲响应率对应的距离行走误差,记为terror(R(i))=f(R(i)),0≤i<n;逐次减小脉冲激光光源的激光脉冲输出能量值为0.1E0至0.15E0之间固定值。
(8)采用数值拟合(计算机处理,为常规技术)方式,得到距离行走误差函数,记为terror(R)=f(R);
其中实时校正方法包括以下步骤:
1)通过对X、Y轴扫描转镜的控制,控制出射激光到达目标的不同位置,达到对目标不同像素点的探测,统计像素位置为(i,j)的光子飞行时间测量值,得到该位置对应的光子飞行时间测量均值
Figure BDA00002655493300053
2)统计像素位置为(i,j)的光子脉冲响应率,得到该位置对应的激光脉冲响应率R(i,j),即:
R(i,j)=Npulse(i,j)/Ntotal(i,j)
其中,Npulse(i,j)表示在像素位置(i,j)处的探测器输出脉冲总数;Ntotal(i,j)表示在像素位置(i,j)处的激光脉冲总数;
3)计算像素位置为(i,j)的距离行走误差,得到距离补偿图像terror(i,j),即:terror(i,j)=f(R(i,j));
4)补偿像素位置为(i,j)的光子飞行时间测量均值并得到校正后的光子飞行时间测量值,即:
t corect ( i , j ) = t me ‾ ( i , j ) + t error ( i , j )
5)统计每个像素点校正后的光子飞行时间测量值,合成目标的三维距离信息图及形状特征的三维距离图像。
实施例
本发明光子计数三维成像激光雷达距离行走误差校正方法中,先验标定方法包括以下步骤:
(1)事先准备高反射率(反射率大于90%)平板一块,放置于光子计数三维成像激光雷达的正前方的10米至20米之间固定距离处;
(2)使得光子计数三维成像激光雷达的扫描机构处于静止状态,设定激光器输出脉冲能量为最大,激光脉冲能量记为E0
(3)使得激光器处于周期工作状态,连续输出多个激光脉冲,激光脉冲重复次数均设定为Ntotal;同时,采用时间相关单光子计数模块分别记录每个激光脉冲中光子飞行时间,记为ttop(0,j)(1≤j≤Ntotal);
(4)计算激光脉冲能量E0对应的激光脉冲响应率为:R(0)=Npusle(0)/Ntotal。其中,表示Npusle(0)在测量持续时间内,探测器响应的脉冲总数;Ntotal表示在在测量持续时间内,激光脉冲的总数;
(5)计算激光脉冲能量E0对应的光子飞行时间测量均值:
t me ‾ ( 0 ) = 1 N total Σ j = 1 N total t tof ( 0 , j ) ;
(6)计算距离行走误差
Figure BDA00002655493300072
其中,ttof为恒定不变量,是光子在光子计数三维激光成像雷达与目标间光子的飞行时间;
(7)逐次减小激光器的激光脉冲输出能量,每次减小的能量值为0.1E0至0.15E0之间固定值,重复步骤(2)至(6),得到不同激光脉冲响应率对应的距离行走误差,记为terror(R(i))=f(R(i)),0≤i<n;
(8)采用最符合的数值拟合方式,得到距离行走误差函数,记为terror(R)=f(R)。
结合图2,图中离散点为实验测得的不同激光脉冲响应率下的距离行走误差,曲线为采用y=axb(a、b为参数)对实验数据进行的拟合。采用terror(R)=aRb对步骤(1)至(8)得到的测试结果进行拟合。采用的时间单位为ns;R为无量纲值,数值变化范围为0至1,得到a=-0.5659,b=3.06,即距离行走误差与脉冲响应率的函数关系表示为:
terror(R)=-0.5659R(3.06)
图3为方法流程图,图4为实验图,本发明实时校正方法包括以下步骤:
(1)通过实验测得单个脉冲,统计像素位置为(i,j)的光子飞行时间测量值,再经过大量脉冲(激光脉冲频率为2.5MHz,积分时间为100ms),取所有脉冲测得该像素位置的光子飞行时间测量值的平均值为该位置对应的光子飞行时间测量均值
Figure BDA00002655493300073
对视场进行三维扫描,可以测得所有像素点的光子飞行时间测量均值,从而得到光子计数三维成像激光雷达的原始距离图像如图5所示;
(2)同时,在步骤(1)中,统计像素位置为(i,j)的测量均值时,可以同时测得该位置的光子脉冲响应率,得到该位置对应的激光脉冲响应率(探测器输出脉冲总数/激光脉冲总数),即:
R(i,j)=Npulse(i,j)/Ntotal(i,j)
其中,Npulse(i,j)表示在像素位置处(i,j)的探测器输出脉冲总数;Ntotal(i,j)表示在像素位置处(i,j)的激光脉冲总数。
得到光子计数三维成像激光雷达输出的原始距离图像对应的激光响应率图像如图6所示;
(3)根据像素点(i,j)的激光响应率,通过之前标定得到的距离行走误差函数,可以计算得到像素位置为(i,j)的距离行走误差,即距离补偿值。得到距离补偿图像,即:
terror(i,j)=f(R(i,j))
其中terror(i,j)为距离补偿值,f为距离行走误差函数,f(x)=-0.5659x(3.06),R(i,j)为激光响应率。
得到光子计数三维成像激光雷达输出的原始距离图像对应的距离行走误差补偿图像如图7所示;
(4)通过原始距离测量均值,再减去距离补偿值,可以得到补偿后的校正值。像素位置为(i,j)的光子飞行时间测量均值经过补偿后得到校正后的光子飞行时间测量值,即:
t corect ( i , j ) = t me ‾ ( i , j ) + t error ( i , j )
其中,tcorect(i,j)为校正后光子飞行时间值,
Figure BDA00002655493300082
为原始光子飞行时间测量均值,terror(i,j)为距离行走误差补偿值。
(5)根据每个像素校正后的光子飞行时间测量值,可以得到目标到光子计数三维成像激光雷达输出的原始距离图像对应的实时校正图像,即目标的三维形状特征,如图8所示。
比较三维成像激光雷达输出的原始距离图像(图5)和校正后的三维距离图像(图8)。图5中红色标出部分,由于背景纸箱贴有胶带,使得这部分的反射率明显高于其他部分,使得三维成像距离图像中此部分存在距离行走误差,与目标三维特征不符合。利用本发明一种适用于光子计数三维成像激光雷达的先验标定距离行走误差校正算法我们得到了校正后的三维成像距离图像(图8),可以从中看出,红色标出部分的距离行走误差被校正去除,得到的三维距离图像与实际目标三维特征更加吻合。
通过上述具体实施例的验证,本发明适用于光子计数三维成像激光雷达的先验标定距离行走误差校正算法可以消除目标(背景)应反射能量差异造成的距离行走误差,避免了距离行走误差的干扰。结果表明,在不增加系统复杂度的前提下,因目标(背景)应反射能量差异造成的距离行走误差是可以有效避免的。

Claims (5)

1.一种光子计数三维成像激光雷达距离行走误差校正系统,其特征在于包括光学系统、单光子探测器、脉冲激光光源、时间相关单光子计数模块和数据处理系统控制模块,光学系统由望远镜镜头(1)、透镜(2)、Y轴扫描转镜(3)、X轴扫描转镜(4)、准直透镜组(5)、四分之一波片(6)、第二半波片(7)、分光棱镜(8)、第一半波片(9)、滤波片(10)、光纤耦合器(11)组成,脉冲激光光源同时发出同步起始信号和脉冲周期激光,同步起始信号输入数据处理系统控制模块,脉冲周期激光依次经过第一半波片(9)、分光棱镜(8)部分反射后在经过第二半波片(7)、四分之一波片(6)、Y轴扫描转镜(3)、准直透镜组(5)、X轴扫描转镜(4)、透镜(2)和望远镜镜头(1),出射到目标,经过目标漫反射被望远镜镜头(1)接收,依次经过透镜(2)、X轴扫描转镜(4)、准直透镜组(5)、Y轴扫描转镜(3)、四分之一波片(6)、第二半波片(7),在分光棱镜(8)处部分透射,经过滤波片(10),达光纤耦合器(11),形成回波信号,被单光子探测器所探测接收,经过时间相关单光子计数模块产生截止信号输入数据处理系统控制模块; 
脉冲激光光源发出周期脉冲激光,同时给予数据处理系统控制模块同步起始信号,经过光学系统后,激光反射信号到达单光子探测器,激励探测器形成响应脉冲,被时间相关单光子计数模块和数据处理系统控制模块接收,结合同步起始信号,即截止信号与同步起始信号的时间差,该时间差乘以光速,形成目标距离信息;通过固定X、Y轴扫描转镜,在数据处理系统控制模块实现单像素的先验标定;通过对X、Y轴扫描转镜的控制,控制出射激光到达目标的不同位置,从而达到对整个视场,即目标的整个平面进行二维扫描,得到每个像素的距离信息及误差信息,合成三维距离图像,对该三维距离图像在数据处理系统控制模块进行实时校正。 
2.根据权利要求1所述的光子计数三维成像激光雷达距离行走误差校正系统,其特征在于单像素数据处理系统控制模块的先验标定过程为: 
(1)事先准备高反射率平板一块作为目标,放置于光子计数三维成像激光雷达的望远镜镜头(1)正前方固定距离处; 
(2)使得光子计数三维成像激光雷达的X轴扫描转镜(4)、Y轴扫描转镜(3)处于静止状态,设定脉冲激光光源输出脉冲能量为最大,激光脉冲能量记为E0; 
(3)使得脉冲激光光源处于周期工作状态,连续输出激光脉冲,激光脉冲重复次数均设定为Ntotal;同时,采用时间相关单光子计数模块分别记录每个激光脉冲中光子飞行时间,记为ttof(0,j),1≤j≤Ntotal
(4)计算激光脉冲能量E0对应的激光脉冲响应率为:R(0)=Npusle(0)/Ntotal,其中,表示Npusle(0)在测量持续时间内,探测器响应的脉冲总数;Ntotal表示在在测量持续时间内,激光脉冲重复次数; 
(5)计算激光脉冲能量E0对应的光子飞行时间测量均值: 
(6)计算距离行走误差
Figure FDA0000489576650000022
其中,ttof为恒定不变量,是光子在光子计数三维激光成像雷达与目标间光子的飞行时间; 
(7)逐次减小脉冲激光光源的激光脉冲输出能量,重复步骤(2)至(6),得到不同激光脉冲响应率对应的距离行走误差,记为terror(R(i))=f(R(i)),0≤i<n; 
(8)采用数值拟合方式,得到距离行走误差函数,记为terror(R)=f(R)。 
3.根据权利要求1所述的光子计数三维成像激光雷达距离行走误差校正系统,其特征在于三维图像数据处理系统控制模块的实时校正过程为: 
1)通过对X、Y轴扫描转镜的控制,控制出射激光到达目标的不同位置,达到对目标不同像素点的探测,统计像素位置为(i,j)的光子飞行时间测量值,根据脉冲个数,得到该位置对应的光子飞行时间测量均值
Figure FDA0000489576650000023
2)统计像素位置为(i,j)的光子脉冲响应率,得到该位置对应的激光脉冲响应率R(i,j),即: 
R(i,j)=Npulse(i,j)/Ntotal(i,j) 
其中,Npulse(i,j)表示在像素位置(i,j)处的探测器输出脉冲总数;Ntotal(i,j)表示在像素位置(i,j)处的激光脉冲总数; 
3)计算像素位置为(i,j)的距离行走误差,得到距离补偿图像terror(i,j),即:terror(i,j)=f(R(i,j)); 
4)补偿像素位置为(i,j)的光子飞行时间测量均值并得到校正后的光子飞行时间测量值,即: 
Figure FDA0000489576650000031
5)统计每个像素点校正后的光子飞行时间测量值,合成目标的三维距离信息图及形状特征。 
4.一种利用权利要求1所述的光子计数三维成像激光雷达距离行走误差校正系统实现校正方法,其特征在于包括先验标定和实时校正方法,其中先验标定方法包括以下步骤: 
(1)事先准备高反射率平板一块作为目标,放置于光子计数三维成像激光雷达的望远镜镜头(1)正前方固定距离处; 
(2)使得光子计数三维成像激光雷达的X轴扫描转镜(4)、Y轴扫描转镜(3)处于静止状态,设定脉冲激光光源输出脉冲能量为最大,激光脉冲能量记为E0; 
(3)使得脉冲激光光源处于周期工作状态,连续输出激光脉冲,激光脉冲重复次数均设定为Ntotal;同时,采用时间相关单光子计数模块分别记录每个激光脉冲中光子飞行时间,记为ttof(0,j)(1≤j≤Ntotal); 
(4)计算激光脉冲能量E0对应的激光脉冲响应率为:R(0)=Npusle(0)/Ntotal,其中,表示Npusle(0)在测量持续时间内,探测器响应的脉冲总数;Ntotal表示在在测量持续时间内,激光脉冲的总数; 
(5)计算激光脉冲能量E0对应的光子飞行时间测量均值: 
(6)计算距离行走误差
Figure FDA0000489576650000033
其中,ttof为恒定不变量,是光子在光子计数三维激光成像雷达与目标间光子的飞行时间; 
(7)逐次减小脉冲激光光源的激光脉冲输出能量,重复步骤(2)至(6),得到不 同激光脉冲响应率对应的距离行走误差,记为terror(R(i))=f(R(i)),0≤i<n; 
(8)采用数值拟合方式,得到距离行走误差函数,记为terror(R)=f(R); 
其中实时校正方法包括以下步骤: 
1)通过对X、Y轴扫描转镜的控制,控制出射激光到达目标的不同位置,达到对目标不同像素点的探测,统计像素位置为(i,j)的光子飞行时间测量值,得到该位置对应的光子飞行时间测量均值
Figure FDA0000489576650000041
2)统计像素位置为(i,j)的光子脉冲响应率,得到该位置对应的激光脉冲响应率R(i,j),即: 
R(i,j)=Npulse(i,j)/Ntotal(i,j) 
其中,Npulse(i,j)表示在像素位置(i,j)处的探测器输出脉冲总数;Ntotal(i,j)表示在像素位置(i,j)处的激光脉冲总数; 
3)计算像素位置为(i,j)的距离行走误差,得到距离补偿图像terror(i,j),即:terror(i,j)=f(R(i,j)); 
4)补偿像素位置为(i,j)的光子飞行时间测量均值并得到校正后的光子飞行时间测量值,即: 
Figure FDA0000489576650000042
5)统计每个像素点校正后的光子飞行时间测量值,合成目标的三维距离信息图及形状特征的三维距离图像。 
5.根据权利要求4所述的光子计数三维成像激光雷达距离行走误差校正系统实现校正方法,其特征在于先验标定的步骤(7)中,逐次减小脉冲激光光源的激光脉冲输出能量值为0.1E0至0.15E0之间固定值。 
CN201210574661.6A 2012-12-26 2012-12-26 光子计数三维成像激光雷达距离行走误差校正系统及方法 Expired - Fee Related CN103064076B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210574661.6A CN103064076B (zh) 2012-12-26 2012-12-26 光子计数三维成像激光雷达距离行走误差校正系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210574661.6A CN103064076B (zh) 2012-12-26 2012-12-26 光子计数三维成像激光雷达距离行走误差校正系统及方法

Publications (2)

Publication Number Publication Date
CN103064076A CN103064076A (zh) 2013-04-24
CN103064076B true CN103064076B (zh) 2014-06-25

Family

ID=48106780

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210574661.6A Expired - Fee Related CN103064076B (zh) 2012-12-26 2012-12-26 光子计数三维成像激光雷达距离行走误差校正系统及方法

Country Status (1)

Country Link
CN (1) CN103064076B (zh)

Families Citing this family (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103576162A (zh) * 2013-10-25 2014-02-12 中国科学院半导体研究所 激光雷达装置及利用该装置测量目标物距离的方法
CN104914446B (zh) * 2015-06-19 2017-06-27 南京理工大学 基于光子计数的三维距离图像时域实时去噪方法
CN105068083B (zh) * 2015-07-08 2018-05-08 中国科学院合肥物质科学研究院 一种实现单像素激光雷达纵向超分辨率成像的方法
JP2018533026A (ja) * 2015-11-05 2018-11-08 ルミナー テクノロジーズ インコーポレイテッド 高分解能のデプスマップを作成するための、改善された走査速度を有するライダーシステム
DK3206045T3 (da) 2016-02-15 2021-06-21 Airborne Hydrography Ab Enkelt foton lidar scanner
CN105807271A (zh) * 2016-05-11 2016-07-27 深圳乐行天下科技有限公司 激光雷达的校准系统及方法
CN107015233B (zh) * 2017-05-10 2023-12-29 南京理工大学紫金学院 一体化光纤式伪随机码幅度调制偏移校正装置
CN110869804B (zh) * 2017-07-13 2023-11-28 苹果公司 用于光发射深度传感器的提前-滞后脉冲计数
CN110608669A (zh) * 2018-06-15 2019-12-24 上海弼智仿生高科技有限公司 三维扫描方法、装置和系统
US11681030B2 (en) 2019-03-05 2023-06-20 Waymo Llc Range calibration of light detectors
CN110018495B (zh) * 2019-04-30 2021-03-26 湖南力研光电科技有限公司 一种条纹管成像激光雷达激光发射随机误差测量及补偿系统
KR102332512B1 (ko) * 2019-11-04 2021-11-29 현대모비스 주식회사 스캐닝 라이다의 오차 보정 방법
US11747453B1 (en) 2019-11-04 2023-09-05 Waymo Llc Calibration system for light detection and ranging (lidar) devices
CN111289101B (zh) * 2019-12-17 2022-08-12 西南技术物理研究所 便携式象限探测器脉冲响应率参数校准装置
CN110940964A (zh) * 2019-12-31 2020-03-31 西安炬光科技股份有限公司 激光雷达及信号识别方法
CN111273260B (zh) * 2020-03-03 2023-10-27 丹阳钒曦光电科技有限公司 一种脉冲激光测距系统晶体管型接收电路距离测量精度提高方法
CN113433563B (zh) * 2020-03-06 2022-06-28 宁波飞芯电子科技有限公司 测距方法及测距装置
CN111708005A (zh) * 2020-08-24 2020-09-25 中国地质大学(武汉) 单光子激光雷达水下光子位移校正、测深方法及装置
CN112904766B (zh) * 2021-01-18 2022-06-21 福州大学 基于像素误差反馈的旋转双棱镜指向控制系统与方法
CN112767284B (zh) * 2021-02-03 2024-03-08 哈尔滨工业大学 基于光子计数熵的激光三维成像云雾后向散射滤除方法及系统
CN113009496B (zh) * 2021-03-05 2024-05-28 南京大学 一种激光测距振镜成像系统及方法
CN113156409B (zh) * 2021-04-07 2022-07-15 南京工业职业技术大学 一种直接序列时间相关光子计数误差补偿方法
CN113466836A (zh) * 2021-06-23 2021-10-01 深圳市欢创科技有限公司 一种测距方法、装置及激光雷达
CN113721228B (zh) * 2021-09-07 2024-05-03 北京航空航天大学 面阵单光子探测系统参数校正及数据处理方法
CN114545428B (zh) * 2022-03-02 2024-02-02 中国科学院光电技术研究所 基于单像素-单光子探测器跟踪测距激光雷达装置及方法
CN118011410B (zh) * 2024-04-09 2024-07-12 深圳市欢创科技股份有限公司 一种测距方法、激光雷达、机器人及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102062861A (zh) * 2010-11-30 2011-05-18 浙江大学 基于单探测器关联成像原理的三维成像方法
CN102608619A (zh) * 2012-03-07 2012-07-25 北京航空航天大学 光子计数压缩采样相控阵激光三维成像方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009029376A1 (de) * 2009-09-11 2011-05-12 Robert Bosch Gmbh Photonendetektor mit paralysierbarem Photonen-empfindlichem Element, insbesondere SPAD, sowie Entfernungsmessgerät mit solchem Photonendetektor

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102062861A (zh) * 2010-11-30 2011-05-18 浙江大学 基于单探测器关联成像原理的三维成像方法
CN102608619A (zh) * 2012-03-07 2012-07-25 北京航空航天大学 光子计数压缩采样相控阵激光三维成像方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A DUAL-THRESHOLD METHOD FOR PHOTON COUNTING IMAGING WITH THE EMCCD;Beibei Zhou 等;《Proceedings of 2010 IEEE 17th International Conference on Image Processing》;20100929;第3333-3336页 *
Beibei Zhou 等.A DUAL-THRESHOLD METHOD FOR PHOTON COUNTING IMAGING WITH THE EMCCD.《Proceedings of 2010 IEEE 17th International Conference on Image Processing》.2010,第3333-3336页.
利用光电经纬仪修正地基红外搜索跟踪系统静态误差的方法;彭晨 等;《红外与激光工程》;20121031;第41卷(第10期);第2791-2794页 *
彭晨 等.利用光电经纬仪修正地基红外搜索跟踪系统静态误差的方法.《红外与激光工程》.2012,第41卷(第10期),第2791-2794页.

Also Published As

Publication number Publication date
CN103064076A (zh) 2013-04-24

Similar Documents

Publication Publication Date Title
CN103064076B (zh) 光子计数三维成像激光雷达距离行走误差校正系统及方法
US7554652B1 (en) Light-integrating rangefinding device and method
CN102540195B (zh) 一种车用五路激光雷达及其控制方法
CN101839981B (zh) 激光成像回波波形和层次特征获取方法及装置
US20090185159A1 (en) Distance measuring method and distance measuring element for detecting the spatial dimension of a target
EP2975428B1 (en) Three-dimensional imaging radar system
CN108594254A (zh) 一种提高tof激光成像雷达测距精度的方法
EP2260325B1 (en) Light-integrating rangefinding device and method
CN110244318B (zh) 基于异步ToF离散点云的3D成像方法
CN111458108B (zh) 一种发射接收光轴平行度测量装置及测量方法
CN105115445A (zh) 基于深度相机与双目视觉复合的三维成像系统及成像方法
CN107907885B (zh) 一种基于单光子计数方法的水下目标探测装置
CN105675633B (zh) 一种x射线分幅相机的标定装置
CN1327240C (zh) 激光测距仪脉冲回波处理方法及装置
CN103994719A (zh) 基于盖革apd阵列的高精度三维成像装置及其使用方法
CN102255655B (zh) 激光通信中跟踪相机兼容实现链路效率的检测方法
CN111796295B (zh) 一种采集器、采集器的制造方法及距离测量系统
WO2021189794A1 (zh) 飞行时间测距系统及其测距方法
CN102944879B (zh) 一种基于mems二维扫描镜的四维成像装置及其成像方法
CN102252828B (zh) 一种监测高反射光学元件在激光辐照下反射率实时变化的方法
EP3862787A1 (en) De-jitter of point cloud data for target recognition
CN108459328A (zh) 一种具有均匀接收光学系统的探测装置
CN104567796A (zh) 一种3d拍摄测距方法
CN211927241U (zh) 一种发射接收光轴平行度测量装置
US11353584B2 (en) Laser rangefinder and method for implementing the same

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140625

Termination date: 20161226