CN104251838B - 违禁物品检测方法和装置 - Google Patents
违禁物品检测方法和装置 Download PDFInfo
- Publication number
- CN104251838B CN104251838B CN201410429717.8A CN201410429717A CN104251838B CN 104251838 B CN104251838 B CN 104251838B CN 201410429717 A CN201410429717 A CN 201410429717A CN 104251838 B CN104251838 B CN 104251838B
- Authority
- CN
- China
- Prior art keywords
- data
- measured object
- standard substance
- absorption spectra
- frequency
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 57
- 239000000126 substance Substances 0.000 claims abstract description 164
- 238000000862 absorption spectrum Methods 0.000 claims abstract description 124
- 238000005259 measurement Methods 0.000 claims description 15
- 230000000875 corresponding Effects 0.000 claims description 12
- 239000000284 extract Substances 0.000 claims description 5
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 235000005035 ginseng Nutrition 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 241000196324 Embryophyta Species 0.000 claims 1
- 230000001131 transforming Effects 0.000 claims 1
- 239000000203 mixture Substances 0.000 description 6
- 238000010521 absorption reaction Methods 0.000 description 4
- 238000000034 method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 229940079593 drugs Drugs 0.000 description 2
- 230000003203 everyday Effects 0.000 description 2
- 239000002360 explosive Substances 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006011 modification reaction Methods 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003595 spectral Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Abstract
本发明公开了一种违禁物品检测方法和装置,所述方法包括:将被测物放置于太赫兹光路中,测量被测物的太赫兹时域波形数据;对太赫兹时域波形数据进行快速傅里叶变换,得到被测物的太赫兹频域波形数据;根据被测物的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出被测物的吸收谱数据;将被测物的吸收谱数据与数据库中预先存储的各种违禁物品的标准品的特征峰数据进行比较,判断被测物是否为其中一种违禁物品;违禁物品的标准品的特征峰数据具体是从该标准品的吸收谱数据的各波峰中选取的波峰的数据。本发明的技术方案中,根据被测物的吸收谱数据与各种违禁物品的标准品的特征峰数据的对比结果,可以更为准确地判断出被测物是否为违禁物品。
Description
技术领域
本发明涉及光谱领域,尤其涉及一种违禁物品检测方法和装置。
背景技术
目前,对违禁物品的检测已经广泛地应用在机场、海关、车站、港口、核电站、政府机关、邮检中心、法院、大使馆等场所,用以检测出人员随身携带、或者行李中的违禁物品,从而保障安全。
现有的一种违禁物品检测方法为:利用X光照射被测物,根据X光透射的原理,显示出被测物的形状轮廓,并对不同密度的被测物分别以不同的颜色进行显示,例如有机物显示为橙色,无机物显示为蓝色,混合物显示为绿色。检测人员根据该违禁物品检测装置显示的形状轮廓和颜色,可以较为容易地识别出枪支、刀具等金属类的违禁物品;但是对于化学类的违禁物品(例如炸药、毒品),该违禁物品的形状通常很难预测,显示的颜色又通常与日用品显示的颜色相同或相近,因此不容易识别出化学类的违禁物品。
现有的另一种违禁物品检测方法为:利用X光照射被测物,根据康普顿散射原理,显示出被测物的形状轮廓,并以不同的灰度分别显示低原子序数、高原子序数的被测物。检测人员根据该违禁物品检测装置显示的形状轮廓和灰度,可以较为容易地识别出枪支、刀具等金属类的违禁物品。但是对于化学类的违禁物品(例如炸药、毒品),该违禁物品的形状通常很难预测;事实上,无论是违禁物品还是日用品都可以只包含低原子序数的成分、只包含高原子序数的成分,或者既包含低原子序数的成分又包含高原子序数的成分;导致该违禁物品显示的灰度通常与日用品显示的灰度相同或相近;造成检测人员根据该违禁物品检测方法不容易识别出化学类的违禁物品。
综上,利用现有的违禁物品检测方法,判断被测物是否为违禁物品的准确性较低;因此,有必要提供一种准确性更高的违禁物品检测方法和装置。
发明内容
针对上述现有技术存在的缺陷,本发明实施例提供了一种违禁物品检测方法和装置,以更为准确地判断出被测物是否为违禁物品。
本发明的技术方案根据一个方面,提供了一种违禁物品检测方法,包括:
将被测物放置于太赫兹光路中,测量所述被测物的太赫兹时域波形数据;
对得到的太赫兹时域波形数据进行快速傅里叶变换,得到所述被测物的太赫兹频域波形数据;
根据所述被测物的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出所述被测物的吸收谱数据;
将所述被测物的吸收谱数据与数据库中预先存储的各种违禁物品的标准品的特征峰数据进行比较,判断所述被测物是否为其中一种违禁物品;其中,所述违禁物品的标准品的特征峰数据具体是从该标准品的吸收谱数据的各波峰中选取的波峰的数据。
较佳地,所述参考太赫兹频域波形数据是根据如下方法预先计算得到的:
在所述太赫兹光路中,测量没有所述被测物的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;
将所述参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据。
较佳地,所述根据所述被测物的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出所述被测物的吸收谱数据,具体包括:
针对所述被测物的太赫兹频域波形数据中的每个频点,根据所述被测物的太赫兹频域波形数据、以及参考太赫兹频域波形数据在该频点的fft值,计算所述被测物在该频点的吸收系数;
由计算得到的所述被测物在各频点的吸收系数组成所述被测物的吸收谱数据。
较佳地,所述将所述被测物的吸收谱数据数据库中预先存储的各种违禁物品的标准品的特征峰数据进行比较,判断所述被测物是否为其中一种违禁物品,具体包括:
对于所述数据库中预先存储的一种违禁物品的标准品的特征峰数据,将该标准品的吸收谱数据的特征峰所对应的频点作为中心频点,将该中心频点两边设定范围的频段作为比对窗口;判断所述被测物的吸收谱数据在所述比对窗口内是否存在波峰;若不存在,则判断所述被测物不是该种违禁物品。
进一步,所述数据库中还存储有该种违禁物品的标准品的吸收谱数据;以及
在所述判断所述被测物的吸收谱数据在所述比对窗口内是否存在波峰后,还包括:
若判断存在波峰,则将所述被测物的吸收谱数据与该种违禁物品的标准品的吸收谱数据进行比较,判断所述被测物是否为该种违禁物品。
较佳地,所述将所述被测物的吸收谱数据与该种违禁物品的标准品的吸收谱数据进行比较,判断所述被测物是否为该种违禁物品,具体包括:
统计所述被测物的吸收谱数据的波峰总数;
对于该标准品的吸收谱数据的每个波峰,在该波峰所对应的频点为中心的设定窗口频段内,若判断所述被测物的吸收谱数据存在波峰,则将该频点记为一个相似点,并计算两个波峰的峰值差,作为该相似点处的峰值差;
根据得到的相似点的个数、每个相似点处的峰值差、所述被测物的吸收谱数据的波峰总数,以及该标准品的吸收谱数据的波峰总数,计算所述被测物与该标准品之间的相似度;
若计算得到的相似度大于设定的相似度阈值,则判断所述被测物是该种违禁物品;否则,判断所述被测物不是该种违禁物品。
较佳地,所述根据得到的相似点的个数、每个相似点处的峰值差、所述被测物的吸收谱数据的波峰总数,以及该标准品的吸收谱数据的波峰总数,计算所述被测物与该标准品之间的相似度,具体为:
根据如下公式2计算得到所述被测物与该标准品之间的相似度S:
其中,nsim、nsam、nstd分别表示所述相似点的个数、所述被测物的吸收谱数据的波峰总数、该标准品的吸收谱数据的波峰总数,i为不小于1且不大于nsim的自然数,pi表示第i个相似点处的峰值差,ω1、ω2、ω3分别表示三个权重系数。
较佳地,所述数据库根据如下方法预先建立:
针对每种违禁物品,提取该种违禁物品的标准品基于太赫兹频段的特征数据:
将该种违禁物品的标准品放置于太赫兹光路中,测量该标准品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据;根据该标准品的太赫兹频域波形数据以及所述参考太赫兹频域波形数据,计算出该标准品的吸收谱数据;从该标准品的吸收谱数据中确定出该标准品的特征峰数据后,将该标准品的特征峰数据作为提取出的该种违禁物品的标准品的特征数据;
将提取出的各种违禁物品的标准品的特征数据存储到所述数据库中。
进一步,在所述将提取出的各种违禁物品的标准品的特征数据存储到所述数据库前,还包括:
对于每种违禁物品,将该种违禁物品的标准品的吸收谱数据也作为该种违禁物品的标准品的特征数据。
本发明的技术方案根据另一个方面,还提供了一种违禁物品检测装置,包括:
时域数据检测模块,用于测量太赫兹光路中经过被测物的太赫兹时域波形数据;
频域数据计算模块,用于对所述时域数据检测模块测量得到的太赫兹时域波形数据进行快速傅里叶变换,得到所述被测物的太赫兹频域波形数据;
吸收谱计算模块,用于根据所述频域数据计算模块得到的所述被测物的太赫兹频域波形数据,以及参考太赫兹频域波形数据,计算出所述被测物的吸收谱数据;
特征比对模块,用于将所述吸收谱计算模块计算出的被测物的吸收谱数据与数据库中预先存储的各种违禁物品的标准品的特征峰数据进行比较,判断所述被测物是否为其中一种违禁物品;其中,所述违禁物品的标准品的特征峰数据具体是从该标准品的吸收谱数据的各波峰中选取的波峰的数据。
进一步,所述违禁物品检测装置,还包括:
参考时域波形测量模块,用于在所述太赫兹光路中,测量没有所述被测物的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;
参考频域波形计算模块,用于将所述参考时域波形测量模块测量得到的参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据进行存储。
进一步,所述违禁物品检测装置,还包括:
标准品数据库建立模块,用于针对每种违禁物品,将该种违禁物品的标准品放置于太赫兹光路中,测量该标准品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据;根据该标准品的太赫兹频域波形数据以及所述参考频域波形计算模块得到的参考太赫兹频域波形数据,计算出该标准品的吸收谱数据;从该标准品的吸收谱数据中确定出该标准品的特征峰数据后,将该种违禁物品的标准品的特征峰数据和吸收谱数据存入到所述数据库中。
本发明的技术方案中,将被测物吸收太赫兹波的吸收谱数据,与每种违禁物品的标准品对太赫兹波的吸收谱数据中的特征数据进行对比,可以根据对比结果更为准确地判断出被测物是否为违禁物品。
而且,还可以根据对比结果计算出相似度的值,并根据相似度的值判断被测物是否为违禁物品,进一步提高了判断的准确性。
此外,被测物吸收太赫兹波的吸收谱数据、以及每种违禁物品的标准品的特征数据都是从各自时域波形数据转换出的频域数据;与时域波形数据相比,数据量大大减少,使得对比时的计算量大大降低;从而可以加快对比速率,提高判断效率。
附图说明
图1为本发明实施例的获取太赫兹光路的参考太赫兹频域波形数据的方法的流程示意图;
图2为本发明实施例的提取每种违禁物品的标准品基于太赫兹频段的特征数据的方法的流程示意图;
图3为本发明实施例的违禁物品检测方法的流程示意图;
图4为本发明实施例的将被测物的吸收谱数据与一种违禁物品的标准品的吸收谱数据进行比较,判断该被测物是否为该种违禁物品的方法的流程示意图;
图5为本发明实施例的被测物的吸收谱数据与一种违禁物品的标准品的吸收谱数据进行比较的实例示意图;
图6为本发明实施例的违禁物品检测装置的内部结构的框架示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举出优选实施例,对本发明进一步详细说明。然而,需要说明的是,说明书中列出的许多细节仅仅是为了使读者对本发明的一个或多个方面有一个透彻的理解,即便没有这些特定的细节也可以实现本发明的这些方面。
本申请使用的“模块”、“系统”等术语旨在包括与计算机相关的实体,例如但不限于硬件、固件、软硬件组合、软件或者执行中的软件。例如,模块可以是,但并不仅限于:处理器上运行的进程、处理器、对象、可执行程序、执行的线程、程序和/或计算机。举例来说,计算设备上运行的应用程序和此计算设备都可以是模块。一个或多个模块可以位于执行中的一个进程和/或线程内。
本发明的发明人考虑到,不同的违禁物品对太赫兹(Thz)频段的电磁波具有不同的吸收特征;采用太赫兹频段的电磁波(本文中简称为太赫兹波)照射被测物,得到被测物吸收太赫兹波的数据;并与预先测得的各种违禁物品吸收太赫兹波的数据进行比较,从而可以更为准确地判断出被测物是否为违禁物品。太赫兹波的波长在30μm-3mm之间,在电磁波谱中介于微波与红外光波之间,属于远红外波段。
而且,被测物吸收太赫兹波的数据、以及违禁物品吸收太赫兹波的数据都可以转换到频域后,进行比较;可以加快比较速率,提高判断的效率。
下面结合附图详细说明本发明的技术方案。
本发明实施例提供的违禁物品检测装置实施违禁物品检测方法之前,需要获取太赫兹光路的参考太赫兹频域波形数据,并建立违禁物品的标准品的数据库。
其中,违禁物品检测装置获取太赫兹光路的参考太赫兹频域波形数据的方法,其流程示意图如图1所示,包括如下步骤:
S101:在太赫兹光路中,测量没有被测物的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据。
具体地,违禁物品检测装置开启太赫兹光路后,周期性采集在该太赫兹光路中未放置被测物的情况下的太赫兹波(简称参考太赫兹波)的信号强度值,得到太赫兹时域波形数据,作为参考太赫兹时域波形数据。
参考太赫兹时域波形数据可以记为k行两列的数组(t,ref),t列中的数据为各采集时刻,ref列中的数据为参考太赫兹波在各采集时刻的信号强度值,k为自然数。
S102:将参考太赫兹时域波形数据进行快速傅里叶变换,得到参考太赫兹频域波形数据。
具体地,违禁物品检测装置对步骤S101中得到的参考太赫兹时域波形数据,进行2m个变换点(即频点)的快速傅里叶变换fft,得到参考太赫兹频域波形数据。m为自然数且2m比2m-1、2m+1更接近于k;较佳地,2m等于1024。
参考太赫兹频域波形数据可以记为2m行两列的数组(ω,fft(ref,2m)),ω列中的数据为各频点,fft(ref,2m)列中的数据为参考太赫兹频域波形数据中的在各频点的fft值。
本发明实施例中,建立违禁物品的标准品的数据库的方法包括:针对每种违禁物品,提取该种违禁物品的标准品基于太赫兹频段的特征数据;并将提取出的各种违禁物品的标准品的特征数据存到该数据库中。
其中,违禁物品检测装置对于每种违禁物品,提取该种违禁物品的标准品基于太赫兹频段的特征数据的方法,其流程示意图如图2所示,包括如下步骤:
S201:将违禁物品的标准品放置于太赫兹光路中,测量该标准品的太赫兹时域波形数据。
具体地,技术人员将违禁物品的标准品放置于太赫兹光路中后,违禁物品检测装置周期性采集太赫兹波的信号强度值,得到该标准品的太赫兹时域波形数据。
违禁物品的标准品的太赫兹时域波形数据可以记为k行两列的数组(t,sam),t列中的数据为各采集时刻,sam列中的数据为该标准品的太赫兹波在各采集时刻的信号强度值。
S202:对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据。
具体地,违禁物品检测装置对步骤S201中得到的违禁物品的标准品的太赫兹时域波形数据,根据与步骤S102中相同的2m个变换点,进行快速傅里叶变换fft,得到该标准品的太赫兹频域波形数据。
违禁物品的标准品的太赫兹频域波形数据可以记为2m行两列的数组(ω,fft(sam,2m)),ω列中的数据为各频点,fft(sam,2m)列中的数据为该标准品的太赫兹频域波形数据中的在各频点的fft值。
S203:根据违禁物品的标准品的太赫兹频域波形数据以及步骤S102中的参考太赫兹频域波形数据,计算出该标准品的吸收谱数据。
具体地,违禁物品检测装置对于步骤S202中得到的违禁物品的标准品的太赫兹频域波形数据中的每个频点ωj,根据该标准品的太赫兹频域波形数据中的在该频点的fft值、参考太赫兹频域波形数据在该频点的fft值,以及如下公式1计算该标准品在该频点的吸收系数αj:
公式1中,j为0到2m-1之间的非负整数;fftj(sam,2m)、fftj(ref,2m)分别表示该标准品的太赫兹频域波形数据、参考太赫兹频域波形数据在频点ωj的fft值。
之后,违禁物品检测装置将由计算得到的违禁物品的标准品在各频点的吸收系数组成该标准品的吸收谱数据。该标准品的吸收谱数据可以记为2m行两列的数组(ω,α),α列中的数据为违禁物品的标准品在各频点的吸收系数。
S204:从违禁物品的标准品的吸收谱数据中确定出该标准品的特征峰数据后,将该标准品的特征峰数据作为提取出的该种违禁物品的标准品的特征数据。
具体地,违禁物品检测装置从步骤S203中得到的违禁物品的标准品的吸收谱数据中,确定出该标准品的吸收谱的各波峰的数据;并从各波峰的数据中确定出特征峰数据后,将该标准品的特征峰数据作为提取出的该种违禁物品的标准品的特征数据。特征峰数据具体包括:特征峰所对应的频点、违禁物品的标准品在该频点的吸收系数。
更优的,违禁物品检测装置还可以将上述步骤S203中得到的违禁物品的标准品的吸收谱数据作为提取出的该种违禁物品的标准品的特征数据。
之后,违禁物品检测装置可以将上述步骤S101-S102中的参考太赫兹时域波形数据、参考太赫兹频域波形数据,以及步骤S201-S204中提取出的各种违禁物品的标准品的特征数据,都存储到预先创建的数据库中。
本发明实施例的违禁物品检测方法的流程示意图,如图3所示,包括如下步骤:
S301:将被测物放置于太赫兹光路中,测量被测物的太赫兹时域波形数据。
具体地,被测物的携带者将被测物放置于太赫兹光路中,违禁物品检测装置周期性采集太赫兹波的信号强度值,得到被测物的太赫兹时域波形数据。
S302:对得到的被测物的太赫兹时域波形数据进行快速傅里叶变换,得到被测物的太赫兹频域波形数据。
具体地,违禁物品检测装置对步骤S301中得到的得到的被测物的太赫兹时域波形数据,根据与步骤S102中相同的2m个变换点,进行快速傅里叶变换,得到被测物的太赫兹频域波形数据。
被测物的太赫兹频域波形数据可以记为2m行两列的数组,一列中的数据为各频点,另一列中的数据为该被测物的太赫兹频域波形数据中的在各频点的fft值。
S303:根据被测物的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出被测物的吸收谱数据。
具体地,违禁物品检测装置对于步骤S302中得到的被测物的太赫兹频域波形数据中的每个频点,根据该被测物的太赫兹频域波形数据、以及参考太赫兹频域波形数据在该频率的fft值,计算出该被测物在该频点的吸收系数。将该被测物的太赫兹频域波形数据在各频点的fft值代替步骤S203中违禁物品的标准品的太赫兹频域波形数据在各频点的fft值后,计算该被测物在每个频点的吸收系数的具体方法,与步骤S203中计算违禁物品的标准品在每个频点的吸收系数的方法相同,此处不再赘述。
之后,违禁物品检测装置将由计算得到的被测物在各频点的吸收系数组成该被测物的吸收谱数据。
S304:将被测物的吸收谱数据与数据库中预先存储的各种违禁物品的标准品的特征峰数据进行比较,判断该被测物是否为其中一种违禁物品。
具体地,违禁物品检测装置对于每种违禁物品,从上述数据库中获取预先存储的该种违禁物品的标准品的特征峰数据,将该标准品的吸收谱数据的特征峰所对应的频点作为中心频点,将该中心频点两边设定范围的频段作为比对窗口;判断被测物的吸收谱数据在该比对窗口内是否存在波峰;若不存在,则判断出该被测物不是该种违禁物品。比对窗口可以根据历史数据、经验或者实际要求等因素确定,不再赘述;用以消除不同时刻、不同测量条件下数据的漂移引起的误差或影响。
违禁物品检测装置若判断出被测物的吸收谱数据在该比对窗口内存在波峰,则从上述数据库中获取该种违禁物品的标准品的吸收谱数据,并将该被测物的吸收谱数据与该种违禁物品的标准品的吸收谱数据进行比较,判断该被测物是否为该种违禁物品,具体方法将在后续进行详细介绍。
下面详细介绍上述步骤S304中,违禁物品检测装置将被测物的吸收谱数据与一种违禁物品的标准品的吸收谱数据进行比较,判断该被测物是否为该种违禁物品的方法,其流程示意图如图4所示,包括如下步骤:
S401:统计被测物的吸收谱数据的波峰总数。
具体地,违禁物品检测装置统计被测物的吸收谱数据的波峰总数;并统计该种违禁物品的标准品的吸收谱数据的波峰总数。
S402:对于该种违禁物品的标准品的吸收谱数据的每个波峰,在该波峰所对应的频点为中心的设定窗口频段内,若判断被测物的吸收谱数据存在波峰,则将该频点记为一个相似点,并计算两个波峰的峰值差,作为该相似点处的峰值差。
具体地,违禁物品检测装置从该种违禁物品的标准品的吸收谱数据中,确定出每个波峰及其所对应的频点;对于确定出的每个波峰,在该波峰所对应的频点为中心的设定窗口频段内,判断被测物的吸收谱数据是否存在波峰;若存在波峰,则将该频点记为一个相似点,并计算两个波峰的峰值差,作为该相似点处的峰值差;若不存在波峰,则忽略该频点相关的数据。例如,如图5所示,横坐标表示频率,纵坐标表示吸收系数;标识为SAMPLE的曲线表示被测物的吸收谱,标识为rdx的曲线表示违禁物品的标准品的吸收谱;对于违禁物品的标准品的吸收谱中的一个波峰,以该波峰对应的频点为中心的设定窗口频段内,被测物的吸收谱存在一个波峰。
之后,违禁物品检测装置计算出相似点的个数。
S403:根据步骤S402中得到的相似点的个数、每个相似点处的峰值差、步骤S401中被测物的吸收谱数据的波峰总数,以及该种违禁物品的标准品的吸收谱数据的波峰总数,计算被测物与该标准品之间的相似度。
具体地,违禁物品检测装置根据步骤S402中得到的相似点的个数nsim、每个相似点处的峰值差、步骤S401中被测物的吸收谱数据的波峰总数nsam、该种违禁物品的标准品的吸收谱数据的波峰总数nstd,以及如下公式2计算出被测物与该标准品之间的相似度S:
公式2中,i为不小于1且不大于nsim的自然数,pi表示第i个相似点处的峰值差,ω1、ω2、ω3分别表示三个权重系数。
S404:若计算得到的相似度大于设定的相似度阈值,则判断被测物是该种违禁物品;否则,判断被测物不是该种违禁物品。
具体地,违禁物品检测装置判断步骤S403中计算得到的相似度S是否大于预先设定的相似度阈值:若是,判断被测物是该种违禁物品,并将该种违禁物品原先存储的名称、所属种类以及计算出的相似度进行显示S,供检查人员参考;否则,判断被测物不是该种违禁物品。
基于上述违禁物品检测方法,本发明实施例的违禁物品检测装置的内部结构的框架示意图,如图6所示,包括:时域数据检测模块601、频域数据计算模块602、吸收谱计算模块603和特征比对模块604。
时域数据检测模块601用于测量太赫兹光路中经过被测物的太赫兹时域波形数据。
频域数据计算模块602用于对时域数据检测模块601测量得到的太赫兹时域波形数据进行快速傅里叶变换,得到被测物的太赫兹频域波形数据。
吸收谱计算模块603用于根据频域数据计算模块602得到的被测物的太赫兹频域波形数据,以及参考太赫兹频域波形数据,计算出被测物的吸收谱数据。
特征比对模块604用于将吸收谱计算模块603计算出的被测物的吸收谱数据与数据库中预先存储的各种违禁物品的标准品的特征峰数据进行比较,判断被测物是否为其中一种违禁物品;其中,违禁物品的标准品的特征峰数据具体是从该标准品的吸收谱数据的各波峰中选取的波峰的数据。
更优的,如图6所示,违禁物品检测装置还包括:参考时域波形测量模块605和参考频域波形计算模块606。
参考时域波形测量模块605用于在太赫兹光路中,测量没有被测物的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据。
参考频域波形计算模块606用于将参考时域波形测量模块605测量得到的参考太赫兹时域波形数据进行快速傅里叶变换,得到参考太赫兹频域波形数据进行存储。
更优的,如图6所示,违禁物品检测装置还包括:标准品数据库建立模块607。
标准品数据库建立模块607用于针对每种违禁物品,将该种违禁物品的标准品放置于太赫兹光路中,测量该标准品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据;根据该标准品的太赫兹频域波形数据以及参考频域波形计算模块606得到的参考太赫兹频域波形数据,计算出该标准品的吸收谱数据;从该标准品的吸收谱数据中确定出该标准品的特征峰数据后,将该种违禁物品的标准品的特征峰数据和吸收谱数据存入到数据库中。
上述时域数据检测模块601、频域数据计算模块602、吸收谱计算模块603、特征比对模块604、参考时域波形测量模块605、参考频域波形计算模块606和标准品数据库建立模块607的功能的具体实现方法,可以参考上述如图1~图4所示的方法流程步骤的具体内容,此处不再赘述。
本发明的技术方案中,将被测物吸收太赫兹波的吸收谱数据,与每种违禁物品的标准品对太赫兹波的吸收谱数据中的特征数据进行对比,可以根据对比结果更为准确地判断出被测物是否为违禁物品。
而且,还可以根据对比结果计算出相似度的值,并根据相似度的值判断被测物是否为违禁物品,进一步提高了判断的准确性。
此外,被测物吸收太赫兹波的吸收谱数据、以及每种违禁物品的标准品的特征数据都是从各自时域波形数据转换出的频域数据;与时域波形数据相比,数据量大大减少,使得对比时的计算量大大降低;从而可以加快对比速率,提高判断效率。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于计算机可读取存储介质中,如:ROM/RAM、磁碟、光盘等。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (9)
1.一种违禁物品检测方法,其特征在于,包括:
将被测物放置于太赫兹光路中,测量所述被测物的太赫兹时域波形数据;
对得到的太赫兹时域波形数据进行快速傅里叶变换,得到所述被测物的太赫兹频域波形数据;
根据所述被测物的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出所述被测物的吸收谱数据;
对于数据库中预先存储的一种违禁物品的标准品的特征峰数据,将该标准品的吸收谱数据的特征峰所对应的频点作为中心频点,将该中心频点两边设定范围的频段作为比对窗口;判断所述被测物的吸收谱数据在所述比对窗口内是否存在波峰;若不存在,则判断所述被测物不是该种违禁物品;其中,所述违禁物品的标准品的特征峰数据具体是从该标准品的吸收谱数据的各波峰中选取的波峰的数据;以及
所述数据库中还存储有该种违禁物品的标准品的吸收谱数据;以及
在所述判断所述被测物的吸收谱数据在所述比对窗口内是否存在波峰后,还包括:
若判断存在波峰,则:
统计所述被测物的吸收谱数据的波峰总数;
对于该标准品的吸收谱数据的每个波峰,在该波峰所对应的频点为中心的设定窗口频段内,若判断所述被测物的吸收谱数据存在波峰,则将该频点记为一个相似点,并计算两个波峰的峰值差,作为该相似点处的峰值差;
根据得到的相似点的个数、每个相似点处的峰值差、所述被测物的吸收谱数据的波峰总数,以及该标准品的吸收谱数据的波峰总数,计算所述被测物与该标准品之间的相似度;
若计算得到的相似度大于设定的相似度阈值,则判断所述被测物是该种违禁物品;否则,判断所述被测物不是该种违禁物品。
2.如权利要求1所述的方法,其特征在于,所述参考太赫兹频域波形数据是根据如下方法预先计算得到的:
在所述太赫兹光路中,测量没有所述被测物的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;
将所述参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据。
3.如权利要求2所述的方法,其特征在于,所述根据所述被测物的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出所述被测物的吸收谱数据,具体包括:
针对所述被测物的太赫兹频域波形数据中的每个频点,根据所述被测物的太赫兹频域波形数据、以及参考太赫兹频域波形数据在该频点的fft值,计算所述被测物在该频点的吸收系数;
由计算得到的所述被测物在各频点的吸收系数组成所述被测物的吸收谱数据。
4.如权利要求3所述的方法,其特征在于,所述根据得到的相似点的个数、每个相似点处的峰值差、所述被测物的吸收谱数据的波峰总数,以及该标准品的吸收谱数据的波峰总数,计算所述被测物与该标准品之间的相似度,具体为:
根据如下公式2计算得到所述被测物与该标准品之间的相似度S:
其中,nsim、nsam、nstd分别表示所述相似点的个数、所述被测物的吸收谱数据的波峰总数、该标准品的吸收谱数据的波峰总数,i为不小于1且不大于nsim的自然数,pi表示第i个相似点处的峰值差,ω1、ω2、ω3分别表示三个权重系数。
5.如权利要求4所述的方法,其特征在于,所述数据库根据如下方法预先建立:
针对每种违禁物品,提取该种违禁物品的标准品基于太赫兹频段的特征数据:
将该种违禁物品的标准品放置于太赫兹光路中,测量该标准品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据;根据该标准品的太赫兹频域波形数据以及所述参考太赫兹频域波形数据,计算出该标准品的吸收谱数据;从该标准品的吸收谱数据中确定出该标准品的特征峰数据后,将该标准品的特征峰数据作为提取出的该种违禁物品的标准品的特征数据;
将提取出的各种违禁物品的标准品的特征数据存储到所述数据库中。
6.如权利要求5所述的方法,其特征在于,在所述将提取出的各种违禁物品的标准品的特征数据存储到所述数据库前,还包括:
对于每种违禁物品,将该种违禁物品的标准品的吸收谱数据也作为该种违禁物品的标准品的特征数据。
7.一种违禁物品检测装置,其特征在于,包括:
时域数据检测模块,用于测量太赫兹光路中经过被测物的太赫兹时域波形数据;
频域数据计算模块,用于对所述时域数据检测模块测量得到的太赫兹时域波形数据进行快速傅里叶变换,得到所述被测物的太赫兹频域波形数据;
吸收谱计算模块,用于根据所述频域数据计算模块得到的所述被测物的太赫兹频域波形数据,以及参考太赫兹频域波形数据,计算出所述被测物的吸收谱数据;
特征比对模块,用于将所述吸收谱计算模块计算出的被测物的吸收谱数据与数据库中预先存储的各种违禁物品的标准品的特征峰数据进行比较,判断所述被测物是否为其中一种违禁物品;其中,所述违禁物品的标准品的特征峰数据具体是从该标准品的吸收谱数据的各波峰中选取的波峰的数据;其中,
将被测物的吸收谱数据与数据库中预先存储的各种违禁物品的标准品的特征峰数据进行比较,判断所述被测物是否为其中一种违禁物品,具体包括:
对于所述数据库中预先存储的一种违禁物品的标准品的特征峰数据,将该标准品的吸收谱数据的特征峰所对应的频点作为中心频点,将该中心频点两边设定范围的频段作为比对窗口;判断所述被测物的吸收谱数据在所述比对窗口内是否存在波峰;若不存在,则判断所述被测物不是该种违禁物品;
所述数据库中还存储有该种违禁物品的标准品的吸收谱数据;以及
在所述判断所述被测物的吸收谱数据在所述比对窗口内是否存在波峰后,还包括:
若判断存在波峰,则:
统计所述被测物的吸收谱数据的波峰总数;
对于该标准品的吸收谱数据的每个波峰,在该波峰所对应的频点为中心的设定窗口频段内,若判断所述被测物的吸收谱数据存在波峰,则将该频点记为一个相似点,并计算两个波峰的峰值差,作为该相似点处的峰值差;
根据得到的相似点的个数、每个相似点处的峰值差、所述被测物的吸收谱数据的波峰总数,以及该标准品的吸收谱数据的波峰总数,计算所述被测物与该标准品之间的相似度;
若计算得到的相似度大于设定的相似度阈值,则判断所述被测物是该种违禁物品;否则,判断所述被测物不是该种违禁物品。
8.如权利要求7所述的装置,其特征在于,还包括:
参考时域波形测量模块,用于在所述太赫兹光路中,测量没有所述被测物的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;
参考频域波形计算模块,用于将所述参考时域波形测量模块测量得到的参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据进行存储。
9.如权利要求8所述的装置,其特征在于,还包括:
标准品数据库建立模块,用于针对每种违禁物品,将该种违禁物品的标准品放置于太赫兹光路中,测量该标准品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据;根据该标准品的太赫兹频域波形数据以及所述参考频域波形计算模块得到的参考太赫兹频域波形数据,计算出该标准品的吸收谱数据;从该标准品的吸收谱数据中确定出该标准品的特征峰数据后,将该种违禁物品的标准品的特征峰数据和吸收谱数据存入到所述数据库中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410429717.8A CN104251838B (zh) | 2014-08-27 | 违禁物品检测方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410429717.8A CN104251838B (zh) | 2014-08-27 | 违禁物品检测方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104251838A CN104251838A (zh) | 2014-12-31 |
CN104251838B true CN104251838B (zh) | 2017-01-04 |
Family
ID=
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1818635A (zh) * | 2006-03-15 | 2006-08-16 | 首都师范大学 | 十二种毒品的THz检测方法及其指纹谱 |
GB2435509A (en) * | 2004-08-26 | 2007-08-29 | Teraview Ltd | Terahertz pulsed spectroscopy method |
CN101620065A (zh) * | 2008-07-01 | 2010-01-06 | 中国科学院上海应用物理研究所 | 麻醉药品的太赫兹检测方法 |
CN102353646A (zh) * | 2011-09-30 | 2012-02-15 | 中国石油大学(北京) | 基于太赫兹时域光谱的干酪根无损检测分析方法 |
CN102749297A (zh) * | 2012-07-10 | 2012-10-24 | 中国计量学院 | 一种基于太赫兹理论模拟光谱的农药鉴别方法 |
CN103969214A (zh) * | 2014-04-17 | 2014-08-06 | 首都师范大学 | 利用太赫兹频段红外光谱技术检测粮食中农药含量的方法 |
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2435509A (en) * | 2004-08-26 | 2007-08-29 | Teraview Ltd | Terahertz pulsed spectroscopy method |
CN1818635A (zh) * | 2006-03-15 | 2006-08-16 | 首都师范大学 | 十二种毒品的THz检测方法及其指纹谱 |
CN101620065A (zh) * | 2008-07-01 | 2010-01-06 | 中国科学院上海应用物理研究所 | 麻醉药品的太赫兹检测方法 |
CN102353646A (zh) * | 2011-09-30 | 2012-02-15 | 中国石油大学(北京) | 基于太赫兹时域光谱的干酪根无损检测分析方法 |
CN102749297A (zh) * | 2012-07-10 | 2012-10-24 | 中国计量学院 | 一种基于太赫兹理论模拟光谱的农药鉴别方法 |
CN103969214A (zh) * | 2014-04-17 | 2014-08-06 | 首都师范大学 | 利用太赫兹频段红外光谱技术检测粮食中农药含量的方法 |
Non-Patent Citations (2)
Title |
---|
基于THz光谱数据库的爆炸物分类识别研究;赵晶晶;《中国优秀硕士学位论文全文数据库 信息科技辑》;20091015(第10期);详见正文第3页倒数第2段,第9-10页"2.2.1 时域数据到频域数据的转换",第42页,第43页图4-3, * |
太赫兹光谱技术在毒品检测中的应用研究;和挺 等;《光谱学与光谱分析》;20130930;第33卷(第9期);2348-2353 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104316487B (zh) | 违禁物品检测方法和装置 | |
Shemer et al. | An experimental study of spatial evolution of statistical parameters in a unidirectional narrow‐banded random wavefield | |
Thomas et al. | Upper bound of 0.28 eV on neutrino masses from the largest photometric redshift survey | |
Soulsby | Selecting record length and digitization rate for near-bed turbulence measurements | |
US11802985B2 (en) | Method and system for analyzing filling for karst reservoir based on spectrum decomposition and machine learning | |
Ellwood et al. | A climate-driven model and development of a floating point time scale for the entire Middle Devonian Givetian Stage: a test using magnetostratigraphy susceptibility as a climate proxy | |
Star et al. | Comparison of ground motions from hybrid simulations to NGA prediction equations | |
Li et al. | Asteroseismology of 36 Kepler subgiants–I. Oscillation frequencies, linewidths, and amplitudes | |
Lumpkin | Global characteristics of coherent vortices from surface drifter trajectories | |
Bouruet-Aubertot et al. | Contrasted turbulence intensities in the Indonesian Throughflow: a challenge for parameterizing energy dissipation rate | |
Cava et al. | Interaction of submeso motions in the Antarctic stable boundary layer | |
CN107255831A (zh) | 一种叠前频散属性的提取方法 | |
CN106814393A (zh) | 一种地层品质因子q的估算方法 | |
CN107632010A (zh) | 一种结合激光诱导击穿光谱对钢铁样品的定量方法 | |
CN106644879A (zh) | 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置 | |
CN103678778A (zh) | 一种放射性物化探信息综合方法 | |
CN106772570A (zh) | 火成岩开启缝地震预测方法 | |
CN105046025A (zh) | 一种核磁共振多相流测量中各相分离的方法 | |
CN107167833A (zh) | 一种γ谱假峰甄别方法、存储介质和系统 | |
Cancelli et al. | Dimensionless criteria for the production‐dissipation equilibrium of scalar fluctuations and their implications for scalar similarity | |
CN105954803A (zh) | 叠后地震反演方法及装置 | |
Drew et al. | The evolution and use of discovery process models at the US Geological Survey | |
Li et al. | Seismic wave field anomaly identification of ultra-deep heterogeneous fractured-vuggy reservoirs: A case study in Tarim Basin, China | |
CN104392084A (zh) | 基于太赫兹频段的标准品特征数据库建立方法和装置 | |
CN107657119A (zh) | 提高伽马能谱数据品质的测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |