CN106645014A - 基于太赫兹光谱的物质识别方法 - Google Patents
基于太赫兹光谱的物质识别方法 Download PDFInfo
- Publication number
- CN106645014A CN106645014A CN201610847955.XA CN201610847955A CN106645014A CN 106645014 A CN106645014 A CN 106645014A CN 201610847955 A CN201610847955 A CN 201610847955A CN 106645014 A CN106645014 A CN 106645014A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- time
- frequency
- sample
- domain
- 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
- 238000000034 method Methods 0.000 title claims abstract description 17
- 239000000463 material Substances 0.000 title abstract description 6
- 238000004611 spectroscopical analysis Methods 0.000 title abstract description 4
- 238000001228 spectrum Methods 0.000 claims abstract description 102
- 238000012549 training Methods 0.000 claims abstract description 11
- 238000013507 mapping Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims abstract description 8
- 230000003595 spectral effect Effects 0.000 claims abstract description 7
- 238000012360 testing method Methods 0.000 claims description 40
- 239000000126 substance Substances 0.000 claims description 24
- 238000001328 terahertz time-domain spectroscopy Methods 0.000 claims description 16
- 238000010521 absorption reaction Methods 0.000 claims description 14
- 238000010606 normalization Methods 0.000 claims description 3
- 238000011426 transformation method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- WPYMKLBDIGXBTP-UHFFFAOYSA-N benzoic acid Chemical compound OC(=O)C1=CC=CC=C1 WPYMKLBDIGXBTP-UHFFFAOYSA-N 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- LPNBBFKOUUSUDB-UHFFFAOYSA-N p-toluic acid Chemical compound CC1=CC=C(C(O)=O)C=C1 LPNBBFKOUUSUDB-UHFFFAOYSA-N 0.000 description 2
- 238000002203 pretreatment Methods 0.000 description 2
- SPSSULHKWOKEEL-UHFFFAOYSA-N 2,4,6-trinitrotoluene Chemical compound CC1=C([N+]([O-])=O)C=C([N+]([O-])=O)C=C1[N+]([O-])=O SPSSULHKWOKEEL-UHFFFAOYSA-N 0.000 description 1
- ALYNCZNDIQEVRV-UHFFFAOYSA-N 4-aminobenzoic acid Chemical compound NC1=CC=C(C(O)=O)C=C1 ALYNCZNDIQEVRV-UHFFFAOYSA-N 0.000 description 1
- 239000005711 Benzoic acid Substances 0.000 description 1
- GUBGYTABKSRVRQ-DCSYEGIMSA-N Beta-Lactose Chemical compound OC[C@H]1O[C@@H](O[C@H]2[C@H](O)[C@@H](O)[C@H](O)O[C@@H]2CO)[C@H](O)[C@@H](O)[C@H]1O GUBGYTABKSRVRQ-DCSYEGIMSA-N 0.000 description 1
- AUNGANRZJHBGPY-UHFFFAOYSA-N D-Lyxoflavin Natural products OCC(O)C(O)C(O)CN1C=2C=C(C)C(C)=CC=2N=C2C1=NC(=O)NC2=O AUNGANRZJHBGPY-UHFFFAOYSA-N 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- WHUUTDBJXJRKMK-UHFFFAOYSA-N Glutamic acid Natural products OC(=O)C(N)CCC(O)=O WHUUTDBJXJRKMK-UHFFFAOYSA-N 0.000 description 1
- WHUUTDBJXJRKMK-VKHMYHEASA-N L-glutamic acid Chemical compound OC(=O)[C@@H](N)CCC(O)=O WHUUTDBJXJRKMK-VKHMYHEASA-N 0.000 description 1
- COLNVLDHVKWLRT-QMMMGPOBSA-N L-phenylalanine Chemical compound OC(=O)[C@@H](N)CC1=CC=CC=C1 COLNVLDHVKWLRT-QMMMGPOBSA-N 0.000 description 1
- AUNGANRZJHBGPY-SCRDCRAPSA-N Riboflavin Chemical compound OC[C@@H](O)[C@@H](O)[C@@H](O)CN1C=2C=C(C)C(C)=CC=2N=C2C1=NC(=O)NC2=O AUNGANRZJHBGPY-SCRDCRAPSA-N 0.000 description 1
- 229930003471 Vitamin B2 Natural products 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- LSQZJLSUYDQPKJ-NJBDSQKTSA-N amoxicillin Chemical compound C1([C@@H](N)C(=O)N[C@H]2[C@H]3SC([C@@H](N3C2=O)C(O)=O)(C)C)=CC=C(O)C=C1 LSQZJLSUYDQPKJ-NJBDSQKTSA-N 0.000 description 1
- 229960003022 amoxicillin Drugs 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 235000010233 benzoic acid Nutrition 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 235000013922 glutamic acid Nutrition 0.000 description 1
- 239000004220 glutamic acid Substances 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- VIKNJXKGJWUCNN-XGXHKTLJSA-N norethisterone Chemical compound O=C1CC[C@@H]2[C@H]3CC[C@](C)([C@](CC4)(O)C#C)[C@@H]4[C@@H]3CCC2=C1 VIKNJXKGJWUCNN-XGXHKTLJSA-N 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- LSQZJLSUYDQPKJ-UHFFFAOYSA-N p-Hydroxyampicillin Natural products O=C1N2C(C(O)=O)C(C)(C)SC2C1NC(=O)C(N)C1=CC=C(O)C=C1 LSQZJLSUYDQPKJ-UHFFFAOYSA-N 0.000 description 1
- FDPIMTJIUBPUKL-UHFFFAOYSA-N pentan-3-one Chemical class CCC(=O)CC FDPIMTJIUBPUKL-UHFFFAOYSA-N 0.000 description 1
- COLNVLDHVKWLRT-UHFFFAOYSA-N phenylalanine Natural products OC(=O)C(N)CC1=CC=CC=C1 COLNVLDHVKWLRT-UHFFFAOYSA-N 0.000 description 1
- ULSIYEODSMZIPX-UHFFFAOYSA-N phenylethanolamine Chemical class NCC(O)C1=CC=CC=C1 ULSIYEODSMZIPX-UHFFFAOYSA-N 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 229960002477 riboflavin Drugs 0.000 description 1
- 239000011716 vitamin B2 Substances 0.000 description 1
- 235000019164 vitamin B2 Nutrition 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/3581—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
- G01N21/3586—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation by Terahertz time domain spectroscopy [THz-TDS]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Data Mining & Analysis (AREA)
- Immunology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Artificial Intelligence (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Toxicology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供了一种基于太赫兹光谱的物质识别方法，用于对待测样品进行识别，包括如下步骤：步骤1，用太赫兹时域光谱系统扫描预定数量的已知样品；步骤2，对每个时域信号图谱进行预处理，得到每个时域信号图谱的归一化频谱；步骤3，将归一化频谱的时间采样频率映射到预定频率得到统一光谱；步骤4，将统一光谱作为特征向量，对SVM分类器进行训练；步骤5，采用任意一台太赫兹时域光谱系统对待测样品进行扫描，对得到的时域信号图谱进行预处理、采样频率映射；步骤6，用训练后的SVM分类器对待测样品的光谱进行分类，得到待测样品的名称。本发明提供的方法消除了仪器差异及主观差异，具有鲁棒性，能够广泛应用。
Description
技术领域
本发明涉及一种物质识别方法，具体涉及一种基于太赫兹光谱的物质识别方法。
背景技术
太赫兹波是频率在0.1THz到10THz范围的电磁波，其波段位于微波与红外之间，具有丰富的科学意义和广阔的应用前景。太赫兹光子的能量只有毫电子伏特，不会因为电离而破坏被测物质，而且大多数极性分子在THz波段内能够产生共振、吸收能量并且被检测得到相应的光谱，因此太赫兹波谱为无损检测及物质识别提供了一种新的有效途径。
近年来，国内外多个研究机构对太赫兹光谱及基于太赫兹光谱的识别进行了研究。由于太赫兹光谱具有指纹谱的特性，物质的不同化学组成与光谱中的波峰位置具有一一对应的关系，因此太赫兹光谱能够用于进行物质种类的识别。
但是，实际应用过程中，由于太赫兹时域光谱系统的机器参数通常不相同，操作人员之间也存在主观差异，使得不同的太赫兹时域光谱系统对同一种样品的检测结果均存在一定差异。因此，现有的基于太赫兹光谱的识别方法往往需要将数据的来源范围限定在同一参数的仪器、甚至是同一台仪器中。这样的识别方法不具有鲁棒性，结果的准确度难以保障，也难以得到广泛的应用。
发明内容
为解决上述问题，提供一种具有鲁棒性的基于太赫兹光谱的识别方法，本发明采用了如下技术方案：
本发明提供了一种基于太赫兹光谱的物质识别方法，用于对待测样品进行识别，其特征在于，包括如下步骤：
步骤1，用太赫兹时域光谱系统扫描预定数量的已知样品，其中每个已知样品扫描预定次数，得到对应的预定数量的时域信号图谱；
步骤2，对每个时域信号图谱进行预处理，得到每个时域信号图谱的归一化频谱；
步骤3，将步骤2中得到的归一化频谱的时间采样频率映射到预定频率得到统一光谱；
步骤4，将步骤3得到的全部统一光谱作为特征向量，对SVM分类器进行训练，使其按照已知样品对统一光谱中的波谱数据进行分类，得到训练后的SVM分类器；
步骤5，采用任意一台太赫兹时域光谱系统对待测样品进行扫描得到待测样品的时域信号图谱，对该时域信号图谱进行步骤2的预处理得到待测样品的归一化频谱，再将该归一化频谱映射至预定频率得到待测样品的统一光谱；
步骤6，用步骤4得到的训练后的SVM分类器对步骤5得到的统一光谱进行分类，得到待测样品的名称。
本发明提供的基于太赫兹光谱的物质识别方法，还可以具有如下技术特征：
其中，步骤2包括如下子步骤：
步骤2.1，根据已知样品的时域参考信号和时域实际信号的反射峰位置差异，去掉时域信号图谱中的反射峰，得到去反射峰的时域信号图谱；
步骤2.2，将步骤2.1中得到的去反射峰的时域信号图谱进行快速傅里叶变换，然后根据布朗-博朗定律计算已知样品的相对吸收系数，使时域信号图谱转换为吸收系数频谱；
步骤2.3，截取吸收系数频谱有效频域范围内的相对吸收系数，得到反映已知样品特征信息的有效频谱；
步骤2.4，采用正交小波变换法去除步骤2.3中得到的有效频谱中的低频基线及高频噪声，得到去噪频谱；
步骤2.5，将步骤2.4中得到的去噪频谱进行归一化处理，得到归一化频谱。
本发明提供的基于太赫兹光谱的物质识别方法，还可以具有如下技术特征：其中，步骤1所采用的太赫兹时域光谱系统与步骤3所采用的太赫兹时域光谱系统的采样频率不同，预定频率与其中的最低采样频率相同。
本发明提供的基于太赫兹光谱的物质识别方法，还可以具有如下技术特征：其中，步骤1所采用的太赫兹时域光谱系统与步骤3所采用的太赫兹时域光谱系统的采样频率不同，预定频率低于其中的最低采样频率。
本发明提供的基于太赫兹光谱的物质识别方法，还可以具有如下技术特征：其中，步骤3中的映射采用样条函数曲线插值法。
发明作用与效果
本发明提供的基于太赫兹光谱的物质识别方法，将预处理后的归一化频谱映射至统一的预定频率，消除了仪器间参数及设置不同所产生的差异；利用SVM分类器进行光谱分类训练，并利用训练完成后的SVM分类器对待测样品光谱进行待测样品的识别，消除了操作人员之间的主观差异。由此，本发明提供的方法消除了仪器差异及主观差异，因而具有鲁棒性，能够广泛应用。
附图说明
图1为本发明的基于太赫兹光谱的物质识别方法的流程图；
图2为正交小波变换前后的太赫兹波谱图；
图3为不同系统测得的同一种样品在预处理后得到的波谱图；
图4为SVM分类器识别测试数据的分类结果图。
具体实施方式
以下结合附图及实施例来说明本发明的具体实施方式。
<实施例>
图1为本发明的基于太赫兹光谱的物质识别方法的流程图。
如图1所示，本实施例中，基于太赫兹光谱的物质识别方法包括如下步骤：
步骤1，采用太赫兹时域光谱系统扫描已知样品，获得每个已知样品在不同采样频率下的时域信号图谱。
在本实施例中，采用来自上海现代光学系统重点实验室(以下简称实验室)的太赫兹时域光谱系统(以下简称THz-TDS)对已知样品进行扫描，采用来自和上海高晶影像科技有限公司(以下简称高晶公司)的THz-TDS对待测样品进行扫描。其中，实验室THz-TDS的采样频率为0.009THz，高晶公司THz-TDS的采样频率为0.01THz。
本实施例中，选择的已知样品为如下13种：4-氨基苯甲酸，A2(1-苯基-2-吡咯烷-1-基乙酮)，A3(2-氨基-1-苯基乙醇)，C5(一种混合物炸药)，D(+)-乳糖水合物，L(+)-酒石酸，TNT，阿莫西林，苯丙氨酸，苯甲酸，对甲基苯甲酸，谷氨酸，维生素B2；每种样品分别用两台THz-TDS各扫描10次，即、每种已知样品得到采样频率0.009THz时域信号图谱10个，扫描结束后，总共得到130张时域信号图谱。
步骤2，对步骤1得到的每个时域信号图谱进行预处理，具体包括如下子步骤：
步骤2.1，以只经过空气的THz信号作为时域参考信号，根据该时域参考信号和时域实际信号的反射峰位置差异，去掉时域信号图谱中的反射峰，得到去反射峰的时域信号图谱。
步骤2.2，将步骤2.1中得到的去反射峰的时域信号图谱进行快速傅里叶变换，然后根据布朗-博朗定律计算已知样品的相对吸收系数，对时域信号图谱进行转换，得到吸收系数频谱。该吸收系数频谱中，纵坐标为吸收系数，横坐标为频率。
步骤2.3，截取0.2THz～1.5THz有效频域范围内的相对吸收系数，得到反映已知样品特征信息的有效频谱。因为在这段频域范围内，物质的吸收峰明显，物质容易简单的被区分；而在1.5THz后有大量的噪声信息，因此不采用。
步骤2.4，采用正交小波变换法去除步骤2.3中得到的有效频谱中的低频基线及高频噪声，得到去噪频谱。
图2为正交小波变换前后的太赫兹波谱图。
图2示出了其中一个图谱在正交小波变换前后的变化。如图2所示，在进行正交小波变换前的波谱(即、图中的原波谱)的吸收系数纵坐标在3～4.5之间，经过变换后变为-0.5～0.5之间。原波谱中，1.2THz～1.5THz区间内的曲线上扬并且具有一些小的波峰，说明此处出现了基线漂移现象和噪声；经过正交小波变换后，该区间的曲线整体波形平缓延伸，并且整体更加平滑，说明正交小波变换消除了基线漂移现象以及噪声。
步骤2.5，将步骤2.4中得到的去噪频谱进行归一化处理，得到归一化频谱。
图3为不同系统测得的同一种样品在预处理后得到的波谱图。
图3示出了样品A2在两种THz-TDS中测得的图谱经过预处理得到的归一化图谱。如图3所示，经过预处理，归一化图谱中的吸收系数纵坐标范围均在0～1之间，并且两个图谱中的特征峰基本一致。
步骤3，将步骤2中得到的归一化频谱的时间采样频率映射到预定频率得到统一光谱，使全部统一光谱的采样频率一致。
在本实施例中，该预定频率采用两台THz-TDS所具有的两个采样频率中最小的频率，即0.009THz。由于采用实验室THz-TDS得到的图谱中，采样频率已经是0.009THz了，因此不进行采样频率映射处理，直接将其归一化图谱作为统一光谱即可。
步骤4，将步骤3得到的各个统一光谱作为特征向量，对SVM分类器进行训练，使其按照已知样品中的波谱数据对统一光谱进行分类，得到训练后的SVM分类器。
本实施例中，采用的SVM分类器为Libsvm工具箱。将步骤3得到的统一光谱在0.2THz～1.5THz的范围内进行等间隔采样，每个统一光谱得到143个波谱数据，将这些波谱数据作为特征值，对SVM分类器进行训练，该130个统一光谱及每个统一光谱对应的143个波谱数据可以看成是一个二维向量。为使分类准确率更高，SVM分类器的训练策略是：采用基于网格搜索法选取最佳惩罚参数和核函数参数的参数寻优算法，利用最佳参数来进行训练。
步骤5，采用任意一台太赫兹时域光谱系统对待测样品进行扫描得到待测样品的时域信号图谱，对该时域信号图谱进行步骤2的预处理得到待测样品的归一化频谱，再将该归一化频谱映射至预定频率得到待测样品的统一光谱。
在本实施例中，采用来自上海高晶影像科技有限公司(以下简称高晶公司)的THz-TDS对待测样品进行扫描，该待测样品选用了与已知样品相同的13种。对每种待测样品进行10次扫描，总共得到130个时域信号图谱。对该130个时域信号图谱进行步骤2及步骤3的处理，由于高晶公司THz-TDS的采样频率与预定频率不同，因此对高晶公司THz-TDS得到的图谱进行采用频率映射处理，将其中的采样频率0.01THz映射到0.009THz得到对应的统一光谱。
步骤6，用步骤4得到的训练后的SVM分类器对步骤5得到的统一光谱进行分类，得到待测样品的名称。
图4为SVM分类器识别测试数据的分类结果图。
将高晶公司THz-TDS得到的130个统一光谱作为测试集，将0.2THz～1.5THz的范围内的全部143维波谱向量作为SVM分类器的特征向量，对训练后的SVM分类器进行分类测试。结果如图4所示，用该训练后的SVM分类器对高晶公司的波谱数据进行分类，其预测测试集分类绝大部分能够落入实际测试集分类，分类的准确率是93.07％。即、对13种待测样品的识别准确率为93.07％。
实施例作用与效果
本实施例提供的基于太赫兹光谱的物质识别方法，将预处理后的归一化频谱映射至统一的预定频率，消除了仪器间参数及设置不同所产生的差异；利用SVM分类器进行光谱分类训练，并利用训练完成后的SVM分类器对待测样品光谱进行待测样品的识别，消除了操作人员之间的主观差异。由此，本发明提供的方法消除了仪器差异及主观差异，因而具有鲁棒性以及广泛应用的潜力。
此外，在本实施例中，已知样品与待测样品各用一台THz-TDS进行扫描，并且这两台THz-TDS的采样频率不同，采样频率映射(即步骤3)所采用的预定频率与其中的最小频率一致，因此可以不用对最小采样频率的图谱进行采用频率映射操作，使得图谱处理过程能够快速完成。
上述实施例仅用于说明本发明的具体实施方式，本发明的方法并不局限于上述实施例的范围。
例如，在实施例中，已知样品与待测样品各用一台THz-TDS进行扫描，并且这两台THz-TDS的采样频率不同，但在本发明中，已知样品的扫描与待测样品的扫描均可以采用任意数量、频率相同或不同的THz-TDS来进行，只要将预定频率设定得不高于所采用的THz-TDS中的最低采样频率即可，如此即可将所有统一光谱中的采样频率统一，使得后续的SVM分类器训练以及SVM识别能够实现。
在实施例中，已知样品及待测样品均为13种，但本发明中，已知样品的数量可以不为13种，待测样品也可以不为13种，只要待测样品属于已知样品中的一种或几种，则本发明的方法均能够实现待测样品的识别。
Claims (5)
1.一种基于太赫兹光谱的物质识别方法，用于对待测样品进行识别，其特征在于，包括如下步骤：
步骤1，用太赫兹时域光谱系统扫描预定数量的已知样品，其中每个所述已知样品扫描预定次数，得到对应的预定数量的时域信号图谱；
步骤2，对每个所述时域信号图谱进行预处理，得到每个所述时域信号图谱的归一化频谱；
步骤3，将步骤2中得到的所述归一化频谱的时间采样频率映射到预定频率得到统一光谱；
步骤4，将步骤3得到的全部所述统一光谱作为特征向量，对SVM分类器进行训练，使其按照所述已知样品对所述统一光谱中的波谱数据进行分类，得到训练后的SVM分类器；
步骤5，采用任意一台太赫兹时域光谱系统对所述待测样品进行扫描得到所述待测样品的时域信号图谱，对该时域信号图谱进行步骤2的预处理得到所述待测样品的归一化频谱，再将该归一化频谱映射至所述预定频率得到所述待测样品的统一光谱；
步骤6，用步骤4得到的所述训练后的SVM分类器对步骤5得到的所述统一光谱进行分类，得到所述待测样品的名称。
2.根据权利要求1所述的基于太赫兹光谱的物质识别方法，其特征在于，
其中，所述步骤2包括如下子步骤：
步骤2.1，根据所述已知样品的时域参考信号和时域实际信号的反射峰位置差异，去掉所述时域信号图谱中的反射峰，得到去反射峰的时域信号图谱；
步骤2.2，将步骤2.1中得到的所述去反射峰的时域信号图谱进行快速傅里叶变换，然后根据布朗-博朗定律计算所述已知样品的相对吸收系数，使所述时域信号图谱转换为吸收系数频谱；
步骤2.3，截取所述吸收系数频谱有效频域范围内的所述相对吸收系数，得到反映所述已知样品特征信息的有效频谱；
步骤2.4，采用正交小波变换法去除步骤2.3中得到的所述有效频谱中的低频基线及高频噪声，得到去噪频谱；
步骤2.5，将步骤2.4中得到的所述去噪频谱进行归一化处理，得到所述归一化频谱。
3.根据权利要求1所述的基于太赫兹光谱的物质识别方法，其特征在于：
其中，步骤1所采用的太赫兹时域光谱系统与步骤3所采用的太赫兹时域光谱系统的采样频率不同，所述预定频率与其中的最低采样频率相同。
4.根据权利要求1所述的基于太赫兹光谱的物质识别方法，其特征在于：
其中，步骤1所采用的太赫兹时域光谱系统与步骤3所采用的太赫兹时域光谱系统的采样频率不同，所述预定频率低于其中的最低采样频率。
5.根据权利要求1所述的基于太赫兹光谱的物质识别方法，其特征在于：
其中，步骤3中的所述映射采用样条函数曲线插值法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610847955.XA CN106645014B (zh) | 2016-09-23 | 2016-09-23 | 基于太赫兹光谱的物质识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610847955.XA CN106645014B (zh) | 2016-09-23 | 2016-09-23 | 基于太赫兹光谱的物质识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106645014A true CN106645014A (zh) | 2017-05-10 |
| CN106645014B CN106645014B (zh) | 2019-04-30 |
Family
ID=58853434
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610847955.XA Active CN106645014B (zh) | 2016-09-23 | 2016-09-23 | 基于太赫兹光谱的物质识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106645014B (zh) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107220625A (zh) * | 2017-06-01 | 2017-09-29 | 中国海洋大学 | 一种对样品的libs光谱自动识别分类方法 |
| CN107561033A (zh) * | 2017-09-21 | 2018-01-09 | 上海理工大学 | 基于太赫兹光谱的混合物中关键物质定性和定量测定方法 |
| CN108088811A (zh) * | 2017-12-27 | 2018-05-29 | 上海理工大学 | 太赫兹波测定混合气体中各组分浓度的方法 |
| CN108279217A (zh) * | 2018-04-28 | 2018-07-13 | 江苏建筑职业技术学院 | 一种基于太赫兹时域光谱的煤岩判别方法 |
| CN108489917A (zh) * | 2018-02-28 | 2018-09-04 | 浙江大学 | 一种去除太赫兹吸收光谱中低频仪器误差和高频振荡噪声的方法 |
| CN108827904A (zh) * | 2018-06-19 | 2018-11-16 | 广东工业大学 | 基于太赫兹光谱的物质识别方法、装置、设备及存储介质 |
| CN108918458A (zh) * | 2018-10-12 | 2018-11-30 | 北京环境特性研究所 | 一种确定材料太赫兹吸收峰的方法 |
| CN109470720A (zh) * | 2018-09-29 | 2019-03-15 | 深圳市华讯方舟太赫兹科技有限公司 | 液体识别方法、向量提取方法、液体识别装置及存储介质 |
| CN110006844A (zh) * | 2019-05-22 | 2019-07-12 | 安徽大学 | 基于函数性主元分析的近红外光谱特征提取方法和系统 |
| WO2020052195A1 (zh) * | 2018-09-14 | 2020-03-19 | 深圳市太赫兹科技创新研究院有限公司 | 玉石鉴别方法、装置、系统和存储介质 |
| WO2020062471A1 (zh) * | 2018-09-29 | 2020-04-02 | 深圳市华讯方舟太赫兹科技有限公司 | 液体识别方法、特征提取方法、液体识别装置及存储装置 |
| CN112129728A (zh) * | 2020-09-27 | 2020-12-25 | 上海理工大学 | 一种药物中咖啡因的定性鉴定和定量测定的方法 |
| CN112485218A (zh) * | 2020-11-05 | 2021-03-12 | 电子科技大学中山学院 | 一种基于人工神经网络的太赫兹危险液体识别方法 |
| CN114088656A (zh) * | 2020-07-31 | 2022-02-25 | 中国科学院上海高等研究院 | 太赫兹光谱物质识别方法及系统、存储介质及终端 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060255277A1 (en) * | 2003-03-21 | 2006-11-16 | Cole Bryan E | Spectroscopy apparatus and associated technique |
| CN103472032A (zh) * | 2013-09-24 | 2013-12-25 | 浙江大学 | 一种利用太赫兹时域光谱技术检测盐酸四环素的方法 |
| CN105527243A (zh) * | 2015-11-30 | 2016-04-27 | 中国科学院重庆绿色智能技术研究院 | 一种基于太赫兹光谱技术测量石墨烯薄膜电导的方法 |
-
2016
- 2016-09-23 CN CN201610847955.XA patent/CN106645014B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060255277A1 (en) * | 2003-03-21 | 2006-11-16 | Cole Bryan E | Spectroscopy apparatus and associated technique |
| CN103472032A (zh) * | 2013-09-24 | 2013-12-25 | 浙江大学 | 一种利用太赫兹时域光谱技术检测盐酸四环素的方法 |
| CN105527243A (zh) * | 2015-11-30 | 2016-04-27 | 中国科学院重庆绿色智能技术研究院 | 一种基于太赫兹光谱技术测量石墨烯薄膜电导的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 赵树森等: "用支持向量机识别毒品的太赫兹吸收光谱", 《中国激光》 * |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107220625A (zh) * | 2017-06-01 | 2017-09-29 | 中国海洋大学 | 一种对样品的libs光谱自动识别分类方法 |
| CN107561033A (zh) * | 2017-09-21 | 2018-01-09 | 上海理工大学 | 基于太赫兹光谱的混合物中关键物质定性和定量测定方法 |
| CN108088811A (zh) * | 2017-12-27 | 2018-05-29 | 上海理工大学 | 太赫兹波测定混合气体中各组分浓度的方法 |
| CN108088811B (zh) * | 2017-12-27 | 2020-09-25 | 上海理工大学 | 太赫兹波测定混合气体中各组分浓度的方法 |
| CN108489917B (zh) * | 2018-02-28 | 2020-09-01 | 浙江大学 | 一种去除太赫兹吸收光谱中低频误差和高频噪声的方法 |
| CN108489917A (zh) * | 2018-02-28 | 2018-09-04 | 浙江大学 | 一种去除太赫兹吸收光谱中低频仪器误差和高频振荡噪声的方法 |
| CN108279217A (zh) * | 2018-04-28 | 2018-07-13 | 江苏建筑职业技术学院 | 一种基于太赫兹时域光谱的煤岩判别方法 |
| CN108279217B (zh) * | 2018-04-28 | 2020-12-22 | 江苏建筑职业技术学院 | 一种基于太赫兹时域光谱的煤岩判别方法 |
| CN108827904A (zh) * | 2018-06-19 | 2018-11-16 | 广东工业大学 | 基于太赫兹光谱的物质识别方法、装置、设备及存储介质 |
| CN108827904B (zh) * | 2018-06-19 | 2021-01-26 | 广东工业大学 | 基于太赫兹光谱的物质识别方法、装置、设备及存储介质 |
| WO2020052195A1 (zh) * | 2018-09-14 | 2020-03-19 | 深圳市太赫兹科技创新研究院有限公司 | 玉石鉴别方法、装置、系统和存储介质 |
| WO2020062471A1 (zh) * | 2018-09-29 | 2020-04-02 | 深圳市华讯方舟太赫兹科技有限公司 | 液体识别方法、特征提取方法、液体识别装置及存储装置 |
| CN109470720A (zh) * | 2018-09-29 | 2019-03-15 | 深圳市华讯方舟太赫兹科技有限公司 | 液体识别方法、向量提取方法、液体识别装置及存储介质 |
| CN109470720B (zh) * | 2018-09-29 | 2023-09-22 | 深圳市重投华讯太赫兹科技有限公司 | 液体识别方法、向量提取方法、液体识别装置及存储介质 |
| CN108918458B (zh) * | 2018-10-12 | 2020-12-04 | 北京环境特性研究所 | 一种确定材料太赫兹吸收峰的方法 |
| CN108918458A (zh) * | 2018-10-12 | 2018-11-30 | 北京环境特性研究所 | 一种确定材料太赫兹吸收峰的方法 |
| CN110006844A (zh) * | 2019-05-22 | 2019-07-12 | 安徽大学 | 基于函数性主元分析的近红外光谱特征提取方法和系统 |
| CN114088656A (zh) * | 2020-07-31 | 2022-02-25 | 中国科学院上海高等研究院 | 太赫兹光谱物质识别方法及系统、存储介质及终端 |
| CN112129728A (zh) * | 2020-09-27 | 2020-12-25 | 上海理工大学 | 一种药物中咖啡因的定性鉴定和定量测定的方法 |
| CN112129728B (zh) * | 2020-09-27 | 2023-08-29 | 上海理工大学 | 一种药物中咖啡因的定性鉴定和定量测定的方法 |
| CN112485218A (zh) * | 2020-11-05 | 2021-03-12 | 电子科技大学中山学院 | 一种基于人工神经网络的太赫兹危险液体识别方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106645014B (zh) | 2019-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106645014A (zh) | 基于太赫兹光谱的物质识别方法 | |
| CN108469560B (zh) | 一种基于快速s变换时频空间模型的电磁干扰客观复杂度评估方法 | |
| CN108279217B (zh) | 一种基于太赫兹时域光谱的煤岩判别方法 | |
| Wang et al. | Comparing the applications of EMD and EEMD on time–frequency analysis of seismic signal | |
| Mahdavi et al. | Application of the local maximum synchrosqueezing transform for seismic data | |
| CN105279379B (zh) | 基于凸组合核函数主成分分析的太赫兹光谱特征提取方法 | |
| Dutta et al. | Image processing based classification of grapes after pesticide exposure | |
| CN107422381B (zh) | 一种基于eemd-ica的地震低频信息流体预测方法 | |
| CN104020135B (zh) | 基于近红外光谱的校正模型建模方法 | |
| CN107561033A (zh) | 基于太赫兹光谱的混合物中关键物质定性和定量测定方法 | |
| Chen et al. | An ameliorated synchroextracting transform based on upgraded local instantaneous frequency approximation | |
| CN113984736B (zh) | 一种基于空间偏移拉曼光谱的带包装食品信号分离方法 | |
| CN104316487A (zh) | 违禁物品检测方法和装置 | |
| CN105447464A (zh) | 一种基于pso的电能质量扰动识别与分类方法 | |
| CN105203495B (zh) | 一种基于希尔伯特-黄变换的光谱信号去噪方法 | |
| Wang et al. | Wood species identification using terahertz time-domain spectroscopy | |
| Cooper et al. | Wavelet based Raman spectra comparison | |
| CN104181508B (zh) | 基于压缩感知的威胁雷达信号检测方法 | |
| CN106770018A (zh) | 近红外光谱快速测定油料种类的方法 | |
| CN113281809A (zh) | 一种地震信号的谱分析方法 | |
| Capobianco et al. | Pigment identification in pictorial layers by HyperSpectral Imaging | |
| Shinzawa et al. | Hydrogen/deuterium (H/D) exchange of gelatinized starch studied by two-dimensional (2D) near-infrared (NIR) correlation spectroscopy | |
| CN105067556A (zh) | 一种快速检测生物液态样品的光学装置及光学方法 | |
| Tao et al. | Discriminant analysis of red wines from different aging ways by information fusion of NIR and MIR spectra | |
| CN109270573B (zh) | 一种快速保频保幅s变换方法 |
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 |