一种去除背景噪声影响的拉曼光谱预处理方法，步骤如下：1）获得信息拉曼光谱及背景光谱；2)分别对信息拉曼光谱和背景拉曼光谱作一阶微分，得到微分差值光谱；3)对微分差值光谱作分段处理；4）对经过步骤3）处理后的微分差值光谱作有效信息保留的处理；5）对经过步骤4）处理后微分差值光谱的有效部分进行分别分段积分；6）对步骤5）获得的积分光谱的突变点作修正，即获得预处理后的拉曼信息光谱。本方法用于拉曼光谱预处理中，能够有效地从复杂背景噪声中提取出有效拉曼信息，方法简单，容易实现。
1.一种去除背景噪声影响的拉曼光谱预处理方法，其特征在于，具体步骤如下：
1)获取信息拉曼光谱Ls及背景拉曼光谱Lb：加入被测物质之前，用拉曼光谱仪对背景进行测量，获得背景拉曼光谱Lb；加入被测物质后，再用拉曼光谱仪测量获得信息拉曼光谱Ls；
2)分别对信息拉曼光谱Ls及背景拉曼光谱Lb作一阶微分处理，并将含有信息拉曼光谱Ls的一阶微分值减去背景拉曼光谱Lb的一阶微分值，得到微分差值光谱CO；
3)将步骤2)得到的微分差值光谱CO作分段处理，将微分差值的绝对值小于阈值σ的分为一段，并用Flag函数标记该段为0，将微分差值大于σ的分为一段，并用Flag函数标记该段为1，将微分差值小于-σ的的分为一段，并用Flag函数标记该段为-1；
4)经过步骤3)处理后的微分差值光谱作有效信息保留处理，保留的有效信息微分差值光谱为C1；
5)对微分差值光谱C1进行分段积分，获得有效拉曼光谱L1；
6)若有效拉曼光谱L1中包含由正值或负值直接变为0的突变点，则对突变点所在区间作线性修正，直至处理后的光谱中没有突变点，即获得预处理后的有效信息拉曼光谱L2。
2.如权利要求1所述的一种去除背景噪声影响的拉曼光谱预处理方法，其特征在于：
步骤3)中的阈值σ＝k*V，其中V为含有有效信息的光谱中选取的比较平衡的一段光谱的散射强度的标准差，标准差计算公式为式中，u为所选取光谱散射强度的平均值，计算分式为 为1，2或3的整
3.如权利要求1所述的一种去除背景噪声影响的拉曼光谱预处理方法，其特征在于，步骤4)中对微分差值光谱作有效信息保留处理的方法为：找到标记函数Flag从1到0的突变点M及标记函数Flag从0到-1的突变点N，若M与N之间的标记函数Flag值全为0，则从M点往前查找微分差值由大于零变为不大于零的突变点M1，直至M1点的微分差值不大于0，从N点往后查找微分差值由小于零到不小于零的突变点N1，直至N1点的微分差值不小于0，则有效信息区间[M1,N1]，保留所有有效信息区间内全部微分差值，其他值设为0，得到处理后的微分差值光谱C1。
4.如权利要求3所述的一种去除背景噪声影响的拉曼光谱预处理方法，其特征在于，步骤5)对微分差值光谱C1进行分段积分的方法为：对微分差值光谱中C1中不为零的有效信息段进行分别分段积分，为零的部分积分值置为0，获得有效拉曼光谱L1。
5.如权利要求1所述的一种去除背景噪声影响的拉曼光谱预处理方法，其特征在于，步骤6)中对突变点进行修改的方法为：设突变点处的散射位移及散射强度分别为x1和y1，往前找到第一个散射强度0的点，记其散射位移为x0，对散射位移在x0到x1之间的点的散射强度进行线性校正，校正公式为 式中xi和yi分别为所修正的点对应的散射位移及散射强度，y'i为修正后的点所对应的散射强度。
一种去除背景噪声影响的拉曼光谱预处理方法
技术领域
[0001] 本发明涉及一种光学领域，特别是一种拉曼光谱去噪处理方法。
背景技术
[0002] 在进行拉曼光谱分析时，被测量物质的拉曼光谱往往会受到背景噪声的影响。这些背景噪声包括机器噪声、荧光干扰、其他物质的拉曼光谱等。背景噪声干扰会影响光谱分析的精度，背景噪声干扰严重会使光谱分析难以有效进行，因此，在进行光谱分析前往往需要找到合适的消除背景噪声影响的光谱预处理方法。
[0003] 常用的去除背景、噪声影响的光谱预处理方法包括 数字平滑滤波：常用的数字平滑滤波方法有窗口移动平均法，窗口移动二项式平滑，通过数字平滑滤波可以提高信噪比，其实质是去除随机误差信号引去的分析信号偏离，对高频噪声有效。(2) 多项式拟合基线校正：利用多项式拟合的基线校正法可以消除荧光背景，这种方法对于荧光背景这样整体平缓的背景消除有效，本质上是对低频噪声的消除。(3)小波变换：小波变换可以通过多分辨率的特性消除噪声，但前提是噪声、背景要与信号有明显的频率差异，才能找到合适的分解方法。(4)FT与HT滤波：低通可以过滤随机信号等高频噪声；高通可以去掉基线漂移等低频信号。以上这些方法在某些情况下能够有效去除机器噪声、荧光等的干扰，但在复杂体系中，这些方法对于其他非探测物质的拉曼光谱干扰则无能为力，面对高强度的荧光背景干扰时也达不到理想的效果。
[0004] 另一方面，由于环境的不稳定，如温度的不稳定，周围环境的变化，表面增强基底的不稳定等，这些因素往往会使为拉曼光谱的散射强度变化很大，直接通过信息拉曼光谱与背景光谱相减的处理方式不能有效去除背景干扰。
发明内容
[0005] 本发明的目的就是提供一种去除背景噪声影响的拉曼光谱预处理方法，它可以有效去除包括高强度荧光的拉曼光谱干扰背景噪声的影响，适用于背景噪声不稳定的情况。
[0006] 本发明的目的是通过这样的技术方案实现的，具体步骤如下：
[0007] 1）获取信息拉曼光谱Ls及背景拉曼光谱Lb，加入被测物质之前，用拉曼光谱仪对背景进行测量，获得背景光谱Lb;加入被测物质后，再用拉曼光谱仪测量获得拉曼信息光谱Ls。
[0008] 2）分别对信息拉曼光谱Ls及背景拉曼光谱Lb作一阶微分处理，并将含有信息光谱Ls的一阶微分值减去背景光谱Lb的一阶微分值，得到微分差值光谱CO；
[0009] 3）将步骤2）得到的微分差值光谱CO作分段处理，将微分差值的绝对值小于阈值σ的分为一段，并用Flag函数标记该段为0，将微分差值大于σ的分为一段，并用Flag函数标记该段为1，将微分差值小于-σ的的分为一段，并用Flag函数标记该段为-1；
[0010] 4）经过步骤3）处理后的微分差值光谱作有效信息保留处理，保留的有效信息微分差值光谱为C1；
[0011] 5）对微分差值光谱C1进行分段积分，获得有效拉曼光谱L1；
[0012] 6）若有效拉曼光谱L1中包含由正值或负值直接变为0的突变点，则对突变点所在区间作线性修正，直至处理后的光谱中没有突变点，即获得预处理后的有效信息拉曼光谱L2。
[0013] 进一步，步骤3）中的阈值σ=k*V，其中V为含有有效信息的光谱中选取的比较平衡的一段光谱的散射强度的标准差，标准差计算公式为
[0014]
[0015] 式中，u为所选取光谱散射强度的平均值，计算分式为 ，k为1，2或3的整数。
[0016] 进一步，步骤4）中对微分差值光谱作有效信息保留处理的方法为：找到标记函数Flag从1到0的突变点1到0的突变点M及标记函数Flag从0到-1的突变点N，若M与N之间的标记函数Flag值全为0，则从M点往前查找微分差值由大于零变为不大于零的突变点M1，直至M1点的微分差值不大于0，从N点往后查找微分差值由小于零到不小于零的突变点N1，直至N1点的微分微分差值不小于0，则有效信息区间[M1,N1]，保留所有有效信息区间内全部微分差值，其他值设为0，得到处理后的微分差值光谱C1。
[0017] 进一步，步骤5）对微分差值光谱C1进行分段积分的方法为：对微分差值光谱中C1中不为零的有效信息段进行分别分段积分，为零的部分积分值置为0，获得有效拉曼光谱L1。
[0018] 进一步，步骤6）中对突变点进行修改的方法为：设突变点处的散射位移及散射强度分别为 和 ，往前找到第一个散射强度0的点，记其散射位移为 ，对散射位移在到 之间的点的散射强度进行线性校正，校正公式为 ，式中 和 分别为所修正的点对应的散射位移及散射强度， 为修正后的点所对应的散射强度。
[0019] 由于采用了上述技术方案，本发明具有如下的优点：
[0020] 本发明采用形状相似性判断的方法，将背景与信息区别分离开来，根据拉曼峰的分布规律对有效的信息进行了保留，因此能够有效地去除背景噪声的影响。本发明是采用基于形状而不是拉曼散射强度的判断方式去除背景噪声、保留有效信息，克服了在背景拉曼光谱散射强度不稳定时，直接采用强度相减的方法容易带来额外噪声的缺点。
[0021] 本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书和权利要求书来实现和获得。
附图说明
[0022] 本发明的附图说明如下。
[0023] 图1为本发明方法的整体流程图；
[0024] 图2 为原始的信息拉曼光谱及背景拉曼光谱的示例图；
[0025] 图3为信息光谱及背景光谱分别一阶微分后相减的处理效果图；
[0026] 图4为分段处理后的效果局部放大图；
[0027] 图5为有效信息保留处理后的效果局部放大图；
[0028] 图6为分段积分处理后的效果局部放大图；
[0029] 图7为突变点修正处理后的效果的局部放大图；
[0030] 图8为经过本发明预处理后得到的信息拉曼光谱实例图。
具体实施方式
[0031] 下面结合附图和实施例对本发明作进一步说明。
[0032] 一种去除背景噪声影响的拉曼光谱预处理方法，如图1所示，包括有对信息光谱及背景光谱分别作一阶微分，然后将信息光谱的微分值减去背景光谱的微分值得到微分差值光谱C0；Flag函数对C0作分段标记处理；有效信息保留处理得到微分差值光谱C1；对微分差值光谱C1中有效信息区间进行分别分段积分，获得有效拉曼光谱L1；进行突变点修正得到本发明预处理后得到最终信息拉曼光谱。
[0033] 未经任何处理的信息拉曼光谱及背景拉曼光谱如图2所示，其中实线表示的是信息拉曼光谱Ls，虚线表示的是背景拉曼光谱Lb。
[0034] 分别对信息拉曼光谱Ls及背景拉曼光谱Lb作一阶微分，并将含有效信息的光谱的一阶微分值减去背景光谱的一阶微分值，得到微分差值光谱C0，如图3所示。
[0035] 接着，C0作分段处理，将微分差值的绝对值小于阈值σ的分为一段，并用Flag函数标记该段为0，将微分差值大于σ的分为一段，并用Flag函数标记该段为1，将微分差值小于-σ的分为一段，并用Flag函数标记该段为-1。处理后，Flag函数分段标记的效果如图4 所示，其中实线部分为微分差值光谱C0的局部放大图，虚线部分表示的是Flag函数分段标记的局部放大效果。
[0036] 其中，阈值σ=k*V,其中V为含有有效信息的光谱中选取的比较平稳的一段光谱的纵坐标散射强度的标准差，标准差的计算公式为 ，u为所选取光谱纵坐标散射强度的平均值,计算公式为 ，k为根据实际应用可灵活定义的阈值系数，确定阈值系数k=1。
[0037] 紧跟着，作有效信息保留处理。其具体步骤为找到标记函数Flag从1到0的突变点M及标记函数Flag从0到-1的突变点N，若M与N之间的标记函数Flag值全为0，则从M点往前查找点M1，直至M1点的微分差值不大于0，从N点往后查找N1点，直至N1点的微分差值不小于0，则有效信息区间为[M1，N1]，保留所有有效信息区间内的全部微分差值，其他值设为0,得到处理后的微分差值光谱C1，如图5所示。
[0038] 然后，进行积分处理。具体积分方法为：对微分差值光谱中C1中不为零的有效信息段进行分别分段积分，为零的部分积分值置为0，获得有效拉曼光谱L1如图6所示，图6中的圆点表示的是有效拉曼光谱L1中由正值或者负值直接变为0的突变点。
[0039] 最后，对L1由正值或者负值直接变为0的突变点进行修正，得到修正后的信息拉曼光谱L2，修正方法的具体方法为：设突变点处的散射位移及散射强度分别为 , ,往前找到第一个散射强度0的点，记其散射位移为 ,对散射位移在 到 之间的点的散射强度进行线性校正，校正公式为 ， ， 分别为所修正的点对应的散射位移及散射强度， 为修正后的点所对应的散射强度。图8为修正后的信息拉曼光谱L2的效果图，图7为信息拉曼光谱L2的局部放大效果图
[0040] 最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。
