一种基于梯度域修正的激光雷达引导滤波去噪方法技术

本发明专利技术公开了一种基于梯度域修正的激光雷达引导滤波去噪方法，输入目标波形G、引导波形I，判断是否存在尖锐起伏，调整参数α；计算目标波形G和引导波形I的协方差，再分别计算目标波形G和引导波形I的方差，逐像素点相乘；根据归一化的G、I波形协方差动态调整滤波器窗宽；对设定窗宽内所有包含像素点p′的线性系数求平均值，得到平均线性系数

【技术实现步骤摘要】
一种基于梯度域修正的激光雷达引导滤波去噪方法
本专利技术属于激光雷达
，具体为一种基于梯度域修正的激光雷达引导滤波去噪方法，是一种通过波形的梯度域信息修正引导滤波的去噪方法，实现对含尖锐凸起、含多脉宽回波分量的激光雷达波形数据动态去噪。
技术介绍
激光雷达LiDAR(LightDetectionandRanging)是一种通过激光脉冲对目标的距离和速度进行探测的主动式遥感技术，原理是向待测目标发射激光脉冲，脉冲在同目标相互作用后，通过采样器对回波信号进行接收采样，对得到的数据进行适当的处理以获得待测目标的相关信息，随着激光成像技术的发展，激光雷达开始具备对目标方向、速度、姿态、形状等信息探测的能力，并被广泛应用于探测领域。在激光雷达系统中，为了获得更加丰富的环境和目标信息，需要对回波波形分析提取如幅值、脉冲宽度、积分强度等信息。对于一些激光雷达的应用场景，如扫描特殊地形和植被覆盖的区域，回波会出现多回波叠加的情况，此时良好的滤波算法可以在去除噪声的同时保留更多有效的信息，是进行波形深入分析的基本保证。激光雷达回波信号的噪声主要由三部分组成：探测器噪声、环境光噪声和读出电路噪声，其中散粒噪声和热噪声占主导地位，通常的表现形式为高斯白噪声。目前在图像领域应用最广泛的滤除高斯噪声的方法有高斯滤波算法、双边滤波算法和引导滤波算法，但是这些算法对波形的尖锐边缘无法有筛选的保留。传统引导滤波在波形中存在尖锐凸起分量的情况下，为了保证大面积平缓分量的滤波效果，设定的参数不适用于尖锐部分，导致尖锐部分经过滤波后趋于平缓；传统引导滤波在无尖锐凸起分量的情况下，窗宽半径固定，无梯度域自适应窗宽的修正，存在峰值处幅值偏低的问题。激光雷达的某些应用场景如海洋水深探测中会出现不同脉宽的回波分量，甚至尖锐凸起分量。引导滤波是一种图像滤波技术，通过引导图I，对目标图像G(输入图像)进行滤波处理，使得最后的输出图像大体上与目标图像G相似，但是纹理部分与引导图I相似，广泛用于图像平滑与边缘保留。激光雷达在一些特定的应用场景全波形呈现出大脉宽变化，甚至出现尖锐分量。传统引导滤波在尖锐存在的情况下需要保证平缓波形的滤波效果，存在对尖锐凸起边缘保留效果差的缺点；在无尖锐凸起的情况存在波峰处幅值偏低，窗宽不能跟随波形梯度域微调与修正等缺点。面对激光雷达不断拓展的应用环境，传统滤波器不能满足激光雷达回波的多种变化，不能适应回波分量大脉宽变化，波形存在尖锐起伏时无法保证平缓分量的滤波效果。因此，亟需一种基于梯度域修正的激光雷达引导滤波去噪方法。
技术实现思路
为解决现有技术存在的缺陷，本专利技术提供一种基于梯度域修正的激光雷达引导滤波去噪方法。为了解决上述技术问题，本专利技术提供了如下的技术方案：本专利技术一种基于梯度域修正的激光雷达引导滤波去噪方法，包括以下步骤：输入目标波形G、引导波形I，判断是否存在尖锐起伏，调整参数α；计算目标波形G和引导波形I的协方差，再分别计算目标波形G和引导波形I的方差，逐像素点相乘；根据归一化的G、I波形协方差，动态调整滤波器窗宽；对设定窗宽内所有包含像素点p′的线性系数求平均值，得到平均线性系数考虑上述尖锐起伏判断参数α、基于梯度域的窗宽修正，得到改进后的引导滤波函数，对输入波形进行滤波处理，输出滤波后的波形。作为本专利技术的一种优选技术方案，本专利技术提供一种基于梯度域修正的激光雷达引导滤波去噪方法，包括以下步骤：步骤一、输入目标波形G，引导波形I，Z与I在以像素k为中心的窗口中存在局部线性关系，则输出波形为：Z(p)＝apG(p)+bp(1)其中，ap、bp为线性系数；去噪后目标波形为：Z(p)＝G(p)-n(p)(2)其中，n(p)为噪声；联立(1)、(2)得到代价函数：其中，p代表输入波形像素点，ap′、bp′分别为对应窗宽Ωδ(p′)的线性系数，ψn为敏感度参数；步骤二、令代价函数E最小，得到：bp′＝μX，δ(p′)-aμG，δ(p′)其中，α为尖锐阈值判断系数，当需要保留输入波形中存在的尖锐边缘时，α＝1；当输入波形中均为平缓高斯分量时，α＝0；ε＝(0.001×(max(G(p))-min(G(p)))2)2步骤三、将引导波形与输入波形的协方差归一化至[0.75，1.25]，得到varG，I(p)，基于梯度域动态变化窗宽；δ(p)＝varG，I(p)×δ0其中，初始窗宽δ0与系统采样率和噪声水平有如下关系：其中，为波形噪声水平参数；步骤四、去设定窗宽内所有包含像素点p′的线性系数平均值；步骤五、得到最终输出结果波形。本专利技术的有益效果是：本专利技术根据波形的梯度域信息改变滤波核，动态调整滤波器窗宽，实现了对尖锐起伏边缘有效保留的同时对平缓区域的良好滤波，基于尖锐阈值判断、梯度域窗宽修正的引导滤波去噪，实现激光雷达大脉宽变化全波形回波的动态去噪。附图说明图1是本专利技术一种基于梯度域修正的激光雷达引导滤波去噪方法的工作流程图。图2是本专利技术一种基于梯度域修正的激光雷达引导滤波去噪方法中基于梯度域信息的动态窗宽对应图。图3是本专利技术一种基于梯度域修正的激光雷达引导滤波去噪方法中基于梯度域修正滤波去噪方法与现有滤波去噪方法的效果对比图。图4是本专利技术一种基于梯度域修正的激光雷达引导滤波去噪方法中基于梯度域修正滤波去噪方法与现有滤波去噪方法在波形存在尖锐起伏情况下的效果对比图。具体实施方式以下结合附图对本专利技术的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本专利技术，并不用于限定本专利技术。为了达到本专利技术的目的，如图1所示，在本专利技术的其中一种实施方式中提供一种基于梯度域修正的激光雷达引导滤波去噪方法，包括以下步骤：输入目标波形G、引导波形I，判断是否存在尖锐起伏，调整参数α；计算目标波形G和引导波形I的协方差，再分别计算目标波形G和引导波形I的方差，逐像素点相乘；根据归一化的G、I波形协方差，动态调整滤波器窗宽；对设定窗宽内所有包含像素点p′的线性系数求平均值，得到平均线性系数考虑上述尖锐起伏判断参数α、基于梯度域的窗宽修正，得到改进后的引导滤波函数，对输入波形进行滤波处理，输出滤波后的波形。为了进一步地优化本专利技术的实施效果，在本专利技术的另一种实施方式中，在前述内容的基础上，本专利技术提供一种基于梯度域修正的激光雷达引导滤波去噪方法，包括以下步骤：步骤一、输入目标波形G，引导波形I，Z与I在以像素k为中心的窗口中存在局部线性关系，则输出波形为：Z(p)＝apG(p)+bp(1)其中，ap、bp为线性系数；去噪后目标波形为：Z(p)＝G(p)-n(p)(2)其中，n本文档来自技高网...

【技术保护点】
1.一种基于梯度域修正的激光雷达引导滤波去噪方法，其特征在于，包括以下步骤：/n输入目标波形G、引导波形I，判断是否存在尖锐起伏，调整参数α；/n计算目标波形G和引导波形I的协方差，再分别计算目标波形G和引导波形I的方差，逐像素点相乘；/n根据归一化的G、I波形协方差，动态调整滤波器窗宽；/n对设定窗宽内所有包含像素点p′的线性系数求平均值，得到平均线性系数

【技术特征摘要】
1.一种基于梯度域修正的激光雷达引导滤波去噪方法，其特征在于，包括以下步骤：
输入目标波形G、引导波形I，判断是否存在尖锐起伏，调整参数α；
计算目标波形G和引导波形I的协方差，再分别计算目标波形G和引导波形I的方差，逐像素点相乘；
根据归一化的G、I波形协方差，动态调整滤波器窗宽；
对设定窗宽内所有包含像素点p′的线性系数求平均值，得到平均线性系数
考虑上述尖锐起伏判断参数α、基于梯度域的窗宽修正，得到改进后的引导滤波函数，对输入波形进行滤波处理，输出滤波后的波形。


2.根据权利要求1所述的基于梯度域修正的激光雷达引导滤波去噪方法，其特征在于，包括以下步骤：
步骤一、输入目标波形G，引导波形I，Z与I在以像素k为中心的窗口中存在局部线性关系，则输出波形为：
Z(p)＝apG(p)+bp(1)
其中，ap、bp为线性系数；去噪后目标波形为：
Z(p)＝G(p)-n(p)(2)
其中，n(p)为噪声；联立(1)、(2)得到代价函数：



其中，p代表输入波形像素点a...

【专利技术属性】
技术研发人员：赵毅强，夏显召，王品权，陈瑞，周意遥，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津;12