基于光子计数的三维距离图像时域实时去噪方法
【技术领域】
[0001] 本发明适用于光子计数三维成像激光雷达系统，特别是一种基于光子计数的三维 距离图像时域实时去噪方法。
【背景技术】
[0002] 距离图像采用三维影像模式描述感兴趣的目标区域，其数据立方中包含了目标的 几何不变特性，可以避免二维图像中的扭曲和混淆，广泛应用于机器视觉、工业模具设计和 军事目标的自动识别与跟踪等领域。目前，已知有若干技术途径可以用于获取目标的距离 图像。其中，采用工作于光子计数模式下的Gm-Aro作为其激光回波信号光子探测器的光子 计数三维成像激光雷达系统，是以光子飞行模式实现对目标的高时间分辨率测量。
[0003] 成像精度和成像时间是光子计数三维成像激光雷达系统的两个重要的性能参数。 传统的基于最大似然估计算法的成像模型，是以单个探测光子为最小单元，通过信号累积 生成光子计数直方图的方式，逐点估计目标的距离信息一对应光子计数直方图的峰值位 置。为了克服单光子探测过程中固有泊松噪声的影响，每个像素点需要足够长的采样积分 时间；而且，每个像素点的采样积分时间为预设固定值。然而，对于未知特性的复杂目标场 景来说，无法直接准确估计每个像素点所需要的采样积分时间，且不同反射率特性的目标 区域所需要的采样积分时间也不相同。因此，在目标特性未知的情况下，距离估计值要么出 现采样不足或是采样饱和的现象，即距离估计精度受噪声影响严重，或是在采样积分时间 增加到某一值后，距离估计精度趋于某一值附近变化，且不依赖至少不明显依赖于采样积 分时间。
[0004] 目前，解决上述问题的方法，已知有两种技术途径：一是Hyunjung Shim和 Seungkyu Lee等人通过采用一种新型的成像模式，即采用多种曝光时间的混合曝光技术， 描述距离估计中出现的采样不足和采样饱和的现象，自适应决定每个像素点的最佳距离 值，提高成像质量（1· Hyunjung Shim and Seungkyu Lee, Hybrid exposure for depth imaging of a time-〇f-flight depth sensor. Optical Express,2014,voI.22 (II),p p:13393-13402. 2. T. Mertens and J. Kautz and F.Van Reeth,Exposure Fusion, IEEE Pacific Conference on Computer Graphics and Applications, 2007,382-390.) 〇 二 是Ahmed Kirmani等人通过建立单光子探测过程的概率统计模型，并结合目标相邻像 素点间的空间相关性，利用第一个探测到回波光子信息，估计目标的三维图像（3. Ahmed kirmani, et al,First-Photon Imaging, Science,2014,vol. 343, pp:58-61. 4. Ahmed kirmani, et al, Spatio-temporal regularization for range imaging with high photon efficiency. SPIE, 2013.)。这两种技术途径均需要较复杂的后续算法处理，不适用于实时 的应用场合；其中，第二种方法在目标场景的低反射率区域，特别是边缘细节部分容易引起 平滑现象。
【发明内容】
[0005] 本发明的目的在于提供一种实现快速提取目标清晰无噪的基于光子计数的三维 距离图像时域实时去噪方法。
[0006] 实现本发明目的的技术解决方案为：一种基于光子计数的三维距离图像时域实时 去噪方法，步骤如下：
[0007] 第一步，估计目标距离的初步先验信息，在目标场景的每个像素点探测过程中，利 用噪声光子和信号光子的不同特性改进传统的基于最大似然估计算法的成像模型，即利用 连续探测到的N个相邻光子为单元，并以该N个光子的飞行时间均值作为目标距离的初步 先验?目息估计；
[0008] 第二步，时域实时去噪，抑制噪声光子的影响，在每个像素点探测过程中，利用上 述第一步得到的目标距离的初步先验信息，作为判断信号光子响应和噪声光子响应的判别 标准，结合噪声光子和信号光子的飞行时间所具有的不同特性，对探测器后续响应到的光 子进行顺序检测，并对这信号光子和噪声光子采取不同的处理方法：忽略被判定为噪声的 光子响应，仅保留被判定为信号的光子响应，并利用其光子飞行时间对目标距离的初步先 验信息进行实时更新，即对目标距离的初步先验信息进行增强，使其逐渐趋近目标距离的 真实值；
[0009] 第三步，逐点估计目标的距离信息，重构出目标的三维距离图像，利用上述第二步 在每个像素点处采集得到的信号光子计数Κ，作为每个像素点停止采样的阈值条件，即在每 个像素点处探测过程中，当探测器响应到K个信号光子响应时，则转移到下一像素点重复 这一过程，直到扫描完整个目标场景为止；并利用该K个信号光子的飞行时间的加权平均 值作为该像素点的距离估计，重构出目标清晰的三维距离图像。
[0010] 本发明与现有技术相比，其显著优点：（1)不仅能够抑制光子计数过程固有泊松 噪声，包括背景噪声光子和暗电流噪声光子的影响，而且避免了生成光子计数直方图的过 程，从而减少了成像时间，提高系统的光子效率；（2)首先利用噪声光子和信号光子的不同 特性，改进传统的基于最大似然估计算法的成像模型，以连续探测到的N个相邻光子为单 元，估计目标距离初步先验信息，并利用后续探测到的信号光子信息对该初步先验信息进 行实时更新，使目标距离估计逐渐收敛于距离真实值；（3)接下来利用目标距离初步先验 信息作为判别标准，并利用采集得到的信号光子计数K作为阈值条件，自适应决定每个像 素点的采样积分时间，能够避免传统的基于最大似然估计算法的成像模型中出现的采样 不足和采样饱和的现象，从而提高成像精度；(4)由于忽略噪声光子计数，仅对信号光子处 理，因此本发明性能较为稳定，在不同强度噪声环境中，都能获得较为稳定的成像质量，且 不增加系统对硬件的要求。
[0011] 下面结合附图对本发明作进一步详细描述。
【附图说明】
[0012] 图1是用于验证本发明的实验目标场景图：图中区域1为高反射率区域（墙），区 域2为低反射率区域（表面覆盖黑布的纸箱）。
[0013] 图2是传统的基于最大似然估计算法的成像方法在采样积分时间为Ims时获取的 目标距离图像：图（b)对应图（a)中的区域1，图（c)对应区域2。
[0014] 图3是传统的基于最大似然估计算法的成像方法在采样积分时间为5ms时获取的 目标距离图像：图（b)对应图（a)中的区域1，图（c)对应区域2。
[0015] 图4是传统的基于最大似然估计算法的成像方法在采样积分时间为20ms时获取 的目标距离图像：图（b)对应图（a)中的区域1，图（c)对应区域2。
[0016] 图5是本发明基于光子计数三维距离图像时域实时去噪方法的流程图。
[0017] 图6是本发明选取N = 3个相邻光子为单元得到的目标距离的初步先验信息。
[0018] 图7是本发明选取阈值K= 10得到的目标三维距离图像。
【具体实施方式】
[0019] 结合图1和图5,本发明基于光子计数的三维距离图像时域实时去噪方法，步骤如 下：
[0020] 第一步，估计目标距离的初步先验信息，在目标场景的每个像素点探测过程中，利 用噪声光子和信号光子的不同特性改进传统的基于最大似然估计算法的成像模型（传统 的成像模型是以单个探测光子为单元，忽略噪声光子和信号光子间的区别，对这两类光子 响应采取同样的处理策略，通过长时间信号累积生成光子计数直方图的方式，估计目标的 距离信息），即利用连续探测到的N个相邻光子为单元（时域N-领域），并以该N个光子的 飞行时间均值作为目标距离的初步先验信息估计；N > 1，K多1。
[0021] 估计目标距离的初步先验信息步骤如下：
[0022] 光子计数过程固有泊松噪声，包括背景噪声和暗电流噪声，其光子飞行时间在整 个探测周期时间区间[0, TJ内服从均匀分布，且相互独立，其中?；为激光脉冲重复周期；而 目标反射激光脉冲回波光子计数，即信号光子计数，其飞行时间集中分布在相对较窄的时 间区间
内，且时间相关，其中^为像素点（i，j)处的目标反射激光 脉冲回波光子的飞行时间的真实值，Tp为激光脉冲宽度，且T ρ〈〈?；;
[0023] 基于噪声光子和信号光子的这一不同特性，在像素点（i，j)处探测过程中，本 发明以连续探测到的N个相邻光子（时域N-领域）为单元，其光子飞行时间集记为
其飞行时间均值： 其顺序绝对偏差（Rank-Ordered Absolute ，
,则根据如下二元假设检验，判断该光子单元 的光子飞行时间均值是否落在区间
内，即可以判断其为噪声光子响 应还是信号光子响应：
[0024] If ROAD (i，j) < (N-1) Tp，Signal-Detection
[0025] If ROA