一种单光子探测时间-强度信息的恢复方法技术

本发明专利技术公开了一种单光子探测时间‑强度信息的恢复方法，包括：对多次探测得到的光子数进行累积平均；对探测距离门内前x个探测时间间隔单元探测到的背景光子噪声数进行平均；计算探测距离门内接收到的回波光子总数的第一估计值；计算探测距离门内接收到的回波光子总数的第二估计值；对第一估计值和第二估计值进行加权求和；得到回波强度分布概率的初始估计值；对探测距离门内各个探测时间间隔单元内接收到的回波光子数进行恢复迭代估计；并计算测距离门内的回波光子总数；判断测距离门内的回波光子总数是否位于第一阈值和第二阈值之间；得到恢复信号的时间‑强度概率分布函数。本发明专利技术能够恢复单光子探测时目标纵深回波信息。

【技术实现步骤摘要】
一种单光子探测时间-强度信息的恢复方法
本专利技术属于激光雷达
，具体涉及一种单光子探测时间-强度信息的恢复方法。
技术介绍
单光子激光雷达有着探测灵敏度高、有效作用距离远等优势，在近年来获得了长足发展。目前，单光子激光雷达在目标信息获取中存在的主要问题在于：（1）在进行目标信息获取时主要采取首光子成像模式，即通过记录目标首个回波光子到达时间来确定目标距离信息，进而实现成像。在对纵深目标探测时，目标回波信号表现为沿时间维度上的强度分布，仅仅通过首光子成像无法获取目标纵深信息。此外，由于噪声和回波光子随机性效应，造成首光子触发时间产生随机抖动，影响对目标距离的准确测定。（2）可以通过多次光子计数信号累积来获得回波光子数随时间变化信息，但是受制于单光子探测器死时间影响，直接累积信号将发生畸变，不能反映准确的目标回波的时间-强度变化信息，也就难以在此基础上准确提前目标纵深信息。
技术实现思路
本专利技术的目的在于提供一种单光子探测时间-强度信息的恢复方法，该恢复方法能够准确反映目标回波的时间-强度变化信息，且能够恢复单光子探测时目标纵深回波信息。为了达到上述目的，本专利技术采用如下技术方案实现：一种单光子探测时间-强度信息的恢复方法，所述恢复方法包括如下步骤：步骤一、对多次探测得到的光子数进行累积平均，得到单个探测时间间隔单元内累积平均信号；步骤二、对探测距离门内前x个探测时间间隔单元探测到的背景光子噪声数进行平均，得到单个探测时间间隔单元内平均光子噪声估计值；步骤三、根据光子雷达方程，计算探测距离门内接收到的回波光子总数的第一估计值；步骤四、根据单个探测时间间隔单元内累积平均信号和平均光子噪声估计值，计算探测距离门内接收到的回波光子总数的第二估计值；步骤五、对第一估计值和第二估计值进行加权求和，得到探测距离门内接收到的回波光子总数的第三估计值；步骤六、根据单个探测时间间隔单元内累积平均信号、平均光子噪声估计值和第三估计值，得到回波强度分布概率的初始估计值；步骤七、根据回波强度分布概率的初始估计值，对探测距离门内各个探测时间间隔单元内接收到的回波光子数进行恢复迭代估计；并根据恢复迭代估计结果，计算测距离门内的回波光子总数；步骤八、判断测距离门内的回波光子总数是否位于第一阈值和第二阈值之间，如是，则更新第三估计值，进入步骤九；如否，则更新第三估计值，返回步骤六；步骤九、根据单个探测时间间隔单元内累积平均信号、平均光子噪声估计值和更新后的第三估计值，得到恢复信号的时间-强度概率分布函数。进一步的，步骤一中，所述单个探测时间间隔单元内累积平均信号为；；其中，Nsum(d)为第d个探测时间间隔单元内累积平均信号；Ni(d)为第i次探测第d个探测时间间隔单元探测到的光子数，d=1,2,3,...,W，i=1,2,3,...,M，M为探测次数，W为探测距离门内探测时间间隔单元数。进一步的，步骤二中，所述单个探测时间间隔单元内平均光子噪声估计值按照如下公式计算：；其中，Nn_even为单个探测时间间隔单元内平均光子噪声估计值；Ni(j)为第i次探测第j个探测时间间隔单元探测的光子数，i=1,2,3,...,M，M为探测次数，j=1,2,3,...,x,，x小于W，W为探测距离门内探测时间间隔单元数。进一步的，步骤四中，所述第二估计值为：；其中，Ns2为第二估计值；Nsum(d)为第d个探测时间间隔单元内累积平均信号；Nn_even为单个探测时间间隔单元内平均光子噪声估计值；d=1,2,3,...,W，W为探测距离门内探测时间间隔单元数。进一步的，步骤五中，所述第三估计值为：；其中，Ns为第三估计值；Ns1和Ns2分别为第一估计值和第二估计值；α为修正系数，取0~1之间。进一步的，步骤六中，所述回波强度分布概率的初始估计值为：；其中，P(d)为探测距离门内第d个探测时间间隔单元内回波强度分布概率的初始估计值；Nsum(d)为第d个探测时间间隔单元内累积平均信号；Nn_even为单个探测时间间隔单元内平均光子噪声估计值；Ns为第三估计值；d=1,2,3,...,W，W为探测距离门内探测时间间隔单元数。进一步的，步骤七中，按照如下公式，对探测距离门内各个探测时间间隔单元内接收到的回波光子数进行恢复迭代估计：；；其中，N(d)为探测距离门内第d个探测时间间隔单元接收到的光子数估计值；N'(d)为探测距离门内前d-1个探测时间间隔单元接收到的总光子数；d=1,2,3,...,W，W为探测距离门内探测时间间隔单元数，N'(0)=0。进一步的，步骤七中，所述测距离门内的回波光子总数为：；其中，N''为测距离门内的回波光子总数。进一步的，步骤八中，所述更新后的第三估计值为：；其中，N's为更新后的第三估计值。进一步的，步骤九的具体实现过程为：步骤91、按照如下公式，计算回波强度分布概率值；；其中，P'(d)为探测距离门内第d个探测时间间隔单元内回波强度分布概率值；步骤92、判断P'(d)和P(d)的差值绝对值是否小于第三阈值，如是，则P'(d)为恢复信号的时间-强度概率分布函数，结束；如否，则令P(d)=P'(d)，返回步骤七。本专利技术的有益效果：本专利技术根据光子雷达方程，实现对探测距离门内接收到的回波光子总数的一次估计；并通过单个探测时间间隔单元内累积平均信号和平均光子噪声估计值，实现对探测距离门内接收到的回波光子总数的二次估计；并通过对两次估计的结果进行加权求和，实现对探测距离门内接收到的回波光子总数的三次估计，避免由于探测器响应死时间效应造成的接收光子数大幅偏离真实目标回波光子数，得到一个较为接近目标真实回波光子数的初始估计值；利用单个探测时间间隔单元内累积平均信号和平均光子噪声估计值以及第三次的估计结果，得到回波强度分布概率的初始估计值，并通过迭代估计，最终恢复目标回波准确的时间-强度分布信息，且避免了噪声、回波光子随机性和探测器死时间效应造成的测距误差，通过恢复出的目标回波信息，保证了目标纵深方向信息提取的准确性；本专利技术通过对探测距离门内各个探测时间间隔单元内接收到的回波光子数进行恢复迭代估计，并利用更新后的测距离门内的回波光子总数对信号恢复进行判断，并在信号恢复判断的基础上，有效地恢复出目标回波的时间-强度分布信息，从而提高了目标纵深方向信息提取的准确性。附图说明图1为本专利技术的单光子探测时间-强度信息恢复方法流程示意图。具体实施方式以下结合附图对本专利技术的具体实施方式作出详细说明。本实施例给出了一种单光子探测时间-强度信息的恢复方法，参考图1，该恢复方法包括如下步骤：步骤一、对多次探测得到的光子数进行累积平均，得到单个探测时间间隔单本文档来自技高网...

【技术保护点】
1.一种单光子探测时间-强度信息的恢复方法，其特征在于，所述恢复方法包括如下步骤：/n步骤一、对多次探测得到的光子数进行累积平均，得到单个探测时间间隔单元内累积平均信号；/n步骤二、对探测距离门内前

【技术特征摘要】
1.一种单光子探测时间-强度信息的恢复方法，其特征在于，所述恢复方法包括如下步骤：
步骤一、对多次探测得到的光子数进行累积平均，得到单个探测时间间隔单元内累积平均信号；
步骤二、对探测距离门内前x个探测时间间隔单元探测到的背景光子噪声数进行平均，得到单个探测时间间隔单元内平均光子噪声估计值；
步骤三、根据光子雷达方程，计算探测距离门内接收到的回波光子总数的第一估计值；
步骤四、根据单个探测时间间隔单元内累积平均信号和平均光子噪声估计值，计算探测距离门内接收到的回波光子总数的第二估计值；
步骤五、对第一估计值和第二估计值进行加权求和，得到探测距离门内接收到的回波光子总数的第三估计值；
步骤六、根据单个探测时间间隔单元内累积平均信号、平均光子噪声估计值和第三估计值，得到回波强度分布概率的初始估计值；
步骤七、根据回波强度分布概率的初始估计值，对探测距离门内各个探测时间间隔单元内接收到的回波光子数进行恢复迭代估计；并根据恢复迭代估计结果，计算测距离门内的回波光子总数；
步骤八、判断测距离门内的回波光子总数是否位于第一阈值和第二阈值之间，如是，则更新第三估计值，进入步骤九；如否，则更新第三估计值，返回步骤六；
步骤九、根据单个探测时间间隔单元内累积平均信号、平均光子噪声估计值和更新后的第三估计值，得到恢复信号的时间-强度概率分布函数。


2.根据权利要求1所述的恢复方法，其特征在于，步骤一中，所述单个探测时间间隔单元内累积平均信号为；

；
其中，Nsum(d)为第d个探测时间间隔单元内累积平均信号；Ni(d)为第i次探测第d个探测时间间隔单元探测到的光子数，d=1,2,3,...,W，i=1,2,3,...,M，M为探测次数，W为探测距离门内探测时间间隔单元数。


3.根据权利要求1所述的恢复方法，其特征在于，步骤二中，所述单个探测时间间隔单元内平均光子噪声估计值按照如下公式计算：

；
其中，Nn_even为单个探测时间间隔单元内平均光子噪声估计值；Ni(j)为第i次探测第j个探测时间间隔单元探测的光子数，i=1,2,3,...,M，M为探测次数，j=1,2,3,...,x,，x小于W，W为探测距离门内探测时间间隔单元数。


4.根据权利要求1~3中任意一项所述的恢复方法，其特征在于，步骤四中，所述第二估计值为...

【专利技术属性】
技术研发人员：胡以华，赵楠翔，徐世龙，方佳节，董骁，侯阿慧，骆盛，雷武虎，李敏乐，魏硕，
申请(专利权)人：中国人民解放军国防科技大学，
类型：发明
国别省市：湖南;43