[发明专利]基于指数软硬择中阈值函数的光谱去噪方法

说明书

技术领域

本发明属于图像处理技术领域，更进一步涉及基于指数软硬择中阈值函数的光谱去噪方法，可作为用于CCD全谱直读光谱仪采集的光谱数据的降噪处理。

背景技术

全谱直读光谱仪是钢铁行业的主要分析检测校准仪器。根据光谱检测器件的不同，目前市场上的全谱直读光谱仪主要有两种：基于光电倍增管(PMT)和基于电荷耦合器件(Charge-coupled Device，CCD)。其中基于CCD的全谱直读光谱仪具有体积小，重量轻，功耗低，使用方便，单次采集光谱范围宽，响应速度快，灵敏度高，数据采集处理快捷方便等特点。因此，采用线阵CCD实现全谱直读光谱分析仪是一个发展趋势。在光谱数据的定量分析中，噪声的存在往往会影响分析结果的准确性，为了提高光谱数据的分析精度，对线阵CCD采集的光谱数据进行降噪处理是十分必要。

根据图像光谱的特点和噪声的统计特性，发展了各种各样的去噪方法。常用的去噪算法有很多种，比如均值滤波、平滑滤波、窗函数法、自适应中值滤波、中值滤波和小波变换滤波。尽管去噪方法多种多样，但降噪研究的工作都是在去噪和保留高频信息之间进行权衡。总体来说减少信号噪声的方法可以在信号空间域和在信号变换域内完成。在全谱直读光谱仪中，线阵CCD采集到的光谱数据含有以下噪声：由杂散光引起的低频背景噪声和由暗电流噪声、转移噪声所组成的高频噪声。这些噪声在光谱图上会形成波动和毛刺，对光谱的峰值位置和波峰的高度会产生一定的影响。根据CCD采集到的光谱特性以及各种去噪方法的特点，采用小波阈值和中值滤波相结合的方法对含噪光谱数据进行处理。

经过小波分解，信号的能量主要分布在低分辨的尺度函数上，而噪声的能量分布不变，仍然均匀分布在所有小波系数上。基于此特性，小波阈值去噪方法就是保留或收缩大的小波系数(低频系数)，剔除小的小波系数(高频系数)。在使用小波阈值法对光谱数据进行去噪处理时，选取不同的阈值函数反映了对小波稀疏模的不同处理策略，通常阈值函数采用Donoho等在文章“Ideal Spatial Adaptation by Wavelet Shrinkage”(Biometrika,Vol.81,No.3,pp.425-455,1994)提出的硬阈值函数(Hard Threshold Function)和软阈值函数(Soft Threshold Function)两种。其中硬阈值函数对于信号的细节特征可以很好地保留，但在阈值±T处不连续，且不可导，同时其剔除了小于阈值的小波系数，不加处理的完全保留了大于阈值的小波系数，因此重构信号容易产生较大的均方差，出现振荡；软阈值函数在阈值±T处是连续的，可以避免去噪后的信号产生振荡现象，但估计的小波系数与实际信号的小波系数之间存在恒定的偏差，使得重构的信号过度平滑，使得边缘相对模糊，同样影响去噪效果。

发明内容

本发明针对上述现有的阈值函数的缺陷，提出了一种指数软硬择中阈值函数，对光谱信号进行去噪处理，以减少了降噪后光谱信号的噪声残留，提高了光谱信号的信噪比。

为了实现上述目标，本发明包括以下步骤：

步骤1：对输入的含噪光谱信号进行多层小波变换。采用db4小波对含噪光谱信号进行5层分解，分别得到小波的低频系数和高频系数；

步骤2：对步骤1得到的各层高频小波系数进行收缩处理。在收缩处理中阈值函数采用了本发明所提出的指数软硬择中阈值函数，即