一种用于微量气体拉曼光谱检测基线校正的包络线迭代方法
【专利摘要】本发明提出一种微量气体拉曼光谱检测基线校正的包络线迭代方法，该方法用diff函数对去噪后的气体拉曼光谱W求数值微分,用find函数和极小值点左边曲线的斜率小于0且右边曲线的斜率大于0的条件求得极小值点，用三次样条曲线插值函数spline和极小值点求得一次包络线L1，对前后谱峰密集度不一致的气体拉曼光谱进行分段处理，并引入包络线迭代思想，使用密集段迭代次数多余稀疏段迭代次数的方式，得到包络线L2(即为气体拉曼光谱基线)，W减去基线L2，完成对W的基线校正。该方法处理结果能够较为准确地校正微量气体拉曼光谱的基线，减少校正后的光谱荧光背景残留，提高气体拉曼光谱定量分析准确度。
【专利说明】
一种用于微量气体拉曼光谱检测基线校正的包络线迭代方法
技术领域
[0001]本发明涉及气体拉曼光谱检测领域，尤其涉及微量气体拉曼光谱的基线校正技术。
【背景技术】
[0002]拉曼光谱(Ramanspectroscopy)是一种基于拉曼散射效应的散射光谱，由印度物理学家Raman所发现，具有非侵入、无损伤和无辐射等优点。随着激光光源的引入、信号检测技术的提高和计算机技术的普及，拉曼光谱技术正日益地应用于化学化工、生物医学检测和材料生产等领域。
[0003]然而受仪器系统、检测环境以及样品等诸多因素的影响，获得的拉曼光谱检测信号总存在不同程度的噪声和背景干扰。由荧光背景等因素使拉曼光谱基线升高的现象称为基线漂移。当基线漂移较大时，将严重影响拉曼特征峰的提取以及被测物质的识别，降低定性定量分析的准确性。因此，基线校正直接关系到光谱分析质量的好坏，对这方面的研究一直备受广大学者关注。目前减小基线漂移的方法主要是两个方面，拉曼光谱分析仪器改良和数据处理方法。拉曼光谱分析仪器改良成本较高且无法完全消除荧光物质的影响，而数据处理方法则成本较低且效果较好。现有的多项式拟合及分段多项式拟合算法是最常用的基线校正方法，但都存在荧光背景残留较多的问题，从而导致基线校正准确度较差，不利于提高气体拉曼光谱检测准确度。而小波变换基线校正法是将拉曼信号分解为一系列正弦波的叠加，能够将信号的低频和高频分开，具有快速、准确等优点，但分解尺度依靠经验，对于不同的实际问题难以达到最佳处理效果，且基线校正易失真或不充分，其计算量和计算复杂度也相对较高。
【发明内容】
[0004]本发明所要解决的问题是针对现有拉曼光谱基线校正方法存在的不足，提出一种微量气体拉曼光谱检测基线校正的包络线迭代方法，该方法基于三次样条曲线插值算法和包络线迭代思想，能够更加准确地表示气体拉曼光谱的基线，减少气体拉曼光谱荧光背景的残留，提高气体浓度的定量分析检测准确度。
[0005]为达到上述目的，本发明采用的技术方案如下:
[0006]—种微量气体拉曼光谱检测基线校正的包络线迭代方法，步骤如下:
[0007]步骤一:对获得的气体拉曼光谱进行去噪处理，定义去噪后的气体拉曼光谱为W;
[0008]步骤二:用dif f函数对W求数值微分，设W共有N+1个点，则数值微分后得到N个值，将其中第I至第N-1值存入矩阵dl中，第2至第N值存入矩阵d2中；
[0009]步骤三:用find函数和极小值点左边曲线的斜率小于O且右边曲线的斜率大于O的条件求得极小值点 indmin，即 indmin = find (dl.*d2〈0&dl〈0)+l ；
[0010]步骤四:用三次样条曲线插值函数spline和步骤三求得的极小值点indmin求得一次包络线L1;[0011 ] 步骤五:判断W谱峰密度，若W前后谱峰密集度一致，重复步骤二至步骤四，对1^进行m次迭代求得包络线L2;若W前后谱峰密集度不同，分离1^谱峰的密集段与稀疏段，分别用步骤二至步骤四对L1稀疏段迭代m次，对密集段迭代m+1次，得到包络线L2,即为气体拉曼光谱检测基线。其中迭代次数m可根据W情况进行增减，，m通常为2?4。
[0012]步骤六:W减去基线1^2，即用迭代的包络线作为基线，原始光谱减去包络线法形成的基线，即完成对W的基线校正。
[0013]在本发明的步骤一中，对测量所得的气体拉曼光谱进行去噪处理，使谱线平滑，可以保证步骤三中极小值点求取的准确度，使步骤五求取的包络线更接近气体拉曼光谱的实际基线，在一定程度上减少了荧光背景的残留。
[0014]在发明的步骤二中，由于计算机只能处理离散数据，故引入diff函数求数值微分，以得到单位横坐标Ax下离散数据的Ay，并以此来代替气体拉曼光谱曲线的斜率。
[0015]在本发明的步骤三中，利用极小值点左边曲线斜率小于O且右边曲线斜率大于O的性质，结合find函数可以准确求得气体拉曼光谱中所有的极小值点横坐标。此外，考虑到步骤二中数值微分后数据总个数的缺失，对极小值点的横坐标进行了加I补偿，保证了基线校正的准确性。
[0016]在本发明的步骤四、五中，对前后谱峰密集度不一致的气体拉曼光谱进行分段处理，并引入包络线迭代思想，使用密集段迭代次数多余稀疏段迭代次数的方式，全面考虑气体拉曼光谱的特点，体现了微量气体拉曼光谱检测基线校正的包络线迭代方法的通用性。
[0017]本发明方法处理结果能够较为准确地校正微量气体拉曼光谱的基线，减少校正后的光谱荧光背景残留，提高气体拉曼光谱检测准确度，对于气体拉曼光谱的定量分析及应用具有非常重要的实际意义。该方法的优点具体如下:
[0018]①使用微量气体拉曼光谱检测基线校正的包络线迭代方法后，可有效减少光谱荧光背景残留，提高气体拉曼光谱检测准确度
[0019]②全面考虑气体拉曼光谱的特点，对气体拉曼光谱进行分段处理，密集段迭代次数多，稀疏段迭代次数少，体现了微量气体拉曼光谱检测基线校正的包络线迭代方法的通用性。
【附图说明】
[0020]图1为微量气体拉曼光谱检测基线校正的包络线迭代方法流程图；
[0021]图2是出拉曼光谱去噪前后的对比图；
[0022]图3是去噪后的出拉曼光谱及迭代包络线L2;
[0023]图4H2拉曼光谱基线校正前后对比图。
【具体实施方式】
[0024]为了更好地理解本发明，下面结合实施例进一步阐述本发明的内容。
[0025]如图1所示为微量气体拉曼光谱检测基线校正的包络线迭代方法流程图，以微量气体H2为例，按该方法对H2拉曼光谱进行基线校正处理，用MATLAB软件编程实现，【具体实施方式】如下:
[0026]①.利用拉曼光谱分析系统对微量H2进行拉曼光谱测量，运用小波模极大值法对获得的出拉曼光谱进行去噪处理，定义去噪后的出拉曼光谱数据为W，其结果如图2所示。
[0027]②.基于diff函数对W求数值微分，W共有2000个点，数值微分后得到1999个值，将第I至1998值存入矩阵dl中，第2至1999值存入矩阵d2中。
[0028]③.基于find函数求极小值点indmin并对其进行加I补偿，即indmin = find(dl.*d2〈0&dl〈0)+l求得的极小值点。
[0029]④.用三次样条曲线插值函数spline和极小值点indmin求得一次包络线Li;
[0030]⑤.在出拉曼光谱数据第1000个点处分段，对左右两边分别重复步骤③和④，进行两次迭代，得到迭代的包络线L2如图3所示。
[0031]⑥.用迭代的包络线L2代替基线，原始光谱W与基线应点的纵坐标相减完成基线校正，基线校正前后对比分析结果如图4所示。
【主权项】
1.一种微量气体拉曼光谱检测基线校正的包络线迭代方法，其特征在于，步骤如下: 步骤一:对获得的气体拉曼光谱进行去噪处理，定义去噪后的气体拉曼光谱为W; 步骤二:用diff函数对W求数值微分，设W共有N+1个点，则数值微分后得到N个值，将其中第I至第N-1值存入矩阵dl中，第2至第N值存入矩阵d2中； 步骤三:用find函数和极小值点左边曲线的斜率小于O且右边曲线的斜率大于O的条件求得极小值点indmin，即 indmin = f ind(dl.*d2〈0&dl〈0)+l ； 步骤四:用三次样条曲线插值函数spline和步骤三求得的极小值点indmin求得一次包络线L1; 步骤五:判断W谱峰密度，若W前后谱峰密集度一致，重复步骤二至步骤四，对1^进行m次迭代求得包络线L2;若W前后谱峰密集度不同，分离1^谱峰的密集段与稀疏段，分别用步骤二至步骤四对L1稀疏段迭代m次，对密集段迭代m+1次，得到包络线L2,即为气体拉曼光谱检测基线； 步骤八:W减去基线L2，即完成对W的基线校正。2.根据权利要求1所述的微量气体拉曼光谱检测基线校正的包络线迭代方法，其特征在于，步骤五的迭代次数m可根据W情况进行增减，，m通常为2?4。
【文档编号】G01N21/65GK106053430SQ201610430372
【公开日】2016年10月26日
【申请日】2016年6月16日
【发明人】陈伟根, 王建新, 杨曼琳, 贺鹏, 王品, 王品一, 万福
【申请人】重庆大学