专利名称：辨识光谱的方法
技术领域：
本发明是有关于一种光谱分析的方法，且特别是有关于一种可自动识别光谱中的化合物的辨识光谱的方法。
背景技术：
傅立叶变换红外光谱(Fourier Transform Infrared Spectroscopy, FTIR)是一种用来获得红外线光谱的重要技术。有些化学物种会吸收特定频率的红外光辐射能量而改变其振动能阶，因而产生不同的光谱。具有FIlR吸收特性的化合物都有独特的光谱图形， 根据这个原理可以利用FIlR检测环境中的化合物。在许多现场量测中，因为有水气、二氧化碳等成分存在，会严重影响光谱图形品质，而且许多光谱图形也会因量测环境变异，造成基线的漂移。这些问题提高光谱图形解析的困难程度，需要更大量的人力与时间来解析光谱图形，严重影响仪器量测结果的时效性。而在传统FIlR分析上，一般须由专业及有经验的光谱图形判读人员来进行分析。 也就是说，FTIR分析对于分析人员的经验与专业程度相当的依赖，一般操作者不易上手使用，据此使得解析光谱图形的技术瓶颈偏高，并且未有客观的判读准则。此外，人力判读的速度有所限制，以致要进行大量光谱图形判读相当耗时。
发明内容
本发明提供一种辨识光谱的方法，提供一种自动化定性分析并且快速判读光谱的功能，据以提供分析者客观的参考。具体而言，本发明提出一种辨识光谱的方法，适于利用电子装置来辨识待测光谱的化合物组合。在光谱辨识方法中，首先对待测光谱执行预处理运算，而获得校正光谱。接着，载入多个化合物的标准光谱，将这些标准光谱对校正光谱执行回归(regression)运算以及假设检定(hypothesis test)，而获得这些化合物的相似度顺序。并且，由相似度顺序的最高者开始，逐一自这些化合物中择一加入至待测化合物组合，以计算待测化合物组合与校正光谱之间的等级相关系数(rank correlation coefficient)。如果目前所选择的待测化合物组合的等级相关系数大于或等于前次所选择的待测化合物组合的等级相关系数， 继续自上述化合物中择一加入至待测化合物组合，以计算待测化合物组合与校正光谱之间的等级相关系数。如果目前所选择的待测化合物组合的等级相关系数小于前次所选择的待测化合物组合的等级相关系数，输出前次所选择的待测化合物组合。在本发明的一实施例中，上述对待测光谱执行预处理运算，而获得校正光谱的步骤包括移除待测光谱的噪声，并且执行基线校正。在本发明的一实施例中，上述移除待测光谱的噪声的步骤包括利用小波算法， 将待测光谱分解为第一近似函数(approximation function)与第一补差函数(detail function)，并且，利用软式门槛值(soft thresholding)算法去除第一补差函数的噪声，而获得噪声移除函数。之后，利用小波还原算法，将第一近似函数与噪声移除函数还原为噪声移除光谱。在本发明的一实施例中，上述执行基线校正的步骤中，依据待测光谱的数据点数， 决定执行小波算法的次数。接着，以执行上述次数的小波算法所获得的第二近似函数作为基线函数。以此基线函数校正噪声移除光谱而获得校正光谱。在本发明的一实施例中，上述在以这些化合物的标准光谱对校正光谱进行回归运算及假设检定，而获得化合物相的相似度顺序中，将各标准光谱对校正光谱进行回归运算， 而获得每一个标准光谱的回归系数。依据假设检定以及这些回归系数，获得每一个标准光谱的估计值，而由各标准光谱的估计值来排序出相似度顺序。在本发明的一实施例中，在上述计算待测化合物组合与校正光谱之间的等级相关系数的步骤中，由待测化合物组合中的待测化合物所对应的标准光谱与校正光谱执行回归运算，而获得组合光谱。利用等级相关系数模型计算组合光谱与校正光谱之间的等级相关系数。在本发明的一实施例中，上述等级相关系数模型为肯德尔(Kendall tau)等级相关系数模型。而待测光谱为傅立叶变换红外光谱。基于上述，本发明依据不同化合物与待测光谱之间的相似度来进行排序，藉以逐一分析不同的化合物组合，以找出待测光谱中可能存在的最佳化合物组合。如此一来，可提供定性分析时较客观的参考，并且有效减少在进行光谱辨识时所需的工时。为让本发明的上述特征和优点能更明显易懂，下文特举实施例，并配合所附图式作详细说明如下。
图1是依照本发明一实施例的一种辨识光谱的方法流程图。图2是依照本发明一实施例的一种利用小波算法后的各解析度的近似函数的示意图。图3是依据本发明一实施例的一种基线校正的示意图。主要元件符号说明S105 S140 本发明一种辨识光谱的方法各步骤。
具体实施例方式在传统傅立叶变换红外光谱(Fourier Transform Infrared Spectroscopy, FTIR)分析上，对于分析人员的经验与专业程度相当的依赖，一般操作者不易上手使用，据此使得解析光谱图形的技术瓶颈偏高，并且未有客观的判读准则。据此，本发明提出一种辨识光谱的方法，可自动识别待测光谱中的化合物，据以提供分析者客观的参考。为了使本发明的内容更为明了，以下特举实施例作为本发明确实能够据以实施的范例。图1是依照本发明一实施例的一种辨识光谱的方法流程图。本实施例是利用电子装置来辨识待测光谱的化合物组合。在此，待测光谱例如为傅立叶变换红外光谱(FTIR)。请参照图1，在步骤S105中，接收待测光谱。接着，对待测光谱执行预处理运算，而获得校正光谱，如步骤SllO及步骤S115所示。由于原始的待测光谱常会有许多来自外界干扰的噪声，而这些噪声往往会影响之后光谱判读的准确性。因此，在接收待测光谱之后，可先对待测光谱执行预处理运算，藉以滤除噪声或进行基线校正等动作。在步骤Sl 10中，移除待测光谱的噪声。例如，利用小波变换(wavelet transform) 算法，将待测光谱分解为近似函数(approximation function)与补差函数(detail function) 0并且，利用软式门槛值算法(soft threshold)来去除补差函数的噪声，而获得噪声移除函数。之后，再利用小波还原算法，将近似函数与噪声移除函数还原为噪声移除光
■i並曰O一般而言，对于信号执行小波变换之后，会获得低频与高频两个部分，其中，高频部分称为近似函数，低频部分称为补差函数。近似函数记录了主要数据结构，乘载大部分数据信息，而补差函数则具有较多的噪声因子。因此，噪声移除主要是在补差函数这一部分进行。在此，利用软式门槛值(soft thresholding)算法来分割噪声与残留信号。而门槛值 (threshold)的选取是让门槛值具有极高的机率大于噪声系数(noise coefficient)的最
大值。软式门槛值算法如下所示
_
权利要求
1.一种辨识光谱的方法，适于利用电子装置来辨识待测光谱的化合物组合，其特征在于，该光谱辨识方法包括对该待测光谱执行预处理运算，而获得校正光谱；载入多个化合物的标准光谱，将所述标准光谱对该校正光谱执行回归运算以及假设检定，而获得所述化合物的相似度顺序；由该相似度顺序的最高者开始，逐一自所述化合物中择一加入至待测化合物组合，以计算该待测化合物组合与该校正光谱之间的等级相关系数；如果目前所选择的该待测化合物组合的所述级相关系数大于或等于前次所选择的该待测化合物组合的所述级相关系数，继续自所述化合物中择一加入至该待测化合物组合， 以计算该待测化合物组合与该校正光谱之间的所述级相关系数；以及如果目前所选择的该待测化合物组合的所述级相关系数小于前次所选择的该待测化合物组合的所述级相关系数，输出前次所选择的该待测化合物组合。
2.根据权利要求1所述的辨识光谱的方法，其特征在于，对该待测光谱执行该预处理运算，而获得该校正光谱的步骤包括移除该待测光谱的噪声；以及执行基线校正。
3.根据权利要求2所述的辨识光谱的方法，其特征在于，移除该待测光谱的噪声的步骤包括利用小波算法，将该待测光谱分解为第一近似函数与第一补差函数； 利用软式门槛值算法，去除该第一补差函数的噪声，而获得噪声移除函数；以及利用小波还原算法，将该第一近似函数与该噪声移除函数还原为噪声移除光谱。
4.根据权利要求3所述的辨识光谱的方法，其特征在于，执行该基线校正的步骤包括 依据该待测光谱的数据点数，决定执行该小波算法的次数；以执行该次数的该小波算法所获得的第二近似函数作为基线函数；以及以该基线函数校正该噪声移除光谱而获得校正光谱。
5.根据权利要求1所述的辨识光谱的方法，其特征在于，在将所述标准光谱对该校正光谱执行该回归运算以及该假设检定，而获得所述化合物的该相似度顺序的步骤包括将每一所述标准光谱对该校正光谱进行该回归运算，而获得每一所述标准光谱的回归系数；以及依据该假设检定以及所述回归系数，获得每一所述标准光谱的估计值，而由每一所述标准光谱的估计值排序出该相似度顺序。
6.根据权利要求1所述的辨识光谱的方法，其特征在于，计算该待测化合物组合与该校正光谱之间的所述级相关系数的步骤包括由该待测化合物组合中的待测化合物所对应的标准光谱与该校正光谱执行该回归运算，而获得组合光谱；以及利用等级相关系数模型计算该组合光谱与该校正光谱之间的所述级相关系数。
7.根据权利要求6所述的辨识光谱的方法，其特征在于，所述等级相关系数模型为肯德尔等级相关系数模型。
8.根据权利要求1所述的辨识光谱的方法，其特征在于，该待测光谱为傅立叶变换红外光谱。
全文摘要
一种辨识光谱的方法，适于利用电子装置来辨识待测光谱所包括的化合物。在光谱辨识方法中，首先对待测光谱执行预处理运算，而获得校正光谱。接着，对多个化合物的标准光谱执行回归运算以及假设检定，而获得这些化合物的相似度顺序。并且，由相似度顺序的最高者开始，逐一自这些化合物中择一加入至待测化合物组合，以计算待测化合物组合与校正光谱之间的等级相关系数。如果目前所选择的待测化合物组合的等级相关系数小于前次所选择的待测化合物组合的等级相关系数，输出目前所选择的待测化合物组合。
文档编号G01N21/35GK102542284SQ20111003952
公开日2012年7月4日 申请日期2011年2月14日 优先权日2010年12月13日
发明者卢鸿兴, 宋隆佑, 杜敬民, 蔡孟原, 谢瑞豪, 陈奕仑 申请人:财团法人工业技术研究院