CN108195817A - 用于去除溶剂干扰的拉曼光谱检测方法 - Google Patents

用于去除溶剂干扰的拉曼光谱检测方法 Download PDF

Info

Publication number
CN108195817A
CN108195817A CN201611126188.XA CN201611126188A CN108195817A CN 108195817 A CN108195817 A CN 108195817A CN 201611126188 A CN201611126188 A CN 201611126188A CN 108195817 A CN108195817 A CN 108195817A
Authority
CN
China
Prior art keywords
spectral signal
raman
solvent
raman spectral
interference
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201611126188.XA
Other languages
English (en)
Other versions
CN108195817B (zh
Inventor
王红球
苟巍
陈卓
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nuctech Co Ltd
Original Assignee
Nuctech Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nuctech Co Ltd filed Critical Nuctech Co Ltd
Priority to CN201611126188.XA priority Critical patent/CN108195817B/zh
Priority to PCT/CN2017/112916 priority patent/WO2018103541A1/zh
Publication of CN108195817A publication Critical patent/CN108195817A/zh
Application granted granted Critical
Publication of CN108195817B publication Critical patent/CN108195817B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/65Raman scattering

Landscapes

  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

本发明公开了一种用于去除溶剂干扰的拉曼光谱检测方法,包括以下步骤:对包含溶剂和样品的溶液的拉曼光谱进行测量,以获得溶液的拉曼光谱信号;对溶剂的拉曼光谱进行测量,以获得溶剂的拉曼光谱信号;从溶液的拉曼光谱信号中逐步减去溶剂的拉曼光谱信号,以获得一系列去除溶剂干扰的拉曼光谱信号;计算一系列去除溶剂干扰的拉曼光谱信号中的每一个的信息熵;选择所计算出的一系列的信息熵中的最大信息熵;和将与该最大信息熵对应的去除溶剂干扰的拉曼光谱信号作为优化的去除溶剂干扰的拉曼光谱信号。

Description

用于去除溶剂干扰的拉曼光谱检测方法
技术领域
本发明涉及拉曼光谱检测领域,尤其涉及一种用于去除溶剂干扰的拉曼光谱检测方法。
背景技术
拉曼光谱是一种分子振动光谱,它可以反映分子的指纹特征,可用于对物质的检测。拉曼光谱检测通过检测待测物对于激发光的拉曼散射效应所产生的拉曼光谱来检测和识别物质。拉曼光谱检测方法已经广泛应用于液体安检、珠宝检测、爆炸物检测、毒品检测、药品检测、农药残留检测等领域。
拉曼光谱检测方法在物质检测中遇到的一个问题是对于溶于溶剂中的样品的检测。由于待测样品有时需要溶于溶剂中,例如在对农药残留进行检测时,因此,不可避免地会使检测受到溶剂的干扰。对于溶于溶剂中的样品,要想采集到它的拉曼光谱,激发光必须同时激发溶剂和样品,而溶剂被光激发后自身会产生一些信号,这会导致采集到的光谱和被测样品本身的拉曼光谱有很大差异,进而导致很多物质无法准确识别。目前对于溶剂的干扰,通常通过主观判断给出一个参数作为去除溶剂干扰的比例系数,任何包含溶剂的光谱信号在去除溶剂干扰时都设定相同的参数,不能很好地适应光谱信号的变化,从而也不能准确地去除溶剂干扰。
因此,如何在拉曼光谱检测中去除溶剂干扰,以准确获得待测样品的拉曼光谱,进而实现对待测物的准确识别,对提高拉曼光谱检测方法在各个应用领域检测的准确性具有重要意义。
发明内容
本发明的目的是提供一种用于去除溶剂干扰的拉曼光谱检测方法,其能够有效地去除溶剂的拉曼光谱对于样品的拉曼光谱信号的干扰,以准确地对样品进行检测和识别。
为了实现上述发明目的,本发明的技术方案通过以下方式来实现:
根据本发明的一个方面,提供一种用于去除溶剂干扰的拉曼光谱检测方法,包括以下步骤:
溶液的拉曼光谱信号获得步骤:对包含溶剂和样品的溶液的拉曼光谱进行测量,以获得溶液的拉曼光谱信号;
溶剂的拉曼光谱信号获得步骤:对溶剂的拉曼光谱进行测量,以获得溶剂的拉曼光谱信号;
去除溶剂干扰的拉曼光谱信号获得步骤:从溶液的拉曼光谱信号中逐步减去溶剂的拉曼光谱信号,以获得一系列去除溶剂干扰的拉曼光谱信号;
信息熵计算步骤:计算一系列去除溶剂干扰的拉曼光谱信号中的每一个的信息熵;
最大信息熵选取步骤:选择所计算出的一系列的信息熵中的最大信息熵;和
优化的去除溶剂干扰的拉曼光谱信号获得步骤:将与该最大信息熵对应的去除溶剂干扰的拉曼光谱信号作为优化的去除溶剂干扰的拉曼光谱信号。
根据一些实施例,所述的去除溶剂干扰的拉曼光谱信号获得步骤包括:
从溶液的拉曼光谱信号中减去比例系数乘以溶剂的拉曼光谱信号的乘积,并且逐步改变比例系数,以从溶液的拉曼光谱信号中逐步减去溶剂的拉曼光谱信号。
根据一些实施例,溶液的拉曼光谱信号为离散数据序列A,溶剂的拉曼光谱信号为离散数据序列B,所述的从溶液的拉曼光谱信号中减去比例系数乘以溶剂的拉曼光谱信号的乘积的步骤用下式表示:
C=A-j*K*B;
其中,C为表示去除溶剂干扰的拉曼光谱信号的离散数据序列,j依次取值为1,2,3,……N且j为自然数,N为预设的计算次数,且K为预设的去除比例。
根据一些实施例,所述拉曼光谱检测方法还包括如下步骤:计算区间确定步骤:确定溶剂的拉曼光谱的特征峰的位置,并且将包括溶剂的拉曼光谱的特征峰的位置的区间作为计算区间;
其中,所述去除溶剂干扰的拉曼光谱信号获得步骤包括:在所述计算区间内,从溶液的拉曼光谱信号中逐步减去溶剂的拉曼光谱信号,以获得一系列去除溶剂干扰的拉曼光谱信号。
根据一些实施例,所述的信息熵计算步骤包括:
根据下面的信息熵计算公式计算一系列去除溶剂干扰的拉曼光谱信号C中的每一个的信息熵:
其中,i表示该拉曼光谱信号的第i个波数,n表示该拉曼光谱信号的信号长度,xi表示第i个波数对应的强度,p(xi)表示该拉曼光谱信号中取强度xi的概率。
根据一些实施例,所述拉曼光谱检测方法还包括如下步骤:在从溶液的拉曼光谱信号中逐步减去溶剂的拉曼光谱信号之前,对溶剂的拉曼光谱信号与溶液的拉曼光谱信号进行归一化处理。
根据一些实施例,K的取值范围为0.005-0.03,且N的取值范围为200-600。
根据本发明的另一方面,还提供一种用于去除溶剂干扰的拉曼光谱检测方法,包括如下步骤:
溶液的拉曼光谱信号获得步骤:对包含溶剂和样品的溶液的拉曼光谱进行测量,以获得溶液的拉曼光谱信号;
溶剂的拉曼光谱信号获得步骤:对溶剂的拉曼光谱进行测量,以获得溶剂的拉曼光谱信号;
归一化步骤:对溶剂的拉曼光谱信号和溶液的拉曼光谱信号作归一化处理;
计算区间确定步骤:确定溶剂的拉曼光谱的特征峰的位置,并且将包括溶剂的拉曼光谱的特征峰的位置的区间作为计算区间;
去除溶剂干扰的区间拉曼光谱信号获得步骤:在计算区间内,从归一化后的溶液的拉曼光谱信号中逐步减去归一化后的溶剂的拉曼光谱信号与比例系数的乘积,通过逐渐改变比例系数,以获得一系列去除溶剂干扰的区间拉曼光谱信号;
信息熵计算步骤:计算一系列去除溶剂干扰的区间拉曼光谱信号中的每一个的信息熵;
优化的比例系数确定步骤:选择所计算出的一系列的信息熵中的最大信息熵,并且选择与最大信息熵对应的比例系数作为优化的比例系数;和
优化的去除溶剂干扰的拉曼光谱信号获得步骤:从归一化后的溶液的拉曼光谱信号中减去归一化后的溶剂的拉曼光谱信号与优化的比例系数的乘积,以获得优化的去除溶剂干扰的拉曼光谱信号。
根据本发明的又一方面,还提供一种电子设备,包括:
存储器,用于存储可执行指令;以及
处理器,用于执行存储器中存储的可执行指令,以执行如上述方面或实施例中任一项所述的方法。
本发明的上述技术方案中的任何一个能够借助计算最大信息熵,从混有溶剂光谱的样品拉曼光谱信号中去除溶剂所导致的干扰信息,从而能够正确获得反映样品自身的属性的拉曼光谱信号。这种方案可以对溶液中的样品的拉曼光谱进行准确的检测,从而有效地识别待测样品。
附图说明
图1示意性地示出根据本发明的实施例的用于去除溶剂干扰的拉曼光谱检测方法的流程图;
图2示意性地示出根据本发明的另一实施例的用于去除溶剂干扰的拉曼光谱检测方法的流程图;
图3示意性地示出将根据本发明的实施例的拉曼光谱检测方法用于去除乙腈干扰的流程图;
图4示意性示出亚胺硫磷乙腈溶液的拉曼光谱;
图5示意性示出乙腈溶剂的拉曼光谱;
图6示意性示出根据本发明实施例的方法计算出的信息熵曲线图;
图7示意性示出通过使用本发明的实施例的拉曼光谱检测方法从亚胺硫磷乙腈溶液中去除乙腈干扰后的溶液的拉曼光谱;和
图8是示出了用于执行根据本发明的实施例的方法的电子设备的示例硬件布置的框图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。在说明书中,相同或相似的附图标号表示相同或相似的部件。下述参照附图对本发明实施方式的说明旨在对本发明的总体发明构思进行解释,而不应当理解为对本发明的一种限制。
在本文中,为了描述方便,使用“第一、第二”、“A、B、C”或“S1/S10、S2/S20”等表述描述方法的步骤,但是,除非有特别说明,这样的表述不应理解为对步骤执行顺序的限制。
在利用拉曼光谱对样品进行检测时,待测样品有时需要溶于溶剂中,例如在检测农药残留时,将待测样品溶于例如乙腈的溶剂中。因而,对于样品的检测可能需要用于检测的激发光同时照射到样品和溶剂上来进行,而溶剂也可能对激发光产生拉曼散射效应,在此情况下,溶剂对样品自身的拉曼光谱信号会产生干扰。而去除这种干扰对于准确有效地检测和识别样品是非常重要的。
图1示意性地示出根据本发明的实施例的用于去除溶剂干扰的拉曼光谱检测方法的流程图。该方法包括:
溶液的拉曼光谱信号获得步骤:对包含溶剂和样品的溶液的拉曼光谱进行测量,以获得溶液的拉曼光谱信号A;
溶剂的拉曼光谱信号获得步骤:对溶剂的拉曼光谱进行测量,以获得溶剂的拉曼光谱信号B;
去除溶剂干扰的拉曼光谱信号获得步骤:从溶液的拉曼光谱信号A中逐步减去溶剂的拉曼光谱信号B,以获得一系列去除溶剂干扰的拉曼光谱信号C;
信息熵计算步骤:计算一系列去除溶剂干扰的拉曼光谱信号C中的每一个的信息熵;
最大信息熵选取步骤:选择所计算出的一系列的信息熵中的最大信息熵;和
优化的去除溶剂干扰的拉曼光谱信号获得步骤:将与该最大信息熵对应的去除溶剂干扰的拉曼光谱信号作为优化的去除溶剂干扰的拉曼光谱信号。
具体地,在所述溶液的拉曼光谱信号获得步骤中,获得的溶液的拉曼光谱信号A是一个离散数据序列,记为Ai,其中,i依次取值为1,2,3,……n且i为自然数,n表示该数据序列的数据点总数。例如,离散数据序列Ai可以是由一组离散的溶液的拉曼光谱信号的峰强度数据构成的向量或矩阵。同样地,在溶剂的拉曼光谱信号获得步骤中,获得的溶剂的拉曼光谱信号B也是一个离散数据序列,记为Bi,例如,由一组离散的溶剂的拉曼光谱信号的峰强度数据构成的向量或矩阵。
在一个示例中,所述的去除溶剂干扰的拉曼光谱信号获得步骤具体包括:从溶液的拉曼光谱信号A中逐步减去j*K*溶剂的拉曼光谱信号B,以获得一系列去除溶剂干扰的拉曼光谱信号Cj。其中,j依次取值为1,2,3,……N且j为自然数,N为预设的计算次数,且K为预设的比例系数。即,Cj=A-j*K*B,j依次取值为1,2,3,……N且j为自然数。相应地,所述一系列去除溶剂干扰的拉曼光谱信号C包括C1、C2、C3、……,CN。应该理解,每一个去除溶剂干扰的拉曼光谱信号Cj也是一个离散数据序列,其也包括n个数据点。
根据本发明的实施例,所述的信息熵计算步骤包括:根据下面的信息熵计算公式计算一系列去除溶剂干扰的拉曼光谱信号中的每一个Cj的信息熵:
其中,i表示该拉曼光谱信号的第i个波数,n表示该拉曼光谱信号的信号长度,xi表示第i个波数对应的强度,p(xi)表示该拉曼光谱信号中取强度xi的概率。
根据信息熵理论,信息熵评价的是随机变量(例如X)等于各个值的平均信息量,即表示随机变量X的不确定度。而对于离散数据序列而言,因为离散情况的均匀分布没有任何的偏向,所以均匀分布应是最不确定的。例如,考虑二元的情况,就是随机变量X只能有X=a或者X=b两种情况,如果p(X=a)概率很大,接近于1,同时p(X=b)概率很小,那么在此情况下,由于比较肯定X=a的发生,所以X的不确定度很小是明显的。因此,应该理解,二元在发生概率相同下应该是最混乱的状态。也就是说,随机变量X在均匀的离散分布下信息熵最大。
具体地,在去除溶剂干扰信息时,由于溶剂本身的拉曼光谱信号具有一些特征峰,所以从包含溶剂和样品的拉曼光谱信号中去除溶剂的拉曼光谱信号的过程是去除一些特征峰的过程,即,使拉曼光谱信号变平缓的过程;进一步地,如果过度去除溶剂干扰信息,那么会使得拉曼光谱增加了一些反向峰,即,使拉曼光谱信号又变得不平缓。而根据上述理论可知,平缓的拉曼光谱信号对应比较均匀的离散分布,也就对应最大信息熵。因此,在逐步去除溶剂干扰的过程中,当获得最大信息熵时,意味着此时恰好去除掉溶剂干扰信息,即,此时对应优化的从包含溶剂和样品的拉曼光谱信号中去除溶剂的拉曼光谱信号的情形。
在一些实施例中,上述的用于去除溶剂干扰的拉曼光谱检测方法还可以包括如下步骤:对溶剂的拉曼光谱信号B和溶液的拉曼光谱信号A进行归一化处理。归一化化是将需要处理的数据通过某种归一化算法处理后限制在预定的范围内,其使得处于两个不同参考系下的数据可以进行比较、计算等,以便于后续数据处理,并且加快后续计算的收敛。例如,在一个示例中,可以将溶剂的拉曼光谱信号B归一化到溶液的拉曼光谱信号A中,以便于后续步骤中进行的从溶液的拉曼光谱信号A中逐步减去溶剂的拉曼光谱信号B的计算。
根据本发明的一些实施例,可以仅在溶剂的拉曼光谱信号的特征峰所在的区间内进行上述的信息熵计算步骤,以减少计算量,从而加快检测速度。图2示意性地示出根据本发明的另一实施例的用于去除溶剂干扰的拉曼光谱检测方法的流程图,这样,该用于去除溶剂干扰的拉曼光谱检测方法可以包括如下步骤:
溶液的拉曼光谱信号获得步骤:对包含溶剂和样品的溶液的拉曼光谱进行测量,以获得溶液的拉曼光谱信号A;
溶剂的拉曼光谱信号获得步骤:对溶剂的拉曼光谱进行测量,以获得溶剂的拉曼光谱信号B;
归一化步骤:归一化溶剂的拉曼光谱信号B和溶液的拉曼光谱信号A;
计算区间确定步骤:确定溶剂的拉曼光谱的特征峰的位置,并且将包括溶剂的拉曼光谱的特征峰的位置的区间作为计算区间;
去除溶剂干扰的区间拉曼光谱信号获得步骤:在计算区间内,从归一化后的溶液的拉曼光谱信号A中逐步减去归一化后的溶剂的拉曼光谱信号B与预设的比例系数的乘积,通过逐渐改变比例系数,以获得一系列去除溶剂干扰的区间拉曼光谱信号C’;
信息熵计算步骤:计算一系列去除溶剂干扰的区间拉曼光谱信号C’中的每一个的信息熵;
优化的比例系数确定步骤:选择所计算出的一系列的信息熵中的最大信息熵,并且选择与最大信息熵对应的比例系数作为优化的比例系数;和
优化的去除溶剂干扰的拉曼光谱信号获得步骤:从归一化后的溶液的拉曼光谱信号A中减去归一化后的溶剂的拉曼光谱信号B与优化的比例系数的乘积,以获得优化的去除溶剂干扰的拉曼光谱信号。
下面,以溶剂为乙腈并且样品为亚胺硫磷为示例,结合图3-7,详细描述根据本发明实施例的用于去除溶剂干扰的拉曼光谱检测方法。如图3所示,该方法可以包括如下步骤:
亚胺硫磷乙腈溶液的拉曼光谱信号获得步骤:对亚胺硫磷乙腈溶液的拉曼光谱进行测量,以获得亚胺硫磷乙腈溶液的拉曼光谱信号A,如图4所示,在图4所示的拉曼光谱图中,横坐标表示拉曼频移或波数(cm-1),纵坐标表示强度(无量纲),这样,拉曼光谱信号A可以是包括若干个离散的图4中示出的纵坐标表示的强度数据的向量,或是包括若干个离散的图4中示出的横坐标表示的拉曼频移和对应的纵坐标表示的强度数据的矩阵。
乙腈的拉曼光谱信号获得步骤:对乙腈本身的拉曼光谱进行测量,以获得乙腈的拉曼光谱信号B,如图5所示,图5中横、纵坐标的含义与图4相同,同样地,拉曼光谱信号B可以是包括若干个离散的图5中示出的纵坐标表示的强度数据的向量,或是包括若干个离散的图5中示出的横坐标表示的拉曼频移和对应的纵坐标表示的强度数据的矩阵。
归一化步骤:归一化乙腈的拉曼光谱信号B和亚胺硫磷乙腈溶液的拉曼光谱信号A。具体地,如果拉曼光谱信号B与拉曼光谱信号A的横坐标不对应,即拉曼光谱信号B与拉曼光谱信号A包括不处于相同拉曼频移位置处的强度数据,那么首先可以根据拉曼光谱信号A、B的原始数据,通过插值计算,计算出其它拉曼频移位置处的强度数据,从而将拉曼光谱信号A、B转换到相同的拉曼频移位置处。例如,在一个示例中,可以设定拉曼频移位置为[350:2:2800]cm-1,其中,350cm-1、2800cm-1分别表示最小、最大拉曼频移位置,2表示间隔步长,即各个拉曼频移位置为350、352、354、356cm-1,以此类推,直至2800cm-1,然后,如果拉曼光谱信号A、B的原始数据中不包括这些拉曼频移位置处的强度数据,那么可以根据拉曼光谱信号A、B的原始数据,通过插值计算,计算出这些拉曼频移位置处的强度数据,从而将拉曼光谱信号A、B转换到相同的拉曼频移位置[350:2:2800]cm-1处,将转换后的拉曼光谱信号A、B分别计为拉曼光谱信号A’、B’。并且,如果拉曼光谱信号B与拉曼光谱信号A的纵坐标不对应,即拉曼光谱信号B与拉曼光谱信号A的强度数据不位于同一参考系中或以不同的量程测量,可以将拉曼光谱信号B’转换到拉曼光谱信号A’中。例如,在一个示例中,可以分别选取拉曼光谱信号A’、B’的纵坐标的最大值max(A’)、max(B’),然后将max(A’)与max(B’)的比值作为纵坐标的归一化系数,使拉曼光谱信号B’的离散数据分别乘以该归一化系数,以将拉曼光谱信号B’转换到拉曼光谱信号A’中,将转换后的拉曼光谱信号B’记为拉曼光谱信号B”。在其它实施例中,也可以保持拉曼光谱信号A不变,直接将拉曼光谱信号B归一化到拉曼光谱信号A中。
计算区间确定步骤:确定乙腈的特征峰的位置,并根据乙腈的特征峰的位置确定包括这些特征峰的位置的区间,以作为计算区间。在图5所示的实施例中,可以确定乙腈的特征峰的位置主要处在900-950cm-1和2200-2300cm-1的区间内,所以可以将区间900-950cm-1和2200-2300cm-1确定为计算区间。
去除乙腈干扰的区间拉曼光谱信号获得步骤:在计算区间900-950cm-1和2200-2300cm-1内,从归一化后的亚胺硫磷乙腈溶液的拉曼光谱信号A’中逐步减去归一化后的乙腈的拉曼光谱信号B”与预设的比例系数的乘积,通过逐渐改变比例系数,以获得一系列去除乙腈干扰的区间拉曼光谱信号C。
信息熵计算步骤:计算一系列去除乙腈干扰的区间拉曼光谱信号C中的每一个的信息熵;
具体地,上述两个步骤的计算过程可以用下列公式表示:
Cj=A-j*K*B
其中,j依次取值为1,2,3,……N且j为自然数,N为预设的计算次数,且K为预设的除去比例,通过逐步改变j,就可以实现j*K的逐步改变,即j*K构成了被逐步改变的比例系数KK。在一些实施例中,K的取值范围为0.005-0.03,且N的取值范围为200-600。优选地,在一个示例中,N设为300,K设为0.01,当N、K取该数值时,能够较快地计算出最大信息熵,又不至于过度增加计算时间。这样,通过依次改变j的取值,可以以0.01为步长,逐步从亚胺硫磷乙腈溶液的拉曼光谱信号A’中减去乙腈的拉曼光谱信号B”,从而获得300个去除乙腈干扰的区间拉曼光谱信号C。
然后,根据下面的信息熵计算公式计算每一个去除乙腈干扰的区间拉曼光谱信号Cj的信息熵,由于存在300个去除乙腈干扰的区间拉曼光谱信号C,所以对应地可以计算出300个信息熵,如图6所示,其示出了以图4、图5所示的拉曼光谱信号为原始数据并通过上述计算方法计算出的信息熵的变化图,在图6中,横坐标表示拉曼光谱信号的序号,纵坐标表示计算出的信息熵。
其中,i表示该拉曼光谱信号的第i个波数,n表示该拉曼光谱信号的信号长度,xi表示第i个波数对应的强度,p(xi)表示该拉曼光谱信号中取强度xi的概率。
优化的比例系数确定步骤:选择所计算出的300个信息熵中的最大信息熵,并且选择与该最大信息熵对应的比例系数作为优化的比例系数,将该优化的比例系数记为KK优化;如图6所示,最大信息熵出现在j=60的位置处,那么KK优化=0.01*60。
在上面的计算过程中,通过将计算区间设定为包括乙腈的特征峰的位置的区间,可以极大减少计算量,从而加快检测速度。
优化的去除乙腈干扰的拉曼光谱信号获得步骤:从亚胺硫磷乙腈溶液的拉曼光谱信号A’中减去乙腈的拉曼光谱信号B”与优化的比例系数的乘积,以获得优化的去除乙腈干扰的拉曼光谱信号C优化。图示的实施例中,C优化=A-60*0.01*B,通过将图4、图5对应的数据代入该公式中,可以得到优化的去除乙腈干扰的拉曼光谱信号C优化的数据,根据该数据绘制的优化的去除乙腈干扰的拉曼谱图如图7所示,比较图4和图7,可以看出,该拉曼光谱信号C优化较好去除了乙腈干扰信号。
根据本发明的又一实施例,还提供一种电子设备,图8是示出了该电子设备的示例硬件布置800的框图。硬件布置800包括处理器806(例如,微处理器(μP)、数字信号处理器(DSP)等)。处理器806可以是用于执行本文描述的方法步骤的不同动作的单一处理单元或者是多个处理单元。布置800还可以包括用于从其他实体接收信号的输入单元802、以及用于向其他实体提供信号的输出单元804。输入单元802和输出单元804可以被布置为单一实体或者是分离的实体。
此外,布置800可以包括具有非易失性或易失性存储器形式的至少一个可读存储介质808,例如是电可擦除可编程只读存储器(EEPROM)、闪存、和/或硬盘驱动器。可读存储介质808包括计算机程序810,该计算机程序810包括代码/计算机可读指令,其在由布置800中的处理器806执行时使得硬件布置800和/或包括硬件布置800在内的设备可以执行例如上面结合图1-3所描述的流程及其任何变形。
计算机程序810可被配置为具有例如计算机程序模块810A~810C架构的计算机程序代码。因此,在例如设备中使用硬件布置800时的示例实施例中,布置800的计算机程序中的代码包括:模块810A,用于…。计算机程序中的代码还包括:模块810B,用于…。计算机程序中的代码还包括:模块810C,用于…。
计算机程序模块实质上可以执行图1-3中所示出的流程中的各个动作,以模拟设备。换言之,当在处理器806中执行不同计算机程序模块时,它们可以对应于设备中的上述不同单元。
尽管上面结合图8所公开的实施例中的代码手段被实现为计算机程序模块,其在处理器806中执行时使得硬件布置800执行上面结合图1~3所描述的动作,然而在备选实施例中,该代码手段中的至少一项可以至少被部分地实现为硬件电路。
处理器可以是单个CPU(中央处理单元),但也可以包括两个或更多个处理单元。例如,处理器可以包括通用微处理器、指令集处理器和/或相关芯片组和/或专用微处理器(例如,专用集成电路(ASIC))。处理器还可以包括用于缓存用途的板载存储器。计算机程序可以由连接到处理器的计算机程序产品来承载。计算机程序产品可以包括其上存储有计算机程序的计算机可读介质。例如,计算机程序产品可以是闪存、随机存取存储器(RAM)、只读存储器(ROM)、EEPROM,且上述计算机程序模块在备选实施例中可以用UE内的存储器的形式被分布到不同计算机程序产品中。
本领域技术人员应当理解,在本发明的实施例中,虽然以亚胺硫磷乙腈溶液为示例详细说明了本发明的技术构思,但是本发明不局限于去除乙腈干扰。
虽然结合附图对本发明进行了说明,但是附图中公开的实施例旨在对本发明优选实施方式进行示例性说明,而不能理解为对本发明的一种限制。
虽然本发明总体构思的一些实施例已被显示和说明,本领域普通技术人员将理解,在不背离本总体发明构思的原则和精神的情况下,可对这些实施例做出改变,本发明的范围以权利要求和它们的等同物限定。

Claims (9)

1.一种用于去除溶剂干扰的拉曼光谱检测方法,包括以下步骤:
溶液的拉曼光谱信号获得步骤:对包含溶剂和样品的溶液的拉曼光谱进行测量,以获得溶液的拉曼光谱信号;
溶剂的拉曼光谱信号获得步骤:对溶剂的拉曼光谱进行测量,以获得溶剂的拉曼光谱信号;
去除溶剂干扰的拉曼光谱信号获得步骤:从溶液的拉曼光谱信号中逐步减去溶剂的拉曼光谱信号,以获得一系列去除溶剂干扰的拉曼光谱信号;
信息熵计算步骤:计算一系列去除溶剂干扰的拉曼光谱信号中的每一个的信息熵;
最大信息熵选取步骤:选择所计算出的一系列的信息熵中的最大信息熵;和
优化的去除溶剂干扰的拉曼光谱信号获得步骤:将与该最大信息熵对应的去除溶剂干扰的拉曼光谱信号作为优化的去除溶剂干扰的拉曼光谱信号。
2.根据权利要求1所述的拉曼光谱检测方法,其特征在于,所述的去除溶剂干扰的拉曼光谱信号获得步骤包括:
从溶液的拉曼光谱信号中减去比例系数乘以溶剂的拉曼光谱信号的乘积,并且逐步改变比例系数,以从溶液的拉曼光谱信号中逐步减去溶剂的拉曼光谱信号。
3.根据权利要求2所述的拉曼光谱检测方法,其特征在于,溶液的拉曼光谱信号为离散数据序列A,溶剂的拉曼光谱信号为离散数据序列B,所述的从溶液的拉曼光谱信号中减去比例系数乘以溶剂的拉曼光谱信号的乘积的步骤用下式表示:
C=A-j*K*B;
其中,C为表示去除溶剂干扰的拉曼光谱信号的离散数据序列,j依次取值为1,2,3,……N且j为自然数,N为预设的计算次数,且K为预设的去除比例。
4.根据权利要求1-3中任一项所述的拉曼光谱检测方法,其特征在于,还包括如下步骤:计算区间确定步骤:确定溶剂的拉曼光谱的特征峰的位置,并且将包括溶剂的拉曼光谱的特征峰的位置的区间作为计算区间;
其中,所述去除溶剂干扰的拉曼光谱信号获得步骤包括:在所述计算区间内,从溶液的拉曼光谱信号中逐步减去溶剂的拉曼光谱信号,以获得一系列去除溶剂干扰的拉曼光谱信号。
5.根据权利要求1-4中任一项所述的拉曼光谱检测方法,其特征在于,所述的信息熵计算步骤包括:
根据下面的信息熵计算公式计算一系列去除溶剂干扰的拉曼光谱信号C中的每一个的信息熵:
其中,i表示该拉曼光谱信号的第i个波数,n表示该拉曼光谱信号的信号长度,xi表示第i个波数对应的强度,p(xi)表示该拉曼光谱信号中取强度xi的概率。
6.根据权利要求1所述的拉曼光谱检测方法,其特征在于,还包括如下步骤:
在从溶液的拉曼光谱信号中逐步减去溶剂的拉曼光谱信号之前,对溶剂的拉曼光谱信号与溶液的拉曼光谱信号进行归一化处理。
7.根据权利要求3所述的拉曼光谱检测方法,其特征在于,K的取值范围为0.005-0.03,且N的取值范围为200-600。
8.一种用于去除溶剂干扰的拉曼光谱检测方法,包括如下步骤:
溶液的拉曼光谱信号获得步骤:对包含溶剂和样品的溶液的拉曼光谱进行测量,以获得溶液的拉曼光谱信号;
溶剂的拉曼光谱信号获得步骤:对溶剂的拉曼光谱进行测量,以获得溶剂的拉曼光谱信号;
归一化步骤:对溶剂的拉曼光谱信号和溶液的拉曼光谱信号作归一化处理;
计算区间确定步骤:确定溶剂的拉曼光谱的特征峰的位置,并且将包括溶剂的拉曼光谱的特征峰的位置的区间作为计算区间;
去除溶剂干扰的区间拉曼光谱信号获得步骤:在计算区间内,从归一化后的溶液的拉曼光谱信号中逐步减去归一化后的溶剂的拉曼光谱信号与比例系数的乘积,通过逐渐改变比例系数,以获得一系列去除溶剂干扰的区间拉曼光谱信号;
信息熵计算步骤:计算一系列去除溶剂干扰的区间拉曼光谱信号中的每一个的信息熵;
优化的比例系数确定步骤:选择所计算出的一系列的信息熵中的最大信息熵,并且选择与最大信息熵对应的比例系数作为优化的比例系数;和
优化的去除溶剂干扰的拉曼光谱信号获得步骤:从归一化后的溶液的拉曼光谱信号中减去归一化后的溶剂的拉曼光谱信号与优化的比例系数的乘积,以获得优化的去除溶剂干扰的拉曼光谱信号。
9.一种电子设备,包括:
存储器,用于存储可执行指令;以及
处理器,用于执行存储器中存储的可执行指令,以执行如权利要求1-8中任一项所述的方法。
CN201611126188.XA 2016-12-08 2016-12-08 用于去除溶剂干扰的拉曼光谱检测方法 Active CN108195817B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201611126188.XA CN108195817B (zh) 2016-12-08 2016-12-08 用于去除溶剂干扰的拉曼光谱检测方法
PCT/CN2017/112916 WO2018103541A1 (zh) 2016-12-08 2017-11-24 用于去除溶剂干扰的拉曼光谱检测方法和电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611126188.XA CN108195817B (zh) 2016-12-08 2016-12-08 用于去除溶剂干扰的拉曼光谱检测方法

Publications (2)

Publication Number Publication Date
CN108195817A true CN108195817A (zh) 2018-06-22
CN108195817B CN108195817B (zh) 2020-09-15

Family

ID=62490664

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611126188.XA Active CN108195817B (zh) 2016-12-08 2016-12-08 用于去除溶剂干扰的拉曼光谱检测方法

Country Status (2)

Country Link
CN (1) CN108195817B (zh)
WO (1) WO2018103541A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109632680A (zh) * 2019-01-29 2019-04-16 吉林大学 一种基于排列熵的水体中磷的检测方法
CN109752364A (zh) * 2019-03-11 2019-05-14 广西科技大学 多组分体系拉曼光谱的不同步褪色荧光消除方法
CN114184600A (zh) * 2021-12-30 2022-03-15 江苏海洋大学 一种基于拉曼光谱背景扣除的水中溶质定量方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115736290B (zh) * 2022-11-07 2023-07-28 南京邦康生物技术有限公司 一种增加骨密度产品的生产工艺控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6867858B2 (en) * 2002-02-15 2005-03-15 Kaiser Optical Systems Raman spectroscopy crystallization analysis method
WO2005090948A3 (en) * 2003-12-29 2006-02-16 Intel Corp Composite organic-inorganic nanoclusters
CN102879378A (zh) * 2011-07-13 2013-01-16 赛默科技便携式分析仪器有限公司 通过拉曼光谱法从含有干扰性荧光污染物的不纯组合物中检测海洛因
CN104749155A (zh) * 2013-12-27 2015-07-01 同方威视技术股份有限公司 用于检测容装体中的样品的拉曼光谱检测方法
CN105699356A (zh) * 2016-02-19 2016-06-22 广西科技大学 通过信息熵判断拉曼光谱的荧光消除程度的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104122246B (zh) * 2013-04-28 2017-03-29 同方威视技术股份有限公司 不同基质乳制品中三聚氰胺含量的拉曼光谱测量方法
CN103439310A (zh) * 2013-09-02 2013-12-11 中国检验检疫科学研究院 一种快速检测三甘醇中水分含量的拉曼光谱方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6867858B2 (en) * 2002-02-15 2005-03-15 Kaiser Optical Systems Raman spectroscopy crystallization analysis method
WO2005090948A3 (en) * 2003-12-29 2006-02-16 Intel Corp Composite organic-inorganic nanoclusters
CN102879378A (zh) * 2011-07-13 2013-01-16 赛默科技便携式分析仪器有限公司 通过拉曼光谱法从含有干扰性荧光污染物的不纯组合物中检测海洛因
CN105842223A (zh) * 2011-07-13 2016-08-10 赛默科技便携式分析仪器有限公司 通过拉曼光谱法从含有干扰性荧光污染物的不纯组合物中检测海洛因
CN104749155A (zh) * 2013-12-27 2015-07-01 同方威视技术股份有限公司 用于检测容装体中的样品的拉曼光谱检测方法
CN105699356A (zh) * 2016-02-19 2016-06-22 广西科技大学 通过信息熵判断拉曼光谱的荧光消除程度的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周广涛等: "基于熵的SAR 干扰效果评估方法", 《航天电子对抗》 *
朱自莹等: "《拉曼光谱在化学中的应用》", 31 December 1998 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109632680A (zh) * 2019-01-29 2019-04-16 吉林大学 一种基于排列熵的水体中磷的检测方法
CN109632680B (zh) * 2019-01-29 2021-07-20 吉林大学 一种基于排列熵的水体中磷的检测方法
CN109752364A (zh) * 2019-03-11 2019-05-14 广西科技大学 多组分体系拉曼光谱的不同步褪色荧光消除方法
CN109752364B (zh) * 2019-03-11 2021-06-08 广西科技大学 多组分体系拉曼光谱的不同步褪色荧光消除方法
CN114184600A (zh) * 2021-12-30 2022-03-15 江苏海洋大学 一种基于拉曼光谱背景扣除的水中溶质定量方法

Also Published As

Publication number Publication date
CN108195817B (zh) 2020-09-15
WO2018103541A1 (zh) 2018-06-14

Similar Documents

Publication Publication Date Title
CN108195817A (zh) 用于去除溶剂干扰的拉曼光谱检测方法
US10408760B2 (en) Raman spectrum detecting method for eliminating package interference and electronic device thereof
CN108241845A (zh) 用于扣除谱图本底的方法和通过拉曼谱图识别物质的方法
JP5071475B2 (ja) 重回帰分析による予測モデルの作成方法、作成装置、作成プログラム
CN106198482B (zh) 基于拉曼光谱的检测保健品中是否添加有西药的方法
Abramo et al. Fourier analysis of multitracer cosmological surveys
Ruan et al. Quantitative detection of harmful elements in alloy steel by LIBS technique and sequential backward selection-random forest (SBS-RF)
EP3088872A1 (en) Raman spectroscopy method used for detecting sample in accommodating body
CN112129741B (zh) 绝缘油老化分析方法、装置、计算机设备和存储介质
SE516343C2 (sv) Metod och anordning för kalibrering av indata
Monakhova et al. Independent components analysis to increase efficiency of discriminant analysis methods (FDA and LDA): Application to NMR fingerprinting of wine
US11610073B2 (en) Device and method for analyzing the state of a system in a noisy context
CN108235733A (zh) 基于拉曼光谱的物质识别方法及云端系统
CN110991004A (zh) 一种轴承故障诊断方法、终端设备及计算机存储介质
CN106703796A (zh) 一种获取油藏的动态储量及水体大小的方法及装置
CN111768287A (zh) 周期识别方法、装置、服务器及可读存储介质
CN113009397A (zh) 预配训练数据使神经网络识别nmr测量中信号的系统和方法
CN108241846A (zh) 用于识别拉曼谱图的方法
CN111157115A (zh) 一种水下布里渊散射光谱获取方法及装置
CN110766402A (zh) 交易顺序依赖漏洞检测方法、系统、电子装置及存储介质
CN113610817B (zh) 一种特征峰识别方法及计算设备、存储介质
CN109246074A (zh) 识别可疑域名的方法、装置、服务器及可读存储介质
CN115541021A (zh) 拉曼光谱特征峰的定位方法、电子设备及存储介质
EP4105672A1 (en) Systems and methods for provisioning training data to enable neural networks to analyze signals in nmr measurements
EP4022268A1 (en) Spectroscopic apparatus and methods for determining components present in a sample

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant