CN106716109A - 用于光谱校正的水干扰估计 - Google Patents
用于光谱校正的水干扰估计 Download PDFInfo
- Publication number
- CN106716109A CN106716109A CN201580041877.8A CN201580041877A CN106716109A CN 106716109 A CN106716109 A CN 106716109A CN 201580041877 A CN201580041877 A CN 201580041877A CN 106716109 A CN106716109 A CN 106716109A
- Authority
- CN
- China
- Prior art keywords
- vector
- principal component
- subregion
- sample spectrum
- coefficient
- 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
Links
- 230000003595 spectral effect Effects 0.000 title claims description 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims description 27
- 238000012937 correction Methods 0.000 title claims description 23
- 239000013598 vector Substances 0.000 claims abstract description 126
- 238000001228 spectrum Methods 0.000 claims abstract description 95
- 239000011159 matrix material Substances 0.000 claims abstract description 41
- 238000000034 method Methods 0.000 claims abstract description 38
- 238000012549 training Methods 0.000 claims abstract description 17
- 239000000126 substance Substances 0.000 claims description 15
- 238000004458 analytical method Methods 0.000 claims description 5
- 230000004438 eyesight Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims 6
- 239000000203 mixture Substances 0.000 description 12
- 230000008569 process Effects 0.000 description 12
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 7
- 238000002835 absorbance Methods 0.000 description 4
- 238000002329 infrared spectrum Methods 0.000 description 4
- 238000000513 principal component analysis Methods 0.000 description 4
- 230000007613 environmental effect Effects 0.000 description 3
- 238000001157 Fourier transform infrared spectrum Methods 0.000 description 2
- 238000000862 absorption spectrum Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 239000012491 analyte Substances 0.000 description 2
- 239000013590 bulk material Substances 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 230000000386 athletic effect Effects 0.000 description 1
- -1 biotoxin Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000002575 chemical warfare agent Substances 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000001320 near-infrared absorption spectroscopy Methods 0.000 description 1
- 238000002610 neuroimaging Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012306 spectroscopic technique Methods 0.000 description 1
- 231100000765 toxin Toxicity 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3504—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/314—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry with comparison of measurements at specific and non-specific wavelengths
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3554—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for determining moisture content
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C99/00—Subject matter not provided for in other groups of this subclass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3504—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
- G01N2021/3545—Disposition for compensating effect of interfering gases
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N2021/3595—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/12—Circuits of general importance; Signal processing
- G01N2201/129—Using chemometrical methods
- G01N2201/1293—Using chemometrical methods resolving multicomponent spectra
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computing Systems (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Crystallography & Structural Chemistry (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Mathematical Physics (AREA)
Abstract
一种方法,该方法包括分解训练集获得含有多个主成分向量的主成分矩阵。该方法还包括通过递增地提供指示所选主成分向量对于所选子区域的权重的系数而可变地拒绝与所选的所述多个主成分向量之一不对应的部分样本光谱向量。校正的光谱向量可通过以下操作而被获取:排除子区域的所述样本光谱向量和相应的主成分向量;将所述样本光谱向量与针对非排除的子区域的主成分矩阵相乘;提供预测干扰向量;以及从所述样本光谱向量减去所述预测干扰向量。
Description
背景技术
一般来说,红外(IR)光谱(这是广义的定义,包括但不限于:近红外(NIR)、中红外(MIR)、和远红外(FIR))是基于分子振动及转动模式,以及组合和泛音频段。因为摩尔吸光系数在电磁频谱的近红外(NIR)区可以非常小,近红外辐射通常可以穿透很远的样本(例如,相比于中红外(MIR)辐射)。因此,红外光谱技术(如漫反射NIR光谱和傅立叶变换红外光谱(FTIR))对于探测散装材料(bulk material)而言是非常有用的。此外,NIR光谱通常需要有限的样本制备或不需要的样本制备。近红外光谱(NIRS)可用于材料表征和分子分析应用,诸如制药、医疗诊断、神经病学、神经影像学、新生儿研究、泌尿学、食品和化学品的质量控制、燃烧产物分析、运动医学、体育科学、体育训练、人体工程学、康复,等等。
发明内容
描述了一种使用例如傅立叶变换红外光谱(FTIR)系统进行样本检测的系统和技术。例如,一种方法包括分解训练集以获得具有多个主成分向量的主成分矩阵。该方法还包括通过递增地执行以下操作而可变地拒绝与所选的所述多个主成分向量之一不对应的部分样本光谱向量:选择子区域的所述样本光谱向量和相应子区域的所选主成分向量;以及将所选子区域的所述样本光谱向量与相应子区域的所选主成分向量相乘,以提供指示所选主成分向量对于所选子区域的权重的系数。该方法还包括基于所述递增地提供的系数,排除子区域的所述样本光谱向量和相应的主成分向量;将所述样本光谱向量与针对非排除的子区域的主成分矩阵相乘,以提供表示所述主成分矩阵的贡献的权重向量;将所述权重向量主乘以所述主成分矩阵,以提供预测干扰向量;以及从所述样本光谱向量减去所述预测干扰向量,以提供校正的光谱向量。
本发明内容还以一种简化的形式提供了选择理念,而该选理念将在在下面的详细描述进行进一步详细描述。本发明内容不旨在定义所要求保护的主题的关键特征或必要特征,也不打算作为帮助确定所要求保护的主题的范围。
附图说明
将参考附图进行详细描述。
图1是水蒸气的傅立叶变换红外光谱(FTIR)汽相谱图。
图2是含有水的样本混合物的光谱图以及针对一定量的水蒸气进行了校正的校正光谱图。
图3是根据本公开示例性实施的第一主成分的滑动窗口估计图。
图4是根据本公开示例性实施的含有水的样本混合物的光谱图以及针对一定量的水蒸气进行了校正的校正光谱图。
图5是用于估计FTIR光谱的水干扰和校正样本混合物中水存在的过程的流程图。
具体实施方式
包括诸如傅里叶变换红外光谱等(FTIR)的红外光谱的光学光谱技术可以在用于用来验证样本中分析物的存在的检测工具中实施,所述分析物包括但不限于危险或违禁品(如炸药、毒品、生物毒素、化学战剂等)。存在或不存在某些分析物(例如,通过光谱分析表示)可能会由于样本混合物中的干扰物质的存在而被掩盖。例如,水蒸汽可以是FTIR汽相光谱的干扰源。图1示出了水蒸气的FTIR汽相谱图。环境因素(如水蒸汽)可以存在于定量分析工具(例如,光谱仪)的气囊(gas cell)中,可以在样本的吸收光谱中呈现正面和负面的频带(bands)。这些因素可能会负面影响识别和量化库化合物(library compounds)的能力的。对这些环境因素进行过度补偿或较低补偿同样会提供有用的置信区间以外的数据。
描述了通过可变的选择来获得FTIR汽相样本光谱中水吸光度的估算方法。可使用数据训练集来获得的主成分矩阵,它可以用来描述与水干扰相关的光谱图案。在窄谱区域的各种范围内确定第一主成分,在感兴趣的光谱范围内提供第一主成分的多个值。假设第一主成分的多个值处于目标光谱范围上,则可在多个窄光谱区域上确定该第一主成分。提供显著区别于其他光谱区域的该第一主成分的值的光谱区域可指示一外围值(outliervalue),该外围值可表示除水之外的化学物质存在于那些光谱区域内。移除这些外围区域可提供对分配给每一主成分的系数或评分(score)的更为精确的估计,而这些系数或评分在被应用至整个光谱时,可用于估计水蒸气在整个光谱上的吸光度。
主成分分析(PCA)是用来识别数据集主要趋势的建模技术。当施加到样本(如水蒸汽)光谱分析的所提供的光谱数据时,该数据的因素/趋势可以被确定以应用到构成其他样本的存在/不存在的因素/趋势。在一个实施例中,可以使用在不同温度及浓度下测量的水蒸气光谱的训练集来创建主成分模型:
Xtrain=TPT (1)
其中Xtrain代表水汽吸收光谱训练集,T代表针对每一因素的系数或评分,且PT表示描述与水干扰相关的光谱图案的因素(例如,主成分载荷向量)。训练集可以包括在不同温度及浓度下测量的水蒸气光谱,以捕捉可在其中寻找到样本混合物的各种各样的环境条件。因此,训练集可以对于使用在此所描述的技术的检测工具所面临的特定场景而言是具有很强鲁棒性的。
将模型应用于样本数据集,可计算每个因素的系数或评分:
Xunknown*P=Tunknown (2)
其中Xunknown代表样本光谱数据集,P代表描述样本数据集的光谱图案的因素,并Tunknown代表所计算的每个因素的系数或评分。这些值可用于通过移除特定水平的水蒸气干扰而校正样本光谱:
Xcorrected=Xunknown-Tunknown*PT (3)
其中Xcorrected表示已通过移除特定水平的水蒸气干扰而被调整的校正后的样本光谱。图2提供了含有水的样本混合物的光谱图以及针对一定量的水蒸气进行了校正的校正光谱图。样本光谱包括存在于水蒸气区域内的其他化学峰值(即,非水峰值)(例如,如图1所示),校正后的数据集可以例如通过指示负峰值而提供对所存在的水量的过估计或误差。图2提供的图形显示了具有负峰的至少两个区域(如在200和202)。剩下的(正的)峰值也可能受到水蒸气量的过估计或误差的影响。
在本公开的实施中,可使用滑动窗口方法来识别除水蒸气之外的特征的光谱区域。这种方法可以用来避免过估计或不能准确提供所存在的水蒸气的量,特别是在样本包括存在于水蒸气区域内的化学峰值(即,非水峰值)的情况下。类似于在此所述的PCA方法,可在窄光谱区域内计算第一主成分(例如,主要为水蒸气的光谱)的系数(T)。该系数可通过使用例如水蒸气光谱训练集(例如,在多种不同浓度及温度下采集的多个光谱)而被导出为评分矩阵。在实施中,可以以10cm-1的步长使用45cm-1(波数)的窗口,从而覆盖约1300cm-1至约2000cm-1的范围,并形成可对应于水蒸气光谱的主峰值区域的约3400cm-1至约4000cm-1的范围(如图1所示)。参见图3,示出了根据本公开示例性实施的第一主成分的滑动窗口估计图,其中包含了针对主成分评分的八十(80)个估计。
评分的阈值可以被确定,其在超过该阈值的评分可以指示包含除水蒸气之外的化学特征的光谱区域。在实施中,评分的中间值(median value)可用作阈值,其中包含显著不同于所述中间值的评分的光谱区域可被标记或归档作为可能包含来自存在于样本混合物内除水蒸气之外的化学物质的特征。
如图3中的示例实施中所示,示出了具有显著不同于评分中间值的值的两个区域的光谱区域(标记为300和302)。这些区域对应于与水以外的化学物质相关的峰值。区域300位于约1300cm-1至约1400cm-1的(对应于窗口1-5)的光谱范围内。区域302位于约1700cm-1至约17750cm-1的(对应于窗口27-30)的光谱范围内。这些区域对应于图2中感兴趣的第一光谱区域内的两个最强峰值(标记为204和206)。由于这些区域对应于非水化学物质,它们可以被从模型因素移除,以提供样本中由于水的贡献的估计。在实施例中,包含显著不同于阈值的评分的区域被从样本混合物光谱及原始模型因素(P)移除,以获取更新的评分估计(例如,更为精确的),类似的方程(2)所述。之后可将更新的评分估计用于(整个光谱范围)原始因素,以估计整个光谱上由于水蒸气所导致的吸光度。然后,从样本混合物光谱中减去由于水蒸汽所导致的吸光度,从而提供无来自水蒸气的实质性干扰的样本光谱。例如,图4显示了样本混合物的光谱图(标记为400),由于水蒸气的存在,该光谱有干扰。图4还包括校正光谱(标记为402),其已根据本公开一示例性实施根据滑动窗口法修改PCA模型而对水汽量进行了校正。在图4中可以看到,相比于根据PCA方法在图2中的区域200及202所发现的对水蒸气量的过估计/误差,校正后的光谱402提供了相同水平的对水蒸气量的过估计/误差。
现参考图5,其示出了根据本公开一示例性实施的用于估计FTIR光谱的水干扰及校正样本混合物内的水存在的过程500的流程图。本文所述的过程的各个方面可以被实施为硬件、软件或其组合。该过程被示为一组方框,该组方框指定了由一个或多个设备(例如,IR光谱仪器仪表、计算机系统控制的IR光谱仪器或红外光谱组件)执行的操作且并不一定限于所显示的执行各个方框的操作的顺序。在实施中,本文所描述的过程的各个方面可通过固定或移动化学检测系统(如汽相检测系统)来实施。
过程500包括分解训练集以获得具有多个主成分向量的主成分矩阵(方框502)。例如,将水蒸气光谱的训练集分解以获得具有多个主成分向量的主成分矩阵。在实施中,训练集包括在各种不同的温度和浓度下测量的水蒸气光谱。
过程500还包括可变地拒绝与多个主成分向量的所选一者(PT)不对应的部分样本光谱向量(X)。方框504包括以下递增的过程:选择子区域的样本光谱向量以及对应区域的所选主成分向量(方框506),以及将所选子区域的该样本光谱向量与对应子区域的所选主成分向量相乘以提供指示所选主成分向量对应所选子区域的权重的系数(方框508)。例如,图3中示出了指示该所选主成分向量的权重的系数的实施,其中示出了主成分向量的滑动窗口估计,其中包含了针对主成分评分的八十(80)个估计。
过程500还包括根据递增提供的系数,排除子区域的样本光谱向量和相应的主成分向量(方框510)。例如,可以确定评分的阈值,其中超过阈值的评分可以指示包括除水蒸气以外的化学特征的光谱区域。在实施中,指示所选主成分向量的权重的系数的中间值可用作所述阈值,其中包含显著不同于所述中间值的评分的光谱子区域可从样本光谱向量及对应的主成分向量中被排除。
过程500进一步包括将样本光谱向量与对于非排除子区域的主成分矩阵相乘,以提供表示主成分矩阵贡献的权重向量(方框512)。例如,在实施中,权重向量(T)被提供如下:
T=X*P
其中X表示样本光谱向量,其中P表示非排除子区域的主成分矩阵。在实施中,T由以下提供:
T=X*P*(PT*P)-1
其中X表示样本光谱向量,P表示非排除子区域的主成分矩阵,PT表示主成分矩阵。
步骤500进一步包括将主成分矩阵乘以权重向量以提供预测干扰向量(方框514)。例如,在实施中,预测干扰向量(X’)可被提供如下:
X′=T*PT
其中T表示权重向量,PT表示主成分矩阵。
步骤500进一步包括从样本光谱向量减去预测干扰向量以提供校正的光谱向量(方框516)。例如,在实施中,根据以下内容提供校正的频谱向量(E):
E=X-X′
其中X表示样本光谱向量,X’表示预测干扰向量。
在实施中,过程500包括分析校正的光谱向量。例如,校正后的光谱可以与感兴趣的各种化学物质的已知光谱库进行比较,以确定校正后的光谱是否与某些感兴趣的化学物质相关联。在分析之后,当校正后的光谱代表感兴趣的化学物质时,过程500可以包括激活一个或多个视觉或听觉标记。例如,当确定校正的频谱对应或与感兴趣的化学物质相关联时,化学检测系统可以激活声音或可见的警报,向用户界面产生可视输出等等。
需要注意的是,对于计算机领域的技术人员而言很显然的是,上述描述的实施例和实施可以使用传统的通用数字计算机按照本说明书的教示来方便地实施。对于软件领域的技术人员而言很显然的是,适当的软件编码可以很容易地由程序员基于本公开的教示而被编写。
可以理解,本文所描述的实施例和实施可以以软件包的形式方便地实施。该软件包可以是具有计算机可执行指令或计算机程序产品的计算机可读存储介质,其采用非短暂的计算机可读存储介质,该存储介质包含存储的计算机代码,该计算机代码可使得计算机执行所公开的功能及在此所公开的过程。计算机可读存储介质可能包括但不限于任何类型的传统软盘、光盘、CD-ROM、磁盘、硬盘、磁光盘、ROM、RAM、EPROM、EEPROM、磁性或光学卡、或任何其他合适的用于存储电子指令的媒介。
此处所使用的术语“约”指的是关于所指定的值或范围的大概值和/或精确值。虽然已通过特定的结构特征和/或方法动作的语言对主题进了描述,但应该理解的是,所附权利要求书中所限定的主题并不限于所述的特定特征或动作。虽然描述了各种不同的配置,但是可在不背离本公开的情况下以各种方式来构建设备、系统、系统组件等。相反,所公开的特定特征及动作均仅作为实施权利要求的示例性形式。
Claims (20)
1.一种计算机实施的方法,包括:
分解训练集以获得具有多个主成分向量的主成分矩阵;
通过递增地执行以下操作而可变地拒绝与所选的所述多个主成分向量之一不对应的部分样本光谱向量:
选择子区域的所述样本光谱向量和相应子区域的所选主成分向量;以及
将所选子区域的所述样本光谱向量与相应子区域的所选主成分向量相乘,以提供指示所选主成分向量对于所选子区域的权重的系数;
基于所述递增地提供的系数,排除子区域的所述样本光谱向量和相应的主成分向量;
将所述样本光谱向量与针对非排除的子区域的主成分矩阵相乘,以提供表示所述主成分矩阵的贡献的权重向量;
将所述权重向量主乘以所述主成分矩阵,以提供预测干扰向量;以及
从所述样本光谱向量减去所述预测干扰向量,以提供校正的光谱向量。
2.根据权利要求1所述的计算机实施的方法,进一步地包括测量水蒸气光谱获得所述训练集。
3.根据权利要求2所述的计算机实施的方法,其中测量水蒸气光谱包括测量水在不同温度或不同压力的至少一个下的光谱。
4.根据权利要求1所述的计算机实施的方法,进一步地包括外围化多数所述递增提供的系数。
5.根据权利要求1所述的计算机实施的方法,其中所述权重向量等于所述样本光谱向量乘以所述非排除子区域的主成分矩阵乘以以下的逆矩阵:所述主成分矩阵乘以所述非排除子区域的主成分矩阵。
6.根据权利要求1所述的计算机实施的方法,其中基于递增提供的系数来排除子区域的样本光谱向量和相应的主成分向量包括:排除大于所述递增提供的系数的阈值的子区域的样本光谱向量和相应的主成分向量。
7.根据权利要求6所述的计算机实施的方法,其中所述递增提供的系数的阈值是所述递增提供的系数的中间值。
8.一种方法,包括:
分解训练集以获得具有多个主成分向量的主成分矩阵;
通过递增地执行以下操作而可变地拒绝与所选的所述多个主成分向量之一不对应的部分样本光谱向量:
选择子区域的所述样本光谱向量和相应子区域的所选主成分向量;以及
将所选子区域的所述样本光谱向量与相应子区域的所选主成分向量相乘,以提供指示所选主成分向量对于所选子区域的权重的系数;
基于所述递增地提供的系数,排除子区域的所述样本光谱向量和相应的主成分向量;
将所述样本光谱向量与针对非排除的子区域的主成分矩阵相乘,以提供表示所述主成分矩阵的贡献的权重向量;
将所述权重向量主乘以所述主成分矩阵,以提供预测干扰向量;
从所述样本光谱向量减去所述预测干扰向量,以提供校正的光谱向量;
分析校正的光谱向量;以及
当所述校正的光谱向量表示感兴趣化学物质时,激活一个或多个视觉或听觉标记。
9.根据权利要求8所述的方法,进一步包括测量水蒸气光谱以获得所述训练集。
10.根据权利要求9所述的方法,其中测量水蒸气光谱包括测量水在不同温度或不同压力的至少一个下的光谱。
11.根据权利要求8所述的方法,进一步地包括外围化多数所述递增提供的系数。
12.根据权利要求8所述的方法,其中所述权重向量等于所述样本光谱向量乘以所述非排除子区域的主成分矩阵乘以以下的逆矩阵:所述主成分矩阵乘以所述非排除子区域的主成分矩阵。
13.根据权利要求8所述的方法,其中基于递增提供的系数来排除子区域的样本光谱向量和相应的主成分向量包括:排除大于所述递增提供的系数的阈值的子区域的样本光谱向量和相应的主成分向量。
14.根据权利要求13所述的方法,其中所述递增提供的系数的阈值是所述递增提供的系数的中间值。
15.具有计算机可执行指令的计算机可读存储介质,该计算机可执行指令用于:
分解训练集以获得具有多个主成分向量的主成分矩阵;
通过递增地执行以下操作而可变地拒绝与所选的所述多个主成分向量之一不对应的部分样本光谱向量:
选择子区域的所述样本光谱向量和相应子区域的所选主成分向量;以及
将所选子区域的所述样本光谱向量与相应子区域的所选主成分向量相乘,以提供指示所选主成分向量对于所选子区域的权重的系数;
基于所述递增地提供的系数,排除子区域的所述样本光谱向量和相应的主成分向量;
将所述样本光谱向量与针对非排除的子区域的主成分矩阵相乘,以提供表示所述主成分矩阵的贡献的权重向量;
将所述权重向量主乘以所述主成分矩阵,以提供预测干扰向量;以及
从所述样本光谱向量减去所述预测干扰向量,以提供校正的光谱向量。
16.根据权利要求15所述的计算机可读存储介质,进一步地包括用于测量水蒸气光谱获得所述训练集的计算机可执行指令。
17.根据权利要求16所述的计算机可读存储介质,其中测量水蒸气光谱包括测量水在不同温度或不同压力的至少一个下的光谱。
18.根据权利要求15所述的计算机可读存储介质,进一步地包括用于外围化多数所述递增提供的系数的计算机可执行指令。
19.根据权利要求15所述的计算机可读存储介质,其中所述权重向量等于所述样本光谱向量乘以所述非排除子区域的主成分矩阵乘以以下的逆矩阵:所述主成分矩阵乘以所述非排除子区域的主成分矩阵。
20.根据权利要求15所述的计算机可读存储介质,其中基于递增提供的系数来排除子区域的样本光谱向量和相应的主成分向量包括:排除大于所述递增提供的系数的中间值的子区域的样本光谱向量和相应的主成分向量。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/446,473 | 2014-07-30 | ||
US14/446,473 US10101269B2 (en) | 2014-07-30 | 2014-07-30 | Estimation of water interference for spectral correction |
PCT/US2015/042928 WO2016019153A1 (en) | 2014-07-30 | 2015-07-30 | Estimation of water interference for spectral correction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106716109A true CN106716109A (zh) | 2017-05-24 |
CN106716109B CN106716109B (zh) | 2019-11-05 |
Family
ID=55179736
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580041877.8A Active CN106716109B (zh) | 2014-07-30 | 2015-07-30 | 用于光谱校正的水干扰估计 |
Country Status (6)
Country | Link |
---|---|
US (2) | US10101269B2 (zh) |
EP (1) | EP3175223B1 (zh) |
CN (1) | CN106716109B (zh) |
CA (1) | CA2956227C (zh) |
PL (1) | PL3175223T3 (zh) |
WO (1) | WO2016019153A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109444072A (zh) * | 2018-10-12 | 2019-03-08 | 南京富岛信息工程有限公司 | 一种近红外光谱模型转移中诡峰的解决方法 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10522252B2 (en) | 2017-06-16 | 2019-12-31 | Carefusion 303, Inc. | Opioid management system |
US11081220B2 (en) | 2018-02-02 | 2021-08-03 | Carefusion 303, Inc. | System and method for dispensing medication |
CN112292733A (zh) | 2018-05-04 | 2021-01-29 | 康尔福盛303公司 | 基于对等团体的异常行为检测 |
WO2020146407A1 (en) | 2019-01-07 | 2020-07-16 | Carefusion 303, Inc. | Machine learning based safety controller |
US11984212B2 (en) | 2019-01-10 | 2024-05-14 | Carefusion 303, Inc. | System for monitoring dose pattern and patient response |
CA3126700A1 (en) | 2019-01-18 | 2020-07-23 | Carefusion 303, Inc. | Medication tracking system |
CN110208248B (zh) * | 2019-06-28 | 2021-11-19 | 南京林业大学 | 一种辨识拉曼光谱异常测量信号的方法 |
CN114166777B (zh) * | 2021-11-30 | 2023-09-01 | 海南电网有限责任公司电力科学研究院 | 一种用于绝缘纸的近红外光谱影响因素校正方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060197957A1 (en) * | 2005-03-07 | 2006-09-07 | Jones Christopher M | Method to reduce background noise in a spectrum |
CN101688832A (zh) * | 2007-07-03 | 2010-03-31 | 皇家飞利浦电子股份有限公司 | 散射组织中物质浓度的光谱法测量 |
CN102313714A (zh) * | 2011-09-21 | 2012-01-11 | 浙江大学 | 一种红花提取液测定方法 |
CN102768195A (zh) * | 2012-06-29 | 2012-11-07 | 杭州中美华东制药有限公司 | 一种虫草菌粉水分含量快速检测方法 |
CN103063605A (zh) * | 2013-01-04 | 2013-04-24 | 云南植物药业有限公司 | 用傅立叶变换近红外光谱仪快速测定三七提取物及其制剂中五种皂苷含量的方法 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0015068B2 (en) * | 1979-02-22 | 1987-03-18 | Beckman Instruments, Inc. | Non-dispersive infrared analyzers |
EP0982582B1 (en) | 1998-08-28 | 2005-06-01 | Perkin-Elmer Limited | Suppression of undesired components in measured spectra |
US6697654B2 (en) | 1999-07-22 | 2004-02-24 | Sensys Medical, Inc. | Targeted interference subtraction applied to near-infrared measurement of analytes |
US7756558B2 (en) * | 2004-05-24 | 2010-07-13 | Trutouch Technologies, Inc. | Apparatus and methods for mitigating the effects of foreign interferents on analyte measurements in spectroscopy |
US6687620B1 (en) * | 2001-08-01 | 2004-02-03 | Sandia Corporation | Augmented classical least squares multivariate spectral analysis |
US6662116B2 (en) * | 2001-11-30 | 2003-12-09 | Exxonmobile Research And Engineering Company | Method for analyzing an unknown material as a blend of known materials calculated so as to match certain analytical data and predicting properties of the unknown based on the calculated blend |
JP2005534012A (ja) | 2002-07-19 | 2005-11-10 | テクニツク・インコーポレーテツド | 工業電解液をリアルタイムモニターする方法と装置 |
US20050285023A1 (en) * | 2004-06-23 | 2005-12-29 | Lambda Solutions, Inc. | Automatic background removal for input data |
US20080112853A1 (en) | 2006-08-15 | 2008-05-15 | Hall W Dale | Method and apparatus for analyte measurements in the presence of interferents |
WO2011051166A1 (en) | 2009-11-02 | 2011-05-05 | Universität Für Bodenkultur Wien | Method for determining polycyclic aromatic hydrocarbon contaminant concentration |
WO2011091280A2 (en) | 2010-01-22 | 2011-07-28 | University Of Massachusetts | Methods and systems for analyte measurement |
US8645082B2 (en) | 2010-09-13 | 2014-02-04 | Mks Instruments, Inc. | Monitoring, detecting and quantifying chemical compounds in a sample |
US8953724B2 (en) * | 2012-06-27 | 2015-02-10 | Andrew Llc | Canceling narrowband interfering signals in a distributed antenna system |
-
2014
- 2014-07-30 US US14/446,473 patent/US10101269B2/en active Active
-
2015
- 2015-07-30 EP EP15827708.7A patent/EP3175223B1/en active Active
- 2015-07-30 CA CA2956227A patent/CA2956227C/en active Active
- 2015-07-30 CN CN201580041877.8A patent/CN106716109B/zh active Active
- 2015-07-30 PL PL15827708T patent/PL3175223T3/pl unknown
- 2015-07-30 WO PCT/US2015/042928 patent/WO2016019153A1/en active Application Filing
-
2018
- 2018-10-15 US US16/160,505 patent/US10942117B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060197957A1 (en) * | 2005-03-07 | 2006-09-07 | Jones Christopher M | Method to reduce background noise in a spectrum |
CN101688832A (zh) * | 2007-07-03 | 2010-03-31 | 皇家飞利浦电子股份有限公司 | 散射组织中物质浓度的光谱法测量 |
CN102313714A (zh) * | 2011-09-21 | 2012-01-11 | 浙江大学 | 一种红花提取液测定方法 |
CN102768195A (zh) * | 2012-06-29 | 2012-11-07 | 杭州中美华东制药有限公司 | 一种虫草菌粉水分含量快速检测方法 |
CN103063605A (zh) * | 2013-01-04 | 2013-04-24 | 云南植物药业有限公司 | 用傅立叶变换近红外光谱仪快速测定三七提取物及其制剂中五种皂苷含量的方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109444072A (zh) * | 2018-10-12 | 2019-03-08 | 南京富岛信息工程有限公司 | 一种近红外光谱模型转移中诡峰的解决方法 |
CN109444072B (zh) * | 2018-10-12 | 2021-02-09 | 南京富岛信息工程有限公司 | 一种近红外光谱模型转移中诡峰的解决方法 |
Also Published As
Publication number | Publication date |
---|---|
PL3175223T3 (pl) | 2021-05-31 |
WO2016019153A1 (en) | 2016-02-04 |
EP3175223B1 (en) | 2020-11-04 |
EP3175223A4 (en) | 2018-04-04 |
US10101269B2 (en) | 2018-10-16 |
CA2956227A1 (en) | 2016-02-04 |
US10942117B2 (en) | 2021-03-09 |
US20190113447A1 (en) | 2019-04-18 |
EP3175223A1 (en) | 2017-06-07 |
CA2956227C (en) | 2024-01-02 |
US20160033394A1 (en) | 2016-02-04 |
CN106716109B (zh) | 2019-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106716109B (zh) | 用于光谱校正的水干扰估计 | |
Galvao et al. | A variable elimination method to improve the parsimony of MLR models using the successive projections algorithm | |
Li et al. | Data fusion strategy in quantitative analysis of spectroscopy relevant to olive oil adulteration | |
He et al. | Baseline correction for Raman spectra using an improved asymmetric least squares method | |
Clairotte et al. | National calibration of soil organic carbon concentration using diffuse infrared reflectance spectroscopy | |
de Araújo Gomes et al. | The successive projections algorithm for interval selection in PLS | |
Han et al. | An ensemble of Monte Carlo uninformative variable elimination for wavelength selection | |
Yun et al. | An efficient method of wavelength interval selection based on random frog for multivariate spectral calibration | |
Duan et al. | Determination of 27 chemical constituents in Chinese southwest tobacco by FT-NIR spectroscopy | |
Jing et al. | Multiblock partial least squares regression based on wavelet transform for quantitative analysis of near infrared spectra | |
CN103558167B (zh) | 一种快速测定腌肉中氯化钠含量的方法 | |
US10345242B2 (en) | Reverse intensity correction for Raman spectral library search | |
Chen et al. | FT-NIR spectroscopy and Whittaker smoother applied to joint analysis of duel-components for corn | |
US10557792B2 (en) | Spectral modeling for complex absorption spectrum interpretation | |
Smith et al. | PRFFECT: A versatile tool for spectroscopists | |
Bertacchini et al. | The impact of chemometrics on food traceability | |
Pan et al. | Optimal partner wavelength combination method with application to near-infrared spectroscopic analysis | |
Nespeca et al. | Multivariate filters combined with interval partial least square method: A strategy for optimizing PLS models developed with near infrared data of multicomponent solutions | |
CN102854151B (zh) | 一种光谱分析中样品集划分的化学计量学方法 | |
Ramírez-Rincón et al. | Separation of optical properties for multicomponent samples and determination of spectral similarity indices based on FEDS0 algorithm | |
Xu et al. | Chemometric methods for evaluation of chromatographic separation quality from two-way data—A review | |
Chen et al. | Investigation of sample partitioning in quantitative near-infrared analysis of soil organic carbon based on parametric LS-SVR modeling | |
Liu et al. | Simultaneous quantitative analysis of three components in mixture samples based on NIR spectra with temperature effect | |
CN115372309A (zh) | 一种关联有效特征谱段选择方法及油品指标含量快检方法 | |
Li et al. | An improved ensemble model for the quantitative analysis of infrared spectra |
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 |