基于光子飞行时间相关性的时域滤波方法
【专利摘要】本发明公开了一种基于光子飞行时间相关性的时域滤波方法。首先利用激光雷达在信噪比比较低的条件下扫描获得光子飞行时间。然后利用光子之间飞行时间相关性，设定滤波条件对光子的飞行时间进行判定，区分出信号光子和噪声光子。接着，去除被判定为噪声光子的探测光子，只保留被判定为信号光子的探测光子。最后，对信号光子使用形心法，提取光子准确的飞行时间，从而得到距离信息。实验表明，在信噪比比较低的情况下，本发明能够有效的滤除三维成像激光雷达扫描得到的信号中的噪声光子，从而达到去噪的目的。
【专利说明】
基于光子飞行时间相关性的时域滤波方法
技术领域
[0001]本发明属于激光雷达三维成像领域，具体涉及一种适用于光子计数激光雷达系统 的基于光子飞行时间相关性的时域滤波方法。 技术背景
[0002] 准确的获取目标场景的三维结构具有广泛的应用。三维成像激光雷达系统以盖革 模式下的雪崩二极管作为探测器，能够高灵敏度和高精度的获取目标场景的深度图像。激 光雷达系统通过对目标场景积分采样，然后累积生成统计直方图，来获取目标场景的三维 信息。
[0003] 为了从统计直方图中提取准确的飞行时间，获得高精度的目标深度估计值，需要 使用回波信号飞行时间提取算法如：峰值法、形心法、拟合法等。另一方面，当噪声较强、目 标反射率较低或积分时间不足时，激光雷达系统累积探测得到的直方图可能堆积不起来， 不会有明显的波峰。此时，使用形心法的话将引入较多的噪声带来误差，由于存在多个峰值 甚至会可能没有峰值，峰值法也不适用，拟合法更不适用。所以在信噪比较低情况下，为了 提高探测概率，就要首先进行滤波，滤除噪声光子，保留信号光子。
[0004] 目前，解决上述问题的方法，已知的技术途径有:HongJin Kong等人采用两个工作 在盖革模式下的Aro对回波信号进行探测。将两个探测器的输出接入与逻辑结构来降低虚 警率，获取清晰的三维距离图像。（[ll]Kong H J，Kim T H，Jo S E，et al.Smart three-dimensional imaging ladar using two Geiger-mode avalanche photodiodes[J] ? Optics express，2011，19(20): 19323-19329 ·)罗韩君等人在理论上进一步研究了这种双 探测器结构的探测性能，指出对于双探测器结构，与逻辑可以很好地抑制虚警率;或逻辑可 提高目标探测概率;使用与逻辑和或逻辑相结合的双探测器结构，可以获得更高的目标探 测概率和更低的虚警率。（[19]罗韩君，詹杰，丰元，等.基于盖革模式AH)的光子计数激光雷 达探测距离研究[J].光电工程，2013,40(12) :80-88.)
[0005] 采用双探测的结构时，回波信号需要经分束器分别进入两个APD，当回波信号的强 度较弱时，这样做可能会降低目标的探测概率;或逻辑虽然可以提高目标的探测概率，但也 提高了虚警率;与逻辑和或逻辑相结合时将需要使用到4个探测器，这无疑会使系统的结构 变得复杂，同时也会引入更多探测器时间抖动带来的误差。
【发明内容】
[0006] 本发明的目的在于，在低信噪比的条件下，提供一种有效滤除噪声光子的时域滤 波方法。
[0007] 实现本发明的技术解决方案为:一种基于光子飞行时间相关性的时域滤波方法步 骤如下：
[0008] 第一步，在低信噪比的情况下，使用三维成像激光雷达系统对目标场景进行扫描， 获得目标场景光子飞行时间分布三维图；
[0009] 第二步，对光子的飞行时间信息进行处理，将当前光子与在它之前到达光子的飞 行时间进行比较来滤波，并且每个像素点中，最早到达的几个光子与最后几个到达的光子 进行比较，这样就能确保每个光子都有其它光子与之比较。计算出每个光子的Ρ%(Μ，Ν)(信 号光子在Μ个光子中和Ν个以上光子相关的概率)和Ρ η。(Μ，Ν)(噪声光子在Μ个光子中和Ν个以 上光子相关的概率），设定:Ρ%(Μ，Ν)彡0.5，？"。(1，《<0.1为滤波条件，对探测到的光子进 行滤波；
[0010] 第三步，滤除噪声光子，并使用形心法提取滤波后光子的飞行时间，从而得到目标 场景的3D图。
[0011] 利用上述三个步骤，本发明能够有效的滤除噪声光子。
[0012] 本发明与现有技术相比，其显著优点：1)本发明利用光子之间的时间相关性来进 行滤波，而不使用双探测器的结构，从而简化了系统的结构，避免了引入更多探测器时间抖 动带来的误差。2)在时间上进行滤波，能快速的达到去噪的目的，提高了三维成像激光雷达 系统的成像时间。3)本发明能有效地提高信号光子的探测概率，减小测距误差，提高成像的 精度。
【附图说明】
[0013]图1是信噪比比较低时得到的直方图。
[0014]图2是在低信噪比条件下，对实验室白墙平面区域进行了 100X 100像素扫描后得 到的光子飞行时间分布三维图。
[0015]图3是对图1的场景不同(Μ，Ν)取值时的滤波效果，（a)(M，N) = (l，l)，（b)(M，N) = (6,2)，（c)(M，N) = (13,3)，（d)(M，N) = (21，4)。
[0016] 图4是实验目标场景，固定在墙上的'Ν'形字母。
[0017] 图5是图3场景滤波前的光子飞行时间分布三维图。
[0018] 图6是图3场景滤波后光子飞行时间分布三维图。
[0019]图7是滤波前用形心法计算得到的距离图。
[0020]图8是滤波后用形心法计算得到的距离图。
【具体实施方式】：
[0021] 当目标场景噪声较强、反射率较低或积分时间不足时，激光雷达系统累积探测得 到的直方图可能堆积不起来，不会有明显的波峰。一些基于直方图的提取飞行时间的算法 如峰值法、形心法、和拟合法将无法使用。因此，在信噪比比较低的情况下，需要首先对回波 信号进行滤波。针对此问题，本发明利用时间相关单光子计数的原理，提出了一种基于光子 之间时间相关性的时域滤波方法。首先利用激光雷达在信噪比比较低的条件下扫描获得光 子飞行时间。然后利用光子之间飞行时间相关性，设定滤波条件对光子的飞行时间进行判 定，区分出信号光子和噪声光子。接着，去除被判定为噪声光子的探测光子，只保留被判定 为信号光子的探测光子。最后，对信号光子使用形心法，提取光子准确的飞行时间，从而得 到距离信息。
[0022] 下面结合附图对本发明作进一步详细描述。
[0023]第一步，在低信噪比的情况下，使用三维成像激光雷达系统对目标场景进行扫描， 获得目标场景的光子飞行时间分布三维图。为了模拟信噪比较低的情况，一方面我们调节 激光器的输出功率设置旋钮来降低激光脉冲的能量，另一方面我们将实验系统的光路暴露 在实验室内的背景光下。图1所示的直方图就是用上述方法对实验室内白墙上某一定点积 分5ms得到的。进一步，在图1的信噪比条件下，我们对该定点周围的平面区域进行了 100X 100像素的扫描，积分时间5ms，将每个像素点的光子飞行时间作图后，得到的光子飞行时间 分布的三维图，如图2所示。
[0024] 第二步，对光子的飞行时间信息进行处理，将当前光子与在它之前到达光子的飞 行时间进行比较来滤波，并且每个像素点中，最早到达的几个光子与最后几个到达的光子 进行比较，这样就能确保每个光子都有其它光子与之比较。计算出每个光子的Psc(M，N)(信 号光子在Μ个光子中和N个以上光子相关的概率)和Pnc(M，N)(噪声光子在Μ个光子中和N个 以上光子相关的概率）。由开关激光测得的数据可知，在距离门内Psig = 0.17，Pnoi = 0.83。 同时，回波脉宽Tf = 800ps，H)L 800-B的激光脉冲波形呈近高斯分布，Tp~3Tf，距离门持续 时间Tgate = 20ns。代入计算得到 Psc = 0 · 16，Pnc = 0 · OSqPsc 和 Pnc也即Psc( 1，1)和 Pnc( 1， 1)，由Psc和Pnc可以进一步计算出在前Μ个光子中探测到N个以上相关光子的概率Psc(M，N) 和Pnc(M，N)。这里为了便于对比滤波后的效果，我们将Pnc(M，N)的值维持在Pnc(l，l)的水 平左右。由此我们分别计算了 (M，N)为(6,2)、（13,3)、（21，4)时的？％(1^)和?此(1^)滤波 后的效果如图3所示。从图3可以看出本发明有效的滤除了噪声光子，提高成像的质量。
[0025] 上述实验是对白墙平面进行扫描，反射率单一，为了查看多反射率时的滤波效果， 我们在白墙平面上固定了一个"N"形字母，如图4所示。字母距离墙面5cm左右，字母左半部 分涂成绿色，反射率由高到低依次为墙面、字母白色部分、字母绿色部分，反射率之比大约 为1.9:1.6:1。同样，我们扫描了100 X 100的像素点，每个像素点积分时间为5ms。
[0026] 实验过程中的噪声水平与第一个实验相同，不开激光积分100ms，在135~155ns的 距离门内接收到的光子数在3900个的水平。打开激光扫描后记录的光子总数为3062290个， 则每个像素点返回光子数的平均值约为306个。以此均值计算得P sig = 0.363，PnQi = 0.637， Psc = 0.253，Pnc = 0.08。若以 PSC(M，N)彡0.5，Pnc(M，N)彡0.1为滤波条件，计算得(M，N)为（11， 3)，P SC(11,3) = 0.5528，Pnc(ll，3) =0.0519。图5、6分别给出了滤波前和滤波后光子飞行时 间分布三维图。
[0027] 第三步，滤除噪声光子，并使用形心法提取滤波后光子的飞行时间，从而得到目标 场景的3D图像。图7、8分别给出了滤波前和滤波后在143~147ns的范围内用形心法计算得 到的距离图。图7和图8相比较可以看出，本发明有效去除了噪声，并简化了系统的结构。
【主权项】
1. 一种基于光子飞行时间相关性的时域滤波方法，其特征在于步骤如下： 第一步，在低信噪比的情况下，使用三维成像激光雷达系统对目标场景进行扫描，获得 目标场景光子飞行时间分布三维图； 第二步，将当前光子与在它之前到达的光子的飞行时间进行比较来滤波，并且每个像 素点中，最早到达的几个光子与最后几个到达的光子进行比较，确保每个光子都有其它光 子与之比较;计算出每个光子的Psc(M，N)和?"。(1^)，？％(1，《是信号光子在1个光子中和~ 个以上光子相关的概率，P nc(M，N)是噪声光子在Μ个光子中和N个以上光子相关的概率，设定 PSC(M，N)彡0.5，？"。(1，《<0.1为滤波条件，对探测到的光子进行滤波； 第三步，滤除噪声光子，并使用形心法提取滤波后光子的飞行时间，从而得到目标场景 的3D距离图像。2. 根据权利要求1所述的基于光子飞行时间相关性的时域滤波方法，其特征在于:第一 步所述低信噪比条件为信噪比小于1，目标回波信号淹没在背景噪声和探测器暗噪声中。3. 根据权利要求1所述的基于光子飞行时间相关性的时域滤波方法，其特征在于:所述 第二步的具体步骤如下： 记时间单元大小为τ，距离门持续时间为Tgate，回波信号相对于距离门开始时的到达时 间为Td，回波脉宽为Tf，回波脉冲持续时间为TP;定义:假如两个光子之间的飞行时间差小于 回波脉宽T f，则称这两个光子为相关光子;同时记积分时间为Tacq，经积分时间Tacq在距离门 内探测到的噪声光子数和光子总数分别为Ν_和N sum，则在距离门内，探测到噪声光子和信 号光子的概率Pnoi和Psig表不为：根据概率论，经推到得到信号光子与下一个光子相关的概率为：噪声光子与下一个光子相关的概率为：信号光子或噪声光子与其它任何一个光子相关的概率都为Ps。和Pn。，在Μ个光子中和N个 光子相关的概率可以看成多重伯努利实验，对于信号光子在Μ个光子中和Ν个以上光子相关 的概率为：对于噪声光子在Μ个光子中和Ν个以上光子相关的概率为：对于某一光子，如果在它前Μ个光子中存在Ν个以上的相关光子，则认为它是信号光子 并保留它，Ρ%(Μ，Ν)和Pnc；(M，N)实际表征的是信号光子和噪声光子被保留的概率;选取的滤 波条件为:P SC(M，N)彡0.5，Pnc(M，NK〇.l，?lM=ll，N=3。4.根据权利要求1所述的基于光子飞行时间相关性的时域滤波方法，其特征在于:第三 步中所述形心法即加权平均法，将直方图中所有点的横纵坐标相乘求和，再除以纵坐标之 和，所求结果作为回波信号的飞行时间。
【文档编号】G01S17/89GK105954733SQ201610437477
【公开日】2016年9月21日
【申请日】2016年6月17日
【发明人】冯振超, 卢斯洋, 何伟基, 方剑, 沈姗姗, 陈钱, 顾国华, 张闻文, 钱惟贤, 任侃, 隋修宝, 邹燕
【申请人】南京理工大学