CN107144829A - 一种高效的激光雷达回波信号去噪方法 - Google Patents
一种高效的激光雷达回波信号去噪方法 Download PDFInfo
- Publication number
- CN107144829A CN107144829A CN201710512546.9A CN201710512546A CN107144829A CN 107144829 A CN107144829 A CN 107144829A CN 201710512546 A CN201710512546 A CN 201710512546A CN 107144829 A CN107144829 A CN 107144829A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mode function
- intrinsic mode
- msup
- msub
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/4802—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种高效激光雷达回波信号去噪方法,将激光雷达回波信号作经验模态分解,得到有限个本征模态函数和残差;将本征模态函数由低阶到高阶依次去除,按顺序将每个剩余分量对应的相关系数与设定的阈值依次进行比较,当某个相关系数开始小于设定的阈值时,确定该相关系数对应的剩余分量中的本征模态函数为低频本征模态函数分量,而剩余分量之前的本征模态函数分量为高频本征模态函数分量;对于高频本征模态函数分量,采用软阈值处理方法对其进行去噪;对于低频本征模态函数分量,采用粗糙惩罚平滑模型进行处理;将处理后的本征模态函数分量与残差进行重构。本发明对于高、低频两个函数组分别进行处理,提高激光雷达回波信号的去噪性能。
Description
技术领域
本发明涉及一种去噪方法,特别是一种高效激光雷达回波信号去噪方法。
背景技术
激光雷达是随着激光器和雷达技术的发展而产生的一种主动式的现代光学遥感设备,具有测量范围广、时空分辨率高和抗干扰能力强等特点。激光雷达的应用十分广泛,在军事领域,激光雷达可用于战场侦察、目标跟踪、障碍回避、空间监视等;在民用领域,激光雷达可用于气象测量、大气研究、遥感遥测、污染监测、医疗诊断等。
然而,在实际应用中,探测器接收到的回波信号往往夹杂着背景光和***中的各类噪声。由于回波信号的强度随距离平方衰减,当探测距离较大时,信号将淹没在较强的噪声之中。所以,如何有效的从强背景噪声中提取出有用信号尤为重要。目前,主流的激光雷达测量***主要采用硬件和软件相结合的方法对回波信号中的背景光噪声进行滤除。其中,硬件去噪主要是通过提高激光器输出功率、在探测器前加入窄带滤光片等方式实现噪声滤除,软件去噪主要是利用滑动平均、傅里叶变换、小波变换、经验模态分解等数字滤波技术实现信号的有效提取。
激光雷达回波信号是一种随距离平方衰减的非线性、非平稳信号,衰减的数量级多达7个数量级或更多。基于该特点,滑动平均以及傅里叶变换不能很好地适用该类信号。小波变换是一种适合处理非线性、非平稳信号的方法,王成等利用小波分析对星载激光雷达波形进行处理,得到更符合原始波形的拟合波形,但该方法存在最优基函数的选择问题(授权号CN20101070853)。经验模态分解是由Huang等人提出的一种具有较强的自适应性的时域信号处理方法,特别适用于处理非线性、非平稳信号,能很好地克服小波变换的基函数选择问题。回波信号经过分解可得到一系列频率从大到小的本征模态函数,噪声往往处于高频部分。因此,目前噪声滤除工作主要针对于高频分量。吴松华等在处理回波信号时,将前五项本征模态函数分量直接去除,与传统的低通滤波和多脉冲平均法相比,所获得的信号平滑性更好、扰动性更小。但是,当回波中存在强突变信号时,分解后的分量会产生模态混叠现象,简单舍弃高频分量,会导致有效信号的损失。为有效提取高频分量中的有用信号,阈值法、Savitzky-Golay滤波等多种技术相继被用于高频分量的平滑处理,提高了***的信噪比。然而,忽略了对低频分量的处理,使得重构的回波信号中依然含有噪声。
发明内容
本发明所要解决的技术问题是提供一种高效激光雷达回波信号去噪方法,提高激光雷达回波信号的信噪比。
为解决上述技术问题,本发明所采用的技术方案是:
一种高效激光雷达回波信号去噪方法,其特征在于包含以下步骤:
步骤一:将激光雷达回波信号作经验模态分解,得到有限个本征模态函数和残差;
步骤二:将本征模态函数由低阶到高阶依次去除,保留此阶数之后的信号作为剩余分量,计算每个剩余分量与激光雷达回波信号间的相关系数,按顺序将每个剩余分量对应的相关系数与设定的阈值依次进行比较,当某个相关系数开始小于设定的阈值时,确定该相关系数对应的剩余分量中的本征模态函数为低频本征模态函数分量,而对应的剩余分量之前的本征模态函数分量为高频本征模态函数分量;
步骤三:对于高频本征模态函数分量,采用软阈值处理方法对其进行去噪;
步骤四:对于低频本征模态函数分量,采用粗糙惩罚平滑模型进行处理;
步骤五:将步骤三和步骤四中处理后的本征模态函数分量与残差进行重构,得到去噪后的激光雷达回波信号。
进一步地,所述步骤一具体步骤为,
1.1识别出激光雷达回波信号f(x)的所有极大值以及极小值点, x为探测距离,形成f(x)的上下包络线U、L,计算上下包络线的均值:M=(U+L)/2;
1.2将f(x)减去上下包络线的均值:h=f(x)-M,判断h是否满足本征模态函数的定义,若满足,则得到第一个本征模态函数分量记为h(1),代表回波信号中的最高频率分量;否则,将h作为新的f(x),重复以上步骤;
1.3将f(x)减去h(1),从f(x)中分离出本征模态函数,得到剩余分量:res=f(x)-h(1),判断res是否满足残差的条件,若满足,保留该分量;否则,将res作为新的f(x),重复以上所有步骤。
进一步地,所述步骤二中剩余信号fm(x)之间的相关系数ρ(m)计算公式为
其中,f(x)为激光雷达回波信号;fm(x)为回波信号去除前m个本征模态函数后的信号;h(i)(x)为第i个本征模态函数分量;N为本征模态函数分量的长度;
所得的相关系数呈单调递减,通过与设定的阈值进行比较,获得最后一个大于阈值的相关系数,从而确定第一个低频本征模态函数分量的位置kth:
其中,last表示在ρ(m)中最后一个大于等于C时,m对应的值,C 的取值范围为[0.75,0.85];L为本征模态函数分量的个数。
进一步地,所述步骤三中对于高频本征模态函数分量采用的软阈值方式,其处理方法为
其中,h(i)(x)为第i个本征模态函数分量;Ti为阈值;
阈值选取为
其中,N为本征模态函数分量的长度;median(|h(i)(x)|)为取第i个本征模态函数分量的中值。
进一步地,所述步骤四中对于低频本征模态函数分量采用的粗糙惩罚平滑模型为
其中,IMFi(n)为第i本征模态函数的第n个点;IMFi'(n)为为第i本征模态函数平滑后的第n个点;λ为惩罚系数;y(x)为待估光滑函数; y(x)”为y(x)的二阶导数。上式右边第一项为基本的最小二乘原理,第二项为增加的粗糙惩罚项;
惩罚系数λ选用交叉验证法进行确定
其中,f(n)为验证样本;fn *(λ)为在一个给定的λ下的估计值;
对于上式,得到的最小值所对应的λ即为所需的光滑参数。
本发明与现有技术相比,具有以下优点和效果:
(1)本发明通过分析激光雷达回波信号与本征模态函数分量之间的相关性,从而有效区分高频与低频本征模态函数分量;
(2)本发明利用软阈值处理方式对高频本征模态函数分量进行去噪,克服硬阈值在重构时产生的波动,提高重构信号的光滑性;另外,利用该方法代替直接去除高频分量的去噪方法,保留信号的完整性,提高信噪比;
(3)本发明通过利用粗糙惩罚光滑模型对低频本征模态函数分量进行平滑处理,可在保持信号光滑性的基础上有效避免过度拟合,从而对噪声进行去除,提高信噪比;
(4)本发明将经验模态分解、粗糙惩罚光滑模型以及软阈值有机结合,与传统的去噪方法进行比较,能在保持信号完整性的基础上,具有更好的光滑性,从而提高了信噪比。
附图说明
图1是本发明的一种高效激光雷达回波信号去噪方法的流程图。
图2是本发明的经验模态分解流程图。
图3是本发明的激光雷达回波信号仿真图。
图4是本发明加噪处理后的激光雷达回波信号。
图5是本发明的实施例信噪比为5dB时,激光雷达回波信号经验模态分解结果。
图6是本发明的实施例信噪比为5dB时,相关系数计算结果图。
图7是本发明的实施例信噪比为5dB时,激光雷达回波信号去噪效果。
图8是本发明的实施例信噪比为5dB时,四种方法的激光雷达回波信号去噪效果。
图9是信噪比为5dB时,各方法去噪后的信噪比表格。
图10是信噪比为5dB时,各方法去噪后的均方差表格。
具体实施方式
下面结合附图并通过实施例对本发明作进一步的详细说明,以下实施例是对本发明的解释而本发明并不局限于以下实施例。
如图1所示,本发明的一种高效激光雷达回波信号去噪方法,包含以下步骤:
步骤一:将激光雷达回波信号作经验模态分解,得到有限个本征模态函数和残差;
如图2所示,经验模态分解具体为:
1.1识别出激光雷达回波信号f(x)(x为探测距离)的所有极大值以及极小值点,形成f(x)的上下包络线U、L,计算上下包络线的均值:M=(U+L)/2;
1.2将f(x)减去上下包络线的均值:h=f(x)-M,判断h是否满足本征模态函数的定义,若满足,则得到第一个本征模态函数分量记为h(1),代表回波信号中的最高频率分量;否则,将h作为新的f(x),重复以上步骤;
1.3将f(x)减去h(1),从f(x)中分离出本征模态函数,得到剩余分量:res=f(x)-h(1),判断res是否满足残差的条件,若满足,保留该分量;否则,将res作为新的f(x),重复以上所有步骤。
步骤二:将本征模态函数由低阶到高阶依次去除,保留此阶数之后的信号作为剩余分量,计算每个剩余分量与激光雷达回波信号间的相关系数,按顺序将每个剩余分量对应的相关系数与设定的阈值依次进行比较,当某个相关系数开始小于设定的阈值时,确定该相关系数对应的剩余分量中的本征模态函数为低频本征模态函数分量,而对应的剩余分量之前的本征模态函数分量为高频本征模态函数分量;
剩余信号fm(x)之间的相关系数ρ(m)计算公式为
其中,f(x)为激光雷达回波信号;fm(x)为回波信号去除前m个本征模态函数后的信号;h(i)(x)为第i个本征模态函数分量;N为本征模态函数分量的长度;
所得的相关系数呈单调递减,通过与设定的阈值进行比较,获得最后一个大于阈值的相关系数,从而确定第一个低频本征模态函数分量的位置kth:
其中,last表示在ρ(m)中最后一个大于等于C时,m对应的值,C 的取值范围为[0.75,0.85];L为本征模态函数分量的个数。
步骤三:对于高频本征模态函数分量,采用软阈值处理方法对其进行去噪;
对于高频本征模态函数分量采用的软阈值方式,其处理方法为
其中,h(i)(x)为第i个本征模态函数分量;Ti为阈值;
阈值选取为
其中,N为本征模态函数分量的长度;median(|h(i)(x)|)为取第i个本征模态函数分量的中值。
步骤四:对于低频本征模态函数分量,采用粗糙惩罚平滑模型进行处理;
对于低频本征模态函数分量采用的粗糙惩罚平滑模型为
其中,IMFi(n)为第i本征模态函数的第n个点;IMFi'(n)为为第i本征模态函数平滑后的第n个点;λ为惩罚系数;y(x)为待估光滑函数; y(x)”为y(x)的二阶导数。上式右边第一项为基本的最小二乘原理,第二项为增加的粗糙惩罚项;
惩罚系数λ选用交叉验证法进行确定
其中,f(n)为验证样本;fn *(λ)为在一个给定的λ下的估计值;
对于上式,得到的最小值所对应的λ即为所需的光滑参数。
步骤五:将步骤三和步骤四中处理后的本征模态函数分量与残差进行重构,得到去噪后的激光雷达回波信号。
下面对本发明的一种高效激光雷达回波信号去噪方法进行仿真分析;
如图3所示,根据激光雷达方程对回波信号进行仿真,其激光雷达方程为
其中,P(r)为距离r处的瞬时接收功率;c为光速;E0为发射的激光脉冲能量;Y(r)为激光雷达几何重叠因子;Ar为接收望远镜的有效接收面积;β(r)为后向散射系数;TtTr为发射和接收光学***的总透射率;α(r)为大气消光系数。
在真实仿真信号中加入高斯白噪声,其信噪比为5dB,此信号作为待处理的激光雷达回波信号,如图4所示。
如图5所示,将回波信号利用经验模态分解,得到8个本征模态函数以及残差。
如图6所示,求解相关系数确定第一个低频本征模态函数的位置: kth=6,其中C选为0.85,则前5项分量为高频分量,后3项为低频分量。
对前5项高频本征模态函数分量,采用软阈值处理方法对其进行去噪;对第6-8项低频本征模态函数分量,采用粗糙惩罚光滑模型进行处理,最终得到去噪后的信号如图7所示。
在相同的信噪比条件下,将本发明提出的方法分别与传统的三种去噪方法进行比较。选用的三种去噪方法为:将低频分量以及残差进行部分重构、对高频分量分量采用硬阈值方法进行处理,以及对回波信号采用小波去噪方法。
其中,小波去噪方法的基函数选为db6,分解层数是3层,采用软阈值进行处理。
分别计算各方法去噪后的输出信噪比(SNRout)和均方差(MSE),公式为:
其中,N为激光雷达回波信号的长度;f(n)为回波信号在第n个点处的值;f*(n)为去噪后第个点的值。
结论:四种方法的去噪效果如图8所示,本发明提出的方法较于其他三种方式,去噪后的信噪比最高,达到16.67dB(图9),均方差最小为2.23×10-22(图10),具有较好的平滑效果,去噪效果更好。
本说明书中所描述的以上内容仅仅是对本发明所作的举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种修改或补充或采用类似的方式替代,只要不偏离本发明说明书的内容或者超越本权利要求书所定义的范围,均应属于本发明的保护范围。
Claims (5)
1.一种高效激光雷达回波信号去噪方法,其特征在于包含以下步骤:
步骤一:将激光雷达回波信号作经验模态分解,得到有限个本征模态函数和残差;
步骤二:将本征模态函数由低阶到高阶依次去除,保留此阶数之后的信号作为剩余分量,计算每个剩余分量与激光雷达回波信号间的相关系数,按顺序将每个剩余分量对应的相关系数与设定的阈值依次进行比较,当某个相关系数开始小于设定的阈值时,确定该相关系数对应的剩余分量中的本征模态函数为低频本征模态函数分量,而对应的剩余分量之前的本征模态函数分量为高频本征模态函数分量;
步骤三:对于高频本征模态函数分量,采用软阈值处理方法对其进行去噪;
步骤四:对于低频本征模态函数分量,采用粗糙惩罚平滑模型进行处理;
步骤五:将步骤三和步骤四中处理后的本征模态函数分量与残差进行重构,得到去噪后的激光雷达回波信号。
2.按照权利要求1所述的一种高效激光雷达回波信号去噪方法,其特征在于:所述步骤一具体步骤为,
1.1识别出激光雷达回波信号f(x)的所有极大值以及极小值点,x为探测距离,形成f(x)的上下包络线U、L,计算上下包络线的均值:M=(U+L)/2;
1.2将f(x)减去上下包络线的均值:h=f(x)-M,判断h是否满足本征模态函数的定义,若满足,则得到第一个本征模态函数分量记为h(1),代表回波信号中的最高频率分量;否则,将h作为新的f(x),重复以上步骤;
1.3将f(x)减去h(1),从f(x)中分离出本征模态函数,得到剩余分量:res=f(x)-h(1),判断res是否满足残差的条件,若满足,保留该分量;否则,将res作为新的f(x),重复以上所有步骤。
3.按照权利要求1所述的一种高效激光雷达回波信号去噪方法,其特征在于:所述步骤二中,剩余信号fm(x)之间的相关系数ρ(m)计算公式为
<mrow>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mi>f</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msup>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,f(x)为激光雷达回波信号;fm(x)为回波信号去除前m个本征模态函数后的信号;h(i)(x)为第i个本征模态函数分量;N为本征模态函数分量的长度;
所得的相关系数呈单调递减,通过与设定的阈值进行比较,获得最后一个大于阈值的相关系数,从而确定第一个低频本征模态函数分量的位置kth:
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>t</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>l</mi>
<mi>a</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
<mrow>
<mn>1</mn>
<mo>&le;</mo>
<mi>m</mi>
<mo>&le;</mo>
<mi>L</mi>
</mrow>
</munder>
<mo>{</mo>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mi>C</mi>
<mo>}</mo>
<mo>+</mo>
<mn>1</mn>
</mrow>
其中,last表示在ρ(m)中最后一个大于等于C时,m对应的值,C的取值范围为[0.75,0.85];L为本征模态函数分量的个数。
4.按照权利要求1所述的一种高效激光雷达回波信号去噪方法,其特征在于:所述步骤三中对于高频本征模态函数分量采用的软阈值方式,其处理方法为
<mrow>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>sgn</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msup>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<msup>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>></mo>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<msup>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>&le;</mo>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,h(i)(x)为第i个本征模态函数分量;Ti为阈值;
阈值选取为
<mrow>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>m</mi>
<mi>e</mi>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msup>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>|</mo>
<mo>)</mo>
</mrow>
</mrow>
<mn>0.6745</mn>
</mfrac>
<msqrt>
<mrow>
<mi>ln</mi>
<mi>N</mi>
</mrow>
</msqrt>
</mrow>
其中,N为本征模态函数分量的长度;median(|h(i)(x)|)为取第i个本征模态函数分量的中值。
5.按照权利要求1所述的一种高效激光雷达回波信号去噪方法,其特征在于:所述步骤四中对于低频本征模态函数分量采用的粗糙惩罚平滑模型为
<mrow>
<mi>S</mi>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>IMF</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
<mo>-</mo>
<msubsup>
<mi>IMF</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>&Integral;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,IMFi(n)为第i本征模态函数的第n个点;IMFi'(n)为为第i本征模态函数平滑后的第n个点;λ为惩罚系数;y(x)为待估光滑函数;y(x)”为y(x)的二阶导数。上式右边第一项为基本的最小二乘原理,第二项为增加的粗糙惩罚项;
惩罚系数λ选用交叉验证法进行确定
<mrow>
<mi>C</mi>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>*</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
其中,f(n)为验证样本;fn *(λ)为在一个给定的λ下的估计值;对于上式,得到的最小值所对应的λ即为所需的光滑参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710512546.9A CN107144829B (zh) | 2017-06-29 | 2017-06-29 | 一种高效的激光雷达回波信号去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710512546.9A CN107144829B (zh) | 2017-06-29 | 2017-06-29 | 一种高效的激光雷达回波信号去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107144829A true CN107144829A (zh) | 2017-09-08 |
CN107144829B CN107144829B (zh) | 2019-11-19 |
Family
ID=59784453
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710512546.9A Active CN107144829B (zh) | 2017-06-29 | 2017-06-29 | 一种高效的激光雷达回波信号去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107144829B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107894586A (zh) * | 2017-10-17 | 2018-04-10 | 南京航空航天大学 | 一种基于同步压缩变换的激光雷达回波信号去噪方法 |
CN108594256A (zh) * | 2018-04-16 | 2018-09-28 | 合肥菲涅尔光电科技有限公司 | 一种基于脉冲编码技术的相干激光雷达 |
CN108845306A (zh) * | 2018-07-05 | 2018-11-20 | 南京信息工程大学 | 基于变分模态分解的激光雷达回波信号去噪方法 |
CN108896456A (zh) * | 2018-04-28 | 2018-11-27 | 南京信息工程大学 | 基于反馈型rbf神经网络的气溶胶消光系数反演方法 |
CN109100009A (zh) * | 2018-06-01 | 2018-12-28 | 国网江苏省电力有限公司南京供电分公司 | 基于经验模态分解emd的分接开关振动信号降噪方法 |
CN109541549A (zh) * | 2018-10-09 | 2019-03-29 | 广东工业大学 | 基于emd和稀疏信号处理的间歇采样转发干扰抑制方法 |
CN109598152A (zh) * | 2018-10-11 | 2019-04-09 | 天津大学 | 基于emd降噪数据预处理的硬件木马检测优化方法 |
CN110564905A (zh) * | 2019-10-08 | 2019-12-13 | 中南大学 | 一种对高炉炉衬冲击回波检测的信号处理方法及*** |
CN111868561A (zh) * | 2018-03-20 | 2020-10-30 | 帕诺森斯有限公司 | 使用本底噪声的自适应识别进行有效信号检测 |
CN112230199A (zh) * | 2019-07-15 | 2021-01-15 | 天津大学 | 一种基于高维特征值分析的激光雷达回波盲去噪方法 |
WO2022061597A1 (zh) * | 2020-09-23 | 2022-03-31 | 深圳市速腾聚创科技有限公司 | 一种信号噪声滤除方法、装置、存储介质及激光雷达 |
WO2022217406A1 (zh) * | 2021-04-12 | 2022-10-20 | 深圳市速腾聚创科技有限公司 | 信号处理方法、装置及可读存储介质 |
CN116343051A (zh) * | 2023-05-29 | 2023-06-27 | 山东景闰工程研究设计有限公司 | 一种基于遥感影像的地质环境监测方法及*** |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012025694A1 (fr) * | 2010-08-27 | 2012-03-01 | France Telecom | Traitement de donnees pour le debruitage d'un signal |
CN105447318A (zh) * | 2015-12-01 | 2016-03-30 | 北京科技大学 | 一种微弱信号去噪方法及装置 |
CN105547465A (zh) * | 2015-12-08 | 2016-05-04 | 华北电力大学(保定) | 一种变压器振动信号绕组状态特征提取方法 |
CN106096242A (zh) * | 2016-06-01 | 2016-11-09 | 浙江浙能北海水力发电有限公司 | 一种基于改进emd分解的尾水管压力脉动综合评价方法 |
US9548189B2 (en) * | 2015-04-23 | 2017-01-17 | Lam Research Corporation | Plasma etching systems and methods using empirical mode decomposition |
CN106706122A (zh) * | 2017-01-24 | 2017-05-24 | 东南大学 | 基于相关系数和emd滤波特性的碰摩声发射信号降噪方法 |
-
2017
- 2017-06-29 CN CN201710512546.9A patent/CN107144829B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012025694A1 (fr) * | 2010-08-27 | 2012-03-01 | France Telecom | Traitement de donnees pour le debruitage d'un signal |
US9548189B2 (en) * | 2015-04-23 | 2017-01-17 | Lam Research Corporation | Plasma etching systems and methods using empirical mode decomposition |
CN105447318A (zh) * | 2015-12-01 | 2016-03-30 | 北京科技大学 | 一种微弱信号去噪方法及装置 |
CN105547465A (zh) * | 2015-12-08 | 2016-05-04 | 华北电力大学(保定) | 一种变压器振动信号绕组状态特征提取方法 |
CN106096242A (zh) * | 2016-06-01 | 2016-11-09 | 浙江浙能北海水力发电有限公司 | 一种基于改进emd分解的尾水管压力脉动综合评价方法 |
CN106706122A (zh) * | 2017-01-24 | 2017-05-24 | 东南大学 | 基于相关系数和emd滤波特性的碰摩声发射信号降噪方法 |
Non-Patent Citations (3)
Title |
---|
BEHNOOD RASTI ET AL.: "Hyperspectral image denoising using first order spectral roughness penalty in wavelet domain", 《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》 * |
王欢雪 等: "一种基于交叉证认和经验模态分解的滤波算法及其在激光雷达回波信号降噪处理中的应用", 《中国激光》 * |
许同乐 等: "基于EMD相关方法的电动机信号降噪的研究", 《船舶力学》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107894586B (zh) * | 2017-10-17 | 2020-01-31 | 南京航空航天大学 | 一种基于同步压缩变换的激光雷达回波信号去噪方法 |
CN107894586A (zh) * | 2017-10-17 | 2018-04-10 | 南京航空航天大学 | 一种基于同步压缩变换的激光雷达回波信号去噪方法 |
CN111868561B (zh) * | 2018-03-20 | 2024-05-14 | 祖克斯有限公司 | 使用本底噪声的自适应识别进行有效信号检测 |
CN111868561A (zh) * | 2018-03-20 | 2020-10-30 | 帕诺森斯有限公司 | 使用本底噪声的自适应识别进行有效信号检测 |
CN108594256A (zh) * | 2018-04-16 | 2018-09-28 | 合肥菲涅尔光电科技有限公司 | 一种基于脉冲编码技术的相干激光雷达 |
CN108896456A (zh) * | 2018-04-28 | 2018-11-27 | 南京信息工程大学 | 基于反馈型rbf神经网络的气溶胶消光系数反演方法 |
CN108896456B (zh) * | 2018-04-28 | 2020-07-31 | 南京信息工程大学 | 基于反馈型rbf神经网络的气溶胶消光系数反演方法 |
CN109100009A (zh) * | 2018-06-01 | 2018-12-28 | 国网江苏省电力有限公司南京供电分公司 | 基于经验模态分解emd的分接开关振动信号降噪方法 |
CN108845306B (zh) * | 2018-07-05 | 2022-04-26 | 南京信息工程大学 | 基于变分模态分解的激光雷达回波信号去噪方法 |
CN108845306A (zh) * | 2018-07-05 | 2018-11-20 | 南京信息工程大学 | 基于变分模态分解的激光雷达回波信号去噪方法 |
CN109541549A (zh) * | 2018-10-09 | 2019-03-29 | 广东工业大学 | 基于emd和稀疏信号处理的间歇采样转发干扰抑制方法 |
CN109541549B (zh) * | 2018-10-09 | 2023-03-07 | 广东工业大学 | 基于emd和稀疏信号处理的间歇采样转发干扰抑制方法 |
CN109598152A (zh) * | 2018-10-11 | 2019-04-09 | 天津大学 | 基于emd降噪数据预处理的硬件木马检测优化方法 |
CN112230199A (zh) * | 2019-07-15 | 2021-01-15 | 天津大学 | 一种基于高维特征值分析的激光雷达回波盲去噪方法 |
CN112230199B (zh) * | 2019-07-15 | 2022-10-25 | 天津大学 | 一种基于高维特征值分析的激光雷达回波盲去噪方法 |
CN110564905A (zh) * | 2019-10-08 | 2019-12-13 | 中南大学 | 一种对高炉炉衬冲击回波检测的信号处理方法及*** |
WO2022061597A1 (zh) * | 2020-09-23 | 2022-03-31 | 深圳市速腾聚创科技有限公司 | 一种信号噪声滤除方法、装置、存储介质及激光雷达 |
WO2022217406A1 (zh) * | 2021-04-12 | 2022-10-20 | 深圳市速腾聚创科技有限公司 | 信号处理方法、装置及可读存储介质 |
CN116343051A (zh) * | 2023-05-29 | 2023-06-27 | 山东景闰工程研究设计有限公司 | 一种基于遥感影像的地质环境监测方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN107144829B (zh) | 2019-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107144829B (zh) | 一种高效的激光雷达回波信号去噪方法 | |
CN108226892B (zh) | 一种基于深度学习的复杂噪声环境下的雷达信号恢复方法 | |
Kaur et al. | A comparative analysis of thresholding and edge detection segmentation techniques | |
CN101661611B (zh) | 基于贝叶斯非局部均值滤波器的实现方法 | |
CN108845306A (zh) | 基于变分模态分解的激光雷达回波信号去噪方法 | |
CN102176001B (zh) | 一种基于透水波段比值因子的水深反演方法 | |
CN104614718B (zh) | 基于粒子群算法的激光雷达波形数据分解的方法 | |
CN109871733A (zh) | 一种自适应海杂波信号去噪方法 | |
CN103245976B (zh) | 基于uwb生物雷达的人体目标与周围环境结构兼容探测方法 | |
CN107167802A (zh) | 一种基于超宽带雷达的呼吸信号检测算法 | |
CN106296655A (zh) | 基于自适应权值和高频阈值的sar图像变化检测方法 | |
CN103456011A (zh) | 一种利用互补信息的改进型高光谱rx异常检测方法 | |
Gomez-Rodriguez et al. | Smoke monitoring and measurement using image processing: application to forest fires | |
CA2719142C (en) | Method for the three-dimensional synthetic reconstruction of objects exposed to an electromagnetic and/or elastic wave | |
CN102506444A (zh) | 基于智能控制的计算机视觉技术的炉膛火焰检测方法 | |
CN101685158B (zh) | 基于隐马尔科夫树模型的sar图像去噪方法 | |
CN103927737A (zh) | 基于非局部均值的sar图像变化检测方法 | |
CN103400383A (zh) | 基于nsct和压缩投影的sar图像变化检测方法 | |
CN104463808A (zh) | 基于空间相关性的高光谱数据降噪方法及*** | |
CN106225681A (zh) | 一种大跨度桥梁健康状态监测装置 | |
CN105158749A (zh) | 高频雷达海杂波幅度统计分布检验方法 | |
CN113887398A (zh) | 一种基于变分模态分解和奇异谱分析的gpr信号去噪方法 | |
Zhang et al. | Denoising method based on CNN-LSTM and CEEMD for LDV signals from accelerometer shock testing | |
CN105118035A (zh) | 基于稀疏表示的自适应光学光斑信号提取方法 | |
CN103169449B (zh) | 呼吸信号检测方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20201228 Address after: 215000 No.8, Xiangjie street, Suzhou high tech Zone, Suzhou City, Jiangsu Province Patentee after: Suzhou aikrypton inno Robot Technology Co.,Ltd. Address before: 210044 No. 219, Ning six road, Nanjing, Jiangsu Patentee before: NANJING University OF INFORMATION SCIENCE & TECHNOLOGY |
|
TR01 | Transfer of patent right |