CN109883547B - 一种基于小波阈值差分的宽波段光谱信号去噪方法 - Google Patents
一种基于小波阈值差分的宽波段光谱信号去噪方法 Download PDFInfo
- Publication number
- CN109883547B CN109883547B CN201910023112.1A CN201910023112A CN109883547B CN 109883547 B CN109883547 B CN 109883547B CN 201910023112 A CN201910023112 A CN 201910023112A CN 109883547 B CN109883547 B CN 109883547B
- Authority
- CN
- China
- Prior art keywords
- wavelet
- layer
- denoised
- matrix
- singular value
- 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.)
- Active
Links
- 230000003595 spectral effect Effects 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000011159 matrix material Substances 0.000 claims abstract description 75
- 238000001228 spectrum Methods 0.000 claims abstract description 40
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 238000012545 processing Methods 0.000 claims abstract description 12
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 230000008569 process Effects 0.000 claims description 15
- 230000000694 effects Effects 0.000 abstract description 6
- 238000005259 measurement Methods 0.000 abstract description 5
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 5
- 229910052753 mercury Inorganic materials 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 3
- YZCKVEUIGOORGS-OUBTZVSYSA-N Deuterium Chemical compound [2H] YZCKVEUIGOORGS-OUBTZVSYSA-N 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 229910052805 deuterium Inorganic materials 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005057 refrigeration Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000005654 stationary process Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
- Complex Calculations (AREA)
Abstract
本发明提出了一种基于小波阈值差分的宽波段光谱信号去噪方法。对含噪光谱信号进行小波分解得到多层小波系数；对多层小波系数行软阈值去噪得到去噪后多层小波系数；通过去噪后每层小波系数构造成Hankel矩阵；根据Hankel矩阵构造奇异值矩阵，通过对奇异值矩阵正交分解得到奇异值；对奇异值进行差分计算构造奇异值差分谱；对奇异值差分谱进行奇异值差分去噪处理；对去噪后的奇异值矩阵进行重构运算得到去噪后Hankel矩阵，以还原为去噪后小波系数，将去噪后小波系数通过Mallat算法进行小波逆变换得到去噪后的光谱信号。本发明去噪效果彻底且能保留光谱有效特征，提高了宽波段光谱仪测量精度和准确度。
Description
技术领域
本发明涉及宽波段光谱信号的去噪方法技术领域，具体地，涉及一种基于小 波阈值差分的宽波段光谱信号去噪方法。
背景技术
现有宽波段微型光谱仪存在两个主要问题，一是光谱范围与分辨率之间存在 限制关系，较难在较宽的光谱范围内取得较为均匀且较高的光谱分辨率；二是光 谱范围包含紫外波段区域，由于CCD探测器的性能会造成仪器紫外响应灵敏度 不高。而合适的光谱数据处理算法可有效改善这两个问题。
光谱数据的处理方法决定了输出光谱的准确性。数据处理过程包括暗信号的 校正、光谱消噪、寻峰、波长和光强定标。其中，暗信号校正是为了去除测量环 境中的暗光噪声等。去噪过程贯穿于光路、电路和软件研发过程的始终。寻峰算 法可获得光谱峰值及其相对应的探测器光敏元序号。此外，采用波长定标算法可 以确定像元序号与光谱波长数值的函数关系，得到光谱仪的波长定标方程；光强 定标是指对不同波段光强的拟合。其中，光谱信号去噪算法的适用性在很大程度 上会决定最终光谱图的精确性。
现有的去噪方法主要有硬件去噪和软件去噪，硬件去噪通常有模拟前端设计 去噪和控温去噪，通过对探测器采取制冷措施和对输出模拟信号进行信号处理实 现噪声抑制。硬件去噪可以在一定程度上抑制噪声，但在实际的工程应用中并不 能彻底消除随机噪声。而且在信号处理的过程中，硬件电路会引入新的随机噪声， 因而，软件去噪是保证光谱仪测量精度必不可少且非常重要的过程。
现有的软件去噪方法主要有傅里叶变换、光谱信号光滑去噪、邻域平均滤波、 中值滤波和小波变换等。其中，小波分析被誉为数学显微镜，对信号有自适应性。 其时间窗口和频率窗口的形状都可改变，具有较好的时频局部化。由于其良好的 时频特性，广泛应用在图像处理和非平稳过程分析、数字滤波等领域中，其表征 信号局部特征的多分辨率分析方法十分适用于光谱信号分析。然而由于小波分解 实质上是对信号进行滤波运算，小波滤波器的性能决定其分解效果，而小波滤波 器两相邻尺度细节信号的频率窗几乎都存在一定的重叠。另外小波分解后各尺度 细节信号频带比较宽，导致有用信号和随机噪声很难准确又完全的分离，去噪存 在信号边缘模糊以及信号特征丢失等现象，且小波阈值去噪结果存在去噪不彻底 或者去噪过度导致信号失真的问题，效果也不尽人意。因而，深入研究光谱仪软 件系统内宽波段光谱去噪方法对于宽波段光谱仪测量精度和准确度的提升具有重要意义。
发明内容
为了解决上述存在的技术问题，本发明的目的是提供一种基于小波阈值差分 的宽波段光谱信号去噪方法，以解决在较宽的光谱范围内CCD对紫外响应灵敏 度不高和光谱仪中多种噪声的去噪的问题。并且通过了相关的对比试验，证明了 该方法去噪彻底且能保留光谱的有效特征，适用于宽波段光谱仪的去噪过程。
本发明的一种基于小波阈值差分的宽波段光谱信号去噪方法，它是通过对原 始光谱信号先应用小波阈值去噪，获得两层小波系数，然后将消除噪声的每层小 波系数均构造成Hankel矩阵，对每个Hankel矩阵均进行奇异值正交分解(SVD)， 选择合适的去噪阶数进行处理，再对处理后结果进行奇异值重构，得到新的 Hankel，再对恢复出去噪后的光谱信号。最终得到理想的去噪结果。
本发明有以下步骤
步骤1：对含噪光谱信号进行小波分解得到多层小波系数；
步骤2：对多层小波系数行软阈值去噪得到去噪后多层小波系数；
步骤3：通过去噪后每层小波系数构造成Hankel矩阵；
步骤4：根据Hankel矩阵构造奇异值矩阵，通过对奇异值矩阵正交分解得 到奇异值；
步骤5：对奇异值进行差分计算构造奇异值差分谱；
步骤6：对奇异值差分谱进行奇异值差分去噪处理；
步骤7：对去噪后的奇异值矩阵进行重构运算得到去噪后Hankel矩阵，以 还原为去噪后小波系数，将去噪后小波系数通过Mallat算法进行小波逆变换得到 去噪后的光谱信号；
作为优选，步骤1中所述对含噪光谱信号进行小波分解得到多层小波系数为：
对含噪光谱信号f(t)(t＝1,2,3,…,N,是离散时间序列号)进行小波分解， 选取小波基为db20，按高频信号到低频信号的顺序，分解为M层小波系数，即： D1，D2，...，DM；
第i层第j个小波系数为di,j，第i层第j个近似信号部分表示为ci,j，其中， i表示分解的层数，i∈[1,M]，j表示第i层第j个系数，j∈[1,N]；H(x)为小波 基函数db20对应的时域小波分解低通滤波器，G(x)为小波基函数db20对应的 时域小波分解高通滤波器；
作为优选，步骤2中所述对多层小波系数行软阈值去噪为：
对每一层小波系数进行软阈值去噪，软阈值的计算公式为：
其中，M为小波系数的层数，N为每层小波系数的数量，δi表示第i层的阈 值，λi表示第i层附加噪声的标准差，其计算公式为：
其中，M为小波系数的层数，N为每层小波系数的数量，将第i层N个小波 系数di,j按从小到大的顺序排序，medium(|di,j|)取第i层Ni个小波系数中间相邻 两项的平均值；
软阈值法的计算公式为：
其中，δi定义为第i层小波系数的阈值，sign(di,j)为符号函数，定义为：
经过阈值去噪后得到第i层第j个新的小波系数d′i,j；
作为优选，步骤3中所述Hankel矩阵的表达式为：
其中，Ai为第i层小波系数构造成的Hankel矩阵，N为含噪信号长度，即 为每层小波系数的数量，mi是矩阵Ai的行数，ni为矩阵Ai的列数，mi+ni＝N+1， ni＜mi，即ni＜(N+1)/2；
作为优选，步骤4中所述奇异值矩阵可以表示为：
其中，Ui为Ai的左奇异矩阵，Vi是Ai的右奇异矩阵，Ui和Vi都是正交矩阵， 且Ui和Vi的列分别是矩阵AiAi T和Ai TAi的特征向量；
步骤4中所述对奇异值矩阵正交分解得到奇异值为：
作为优选，步骤5中所述奇异值差分谱定义为
bi,k＝σi,k-σi,k+1i＝1,2,…,M k＝1,2,…,pi-1
其中，则序列Bi＝(bi,1，bi,2，...，bi,p-1)称为矩阵Ai的奇异值差分谱，bi,k为矩 阵Ai的第k个奇异值差分谱点；
作为优选，步骤6中所述对奇异值差分谱进行奇异值差分去噪处理为：
通过比较找到差分谱Bi＝(bi,1，bi,2，...，bi,p-1)的峰值，即满足bi,r＞bi,r-1且 bi,r＞bi,r+1，将所有的bi,r对应的奇异值σi,r作为噪声信号的奇异值，进行置0处理， 即令σi,r＝0；
将峰值全部去掉即可完成奇异值去噪过程，得到去噪后的奇异值矩阵S′i：
S′i＝diag(σ′i，1，σ′i，2，…，σ′i，p)i＝1，2，…，M，k＝1，2，…，pi-1
其中，σ′i,1,σ′i,2,…,σ′i,p为去噪后的奇异值，M为小波系数的层数；
作为优选，步骤7中所述进行重构运算得到去噪后的Hankel矩阵为：
A′i＝UiS′iVi T
其中，A′i为第i层小波系数去噪后的Hankel矩阵，则矩阵A′i的每个元素可 以看做是第i层去噪后小波系数：
其中，i表示分解的层数，i∈[1,M]，j表示第i层第j个系数，j∈[1,N],第 i+1层第j个去噪后小波系数为d″i+1,j，第i+1层第j个近似信号部分表示为ci+1,j， 其中，i表示分解的层数，i∈[1,M]，j表示第i层第j个系数，j∈[1,N]；
对每一层A′i矩阵的元素d″i,j通过Mallat算法进行小波逆变换，即可得到去噪 后的光谱信号；
Mallat算法的小波逆变换可以表示为：
其中，d″i+1,j为第i+1层第j个去噪后小波系数，ci+1,j为第i+1层第j个近似 信号部分表示，c′i,j为第i层第j个去噪后近似信号部分表示，i表示分解的层数， i∈[1,M]，j表示第i层第j个系数，j∈[1,N],H(x)为小波基函数db20对应的 时域小波分解低通滤波器，G(x)为小波基函数db20对应的时域小波分解高通滤 波器，f′(j)为去噪后的光谱信号。
与现有技术相比，本发明具有以下优点：
本发明基于小波阈值—奇异值差分的宽波段光谱信号去噪方法，使用小波阈 值对高频信号的低层小波系数进行消噪处理，对高频信号的全部小波系数进行奇 异值分解，并采用差分法选择奇异值去噪阶次，实现光谱的去噪过程。相较于单 纯的小波默认阈值去噪，谱线峰值均存在，且削减程度很低，不存在明显的信号 失真；相较于单纯的奇异值正交分解算法，去噪效果彻底，信噪比高，且峰值谱 线削减较少，未出现振铃等信号失真。该方法去噪效果彻底且能保留光谱有效特 征，基本适用于宽波段光谱仪的去噪过程，提高了宽波段光谱仪测量精度和准确 度。
附图说明
图1：为本发明方法流程图；
图2：为对含噪氘灯信号的去噪结果与原始信号对比；
图3：为对汞灯光谱信号的去噪结果与标准光谱信号对比。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明，下面结合附图及实施例对 本发明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解 释本发明，并不用于限定本发明。
下面结合图1至图3介绍本发明的实施方式，具体为：
步骤1：对含噪光谱信号进行小波分解得到多层小波系数；
步骤1中所述对含噪光谱信号进行小波分解得到多层小波系数为：
对含噪光谱信号f(t)(t＝1,2,3,…,N,是离散时间序列号)进行小波分解， 选取小波基为db20，按高频信号到低频信号的顺序，分解为M层小波系数，即： D1，D2，...，DM；
第i层第j个小波系数为di,j，第i层第j个近似信号部分表示为ci,j，其中， i表示分解的层数，i∈[1,M]，j表示第i层第j个系数，j∈[1,N]；H(x)为小波 基函数db20对应的时域小波分解低通滤波器，G(x)为小波基函数db20对应的 时域小波分解高通滤波器，M＝5,N＝32；
步骤2：对多层小波系数行软阈值去噪得到去噪后多层小波系数；
步骤2中所述对多层小波系数行软阈值去噪为：
对每一层小波系数进行软阈值去噪，软阈值的计算公式为：
其中，M为小波系数的层数，N为每层小波系数的数量，M＝5,N＝32，δi表 示第i层的阈值，λi表示第i层附加噪声的标准差，其计算公式为：
其中，M为小波系数的层数，N为每层小波系数的数量，将第i层N个小波 系数di,j按从小到大的顺序排序，medium(|di,j|)取第i层Ni个小波系数中间相邻 两项的平均值；
软阈值法的计算公式为：
其中，δi定义为第i层小波系数的阈值，sign(di,j)为符号函数，定义为：
经过阈值去噪后得到第i层第j个新的小波系数d′i,j；
步骤3：通过去噪后每层小波系数构造成Hankel矩阵；
步骤3中所述Hankel矩阵的表达式为：
其中，Ai为第i层小波系数构造成的Hankel矩阵，N为含噪信号长度，N＝32 即为每层小波系数的数量，mi是矩阵Ai的行数，ni为矩阵Ai的列数，mi+ni＝N+1，ni＜mi，即ni＜(N+1)/2；
步骤4：根据Hankel矩阵构造奇异值矩阵，通过对奇异值矩阵正交分解得 到奇异值；
步骤4中所述奇异值矩阵可以表示为：
其中，Ui为Ai的左奇异矩阵，Vi是Ai的右奇异矩阵，Ui和Vi都是正交矩阵， 且Ui和Vi的列分别是矩阵AiAi T和Ai TAi的特征向量；
步骤4中所述对奇异值矩阵正交分解得到奇异值为：
步骤5：对奇异值进行差分计算构造奇异值差分谱；
作为优选，步骤5中所述奇异值差分谱定义为
bi,k＝σi,k-σi,k+1i＝1,2,…,M k＝1,2,…,pi-1
其中，则序列Bi＝(bi,1，bi,2，...，bi,p-1)称为矩阵Ai的奇异值差分谱，bi,k为矩 阵Ai的第k个奇异值差分谱点，M＝5；
步骤6：对奇异值差分谱进行奇异值差分去噪处理；
作为优选，步骤6中所述对奇异值差分谱进行奇异值差分去噪处理为：
通过比较找到差分谱Bi＝(bi,1，bi,2，...，bi,p-1)的峰值，即满足bi,r＞bi,r-1且 bi,r＞bi,r+1，将所有的bi,r对应的奇异值σi,r作为噪声信号的奇异值，进行置0处理， 即令σi,r＝0；
将峰值全部去掉即可完成奇异值去噪过程，得到去噪后的奇异值矩阵S′i：
S′i＝diag(σ′i，1，σ′i，2，…，σ′i，p)i＝1，2，…，M，k＝1，2，…，pi-1
其中，σ′i,1,σ′i,2,…,σ′i,p为去噪后的奇异值，M为小波系数的层数；
步骤7：对去噪后的奇异值矩阵进行重构运算得到去噪后Hankel矩阵，以 还原为去噪后小波系数，将去噪后小波系数通过Mallat算法进行小波逆变换得到 去噪后的光谱信号；
作为优选，步骤7中所述进行重构运算得到去噪后的Hankel矩阵为：
A′i＝UiS′iVi T
其中，A′i为第i层小波系数去噪后的Hankel矩阵，则矩阵A′i的每个元素可 以看做是第i层去噪后小波系数：
其中，M＝5,N＝32，i表示分解的层数，i∈[1,M]，j表示第i层第j个系数， j∈[1,N],第i+1层第j个去噪后小波系数为d″i+1,j，第i+1层第j个近似信号部分 表示为ci+1,j，其中，i表示分解的层数，i∈[1,M]，j表示第i层第j个系数， j∈[1,N]；
对每一层A′i矩阵的元素d″i,j通过Mallat算法进行小波逆变换，即可得到去噪 后的光谱信号；
Mallat算法的小波逆变换可以表示为：
其中，d″i+1,j为第i+1层第j个去噪后小波系数，ci+1,j为第i+1层第j个近似 信号部分表示，c′i,j为第i层第j个去噪后近似信号部分表示，i表示分解的层数， i∈[1,M]，j表示第i层第j个系数，j∈[1,N],H(x)为小波基函数db20对应的 时域小波分解低通滤波器，G(x)为小波基函数db20对应的时域小波分解高通滤 波器，f′(j)为去噪后的光谱信号。
图1是本发明方法流程图，具体步骤如上面所述。
图2是对含噪氘灯信号的去噪结果与原始信号对比，487.50nm、582.23nm 和656.05nm处谱线峰值均存在，且削减程度很低，不存在明显的信号失真，去 噪效果彻底，信噪比高，且峰值谱线削减较少，未出现振铃等信号失真，实验结 果基本能反映光谱信号的基本特征。
图3是在实验室条件下，使用宽波段微型光谱仪对汞灯信号进行采集，使用 小波阈值-奇异值差分去噪算法对其进行去噪处理，结果如图3所示。实验结果 中，低压汞灯光谱图的6条特征谱线：253.65nm、365.01nm、404.66nm、435.84nm、 546.07nm、576.96nm均可清晰识别，与汞灯标准光谱对比，差别较小，且噪声 含量很少，即能够较准确得到汞灯光谱谱线。
应当理解的是，上述针对较佳实施例的描述较为详细，并不能因此而认为是 对本发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不 脱离本发明权利要求所保护的范围情况下，还可以做出替换或变形，均落入本发 明的保护范围之内，本发明的请求保护范围应以所附权利要求为准。
Claims (1)
1.一种基于小波阈值差分的宽波段光谱信号去噪方法，其特征在于，包括：
步骤1：对含噪光谱信号进行小波分解得到多层小波系数；
步骤2：对多层小波系数进行软阈值去噪得到去噪后多层小波系数；
步骤3：通过去噪后每层小波系数构造成Hankel矩阵；
步骤4：根据Hankel矩阵构造奇异值矩阵，通过对奇异值矩阵正交分解得到奇异值；
步骤5：对奇异值进行差分计算构造奇异值差分谱；
步骤6：对奇异值差分谱进行奇异值差分去噪处理；
步骤7：对去噪后的奇异值矩阵进行重构运算得到去噪后Hankel矩阵，以还原为去噪后小波系数，将去噪后小波系数通过Mallat算法进行小波逆变换得到去噪后的光谱信号；
步骤1中所述对含噪光谱信号进行小波分解得到多层小波系数为：
对含噪光谱信号f(t)(t＝1,2,3,…,N,是离散时间序列号)进行小波分解，选取小波基为db20，按高频信号到低频信号的顺序，分解为M层小波系数，即：D1，D2，...，DM；
第i层第j个小波系数为di,j，第i层第j个近似信号部分表示为ci,j，其中，i表示分解的层数，i∈[1,M]，j表示第i层第j个系数，j∈[1,N]；H(x)为小波基函数db20对应的时域小波分解低通滤波器，G(x)为小波基函数db20对应的时域小波分解高通滤波器；
步骤2中所述对多层小波系数进行软阈值去噪为：
对每一层小波系数进行软阈值去噪，软阈值的计算公式为：
其中，M为小波系数的层数，N为每层小波系数的数量，δi表示第i层的阈值，λi表示第i层附加噪声的标准差，其计算公式为：
其中，M为小波系数的层数，N为每层小波系数的数量，将第i层N个小波系数di,j按从小到大的顺序排序，medium(|di,j|)取第i层Ni个小波系数中间相邻两项的平均值；
软阈值法的计算公式为：
其中，δi定义为第i层小波系数的阈值，sign(di,j)为符号函数，定义为：
经过软阈值法去噪后得到第i层第j个新的小波系数d′i,j；
步骤3中所述Hankel矩阵的表达式为：
其中，Ai为第i层小波系数构造成的Hankel矩阵，N为含噪信号长度，即为每层小波系数的数量，mi是矩阵Ai的行数，ni为矩阵Ai的列数，mi+ni＝N+1，ni＜mi，即ni＜(N+1)/2；
步骤4中所述奇异值矩阵可以表示为：
其中，Ui为Ai的左奇异矩阵，Vi是Ai的右奇异矩阵，Ui和Vi都是正交矩阵，且Ui和Vi的列分别是矩阵AiAi T和Ai TAi的特征向量；
步骤4中所述对奇异值矩阵正交分解得到奇异值为：
步骤5中所述奇异值差分谱定义为
bi,k＝σi,k-σi,k+1
其中，i＝1,2,...,M,k＝1,2,...,pi-1，则序列
称为矩阵Ai的奇异值差分谱，bi,k为矩阵Ai的第k个奇异值差分谱点，pi为Ai的奇异值数量；
步骤6中所述对奇异值差分谱进行奇异值差分去噪处理为：
通过比较找到差分谱的峰值，即满足bi,r＞bi,r-1且bi,r＞bi,r+1，将所有的bi,r对应的奇异值σi,r作为噪声信号的奇异值，进行置0处理，即令σi,r＝0；
将峰值全部去掉即可完成奇异值去噪过程，得到去噪后的奇异值矩阵S′i：
其中，
为去噪后的奇异值，M为小波系数的层数，pi为Ai的奇异值数量；
步骤7中所述进行重构运算得到去噪后的Hankel矩阵为：
A′i＝UiS′iVi T
其中，A′i为第i层小波系数去噪后的Hankel矩阵，则矩阵A′i的每个元素可以看做是第i层去噪后小波系数：
其中，i表示分解的层数，i∈[1,M]，j表示第i层第j个系数，j∈[1,N],第i+1层第j个去噪后小波系数为d″i+1,j，第i+1层第j个近似信号部分表示为ci+1,j，其中，i表示分解的层数，i∈[1,M]，j表示第i层第j个系数，j∈[1,N]；
对每一层A′i矩阵的元素d″i,j通过Mallat算法进行小波逆变换，即可得到去噪后的光谱信号；
Mallat算法的小波逆变换可以表示为：
其中，d″i+1,j为第i+1层第j个去噪后小波系数，ci+1,j为第i+1层第j个近似信号部分表示，c′i,j为第i层第j个去噪后近似信号部分表示，i表示分解的层数，i∈[1,M]，j表示第i层第j个系数，j∈[1,N],H(x)为小波基函数db20对应的时域小波分解低通滤波器，G(x)为小波基函数db20对应的时域小波分解高通滤波器，f′(j)为去噪后的光谱信号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910023112.1A CN109883547B (zh) | 2019-01-10 | 2019-01-10 | 一种基于小波阈值差分的宽波段光谱信号去噪方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910023112.1A CN109883547B (zh) | 2019-01-10 | 2019-01-10 | 一种基于小波阈值差分的宽波段光谱信号去噪方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109883547A CN109883547A (zh) | 2019-06-14 |
| CN109883547B true CN109883547B (zh) | 2021-03-16 |
Family
ID=66925873
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910023112.1A Active CN109883547B (zh) | 2019-01-10 | 2019-01-10 | 一种基于小波阈值差分的宽波段光谱信号去噪方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109883547B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110503060B (zh) * | 2019-08-28 | 2022-05-10 | 中南大学 | 一种光谱信号去噪方法及其系统 |
| CN111949919B (zh) * | 2020-06-14 | 2022-04-01 | 杭州谱育科技发展有限公司 | 降低aes检测中噪声的方法 |
| CN112098774A (zh) * | 2020-09-10 | 2020-12-18 | 厦门理工学院 | 一种剩余电流故障判定方法、终端设备及存储介质 |
| CN113324918B (zh) * | 2021-01-26 | 2023-01-24 | 核工业北京地质研究院 | 一种岩石光谱去噪方法 |
| CN112998704A (zh) * | 2021-03-04 | 2021-06-22 | 研和智能科技(杭州)有限公司 | 一种可穿戴设备血氧饱和度计算方法 |
| CN113569731A (zh) * | 2021-07-27 | 2021-10-29 | 广东电网有限责任公司 | 卫星影像杆塔智能识别方法、系统、设备及介质 |
| CN113533238B (zh) * | 2021-09-15 | 2021-12-28 | 武汉敢为科技有限公司 | 一种基于吸收光谱的六氟化硫分解气体的检测方法及系统 |
| CN114004254A (zh) * | 2021-10-27 | 2022-02-01 | 云南电网有限责任公司昆明供电局 | 电力电缆局部放电信号白噪声及脉冲型干扰滤除算法 |
| CN116541657A (zh) * | 2023-04-26 | 2023-08-04 | 湖北工业大学 | 基于小波变换和优化奇异值分解的紫外光谱去噪方法、系统、设备及介质 |
| CN116499469B (zh) * | 2023-06-28 | 2023-09-08 | 北京航空航天大学 | 利用神经网络模型在线学习和补偿的gnss/ins组合导航方法 |
| CN117072155A (zh) * | 2023-10-16 | 2023-11-17 | 克拉玛依市红都有限责任公司 | 一种应用于石油开采的温度监测方法、系统、设备及介质 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8759778B2 (en) * | 2007-09-27 | 2014-06-24 | Anis Rahman | Terahertz time domain and frequency domain spectroscopy |
| CN103217409B (zh) * | 2013-03-22 | 2015-02-18 | 中国科学院重庆绿色智能技术研究院 | 一种拉曼光谱预处理方法 |
| CN106404178B (zh) * | 2016-08-30 | 2018-12-04 | 中国科学院地理科学与资源研究 | 基于小波变换的高光谱热红外地表温度与发射率分离方法 |
| CN106885978B (zh) * | 2017-04-20 | 2019-04-16 | 重庆大学 | 一种基于绝缘油拉曼光谱小波包能量熵的油纸绝缘老化诊断方法 |
-
2019
- 2019-01-10 CN CN201910023112.1A patent/CN109883547B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109883547A (zh) | 2019-06-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109883547B (zh) | 一种基于小波阈值差分的宽波段光谱信号去噪方法 | |
| CN112906578B (zh) | 一种桥梁时间序列位移信号去噪的方法 | |
| CN104182946A (zh) | 基于指数软硬择中阈值函数的光谱去噪方法 | |
| CN107918710B (zh) | 基于凸优化的非下采样图滤波器组的设计方法 | |
| CN113238193B (zh) | 一种多分量联合重构的sar回波宽带干扰抑制方法 | |
| Kober | Robust and efficient algorithm of image enhancement | |
| CN111242854A (zh) | 一种图像去噪方法 | |
| Zhu et al. | Deconvolution methods based on convex regularization for spectral resolution enhancement | |
| CN110032968B (zh) | 基于双树复小波和自适应半软阈值法的去噪方法 | |
| Duijster et al. | Wavelet-based EM algorithm for multispectral-image restoration | |
| CN109724693B (zh) | 一种基于平稳小波的融合光谱去噪方法 | |
| CN112730712B (zh) | 提高lc-ms数据信噪比的方法 | |
| CN105187341A (zh) | 一种基于交叉验证的平稳小波变换去噪方法 | |
| Özaydın et al. | Speech enhancement using maximal overlap discrete wavelet transform | |
| CN108133203B (zh) | 一种基于稀疏表示的拉曼光谱特征提取方法 | |
| CN113095113A (zh) | 一种用于水下目标识别的小波线谱特征提取方法及系统 | |
| CN113109873B (zh) | 一种基于秩残差约束的沙漠地区地震信号噪声抑制方法 | |
| CN115795272A (zh) | 一种基于分数阶迭代离散小波变换的谱线去噪方法 | |
| Shi et al. | Restoration of atmospheric turbulence degraded images | |
| CN111161128A (zh) | 基于频域方向滤波的图像变换及其在稀疏分解中的应用 | |
| Lou et al. | An approach based on simplified KLT and wavelet transform for enhancing speech degraded by non-stationary wideband noise | |
| Xu et al. | Random seismic noise attenuation by learning-type overcomplete dictionary based on K-singular value decomposition algorithm | |
| CN109346097B (zh) | 一种基于Kullback-Leibler差异的语音增强方法 | |
| Panisetti et al. | An unbiased risk estimator for multiplicative noise—Application to 1-D signal denoising | |
| Tan et al. | The Research on Denoising of SAR Image Based on Improved K-SVD Algorithm |
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 |