基于光子计数熵的激光三维成像云雾后向散射滤除方法及系统

摘要

基于光子计数熵的激光三维成像云雾后向散射滤除方法及系统，涉及激光三维成像技术领域。本发明是为了解决现有云雾环境激光三维成像目标信号峰淹没于后向散射信号，提取目标信号峰需要长时间统计探测，测量实时性差，且云雾多径效应导致目标信号峰展宽，距离分辨率低的问题。本发明利用两步滤噪处理滤除后向散射噪声，提高信噪比，能够获取高精度目标距离像。利用滑窗处理结合光子计数熵滤除后向散射噪声，同时增强目标信号的信噪比，实现高精度实时距离像获取。

说明书

技术领域

本发明属于激光三维成像技术领域。

背景技术

云雾环境激光三维成像对科学和技术领域都具有重要的意义，无论在地形测绘、自动驾驶等民用方面，还是战场环境预警、导航等军用方面。云雾环境激光三维成像目标信号峰淹没于后向散射信号，提取目标信号峰需要长时间统计探测，测量实时性差，且云雾多径效应导致目标信号峰展宽，距离分辨率低。

发明内容

本发明是为了解决现有云雾环境激光三维成像目标信号峰淹没于后向散射信号，提取目标信号峰需要长时间统计探测，测量实时性差，且云雾多径效应导致目标信号峰展宽，距离分辨率低的问题，现提供基于光子计数熵的激光三维成像云雾后向散射滤除方法及系统。

基于光子计数熵的激光三维成像云雾后向散射滤除方法，该方法基于以下装置实现，所述装置包括：扫描装置、发射/接收系统、盖格模式雪崩光电二极管；

光源发射的激光脉冲能够入射至扫描装置，扫描装置能够使得其接收到的光束以任意角度入射至视场内，扫描装置的出射光束能够通过发射/接收系统照射至被测目标，所述被测目标位于云雾环境中，经被测目标反射后的光束能够被发射/接收系统接收，发射/接收系统能够将反射后的光束照射至盖格模式雪崩光电二极管的光敏面；

设K为光源发射激光脉冲的总次数，j＝1,2,...,L+M，L为包含目标的回波信号的总数，M为不包含目标的回波信号的总数，所述基于光子计数熵的激光三维成像云雾后向散射滤除方法包括以下步骤：

初始化：启动光源发出激光脉冲，使j＝1，

步骤一：调整扫描装置使其接收到的光束以角度α

步骤二：将光源发射激光脉冲的时刻标记为计时起点，将盖格模式雪崩光电二极管发生响应的时刻标记为计时终点，将计时起点与计时终点作差获得飞行时间，记录该飞行时间内的回波信号，所述回波信号为光子计数次数，

步骤三：判断j是否等于L+M，是则执行步骤四，否则使j＝j+1，然后返回步骤一，

步骤四：分别计算每个回波信号对应的光子计数熵，并依据光子计数熵判断每个回波信号是否包含目标，

步骤五：利用不包含目标的回波信号特征拟合视场内后向散射峰的分布，将该分布作为估计后的后向散射信号，

步骤六：将估计后的后向散射信号分别与每个包含目标的回波信号作差，使得后向散射噪声被滤除。

进一步的，上述步骤四中根据下式计算光子计数熵：

其中，Y(k)为{y(n)}的离散傅里叶变换结果，|Y(k)|

进一步的，上述步骤四中判断每个回波信号是否包含目标的具体方法为：

预设光子计数熵阈值为H

若H

若H

进一步的，上述步骤五中，将所有不包含目标的回波信号进行统计平均，获得统计平均值

其中，

利用gamma分布对

进一步的，上述基于光子计数熵的激光三维成像云雾后向散射滤除方法还包括以下步骤：

步骤七：利用滑窗将步骤六作差后所有差值向量均匀划分为Q个计数子区间，分别对每个计数子区间进行离散傅里叶变换，并计算差值向量的光子计数熵，

步骤八：选取最小光子计数熵对应计数子区间的时延量作为目标距离。

进一步的，上述步骤七中，滑窗划分的第q个计数子区间y

其中，q＝0,1,…,Q-1，

进一步的，根据下式对第q个计数子区间y

其中，Y

基于光子计数熵的激光三维成像云雾后向散射滤除系统，包括：扫描装置、发射/接收系统、计时模块、盖格模式雪崩光电二极管、目标点检测模块、后向散射估计模块和后向散射滤除模块；

光源发射的激光分成两部分、并分别入射至扫描装置和计时模块，

扫描装置使得其接收到的一部分光束能够以任意角度入射至视场内，扫描装置的出射光束通过发射/接收系统照射至被测目标，所述被测目标位于云雾环境中，经被测目标反射后的光束被发射/接收系统接收，发射/接收系统将反射后的光束照射至盖格模式雪崩光电二极管的光敏面，盖格模式雪崩光电二极管的响应信号输出端连接计时模块的响应信号输入端，

计时模块将接收到响应信号时刻标记为计时终点，计时模块将接收到光源发射激光的时刻标记为计时起点，计时模块还用于将计时起点与计时终点作差获得飞行时间及该飞行时间内的回波信号，

目标点检测模块用于根据每个回波信号计算光子计数熵，并依据光子计数熵判断每个回波信号是否包含目标，

后向散射估计模块用于利用不包含目标的回波信号特征拟合视场内后向散射峰的分布，将该分布作为估计后的后向散射信号，

后向散射滤除模块用于将估计后的后向散射信号分别与每个包含目标的回波信号作差，使得后向散射噪声被滤除。

进一步的，上述基于光子计数熵的激光三维成像云雾后向散射滤除系统还包括随机噪声滤除模块，随机噪声滤除模块用于利用滑窗将作差后所有差值向量均匀划分为Q个计数子区间，分别对每个计数子区间进行离散傅里叶变换，并计算差值向量的光子计数熵，选取最小光子计数熵对应计数子区间的时延量作为目标距离。

与现有技术相比，本发明的有益效果是：

本发明无需现有方法依靠长时间统计探测累积目标信号峰，而是在不增加系统复杂度的前提下，利用两步滤噪处理滤除后向散射噪声，提高信噪比，能够获取高精度目标距离像。利用滑窗处理结合光子计数熵滤除后向散射噪声，同时增强目标信号的信噪比，实现高精度实时距离像获取。

附图说明

图1为基于光子计数熵的激光三维成像云雾后向散射滤除系统的结构框图；

图2为光子计数熵滤除后向散射噪声流程示意图；

图3为光子计数直方图，其中A表示介质峰，B表示目标峰；

图4为后向散射噪声粗滤除结果示意图；

图5为基于光子计数熵的随机噪声滤除结果示意图。

具体实施方式

具体实施方式一：参照图2至图5具体说明本实施方式，本实施方式所述的基于光子计数熵的激光三维成像云雾后向散射滤除方法，该方法基于以下装置实现。

所述装置包括：扫描装置2、发射/接收系统3、盖格模式雪崩光电二极管5。

光源发射的激光脉冲能够入射至扫描装置2，扫描装置2能够使得其接收到的光束以任意角度入射至视场内，扫描装置2的出射光束能够通过发射/接收系统3照射至被测目标10，所述被测目标10位于云雾环境中，经被测目标10反射后的光束能够被发射/接收系统3接收，发射/接收系统3能够将反射后的光束照射至盖格模式雪崩光电二极管5的光敏面。

设K为光源发射激光脉冲的总次数，j＝1,2,...,L+M，L为包含目标的回波信号的总数，M为不包含目标的回波信号的总数，所述基于光子计数熵的激光三维成像云雾后向散射滤除方法包括以下步骤：

初始化：启动光源发出激光脉冲，使j＝1，

步骤一：调整扫描装置2使其接收到的光束以角度α

步骤二：将光源发射激光脉冲的时刻标记为计时起点，将盖格模式雪崩光电二极管5发生响应的时刻标记为计时终点，将计时起点与计时终点作差获得飞行时间，记录该飞行时间内的回波信号，所述回波信号为光子计数次数，

步骤三：判断j是否等于L+M，是则执行步骤四，否则使j＝j+1，然后返回步骤一。

步骤四：依据光子计数直方图中的数据，分别计算每个回波信号对应的光子计数熵，具体的，根据下式计算光子计数熵：

其中，Y(k)为{y(n)}的离散傅里叶变换结果，|Y(k)|

之后，预设光子计数熵阈值为H

以下为滤噪过程，如图2所示。

步骤五：利用不包含目标的回波信号特征拟合视场内后向散射峰的分布，将该分布作为估计后的后向散射信号。具体如下：

将所有不包含目标的回波信号进行统计平均，获得统计平均值

其中，

利用gamma分布对

步骤六：将估计后的后向散射信号分别与每个包含目标的回波信号作差，使得后向散射噪声被滤除。

步骤七：利用滑窗将步骤六作差后所有差值向量均匀划分为Q个计数子区间，设滑窗划分的第q个计数子区间y

其中，q＝0,1,…,Q-1，

根据下式对第q个计数子区间y

其中，Y

之后计算差值向量的光子计数熵。

步骤八：选取最小光子计数熵对应计数子区间的时延量作为目标距离。

本实施方式利用云雾后向散射信号和目标回波信号的光子计数熵特征差异，区分目标像素点和非目标像素点，统计全部非目标像素点的回波信号特征并估计后向散射峰，然后在目标像素点的回波信号中减除后向散射峰，并利用滑窗处理结合光子计数熵滤除剩余的随机噪声，提高信噪比，获得高精度目标距离像。

具体实施方式二：参照图1、3、4和5具体说明本实施方式，本实施方式所述的基于光子计数熵的激光三维成像云雾后向散射滤除系统，包括：扫描装置2、发射/接收系统3、计时模块4、盖格模式雪崩光电二极管5、目标点检测模块6、后向散射估计模块7、后向散射滤除模块8和随机噪声滤除模块9。

由激光器1发射的激光分成两部分、并分别入射至扫描装置2和计时模块4，

扫描装置2使得其接收到的一部分光束能够以任意角度入射至视场内，扫描装置2的出射光束通过发射/接收系统3照射至被测目标10，所述被测目标10位于云雾环境中，经被测目标10反射后的光束被发射/接收系统3接收，发射/接收系统3将反射后的光束照射至盖格模式雪崩光电二极管5的光敏面，盖格模式雪崩光电二极管5的响应信号输出端连接计时模块4的响应信号输入端，

计时模块4将接收到响应信号时刻标记为计时终点，计时模块4将接收到光源发射激光的时刻标记为计时起点，计时模块4还用于将计时起点与计时终点作差获得飞行时间及该飞行时间内的回波信号。累积K个发射激光脉冲后，获得光子到达时间计数直方图，如图3所示。

目标点检测模块6用于根据每个回波信号计算光子计数熵，并依据光子计数熵判断每个回波信号是否包含目标。其中，根据下式计算光子计数熵：

其中，Y(k)为{y(n)}的离散傅里叶变换结果，|Y(k)|

预设光子计数熵阈值为H

后向散射估计模块7用于利用不包含目标的回波信号特征拟合视场内后向散射峰的分布，将该分布作为估计后的后向散射信号。其中，将所有不包含目标的回波信号进行统计平均，获得统计平均值

其中，

利用gamma分布对

后向散射滤除模块8用于将估计后的后向散射信号分别与每个包含目标的回波信号作差，使得后向散射噪声被粗略滤除。

由于GM-APD概率性响应回波光子，光子计数存在随机起伏，后向散射滤除模块8处理之后光子计数直方图中仍然存在大量的随机噪声，且减除处理没有增强目标信号峰，信号峰仍然淹没于噪声中。因此，需要将后向散射滤除模块8处理后含有大量随机噪声的信号输入随机噪声滤除模块9，利用滑窗处理结合光子计数熵的滤噪方法滤除随机噪声，增强目标信号峰。具体的，如图2所示，随机噪声滤除模块9用于利用滑窗将步骤六作差后所有差值向量均匀划分为Q个计数子区间，分别对每个计数子区间进行离散傅里叶变换，并计算差值向量的光子计数熵。由于目标信号峰随机性低，光子计数熵远低于随机噪声，处理后目标信号峰得到增强，信噪比提高。

最后在滤噪处理后的数据中搜索光子计数熵最小的子区间，该子区间的时延量作为目标距离的估计，从而获取高精度目标距离像。

本实施方式以云雾后向散射信号和目标回波信号的光子计数熵特征差异为基础，通过光子计数熵准确区分目标像素点和非目标像素点，利用非目标像素点获得的接收信号估计云雾后向散射峰，然后在目标像素点减除后向散射信号，并利用光子计数熵结合滑窗滤除剩余的随机噪声，提高信噪比，获得高精度目标距离像。