块进行光谱预处理；
[0041 ]系统根据设定取一测量区间；
[0042]收集测量区间内的光谱数据，进行区间相似度计算；
[0043]将区间移动到下一区间单位重复计算；
[0044]获取全部的区间相似性。
[0045]结合实施例对本发明的方法分析过程进行详述:
[0046]如图2-1，2-2，2-3光谱所示，本实施例将三种不同名称的塑料样品分别放在Xplora共焦拉曼光谱仪的扫描平台上，在LabSpec软件上选择激发波长为785nm，设置仪器的积分时间为10s，通光孔径为300mm，光谱仪狭缝大小100mm。对每种样品在不同部位分别测量三次，得到一种样品的三组拉曼光谱数据，将其绘成光谱曲线。图中编号“样品1-1、1_
2、1_3”分别表示样品I的第一次、第二次、第三次测量数据，以此类推。从图2-1?2-3可见，同一样品三次测量的拉曼光谱曲线差别很小，而不同样品之间的拉曼光谱有较大差异。
[0047]由于拉曼光谱仪采集到的样品原始光谱包含荧光背景、检测器(CCD)噪声、激光器功率波动等干扰信息，会影响后续的比对和信号处理。因此，需要对测量得到的原始光谱数据预处理，以利于后续有效信息的提取。
[0048]本发明中对数据的预处理分如下两个步骤:I)对光谱进行平滑滤波;2)荧光背景扣除和归一化。图3是样品I第一次测量的原始拉曼光谱和平滑去背景后的光谱比较图，从图中可以看出处理后的光谱基线基本处在O附近.
[0049]两种物质的拉曼光谱相似性可以用相似性测度(以下简称“相似度”)这一参数定量地表示。每种物质的拉曼光谱是组成该物质分子结构的反映，具有其独特的结构和模式特性。令光谱数据点数为模式空间的维数，这样一张图谱在模式空间中即可表达为一个模式向量，N张图谱间相似性的分析就转化为计算模式空间中N个模式向量的相似度。采用相似性函数计算相似度的方法较为普遍，本发明采用夹角余弦法、广义Jaccard系数法这两种较有代表性的方法进行相似度的计算，并对计算结果进行讨论。
[0050]从图4-1和图4-2均能看出，同一样品三次测量的相似度均在0.92以上，不同样品之间的相似度均在0.32以下，可以从相似度明显地区分相同和不同样品。用夹角余弦法对同一样品在三个不同部位分别进行三次测量时得到的拉曼光谱数据进行计算时，得到的相似度结果比用广义Jaccard系数法得到的结果更大;而广义Jaccard系数法在计算不同样品的相似度时将给出比夹角余弦法更低的结果。因此，为了更准确地判别两种样品的相似性，可以将夹角余弦和广义Jaccard系数法结合使用。
[0051]从上述分析可知，由于第一类样品之间拉曼光谱峰值出现的位置和幅度有较大的差异，通过夹角余弦和Jaccard系数法计算相似度可以比较容易地区分出两种样品是否为同一种材料。但如果我们对具有相同名称、仅在有无添加阻燃剂方面有区别的第二类塑料样品做拉曼光谱(如图5所示)分析时，就可以看出有/无阻燃剂的ABS样品拉曼光谱基本相同，添加阻燃剂的样品在250cm-1和825cm-1附近出现新的尖峰，而未添加阻燃剂的样品在这两个波数位置没有波峰。对图5光谱采用广义Jaccard系数法进行全区段(全局)相似度分析，这两条光谱的相似度可达0.85。如果样品中添加物含量进一步减小，两种相同基质的样品拉曼光谱将仅在某个局部出现很小的差异。因此，对两种成份相差很小的样品采用全局的拉曼光谱相似度分析，很难从相似度计算结果将样品进行区分。
[0052]为解决在采用全局光谱数据进行相似度分析时遇到的:I)两种样品拉曼光谱差别出现位置事先未知；2)光谱差别很小，导致从相似度计算结果难以区分两种不同样品的难题，本发明根据两种差异很小的样品在拉曼光谱结构上的特点，创新地提出在拉曼光谱的全部区间，以一定的光谱间隔分别计算不同波数位置的局部光谱区间相似度，得到局部光谱相似度随拉曼光谱波数的变化曲线。因此，如果两种样品的拉曼光谱仅在某个波数附近的局部区间存在差异，则在此波数计算得到的局部相似度将出现较大的波动。图6是当光谱测量区间宽度设为5cm-1，对图5光谱数据进行移动局部相似度分析得到的结果。
[0053]从图6可以看出，当窗口移动到添加阻燃剂产生新的波峰处时，两种样品的相似度明显的下降，而在其他位置相似度基本为I。相比之下，对同一样品不同次数的测试数据进行如上运算时，在不同波数位置的相似度基本都为I，结果如图7所示。
[0054]以上所述，仅为本发明较佳的【具体实施方式】，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
【主权项】
1.一种通过拉曼光谱检测物品一致性的方法，其特征在于包括以下步骤: 1)使用光谱仪获取不同样品的光谱数据； 2)对获取的光谱数据进行预处理； 3)设定一段光谱间隔区间，分别计算不同波数位置的局部光谱间隔区间相似度； 4)得到局部光谱相似度随拉曼光谱波数的变化曲线。 5)根据曲线得出物品的一致性。2.如权利要求1所述的通过拉曼光谱检测物品一致性的方法，其特征在于:所述样品为塑料制品。3.如权利要求1所述的通过拉曼光谱检测物品一致性的方法，其特征在于:在步骤I)时，对同一样品测量三次。4.如权利要求1所述的通过拉曼光谱检测物品一致性的方法，其特征在于:在步骤2)时，所述预处理包括对光谱进行平滑滤波、荧光背景扣除、归一化。5.如权利要求1所述的通过拉曼光谱检测物品一致性的方法，其特征在于:在步骤3)设定的光谱间隔为5cm—、6.—种通过拉曼光谱检测物品一致性的系统，包括拉曼光谱测量系统，所述测量系统包括激发光源、光纤或光束传输光路、含波长选择滤光镜的检测探头、光谱检测单元、信号处理与显示，其特征在于:还包含一数据处理模块，所述数据处理模块包括数据收集模块用于收集所述光谱测量系统反馈的光谱数据，数据计算模块用于进行相似度计算，数据结果显示模块，所述数据收集模块具有一可移动的测量区间。7.如权利要求6所述的通过拉曼光谱检测物品一致性的系统，其特征在于:所述测量区间收集区间内光谱数据，传输至数据计算模块计算，之后所述测量区间收集下一单位区间内的光谱数据。
【专利摘要】一种通过拉曼光谱检测物品一致性的方法，包括以下步骤：1)使用光谱仪获取不同样品的光谱数据；2)对获取的光谱数据进行预处理；3)设定一段光谱间隔区间，分别计算不同波数位置的局部光谱间隔区间相似度；4)得到局部光谱相似度随拉曼光谱波数的变化曲线。5)根据曲线得出物品的一致性。
【IPC分类】G01N21/65
【公开号】CN105606584
【申请号】CN201510932814
【发明人】吴兆喜, 董小鹏, 赖博文, 陆豆
【申请人】厦门出入境检验检疫局检验检疫技术中心, 厦门大学
【公开日】2016年5月25日
【申请日】2015年12月15日