专利名称：光子计数三维成像激光雷达距离行走误差校正系统及方法
技术领域：
本发明属于光学、激光雷达和图像处理技术，特别是一种光子计数三维成像激光雷达距离行走误差校正系统及方法。
背景技术：
距离图采用三维影像模式描述感兴趣区域，数据立方中包含了目标的几何不变特性，可以避免二维图像中的扭曲和混淆，广泛应用于工业模具设计、军事目标自动识别等领域。目前已知有若干技术途径可以用于获取距离图像。其中，光子计数三维激光成像雷达采用具有单光子灵敏度、工作于盖革模式下的雪崩二极管作为激光回波信号光子探测器，以光子飞行模式实现高时间分辨率测量。对目标进行多脉冲重复测量，再利用统计原理、光子计数技术得到目标单一像素的距离信息，通过扫描整个视场，得到目标的完整三维距离信息。光子计数三维激光成像雷达的测量精度是最为重要的参数。精度定义为测量结果与实际距离间的差异，是探测器响应特性的随工作时间、工作状态和外部输入条件等因素改变而产生测量漂移，也称之为距离行走误差，以光子飞行时间测量误差或测量长度误差单位表示。为了克服距离行走误差，已有两种技术途径。(I)假定造成距离行走误差的原因是GmAro输出电流脉冲上升速率因激光脉冲回波光子数多寡存在差异，提出一种采用多级阈值比较电路的校正方法[G.Kirchner，F. Koidl, et al，Proc. SPIE. 3218，106-112(1997).]。(2)认为造成距离行走误差是激光脉冲回波光子数差异造成的GmAPD探测概率变化，进而提出了基于GmAro探测概率模型的查表式校正方法[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轴扫描转镜、四分之一波片、第二半波片，在分光棱镜处部分透射，经过滤波片，达光纤耦合器，形成回波信号，被单光子探测器所探测接收，经过时间相关单光子计数�？椴刂剐藕攀淙胧荽硐低晨刂颇？椋宦龀寮す夤庠捶⒊鲋芷诼龀寮す猓备枋荽硐低晨刂颇？橥狡鹗夹藕牛庋低澈螅す夥瓷湫藕诺酱锏ス庾犹讲馄鳎だ讲馄餍纬上煊β龀澹皇奔湎喙氐ス庾蛹剖？楹褪荽硐低晨刂颇？榻邮眨岷贤狡鹗夹藕牛唇刂剐藕庞胪狡鹗夹藕诺氖奔洳睿檬奔洳畛艘怨馑伲纬赡勘昃嗬胄畔ⅲ煌ü潭╔、Y轴扫描转镜，在数据处理系统控制�？槭迪值ハ袼氐南妊楸甓ǎ煌ü訶、Y轴扫描转镜的控制，控制出射激光到达目标的不同位置，从而达到对整个视场，即目标的整个平面进行二维扫描，得到每个像素的距离信息及误差信息，合成三维距离图像，对该三维距离图像在数据处理系统控制�？榻惺凳毙Ｕ�。本发明与现有技术相比，其显著优点(I)同时考虑GmAro上升电流和探测率的影响，更符合实际情况，更准确，能够消除距离行走误差；(2)测量的同时进行误差校正，实时性好；(3)实施容易，无需额外增加硬件系统复杂性；(4)避免了测量误差对距离信息的干扰，提闻成像质量。下面结合附图对本发明作进一步详细描述。
图1是本发明光子计数三维成像激光雷达距离行走误差校正系统的示意图。图2是本发明采用的光子计数三维激光成像雷达的先验标定实验测量结果及数值拟合结果。图3是本发明适用于光子计数三维成像激光雷达的先验标定距离行走误差校正方法的实时校正方法流程示意图。图4是本发明实施效果验证中所采用的目标及实验场景实物图。图5是本发明实施效果验证中光子计数三维成像激光雷达输出的原始距离图像。图6是本发明实施效果验证中光子计数三维成像激光雷达输出的原始距离图像对应的响应率图像。图7是本发明实施效果验证中光子计数三维成像激光雷达输出的原始距离图像对应的距离行走误差补偿图像。图8是本发明实施效果验证中光子计数三维成像激光雷达输出的原始距离图像对应的实时校正图像。
具体实施例方式本发明适用于光子计数三维成像激光雷达的先验标定距离行走误差校正方法适用装置为光子计数三维成像激光雷达。首先通过多组实验，实验中改变出射激光的能量，从而得到不同的回波光子数，定义回波光子数与出射脉冲总数的比值为激光脉冲响应率。每次实验可测得相应脉冲响应率下的距离行走误差。然后采用数学拟合的方法，获取距离行走误差关于激光脉冲响应率的函数，记为距离行走误差函数。通过此次标定，就能得到该实验系统的距离行走误差校正函数，即校正方法。实验测量中，得到含有距离行走误差的原始三维距离图像，同时可以得到三维距离图像中每个像素点的激光脉冲响应率，根据此脉冲响应率，结合误差行走函数，能够得到该像素点的距离行走误差，然后对其进行补偿校正，就能够得到更接近真实目标的三维距离图像。结合图1，光子计数三维成像激光雷达距离行走误差校正(实验)系统，包括光学系统、单光子探测器、脉冲激光光源、时间相关单光子计数�？�(市购产品，如PicoHarp300)和数据处理系统控制�？�,光学系统由望远镜镜头1、透镜2, Y轴扫描转镜3、X轴扫描转镜4、准直透镜组5、四分之一波片6、第二半波片7、分光棱镜8、第一半波片9、滤波片10、光纤耦合器11组成，脉冲激光光源同时发出同步起始信号和脉冲周期激光，同步起始信号输入数据处理系统控制�？椋龀逯芷诩す庖来尉谝话氩ㄆ�9、分光棱镜8部分反射后再经过第二半波片7、四分之一波片6、Y轴扫描转镜3、准直透镜组5、Χ轴扫描转镜4、透镜2和望远镜镜头1，出射到目标，经过目标漫反射被望远镜镜头I接收，依次经过透镜2、Χ轴扫描转镜4、准直透镜组5、Υ轴扫描转镜3、四分之一波片6、第二半波片7，在分光棱镜8处部分透射，经过滤波片10，达光纤耦合器11，形成回波信号，被单光子探测器所探测接收，经过时间相关单光子计数�？椴刂剐藕攀淙胧荽硐低晨刂颇？椋宦龀寮す夤庠捶⒊鲋芷诼龀寮す猓备枋荽硐低晨刂颇？橥狡鹗夹藕牛庋低澈螅す夥瓷湫藕诺酱锏ス庾犹讲馄鳎だ讲馄餍纬上煊β龀澹皇奔湎喙氐ス庾蛹剖？楹褪荽硐低晨刂颇？榻邮眨岷贤狡鹗夹藕牛唇刂剐藕庞胪狡鹗夹藕诺氖奔洳睿檬奔洳畛艘怨馑伲纬赡勘昃嗬胄畔ⅲ煌ü潭╔、Y轴扫描转镜，在数据处理系统控制模 块实现单像素的先验标定；通过对X、Y轴扫描转镜的控制，控制出射激光到达目标的不同位置，从而达到对整个视场，即目标的整个平面进行二维扫描，得到每个像素的距离信息及误差信息，合成三维距离图像，对该三维距离图像在数据处理系统控制模块进行实时校正。上述单像素数据处理系统控制模块的先验标定过程和实时校正过程同方法中步骤。结合图3，本发明利用上述的光子计数三维成像激光雷达距离行走误差校正系统实现校正方法，包括先验标定和实时校正方法，其中先验标定方法包括以下步骤(I)事先准备高反射率平板一块作为目标，放置于光子计数三维成像激光雷达的望远镜镜头I正前方固定距离处；(2)使得光子计数三维成像激光雷达的X轴扫描转镜4、Υ轴扫描转镜3处于静止状态，设定脉冲激光光源输出脉冲能量为最大，激光脉冲能量记为Etl ；(3)使得脉冲激光光源处于周期工作状态，连续输出激光脉冲，激光脉冲重复次数均设定为Nttrtal ;同时，采用时间相关单光子计数�？榉直鸺锹济扛黾す饴龀逯泄庾臃尚惺奔洌俏� ttQf(0，j)(l 彡 j 彡Ntrtal)；(4)计算激光脉冲能量Etl对应的激光脉冲响应率为R(0) = Npusle (O)/Nttrtal，其中，表示Npusle (O)在测量持续时间内，探测器响应的脉冲总数；Ntotal表示在在测量持续时间内，激光脉冲的总数；(5)计算激光脉冲能量Etl对应的光子飞行时间测量均值
权利要求
1.一种光子计数三维成像激光雷达距离行走误差校正系统，其特征在于包括光学系统、单光子探测器、脉冲激光光源、时间相关单光子计数模块和数据处理系统控制�？椋庋低秤赏毒稻低�(I)、透镜(2)、Y轴扫描转镜(3)、Χ轴扫描转镜(4)、准直透镜组(5)、四分之一波片￠)、第二半波片(7)、分光棱镜(8)、第一半波片(9)、滤波片(10)、光纤耦合器(11)组成，脉冲激光光源同时发出同步起始信号和脉冲周期激光，同步起始信号输入数据处理系统控制�？椋龀逯芷诩す庖来尉谝话氩ㄆ�(9)、分光棱镜(8)部分反射后在经过第二半波片(7)、四分之一波片(6)、Y轴扫描转镜(3)、准直透镜组(5)、X轴扫描转镜(4)、透镜(2)和望远镜镜头(I)，出射到目标，经过目标漫反射被望远镜镜头(I)接收，依次经过透镜(2)、X轴扫描转镜(4)、准直透镜组(5)、Y轴扫描转镜(3)、四分之一波片￠)、第二半波片(7)，在分光棱镜(8)处部分透射，经过滤波片(10)，达光纤耦合器(11)，形成回波信号，被单光子探测器所探测接收，经过时间相关单光子计数�？椴刂剐藕攀淙胧荽硐低晨刂颇？椋� 脉冲激光光源发出周期脉冲激光，同时给予数据处理系统控制�？橥狡鹗夹藕牛庋低澈螅す夥瓷湫藕诺酱锏ス庾犹讲馄鳎だ讲馄餍纬上煊β龀澹皇奔湎喙氐ス庾蛹剖？楹褪荽硐低晨刂颇？榻邮眨岷贤狡鹗夹藕牛唇刂剐藕庞胪狡鹗夹藕诺氖奔洳睿檬奔洳畛艘怨馑伲纬赡勘昃嗬胄畔ⅲ煌ü潭╔、Y轴扫描转镜，在数据处理系统控制�？槭迪值ハ袼氐南妊楸甓ǎ煌ü訶、Y轴扫描转镜的控制，控制出射激光到达目标的不同位置，从而达到对整个视�。茨勘甑恼銎矫娼卸瑁玫矫扛鱿袼氐木嗬胄畔⒓拔蟛钚畔ⅲ铣扇嗬胪枷瘢愿萌嗬胪枷裨谑荽硐低晨刂颇？榻惺凳毙Ｕ�
2.根据权利要求1所述的光子计数三维成像激光雷达距离行走误差校正系统，其特征 在于单像素数据处理系统控制�？榈南妊楸甓ü涛� (1)事先准备高反射率平板一块作为目标，放置于光子计数三维成像激光雷达的望远镜镜头(I)正前方固定距离处； (2)使得光子计数三维成像激光雷达的X轴扫描转镜(4)、Y轴扫描转镜(3)处于静止状态，设定脉冲激光光源输出脉冲能量为最大，激光脉冲能量记为Etl ； (3)使得脉冲激光光源处于周期工作状态，连续输出激光脉冲，激光脉冲重复次数均设定为Ntotal ;同时，采用时间相关单光子计数�？榉直鸺锹济扛黾す饴龀逯泄庾臃尚惺奔洌俏� ttof (0，j) (I ^ j ^ Ntotal)； (4)计算激光脉冲能量Etl对应的激光脉冲响应率为R(0)=Npusle(O)/Nttrtal，其中，表示Npusle (O)在测量持续时间内，探测器响应的脉冲总数；Ntotal表示在在测量持续时间内，激光脉冲的总数； (5)计算激光脉冲能量Etl对应的光子飞行时间测量均值
3.根据权利要求1所述的光子计数三维成像激光雷达距离行走误差校正系统，其特征在于三维图像数据处理系统控制�？榈氖凳毙Ｕ涛� 1)通过对X、Y轴扫描转镜的控制，控制出射激光到达目标的不同位置，达到对目标不同像素点的探测，统计像素位置为(i，j)的光子飞行时间测量值，根据脉冲个数，得到该位置对应的光子飞行时间测量均值; 2)统计像素位置为(i，j)的光子脉冲响应率，得到该位置对应的激光脉冲响应率R(i，j)，即 RU，j) = Npulse (i, j)/Ntotal (i，j) 其中，Npulse(i，j)表示在像素位置(i，j)处的探测器输出脉冲总数；Nt()tal(i，j)表示在像素位置(i，j)处的激光脉冲总数； 3)计算像素位置为(i，j)的距离行走误差，得到距离补偿图像即te OT(i，j) = f (R(i, j))； 4)补偿像素位置为(i，j)的光子飞行时间测量均值并得到校正后的光子飞行时间测量值，即 Korect (D) = ^me (β，I) + Kmr I) 5)统计每个像素点校正后的光子飞行时间测量值，合成目标的三维距离信息图及形状特征。
4.一种利用权利要求1所述的光子计数三维成像激光雷达距离行走误差校正系统实现校正方法，其特征在于包括先验标定和实时校正方法，其中先验标定方法包括以下步骤 (1)事先准备高反射率平板一块作为目标，放置于光子计数三维成像激光雷达的望远镜镜头(I)正前方固定距离处； (2)使得光子计数三维成像激光雷达的X轴扫描转镜(4)、Υ轴扫描转镜(3)处于静止状态，设定脉冲激光光源输出脉冲能量为最大，激光脉冲能量记为Etl ； (3)使得脉冲激光光源处于周期工作状态，连续输出激光脉冲，激光脉冲重复次数均设定为Ntotal ;同时，采用时间相关单光子计数�？榉直鸺锹济扛黾す饴龀逯泄庾臃尚惺奔洌俏� ttof (O，j) (I ^ j ^ Ntotal)； (4)计算激光脉冲能量Etl对应的激光脉冲响应率为R(0)=Npusle(O)/Nttrtal，其中，表示Npusle (O)在测量持续时间内，探测器响应的脉冲总数；Ntotal表示在在测量持续时间内，激光脉冲的总数； (5)计算激光脉冲能量Etl对应的光子飞行时间测量均值(6)计算距离行走误差Imw(O)= ^(O)-Q，其中，ttof为恒定不变量，是光子在光子计数三维激光成像雷达与目标间光子的飞行时间； (7)逐次减小脉冲激光光源的激光脉冲输出能量，重复步骤(2)至￠)，得到不同激光脉冲响应率对应的距离行走误差，记为terror(R⑴)=f(R(i))，0彡i < η ;(8)采用数值拟合方式，得到距离行走误差函数，记为(R) =f(R)； 其中实时校正方法包括以下步骤 1)通过对X、Y轴扫描转镜的控制，控制出射激光到达目标的不同位置，达到对目标不同像素点的探测，统计像素位置为(i，j)的光子飞行时间测量值，得到该位置对应的光子飞行时间测量均植； 2)统计像素位置为(i，j)的光子脉冲响应率，得到该位置对应的激光脉冲响应率R(i，j)，即 RU，j) = Npulse (i, j)/Ntotal (i，j) 其中，Npulse(i，j)表示在像素位置(i，j)处的探测器输出脉冲总数；Nt()tal(i，j)表示在像素位置(i，j)处的激光脉冲总数； 3)计算像素位置为(i，j)的距离行走误差，得到距离补偿图像即te OT(i，j) = f (R(i, j))； 4)补偿像素位置为(i，j)的光子飞行时间测量均值并得到校正后的光子飞行时间测量值，即 ^cored /) — ^me (β, j) + ^error (β，j) 5)统计每个像素点校正后的光子飞行时间测量值，合成目标的三维距离信息图及形状特征的三维距离图像。
5.根据权利要求4所述的光子计数三维成像激光雷达距离行走误差校正方法，其特征在于先验标定的步骤(7)中，逐次减小脉冲激光光源的激光脉冲输出能量值为O.1Etl至O.15E。之间固定值。
全文摘要
本发明公开了一种光子计数三维成像激光雷达距离行走误差校正系统及方法，首先对光子计数三维成像激光雷达的距离行走误差进行多点标定，得到不同激光脉冲响应率下的距离行走误差；然后采用数学拟合的方法，获取距离行走误差关于激光脉冲响应率的函数，记为距离行走误差函数；最后，在光子计数三维成像激光雷达的工作过程中，统计原始三维距离图像的激光脉冲响应率分布，采用距离行走误差函数对距离行走误差进行预测，对原始三维距离图像应用补偿进行校正，得到校正后的三维距离图像。本发明能够消除距离行走误差。
文档编号G01S7/497GK103064076SQ20121057466
公开日2013年4月24日 申请日期2012年12月26日 优先权日2012年12月26日
发明者何伟基, 陈钱, 顾国华, 司马博羽, 陈云飞, 张闻文, 钱惟贤, 隋修宝, 于雪莲, 路东明 申请人:南京理工大学
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
