一种基于特征峰的光谱重构方法与流程

1.本发明涉及光谱重构技术领域，尤其涉及一种基于特征峰的光谱重构方法。


背景技术：

2.近年来，近红外光谱分析技术发展十分迅速，已在化工，制药，军工，食品等多个领域都获得了应用。近红外光谱技术属于分子光谱技术，可以在分子水平上表明物质成分和性质信息，不论对经济还是社会影响来说，都取得了非常高的效益，极具发展潜力。
3.然而，目前大多数物质成分和性质信息检测主要使用大型实验室近红外光谱仪器进行，这些方法虽然定量准确灵敏度高，但所需设备体积庞大，设备费用昂贵，样品制备时间长且样品制作方法严格，检测设备和样品制备需要专业人员操作，检测环境固定，且分析时间长，不适用于现场检测，不便于推广使用。
4.伴随着便携式近红外光谱技术的发展，市场主流的大型近红外光谱仪设备都朝着体积小巧，价格低廉的便携式方向发展。但是便携式近红外光谱仪易于受光源、检测器、使用方法、环境条件等影响，其光谱分析能力偏弱，尤其是在样品数量十分有限的情况下，便携式近红外光谱分析技术容易受到异常样本的影响进而降低其预测分析能力。在便携式近红外光谱分析技术的实际应用中，往往样品数量十分有限，且便携式近红外光谱设备采集到的数据也较为冗余，包含了过多与样品相关性较小的数据信息，对建模分析工作带来了较大的工作量与难度，由此如何获取一种既能缩减光谱自身数据量，又可以最大程度的保证样品特性，且增大样品光谱数据库的方法成为了必须解决的问题。


技术实现要素：

5.本发明的目的是克服现有技术中存在的不足，提供一种基于特征峰的光谱重构方法。该方法首先采集待测样品的光谱数据，然后对待测样品光谱数据进行二阶求导，得到待测样品特征峰，根据便携式近红外光谱波长点权重特性选取一定数目的特征峰光强点替代原始光谱数据，然后根据样品组分标定值进行插值运算，计算虚拟插值点，再根据样品替换后的光谱数据及虚拟插值点数目计算虚拟光谱数据，最后将虚拟光谱数据与虚拟组分标定值进行一一对应完成光谱数据库的扩充。该方法通过特征峰信息选取权重系数较高的部分光强点代替原始光谱数据，进而在保留光谱数据特征信息前提下极大程度减小光谱数据量，提高光谱重构效率，有效扩充样品光谱数据库，解决便携式近红外光谱分析技术需要大量实体样品才能有较好分析预测效果的问题。
6.为了实现上述目的，本发明采用的技术方案是：一种基于特征峰的光谱重构方法，包括如下步骤：
7.a.采集待测样品的光谱数据；
8.b.对待测样品光谱数据进行二阶求导，得到待测样品特征峰；
9.c.根据便携式近红外光谱波长点权重特性选取一定数目的特征峰光强点替代原始光谱数据；
10.d.根据样品组分标定值进行插值运算，计算虚拟插值点，再根据样品替换后的光谱数据及虚拟插值点数目计算虚拟光谱数据；
11.e.将虚拟光谱数据与虚拟组分标定值进行一一对应完成光谱数据库的扩充。
12.进一步的技术方案为，步骤a中所述采集待测样品的光谱数据，采用波长均分式便携近红外光谱仪对待测样品光谱数据进行采集，其波长范围为,分辨率为，则每个待测样品的实际光谱数据包含50个光强值点。
13.进一步的技术方案为，步骤b中二阶求导光谱的半高宽只有原始谱图半峰宽的，能够简单的分辨出强峰两侧的小肩峰，在正确测定峰位及确定肩峰位置极为有效，通过二阶求导可以清晰的分辨出待测样品谱图的峰值，即特征峰波长点位置。
14.进一步的技术方案为，步骤c具体为根据待测样品光谱数据的二阶求导光谱图形，选取特征峰周围权重系数较高的波长点替代原始光谱数据。
15.进一步的技术方案为，步骤d具体为将待测样品的组分标定值进行顺序排布，每两个相邻标定值之间进行插值运算，计算得到虚拟插值点，再根据具体的虚拟插值点数目以及相邻标定值样品替换后的光谱数据计算得到虚拟光谱数据。
16.进一步的技术方案为，步骤e具体为将虚拟组分标定值按照顺序放入实际待测样品组分标定值之间，再将得到的虚拟光谱数据按照顺序放入替换后的实际待测样品光谱数据中，虚拟组分标定值与虚拟光谱数据进行一一对应，进而有效扩充样品光谱数据库。
17.与现有技术相比，本发明具有如下有益效果：提供一种基于特征峰的光谱重构方法，该方法通过特征峰信息选取权重系数较高的部分光强点代替原始光谱数据，进而在保留光谱数据特征信息前提下极大程度减小光谱数据量，提高光谱重构效率，有效扩充样品光谱数据库，解决便携式近红外光谱分析技术需要大量实体样品才能有较好分析预测效果的问题。
附图说明
18.图1是本发明一种基于特征峰的光谱重构方法的示意图；
19.图2是本发明待测样品二阶求导后的光谱数据图。
具体实施方式
20.下面结合附图和具体的实施方式对本发明作进一步详细的描述。
21.实施例1
22.如图1所示，图1中101是采集待测样品的光谱数据。采用波长均分式便携近红外光谱仪对待测样品光谱数据进行采集，其波长范围为1758nm～2150nm,分辨率为8nm，则每个待测样品的实际光谱数据包含50个光强值点。
23.在本实施例中，不同待测样品其光谱特征峰并不相同，采用波长均分式便携式近红外光谱仪对待测样品进行光谱数据采集可以最大均匀程度的采集到样品光谱数据，根据本专利的实际待测样品，选用波长范围为1758nm～2150nm,分辨率为8nm的便携式近红外光谱仪进行数据采集，可以计算出每一条光谱数据包含的光强点为n＝1+(2150-1758)/8＝50个，第1～50个波长点对应的波长范围为(1758nm，1766nm，
……
2150nm)，每个待测样品采集获得的光谱数据实际表示为50个波长点上光强值的集合矩阵。为了尽可能的减少测量误
差，每个待测样品进行四次光谱数据采集取平均值作为最终光谱数据值。
24.图1中102是对待测样品光谱数据进行二阶求导，得到待测样品特征峰。二阶求导光谱的半高宽只有原始谱图半峰宽的1/3左右，能够简单的分辨出强峰两侧的小肩峰，在正确测定峰位及确定肩峰位置极为有效，通过二阶求导可以清晰的分辨出待测样品谱图的峰值，即特征峰波长点位置。
25.在本实施例中，如图2所示，对采集得到的待测样品光谱数据进行二阶求导，在波长点为1790nm的位置出现最为明显的谱峰，该峰值点即为该待测样品的特征峰。
26.图1中103是根据便携式近红外光谱波长点权重特性选取一定数目的特征峰光强点替代原始光谱数据。根据待测样品光谱数据的二阶求导光谱图形，选取特征峰周围权重系数较高的波长点替代原始光谱数据。
27.在本实施例中，如图2所示，待测样品光谱数据的二阶求导光谱图形中，其特征峰在波长点1790nm的位置，其中特征峰周围光强点均有较高的权重系数，即在权重系数较高的波段范围为1758nm～1838nm。由以上可知，替代原始50个光强数值点光谱数据的重组光谱数据包含的光强值点数目为m＝1+(1838-1758)/8＝11个，第1～11个波长点对应的波长范围为(1758nm，1766nm，
……
1838nm)，重构之后每个待测样品光谱数据实际表示为11个波长点上光强值的集合矩阵。相较于待测样品的原始光谱，最大程度的保留特征信息之外，每条光谱数据的光强值点数目得到了极大的缩减，有效的提高光谱分析效率。
28.图1中104是根据样品组分标定值进行插值运算，计算虚拟插值点，再根据样品替换后的光谱数据及虚拟插值点数目计算虚拟光谱数据。将待测样品的组分标定值进行顺序排布，每两个相邻标定值之间进行插值运算，计算得到虚拟插值点。再根据具体的虚拟插值点数目以及相邻标定值样品替换后的光谱数据计算得到虚拟光谱数据。
29.在本实施例中，如若需要增大四倍光谱数据库，将待测样品的组分标定值进行顺序排布，每两个相邻标定值之间插入8个虚拟标定值，假定该两个相邻的标定值分别是d1，d2，则其标定值差值范围为d＝d
1-d2，继而可以得到插入的第一个虚拟标定值x1＝d/9+d2，第二个虚拟标定值x2＝2d/9+d2，依次类推可以得到8个虚拟标定值。同理，假定该两个待测样品的重组光谱数据集合矩阵分别为p1＝(k1，k2……k11
)，p2＝(s1，s2，
……s11
)，则其光谱数据集合矩阵差值范围为p＝p
1-p2＝(k
1-s1，k
2-s2，
……k11-s
11
)，继而得到第一个虚拟光谱数据y1＝p/9+p2＝((k
1-s1)/9+s1，(k
2-s2)/9+s2，
……
(k
11-s
11
)/9+s
11
)，第二个虚拟光谱数据y2＝2p/9+p2＝(2(k
1-s1)/9+s1，2(k
2-s2)/9+s2，
……
2(k
11-s
11
)/9+s
11
)，依次类推可以得到8个虚拟光谱数据。
30.图1中105是将虚拟光谱数据与虚拟组分标定值进行一一对应完成光谱数据库的扩充。将虚拟组分标定值按照顺序放入实际待测样品组分标定值之间，再将得到的虚拟光谱数据按照顺序放入替换后的实际待测样品光谱数据中，虚拟组分标定值与虚拟光谱数据进行一一对应，进而有效扩充样品光谱数据库。
31.在本实施例中，将步骤104中的实际标定值d1，d2以及通过计算得到的8个虚拟标定值x1，x2，
……
x8进行顺序排布，形成新的标定值集合d1，x1，x2，
……
x8，d2。再将步骤104中的实际重组后光谱数据p1，p2以及通过计算得到的8个虚拟光谱数据y1，y2，
……
y8进行顺序排布，形成新的光谱数据集合p1，y1，y2，
……
y8，d2。最后将新的标定值集合与新的光谱数据集合进行一一对应，即d1对应p1，x1对应y1，x8对应y8，并将对应完好之后的样品数据放入光谱
数据库中，进而有效扩充样品光谱数据库。
32.尽管这里参照本发明的解释性实施例对本发明进行了描述，上述实施例仅为本发明较佳的实施方式，本发明的实施方式并不受上述实施例的限制，应该理解，本领域技术人员可以设计出很多其他的修改和实施方式，这些修改和实施方式将落在本技术公开的原则范围和精神之内。