CN104345049A - 激光诱导击穿光谱的小波阈值降噪的阈值校正方法 - Google Patents

Abstract

本发明公开了一种用于激光诱导击穿光谱的小波阈值降噪的阈值校正方法，基于合理的阈值校正数学模型，通过灰色理论和模糊理论计算得到校正后的阈值并用于激光诱导击穿光谱的小波阈值降噪。通过本方法取得了很好的激光诱导击穿光谱的降噪效果，在提高了信噪比的同时降低了检出限。

Description

激光诱导击穿光谱的小波阈值降噪的阈值校正方法

技术领域

本发明涉及激光诱导击穿光谱的信号预处理分析领域，具体是一种基于合理的阈值校正数学模型，通过灰色理论和模糊理论计算得到校正后的阈值并用于激光诱导击穿光谱的小波阈值降噪。

背景技术

激光诱导击穿光谱（LIBS）是一种典型的原子发射光谱测量技术。它利用聚焦的强脉冲激光将待测样品激发成等离子体而形成发射光谱，通过分析等离子体中原子或离子光谱来实现对样品的元素分析。LIBS发射强度的严重波动受到以冲击波形式的等离子传播的影响，并且受到来自多种噪声源，比如光电耦合器的暗电流，电子电路的热噪声以及光谱仪的杂散光和等离子体的连续辐射等。这些噪声夹杂在LIBS光谱信号中，不利于后续的样品元素分析。

为了解决LIBS光谱信号的降噪问题，人们通常采用的手段有两类：直接影响测量采集的硬件优化方式（设计滤波器、隔离器和探测器）以及LIBS光谱信号采集之后使强度增强的软件方法（设计数字滤波器、平滑滤波器和小波滤波器）。硬件优化方式由于其研发周期长和耗费大，造成了硬件优化手段的应用存在很大的局限性。软件方法不仅具有比硬件优化方式更高的精度而且还具有硬件优化方式不能比拟的可靠性，甚至能够实现硬件优化方式在理论上也无法达到的性能。软件滤波方法主要包括平滑滤波器、数字滤波器和小波滤波方法。平滑滤波器虽然能够保留LIBS信号的形状特征但是会导致LIBS信号在幅值上的损失较大，这不利于后续的LIBS定量分析。数字滤波器由于需要选择的合理滤波器参数较多会增加算法的复杂性并存在计算量大的不足。小波滤波方法主要包括三种：系数相关性方法、模极大值方法和阈值方法。系数相关性方法不仅需要定义直接影响到降噪效果的相关系数，而且需要迭代，计算量较大。模极大值方法存在一个由模极大值重构小波系数的问题，由于其重构算法的不同会导致降噪结果的不理想，此外还有算法复杂、计算量较大、收敛缓慢和稳定性较差的不足。阈值方法不仅具有需要选择的滤波参数（阈值和阈值函数）较少的优点，而且实现最简单、计算量最小。

发明内容

为了解决现有小波阈值降噪中选择阈值的困难所引起的降噪效果不理想的不足，本发明的目的在于提出了一种基于合理的阈值校正数学模型，通过灰色理论和模糊理论计算得到校正后的阈值并用于激光诱导击穿光谱的小波阈值降噪。

本发明为实现上述目的所采用的技术方案是：一种激光诱导击穿光谱的小波阈值降噪的阈值校正方法，包括以下步骤：

输入整个激光诱导击穿光谱信号，利用Shannon熵来选择小波基函数利用白噪声检验方法来确定小波阈值降噪的分解层数J；

根据所述小波基函数以及小波阈值降噪的分解层数J对激光诱导击穿光谱信号进行离散小波变换，得到每一个分解层j上的近似系数cA_j和细节系数cD_j，j=1,…,J；

在每一个分解层j上，基于Donoho阈值λ_j，根据灰色系统理论计算得到上阈值的修正参数γ_j；根据模糊系统理论计算小波系数模糊集合的隶属度：近似系数模糊集合的隶属度和细节系数模糊集合的隶属度

在每一个分解层j上，根据和计算得到下阈值的修正参数C_j；

在每一个分解层j上，根据上阈值的修正参数γ_j和下阈值λ_1j的修正参数C_j，计算得到新的上阈值λ_2j=γ_j·λ_j和新的下阈值λ_1j=C_j·λ_2j；

在每一个分解层j上，将新的上阈值λ_2j和新的下阈值λ_1j代入半软阈值函数T_j，计算得到校正后的细节系数

通过Mallet方法，利用cA_J和校正后的细节系数对激光诱导击穿光谱信号进行小波重构，完成激光诱导击穿光谱信号降噪过程。

所述的每一个分解层j上的近似系数cA_j={a_j,1,…,a_j,i,…,d_j,k}和细节系数cD_j={d_j,1,…,d_j,i,…,d_j,k},其中，i代表小波系数的序号，k代表小波系数的个数。

所述Donoho阈值的计算公式其中，σ_j是噪声方差的估计值Median(cD_j)表示在每一个分解层j上细节系数cD_j的中位数，L_j是每一个分解层j上细节系数cD_j的个数。

所述的上阈值的修正参数γ_j的计算公式：

${\mathrm{\γ}}_j=\frac{\underset{j}{\min }\underset{i}{\min }|a_{J,i}-d_{j,i}|+\mathrm{\ξ}\underset{j}{\max }\underset{i}{\max }|a_{J,i}-d_{j,i}|}{|a_{J,i}-d_{j,i}|+\mathrm{\ξ}\underset{j}{\max }\underset{i}{\max }|a_{J,i}-d_{j,i}|}$

其中，ξ是辨识系数，取ξ=0.5；当j≠J时，对a_J,k进行上抽样插零，其中，j代表分解层数，j=1,…,J；i代表每一分解层数上小波系数的序号，max代表极大值，min代表极小值。

所述小波系数模糊集合的隶属度的计算公式：

$\mathrm{\&mu;}\left(x\right)=\left\{\begin{array}{c}0,x\&le;p\\ \frac{1}{2}+\frac{1}{2}\sin \left[\frac{\mathrm{\&pi;}}{q-p}(x-\frac{p+q}{2})\right],p<x\&le;q\\ 1,x>q\end{array}\right.$

其中，p是隶属度函数μ(x)等于0时自变量x的取值，q是隶属度函数μ(x)等于1时自变量x的取值；

所述p和q的计算公式是：

$\left\{\begin{array}{c}p=\frac{t_1+t_2}{2}+\frac{\mathrm{\&pi;}(t_1-t_2)}{4\arcsin (2\mathrm{\&Delta;}-1)}\\ q=\frac{t_1+t_2}{2}-\frac{\mathrm{\&pi;}(t_1-t_2)}{4\arcsin (2\mathrm{\&Delta;}-1)}\end{array}\right.$

其中，t₁是小波系数模糊集合的绝对极小值，t₂是小波系数模糊集合的绝对极大值；⊿→0；

将|a_j,i|和|d_j,i|代入到μ(x)中就可以分别得到和的隶属度。

所述下阈值的修正参数C_j=1-N_j；其中，N_j是欧几里德贴近度，N_j的计算公式： $N_j(\stackrel{\&OverBar;}{A},\stackrel{\&OverBar;}{D})=1-\sqrt{\frac{1}{k}\underset{i=1}{\overset{k}{\mathrm{\&Sigma;}}}{[{\mathrm{\&mu;}}_{\stackrel{\&OverBar;}{A}}\left(\right|a_{j,i}\left|\right)-{\mathrm{\&mu;}}_{\stackrel{\&OverBar;}{A}}\left(\right|d_{j,i}\left|\right)]}^2};$ 其中，i代表小波系数的序号，k代表小波系数的个数。

所述半软阈值函数的计算公式：

其中，sgn（.）代表符号函数，j代表分解层数，j=1,…,J；λ_1j代表每一个分解层上的下阈值，λ_2j代表每一个分解层上的上阈值，cD_j代表每个分解层上的细节小波系数；代表每一个分解层上修正后的小波细节系数，T_j代表半软阈值函数。

所述小波重构是根据Mallet方法，利用cA_J和校正后的细节系数进行小波逆变换完成了小波信号重构。

本发明具有以下优点：

1、本发明激光诱导击穿光谱的小波阈值降噪的阈值校正方法建立了上阈值和下阈值校正的数学模型，不仅根据不同分解层数上的小波系数特征进行校正，而且还根据半软阈值函数的特点进行调节，提高了小波阈值选择的自适应性、准确性和可靠性。

2、本发明的方法不需要人为参与小波阈值的设定，极大地减小了人为因素的干扰，而且本方法算法简单、计算速度快，易于实施。

3、本发明的方法可以提高LIBS光谱信号的信噪比，降低了最低检出限。

附图说明

图1本发明工作流程图；

图2为本发明铜合金样品未进行降噪处理的局部信号谱图；

图3为本发明对铜合金样品进行小波阈值降噪的小波细节系数的结果图；

图4为本发明对铜合金样品进行小波阈值降噪前后的局部信号对比谱图。

具体实施方式

下面结合附图及实施例对本发明做进一步的详细说明。

如图1所示，本发明基于合理的阈值校正数学模型，通过灰色理论和模糊理论计算得到校正后的阈值并用于激光诱导击穿光谱的小波阈值降噪的工作流程为：

基于合理的阈值校正数学模型，通过灰色理论和模糊理论计算得到校正后的阈值并用于激光诱导击穿光谱的小波阈值降噪，具体步骤如下：

步骤1）输入整个LIBS光谱信息；

步骤2）利用Shannon熵来选择小波基函数利用白噪声检验方法来确定小波阈值降噪的分解层数J；

步骤3）根据2）中确定的小波基函数以及小波阈值降噪的分解层数J对1）中光谱信号进行离散小波变换得到每一个分解层j上的近似系数cA_j和细节系数cD_j（j=1,…,J）；

所述的每一个分解层j上的近似系数cA_j={a_j,1,…,a_j,i,…,d_j,k}和细节系数cD_j={d_j,1,…,d_j,i,…,d_j,k},其中，i代表小波系数的序号，k代表小波系数的个数。

步骤4）在每一个分解层j上，基于Donoho阈值λ_j，根据灰色系统理论计算得到上阈值的修正参数γ_j；

Donoho阈值的计算公式σ_j是噪声方差的估计值Median(cD_j)表示在每一个分解层j上细节系数cD_j的中位数，L_j是每一个分解层j上细节系数cD_j的个数。根据灰色系统理论，通过计算每一分解层j上的灰色关联度得到所述的上阈值的修正参数γ_j。所述的上阈值的修正参数γ_j的计算公式： ${\mathrm{\&gamma;}}_j=\frac{\underset{j}{\min }\underset{i}{\min }|a_{J,i}-d_{j,i}|+\mathrm{\&xi;}\underset{j}{\max }\underset{i}{\max }|a_{J,i}-d_{j,i}|}{|a_{J,i}-d_{j,i}|+\mathrm{\&xi;}\underset{j}{\max }\underset{i}{\max }|a_{J,i}-d_{j,i}|}.$ 其中，ξ是辨识系数，取ξ=0.5；当j≠J时，对a_J,k进行上抽样插零。

步骤5）在每一个分解层j上，根据模糊系统理论计算小波系数模糊集合的隶属度：近似系数模糊集合的隶属度μ_A(|a_j,i|)和细节系数模糊集合的隶属度

根据模糊系统理论，建立每一分解层j上的小波系数模糊集合：近似系数模糊集合和细节系数模糊集合计算近似系数模糊集合的隶属度μ_A(|a_j,i|)和细节系数模糊集合的隶属度的计算公式：

$\mathrm{\&mu;}\left(x\right)=\left\{\begin{array}{c}0,x\&le;p\\ \frac{1}{2}+\frac{1}{2}\sin \left[\frac{\mathrm{\&pi;}}{q-p}(x-\frac{p+q}{2})\right],p<x\&le;q\\ 1,x>q\end{array}\right..$ 其中，p是隶属度函数μ(x)等于0时自变量x的取值，q是隶属度函数μ(x)等于1时自变量x的取值。p和q的计算公式是： $\left\{\begin{array}{c}p=\frac{t_1+t_2}{2}+\frac{\mathrm{\&pi;}(t_1-t_2)}{4\arcsin (2\mathrm{\&Delta;}-1)}\\ q=\frac{t_1+t_2}{2}-\frac{\mathrm{\&pi;}(t_1-t_2)}{4\arcsin (2\mathrm{\&Delta;}-1)}\end{array}\right..$ 其中，t₁是小波系数模糊集合的绝对极小值，t₂是小波系数模糊集合的绝对极大值；⊿→0。

步骤6）在每一个分解层j上，根据5）中得到的和计算得到下阈值的修正参数C_j；

所述的下阈值的修正参数C_j=1-N_j。其中，N_j是欧几里德贴近度，N_j的计算公式：其中，i代表小波系数的序号，k代表小波系数的个数。

步骤7）在每一个分解层j上，根据4）中得到上阈值的修正参数γ_j和根据6）下阈值λ_1j的修正参数C_j，计算得到新的上阈值λ_2j=γ_j·λ_j和新的下阈值λ_1j=C_j·λ_2j；

步骤8）在每一个分解层j上，将新的上阈值λ_2j和新的下阈值λ_1j代入半软阈值函数T_j，计算得到校正后的细节系数

所述的半软阈值函数的计算公式：其中，sgn（.）代表符号函数。

步骤9）通过Mallet方法，利用3）中得到的cA_J和8）中得到的校正后的细节系数对激光诱导击穿光谱信号进行小波重构，完成激光诱导击穿光谱信号降噪过程。

图2、图3及图4给出了本实施例针对铜合金样品的测定过程。

请参阅图1，本发明的具体分析铜合金实施例步骤为：

步骤1）输入整个光谱信息；

本实施例中选择波长在220～375nm之间的光谱数据。图2中黑色实线是所选择的光谱局部信号。

步骤2）步骤2）利用Shannon熵来选择小波基函数利用白噪声检验方法来确定小波阈值降噪的分解层数J；

本实施例中确定的小波基函数是‘db5’，确定小波阈值降噪的分解层数J=4。

步骤3）根据2）中确定的小波基函数以及小波阈值降噪的分解层数J对1）中光谱信号进行离散小波变换得到每一个分解层j上的近似系数cA_j和细节系数cD_j（j=1,…,J）；

步骤4）在每一个分解层j上，基于Donoho阈值λ_j，根据灰色系统理论计算得到上阈值的修正参数γ_j；

本实施例中Donoho阈值λ₁=2.3982,λ₂=3.5038,λ₃=6.0059,λ₄=14.2357。

上阈值的修正参数γ₁=0.9996，γ₂=0.8878，γ₃=0.8179，γ₄=0.7948。

步骤5）在每一个分解层j上，根据模糊系统理论计算小波系数模糊集合的隶属度；

步骤6）在每一个分解层j上，根据5）中得到的μ_A(|a_j,i|)和计算得到下阈值的修正参数C_j；

本实施例中下阈值的修正参数C₁=0.4938,C₂=0.5625,C₃=0.5443,C₄=0.9342.

步骤7）在每一个分解层j上，根据4）中得到上阈值的修正参数γ_j和根据6）下阈值λ_1j的修正参数C_j，计算得到新的上阈值λ_2j=γ_j·λ_j和新的下阈值λ_1j=C_j·λ_2j；

本实施例中新的上阈值λ₂₁=2.3972，λ₂₂=3.1106，λ₂₃=4.9122，λ₂₄=11.3145；新的下阈值λ₁₁=1.1837，λ₁₂=1.7497，λ₁₃=2.6737，λ₁₄=10.5700。

步骤8）在每一个分解层j上，将新的上阈值λ_2j和新的下阈值λ_1j代入半软阈值函数T_j，计算得到校正后的细节系数

图3是对铜合金样品进行小波阈值降噪的小波细节系数的结果图；

步骤9）通过Mallet方法，利用3）中得到的cA_J和8）中得到的校正后的细节系数对激光诱导击穿光谱信号进行小波重构，完成激光诱导击穿光谱信号降噪过程。

图4中的红色实线是LIBS光谱信号降噪的最终结果。从图4的对比降噪前后的效果可以看出，本方法的选择小波阈值的效果很好，不仅很好的保留了原来的LIBS光谱信号的特征，而且提高了信噪比，降低了最低检出限。

Claims (8)

1.一种激光诱导击穿光谱的小波阈值降噪的阈值校正方法，其特征在于，包括以下步骤：

输入整个激光诱导击穿光谱信号，利用Shannon熵来选择小波基函数利用白噪声检验方法来确定小波阈值降噪的分解层数J；

根据所述小波基函数以及小波阈值降噪的分解层数J对激光诱导击穿光谱信号进行离散小波变换，得到每一个分解层j上的近似系数cA_j和细节系数cD_j，j=1,…,J；

在每一个分解层j上，基于Donoho阈值λ_j，根据灰色系统理论计算得到上阈值的修正参数γ_j；根据模糊系统理论计算小波系数模糊集合的隶属度：近似系数模糊集合的隶属度和细节系数模糊集合的隶属度

在每一个分解层j上，根据和计算得到下阈值的修正参数C_j；

在每一个分解层j上，根据上阈值的修正参数γ_j和下阈值λ_1j的修正参数C_j，计算得到新的上阈值λ_2j=γ_j·λ_j和新的下阈值λ_1j=C_j·λ_2j；

在每一个分解层j上，将新的上阈值λ_2j和新的下阈值λ_1j代入半软阈值函数T_j，计算得到校正后的细节系数

通过Mallet方法，利用cA_J和校正后的细节系数对激光诱导击穿光谱信号进行小波重构，完成激光诱导击穿光谱信号降噪过程。

2.根据权利要求1所述的激光诱导击穿光谱的小波阈值降噪的阈值校正方法，其特征在于，所述的每一个分解层j上的近似系数cA_j={a_j,1,…,a_j,i,…,d_j,k}和细节系数cD_j={d_j,1,…,d_j,i,…,d_j,k},其中，i代表小波系数的序号，k代表小波系数的个数。

3.根据权利要求1所述的激光诱导击穿光谱的小波阈值降噪的阈值校正方法，其特征在于，所述Donoho阈值的计算公式其中，σ_j是噪声方差的估计值Median(cD_j)表示在每一个分解层j上细节系数cD_j的中位数，L_j是每一个分解层j上细节系数cD_j的个数。

4.根据权利要求1所述的激光诱导击穿光谱的小波阈值降噪的阈值校正方法，其特征在于，所述的上阈值的修正参数γ_j的计算公式：

${\mathrm{\&gamma;}}_j=\frac{\underset{j}{\min }\underset{i}{\min }|a_{J,i}-d_{j,i}|+\mathrm{\&xi;}\underset{j}{\max }\underset{i}{\max }|a_{J,i}-d_{j,i}|}{|a_{J,i}-d_{j,i}|+\mathrm{\&xi;}\underset{j}{\max }\underset{i}{\max }|a_{J,i}-d_{j,i}|}$

其中，ξ是辨识系数，取ξ=0.5；当j≠J时，对a_J,k进行上抽样插零，其中，j代表分解层数，j=1,…,J；i代表每一分解层数上小波系数的序号，max代表极大值，min代表极小值。

5.根据权利要求1所述的激光诱导击穿光谱的小波阈值降噪的阈值校正方法，其特征在于，所述小波系数模糊集合的隶属度的计算公式：

$\mathrm{\&mu;}\left(x\right)=\left\{\begin{array}{c}0,x\&le;p\\ \frac{1}{2}+\frac{1}{2}\sin \left[\frac{\mathrm{\&pi;}}{q-p}(x-\frac{p+q}{2})\right],p<x\&le;q\\ 1,x>q\end{array}\right.$

其中，p是隶属度函数μ(x)等于0时自变量x的取值，q是隶属度函数μ(x)等于1时自变量x的取值；

所述p和q的计算公式是：

$\left\{\begin{array}{c}p=\frac{t_1+t_2}{2}+\frac{\mathrm{\&pi;}(t_1-t_2)}{4\arcsin (2\mathrm{\&Delta;}-1)}\\ q=\frac{t_1+t_2}{2}-\frac{\mathrm{\&pi;}(t_1-t_2)}{4\arcsin (2\mathrm{\&Delta;}-1)}\end{array}\right.$

其中，t₁是小波系数模糊集合的绝对极小值，t₂是小波系数模糊集合的绝对极大值；⊿→0；

将|a_j,i|和|d_j,i|代入到μ(x)中就可以分别得到和的隶属度。

6.根据权利要求1所述的激光诱导击穿光谱的小波阈值降噪的阈值校正方法，其特征在于，所述下阈值的修正参数C_j=1-N_j；其中，N_j是欧几里德贴近度，N_j的计算公式： $N_j(\stackrel{\&OverBar;}{A},\stackrel{\&OverBar;}{D})=1-\sqrt{\frac{1}{k}\underset{i=1}{\overset{k}{\mathrm{\&Sigma;}}}{[{\mathrm{\&mu;}}_{\stackrel{\&OverBar;}{A}}\left(\right|a_{j,i}\left|\right)-{\mathrm{\&mu;}}_{\stackrel{\&OverBar;}{A}}\left(\right|d_{j,i}\left|\right)]}^2};$ 其中，i代表小波系数的序号，k代表小波系数的个数。

7.根据权利要求1所述的激光诱导击穿光谱的小波阈值降噪的阈值校正方法，其特征在于，所述半软阈值函数的计算公式：

其中，sgn（.）代表符号函数，j代表分解层数，j=1,…,J；λ_1j代表每一个分解层上的下阈值，λ_2j代表每一个分解层上的上阈值，cD_j代表每个分解层上的细节小波系数；代表每一个分解层上修正后的小波细节系数，T_j代表半软阈值函数。

8.根据权利要求1所述的激光诱导击穿光谱的小波阈值降噪的阈值校正方法，其特征在于，所述小波重构是根据Mallet方法，利用cA_J和校正后的细节系数进行小波逆变换完成了小波信号重构。