CN105784637B - 标识光谱差异性的方法 - Google Patents
标识光谱差异性的方法 Download PDFInfo
- Publication number
- CN105784637B CN105784637B CN201610192572.3A CN201610192572A CN105784637B CN 105784637 B CN105784637 B CN 105784637B CN 201610192572 A CN201610192572 A CN 201610192572A CN 105784637 B CN105784637 B CN 105784637B
- Authority
- CN
- China
- Prior art keywords
- moving window
- scale
- spectrogram
- spectrum
- window
- 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.)
- Active
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000003595 spectral effect Effects 0.000 claims abstract description 55
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 238000002329 infrared spectrum Methods 0.000 abstract description 4
- 239000000126 substance Substances 0.000 description 2
- 238000001069 Raman spectroscopy Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000004611 spectroscopical analysis Methods 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/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
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (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)
- Spectrometry And Color Measurement (AREA)
Abstract
一种标识光谱差异性的方法,该方法为将谱图定义多个不同尺度的移动窗口;在每个尺度下,计算窗口内互相比对的光谱的向量夹角值,通过移动窗口,得到整个光谱范围的系列夹角值;再调整窗口大小,从而得到各个尺度下的完整谱线差异;然后,考察谱线不同位置和尺度的差异分布,由差异分布判断谱线的一致性。本发明通过判断测量光谱与标准光谱的一致性,可以实现近红外光谱的定性判断,所用设备简单,成本低,而且增强了辨识的合理性,方法简便直观,便于理解。
Description
技术领域
本发明涉及光谱学及仪器分析领域,特别是一种标识光谱差异性的方法。
背景技术
光谱是一种表达物质特性的手段,通常中红外、拉曼等包含的清晰特征峰会被用来进行物质定性。但是近红外光谱通常缺乏清晰的特征峰,影响到定性的准确性。目前,近红外定性多采用模式识别方法,通过较为笼统的相关系数、夹角余弦值等相似系数判别。
由于近红外光谱特征不显著,仪器响应差异也会引起既有相似系数缺乏特异性。
以图中所示谱线为例,图1是样本1在低光度值的谱线,图2是样本1叠加了一个其他响应峰(即样本2)的谱线,图3是仪器在大响应值下样本1的谱线,图4是将三张谱线放在一起进行对比,为了对比将图3的响应值按比例做了缩小,在该图4中,y1所示的谱线为样本1的谱线,y2所示的谱线为样本1的谱线,y3所示的谱线为样本1大响应值按比例缩小的谱线。由图中可以看出三条谱线不一致,图2与图1的差异由不同样本导致,而图3与图1是同一样本,由于仪器的非线性响应导致。
图1与图2的相关系数为0.9991,图1与图3的相关系数为0.9981,同一台光谱仪对同一样本的响应差异大于不同样本的响应差异。如果以相关系数度量光谱一致性,根据相关系数判断,将得到不合理的结论。
发明内容
本发明要解决的技术问题是:提供一种标识光谱差异性的方法,以实现近红外光谱的定性判断。
解决上述技术问题的技术方案是:一种标识光谱差异性的方法,该方法为将谱图定义多个不同尺度的移动窗口;在每个尺度下,计算窗口内互相比对的光谱的向量夹角值,通过移动窗口,得到整个光谱范围的系列夹角值;再调整窗口大小,从而得到各个尺度下的完整谱线差异;然后,考察谱线不同位置和尺度的差异分布,由差异分布判断谱线的一致性。
本发明的进一步技术方案是:该方法包括以下步骤:
S1. 选择被比较的两个谱图;
S2.根据谱图所包含的数据点数,分别以总数据点数的1/2、1/4、1/8、……、1/2n分割,最小尺度的移动窗口内包含8至10个数据点,定义不同尺度下的移动窗口大小;
S3.从最小尺度的移动窗口开始,从谱图左端移动,计算移动窗口内两条谱图的向量夹角;
S4.将移动窗口移动至下一点,计算向量夹角,直至移动窗口移至谱图右端;
S5. 将移动窗口移动计算得到的向量夹角保存至数据序列Stemp;
S6.选择下一尺度的移动窗口,从谱图左端移动,计算移动窗口内两条谱图的向量夹角,将移动窗口移动至下一点,计算向量夹角,直至移动窗口移至谱图右端,将移动窗口移动计算得到的向量夹角保存至数据序列Stemp,重复步骤S6的以上操作,直至所有大小的移动窗口都完成计算,得到一个完整光谱在不同尺度下的向量夹角数据序列Sall;
S7.计算数据序列Sall的方差值D;
S8.将方差值D除以π/2,得到光谱差异系数τ ;
S9.根据光谱差异系数τ判断两个谱图的相似度:光谱差异系数τ的值介于0和1之间,它越趋近于0,代表被比较的两个谱图越相似;光谱差异系数τ趋近于1,表明两个谱图的差异越大。
本发明的再进一步技术方案是:在步骤S1中,选择的谱图要求谱图的横坐标一致,也就是对应的数据点处于相同的波长。
由于采用上述结构,本发明之标识光谱差异性的方法与现有技术相比,具有以下有益效果:
1.可实现近红外光谱的定性判断:
由于本发明是将谱图定义多个不同尺度的移动窗口;在每个尺度下,计算窗口内互相比对的光谱的向量夹角值,通过移动窗口,得到整个光谱范围的系列夹角值;再调整窗口大小,从而得到各个尺度下的完整谱线差异;然后,考察谱线不同位置和尺度的差异分布,由差异分布判断谱线的一致性。因此,本发明通过判断测量光谱与标准光谱的一致性,可以实现近红外光谱的定性判断。
2.所用设备简单:
本发明无需复杂的设备,成本较低。
3. 增强了辨识的合理性:
从图5-图7的角度序列图可以看出,本发明加入了细节辨识,将原来采用粗略的相关系数整体辨识转化增加了细节辨识,有效增强了辨识的合理性。
4. 方法简便直观,便于理解:
由于向量间的夹角越大,相似性越低,其最大值是π/2,最小值为0。本发明将角度序列方差值归一至光谱差异系数τ,根据光谱差异系数τ判断两个谱图的相似度:光谱差异系数τ的值介于0和1之间,它越趋近于0,代表被比较的两个谱图越相似;光谱差异系数τ趋近于1,表明两个谱图的差异越大。因此,本方法比较简便直观,也便于理解。
下面,结合附图和实施例对本发明之标识光谱差异性的方法的技术特征作进一步的说明。
附图说明
图1:样本1在低光度值的谱线,
图2:样本1叠加了一个其他响应峰(即样本2)的谱线,
图3:仪器在大响应值下样本1的谱线,
图4:三条谱线的对比图,
图5:图1和图2的移动窗口夹角序列值,
图6:图1和图3的移动窗口夹角序列值,
图7:图2和图3的移动窗口夹角序列值。
具体实施方式
一种标识光谱差异性的方法,该方法为将谱图定义多个不同尺度的移动窗口;在每个尺度下,计算窗口内互相比对的光谱的向量夹角值,通过移动窗口,得到整个光谱范围的系列夹角值;再调整窗口大小,从而得到各个尺度下的完整谱线差异;然后,考察谱线不同位置和尺度的差异分布,由差异分布判断谱线的一致性。
该方法包括以下步骤:
S1. 选择被比较的两个谱图,选择的谱图要求谱图的横坐标一致,也就是对应的数据点处于相同的波长;
S2.根据谱图所包含的数据点数,分别以总数据点数的1/2、1/4、1/8、……、1/2n分割,最小尺度的移动窗口内包含8至10个数据点,定义不同尺度下的移动窗口大小;
S3.从最小尺度的移动窗口开始,从谱图左端移动,计算移动窗口内两条谱图的向量夹角;
S4.将移动窗口移动至下一点,计算向量夹角,直至移动窗口移至谱图右端;
S5. 将移动窗口移动计算得到的向量夹角保存至数据序列Stemp;
S6.选择下一尺度的移动窗口,从谱图左端移动,计算移动窗口内两条谱图的向量夹角,将移动窗口移动至下一点,计算向量夹角,直至移动窗口移至谱图右端,将移动窗口移动计算得到的向量夹角保存至数据序列Stemp,重复步骤S6的以上操作,直至所有大小的移动窗口都完成计算,得到一个完整光谱在不同尺度下的向量夹角数据序列Sall;
S7.计算数据序列Sall的方差值D;
S8.将方差值D除以π/2,得到光谱差异系数τ ;
S9.根据光谱差异系数τ判断两个谱图的相似度:光谱差异系数τ的值介于0和1之间,它越趋近于0,代表被比较的两个谱图越相似;光谱差异系数τ趋近于1,表明两个谱图的差异越大。
以下是具体的实施案例:
实施例一:
S1.先择图1、图2、图3所示的三张谱图,分别比较图1与图2,图1与图3,图2与图3的差异;
S2. 图1、图2、图3中分别含有1026个数据点,分别取8、16、32、64、128、256、512个数据点的移动窗口;
S3. 自谱图的左端,计算8个数据点移动窗口内的两条谱线的向量夹角;
S4.移动窗口移至下一点,计算向量夹角,直至移动窗口到达谱图的右端;
S5. 将移动窗口移动计算得到的向量夹角保存至数据序列Stemp;
S6.增大数据窗口,从谱图左端移动,计算移动窗口内两条谱图的向量夹角,将移动窗口移动至下一点,计算向量夹角,直至移动窗口移至谱图右端,将移动窗口移动计算得到的向量夹角保存至数据序列Stemp,重复步骤S6的以上操作;直至所有大小的数据窗口计算完成;将Stemp输出为Sall,得到一个有关向量夹角的数据序列Sall;
S7.计算数据序列Sall的方差值D;
S8.将方差值D除以π/2,得到光谱差异系数分别为:τ13=0.0177 ,τ12=0.1257,τ32=0.1227;
S9. 根据光谱差异系数τ判断相似度:
由光谱差异系数τ13=0.0177 ,τ12=0.1257,τ32=0.1227可看出,同一样本谱图1和图3之间的差异,显著小于与不同样本的谱图2的差异,推断出图1与图3为同一样本光谱。
Claims (2)
1.一种标识光谱差异性的方法,其特征在于:该方法为将谱图定义多个不同尺度的移动窗口;在每个尺度下,计算窗口内互相比对的光谱的向量夹角值,通过移动窗口,得到整个光谱范围的系列夹角值;再调整窗口大小,从而得到各个尺度下的完整谱线差异;然后,考察谱线不同位置和尺度的差异分布,由差异分布判断谱线的一致性;该方法包括以下步骤:
S1. 选择被比较的两个谱图;
S2.根据谱图所包含的数据点数,分别以总数据点数的1/2、1/4、1/8、……、1/2n分割,最小尺度的移动窗口内包含8至10个数据点,定义不同尺度下的移动窗口大小;
S3.从最小尺度的移动窗口开始,从谱图左端移动,计算移动窗口内两条谱图的向量夹角;
S4.将移动窗口移动至下一点,计算向量夹角,直至移动窗口移至谱图右端;
S5. 将移动窗口移动计算得到的向量夹角保存至数据序列Stemp;
S6. 选择下一尺度的移动窗口,从谱图左端移动,计算移动窗口内两条谱图的向量夹角,将移动窗口移动至下一点,计算向量夹角,直至移动窗口移至谱图右端,将移动窗口移动计算得到的向量夹角保存至数据序列Stemp,重复步骤S6的以上操作,直至所有大小的移动窗口都完成计算,得到一个完整光谱在不同尺度下的向量夹角数据序列Sall;
S7.计算数据序列Sall的方差值D;
S8.将方差值D除以π/2,得到光谱差异系数τ ;
S9.根据光谱差异系数τ判断两个谱图的相似度:光谱差异系数τ的值介于0和1之间,它越趋近于0,代表被比较的两个谱图越相似;光谱差异系数τ趋近于1,表明两个谱图的差异越大。
2.根据权利要求1所述的标识光谱差异性的方法,其特征在于:在步骤S1中,选择的谱图要求谱图的横坐标一致,也就是对应的数据点处于相同的波长。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610192572.3A CN105784637B (zh) | 2016-03-30 | 2016-03-30 | 标识光谱差异性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610192572.3A CN105784637B (zh) | 2016-03-30 | 2016-03-30 | 标识光谱差异性的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105784637A CN105784637A (zh) | 2016-07-20 |
CN105784637B true CN105784637B (zh) | 2018-08-17 |
Family
ID=56392363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610192572.3A Active CN105784637B (zh) | 2016-03-30 | 2016-03-30 | 标识光谱差异性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105784637B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106323898B (zh) * | 2016-08-23 | 2019-01-15 | 广西科技大学 | 混合体系中添加物直接光谱定量的背景库扩充方法 |
CN106383091B (zh) * | 2016-08-23 | 2019-03-01 | 广西科技大学 | 通过光谱直接定量添加物含量的方法 |
CN106802283B (zh) * | 2016-12-31 | 2018-11-30 | 华中科技大学 | 一种获取谱差异的方法及系统 |
CN109975232B (zh) * | 2017-12-28 | 2023-08-01 | 交通运输部科学研究院 | 一种沥青及沥青改性添加剂的检测方法 |
CN110208666B (zh) * | 2019-07-03 | 2021-07-16 | 云南电网有限责任公司电力科学研究院 | 局部放电特征光谱的选取方法 |
CN111426648B (zh) * | 2020-03-19 | 2023-04-07 | 甘肃省交通规划勘察设计院股份有限公司 | 一种红外谱图的相似度确定方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102880861A (zh) * | 2012-09-05 | 2013-01-16 | 西安电子科技大学 | 基于线性预测倒谱系数的高光谱图像分类方法 |
CN104502302A (zh) * | 2014-11-04 | 2015-04-08 | 中国航天科工集团第三研究院第八三五七研究所 | 混油的太赫兹时域波形多参数组合定量分析方法 |
CN104777143A (zh) * | 2014-01-15 | 2015-07-15 | 中国人民解放军第二军医大学 | 一种用于过期药的、基于拉曼光谱的相似度鉴定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9274052B2 (en) * | 2013-07-10 | 2016-03-01 | Canon Kabushiki Kaisha | Feature vector for classifying specular objects based on material type |
-
2016
- 2016-03-30 CN CN201610192572.3A patent/CN105784637B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102880861A (zh) * | 2012-09-05 | 2013-01-16 | 西安电子科技大学 | 基于线性预测倒谱系数的高光谱图像分类方法 |
CN104777143A (zh) * | 2014-01-15 | 2015-07-15 | 中国人民解放军第二军医大学 | 一种用于过期药的、基于拉曼光谱的相似度鉴定方法 |
CN104502302A (zh) * | 2014-11-04 | 2015-04-08 | 中国航天科工集团第三研究院第八三五七研究所 | 混油的太赫兹时域波形多参数组合定量分析方法 |
Non-Patent Citations (4)
Title |
---|
Comparing Similar Spectra: From Similarity Index to Spectral Contrast Angle;Katty X. Wan et al;《J Am Soc Mass Spectrom》;20021231;第13卷;第85-88页 * |
中药指纹图谱相似度计算方法探析;关洪月 等;《中国实验方剂学杂志》;20110930;第17卷(第18期);第282-287页 * |
中药材红外光谱阵列相关系数比对程序的设计与检验;徐永群 等;《计算机与应用化学》;20020531;第19卷(第3期);第223-224页第1节 * |
中药金礞石红外指纹图谱相似度分析;王栋 等;《光谱学与光谱分析》;20111031;第31卷(第10期);第2716-2717页第3-4节 * |
Also Published As
Publication number | Publication date |
---|---|
CN105784637A (zh) | 2016-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105784637B (zh) | 标识光谱差异性的方法 | |
CN105352895B (zh) | 高光谱遥感数据植被信息提取方法 | |
US11493447B2 (en) | Method for removing background from spectrogram, method of identifying substances through Raman spectrogram, and electronic apparatus | |
CN104458693B (zh) | 用于毒品检测的拉曼光谱测量方法 | |
US10969952B2 (en) | Color and texture match ratings for optimal match selection | |
JP3901137B2 (ja) | におい識別装置 | |
Esbensen et al. | The RPD myth… | |
Cui et al. | Chemometric algorithms for analyzing high dimensional temperature dependent near infrared spectra | |
CN106918566A (zh) | 基于高光谱数据的土壤水分、盐分信息联合提取方法 | |
CN105628581A (zh) | 一种基于高光谱技术的致密砂岩储层露头孔隙度表征方法 | |
CN109564199A (zh) | 分析数据处理方法和分析数据处理装置 | |
Salvatore et al. | Classification methods of multiway arrays as a basic tool for food PDO authentication | |
CN105222892A (zh) | 一种紫外光谱仪的波长定标方法 | |
CN107632010A (zh) | 一种结合激光诱导击穿光谱对钢铁样品的定量方法 | |
Petrova et al. | Measurements and calculations of Ar-broadening parameters of water vapour transitions in a wide spectral region | |
CN107402192A (zh) | 一种快速分析香精香料质量稳定性的方法 | |
Ryu et al. | Two-dimensional correlation analysis and waterfall plots for detecting positional fluctuations of spectral changes | |
CN110836878B (zh) | 一种卷积插值耦合高斯混合模型快速三维荧光寻峰方法 | |
CN101144736A (zh) | 不同类型光谱仪的对比方法 | |
CN102103079B (zh) | 一种光谱分析方法 | |
Winterauer et al. | Nanoscale spatial resolution in far-field Raman imaging using hyperspectral unmixing in combination with positivity constrained super-resolution | |
Wallace et al. | Data analysis methods for synthetic polymer mass spectrometry: Autocorrelation | |
Bennedsen et al. | Discretization of L\'evy semistationary processes with application to estimation | |
CN113496218B (zh) | 一种高光谱遥感敏感波段选择方式的评价方法和系统 | |
CN111426648B (zh) | 一种红外谱图的相似度确定方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20201103 Address after: 201201 room 413, building 1, 111 home business road, Pudong New Area, Shanghai. Patentee after: SHANGHAI OCEANHOOD OPTO-ELECTRONICS TECH Co.,Ltd. Address before: 545006 the Guangxi Zhuang Autonomous Region East Road, Liuzhou, No. 268 Patentee before: GUANGXI University OF SCIENCE AND TECHNOLOGY |