CN104316487A - 违禁物品检测方法和装置 - Google Patents
违禁物品检测方法和装置 Download PDFInfo
- Publication number
- CN104316487A CN104316487A CN201410482295.0A CN201410482295A CN104316487A CN 104316487 A CN104316487 A CN 104316487A CN 201410482295 A CN201410482295 A CN 201410482295A CN 104316487 A CN104316487 A CN 104316487A
- Authority
- CN
- China
- Prior art keywords
- data
- checked
- absorption spectra
- terahertz
- crest
- 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
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种违禁物品检测方法和装置,该方法包括:采集待检物品的吸收谱数据;确定出采集的吸收谱数据中的各波峰;针对每个波峰,根据该波峰对应的太赫兹吸收值、该波峰的相邻点的太赫兹吸收值,计算出该波峰的峰值明显度;将计算出的波峰明显度与预设的干扰阈值进行比较,根据比较结果判断该波峰是否为干扰点;利用滑动平均方法对判断出的待检物品的吸收谱数据中的干扰点进行数据修复,得到待检物品的修复后吸收谱数据;将待检物品的修复后吸收谱数据与违禁物品指纹谱库中的各种违禁物品的标准品的吸收谱数据进行比对,根据比对结果判断待检物品是否为其中一种违禁物品。应用本发明,可以提高检测结果的准确度。
Description
技术领域
本发明涉及光谱领域,尤其涉及一种违禁物品检测方法和装置。
背景技术
目前,对违禁物品的检测已经广泛地应用在机场、海关、车站、港口、核电站、政府机关、邮检中心、法院、大使馆等场所,用以检测出人员随身携带、或者行李中的违禁物品,从而保障安全。
现有的一种违禁物品检测方法为:利用X光照射待检物品,根据X光透射的原理,显示出待检物品的形状轮廓,并对不同密度的待检物品分别以不同的颜色进行显示,例如有机物显示为橙色,无机物显示为蓝色,混合物显示为绿色。检测人员根据现有的违禁物品检测装置显示的形状轮廓和颜色,可以较为容易地识别出枪支、刀具等金属类的违禁物品;但是对于化学类的违禁物品(例如炸药、毒品),该违禁物品的形状通常很难预测,显示的颜色又通常与日用品显示的颜色相同或相近,因此不容易识别出化学类的违禁物品。
现有的另一种违禁物品检测方法为:利用X光照射待检物品,根据康普顿散射原理,显示出待检物品的形状轮廓,并以不同的灰度分别显示低原子序数、高原子序数的待检物品。检测人员根据现有的违禁物品检测装置显示的形状轮廓和灰度,可以较为容易地识别出枪支、刀具等金属类的违禁物品。但是对于化学类的违禁物品(例如炸药、毒品),该违禁物品的形状通常很难预测;事实上,无论是违禁物品还是日用品都可以只包含低原子序数的成分、只包含高原子序数的成分,或者既包含低原子序数的成分又包含高原子序数的成分;导致该违禁物品显示的灰度通常与日用品显示的灰度相同或相近;造成检测人员根据该违禁物品检测方法不容易识别出化学类的违禁物品。
综上,利用现有的违禁物品检测方法,判断待检物品是否为违禁物品的准确性较低;因此,有必要提供一种准确性更高的违禁物品检测方法和装置。
发明内容
针对上述现有技术存在的缺陷,本发明提供了一种违禁物品检测方法和装置,用以提高检测结果的准确性。
根据本发明的一个方面,提供了一种违禁物品检测方法,包括:
将待检物品放置于预设的太赫兹光路中,采集所述待检物品的吸收谱数据;
确定出采集的所述待检物品的吸收谱数据中的各波峰;
针对确定出的每个波峰,根据该波峰对应的太赫兹吸收值、该波峰的相邻点的太赫兹吸收值,计算出该波峰的峰值明显度;并将计算出的波峰明显度与预设的干扰阈值进行比较,根据比较结果判断该波峰是否为干扰点;
利用预设的滑动平均方法对判断出的所述待检物品的吸收谱数据中的干扰点进行数据修复,得到所述待检物品的修复后吸收谱数据;
将所述待检物品的修复后吸收谱数据与违禁物品指纹谱库中的各种违禁物品的标准品的吸收谱数据进行比对,根据比对结果判断所述待检物品是否为其中一种违禁物品。
较佳地,所述将待检物品放置于预设的太赫兹光路中,采集所述待检物品的吸收谱数据,具体包括:
将待检物品放置于预设的太赫兹光路中,测量所述待检物品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到所述待检物品的太赫兹频域波形数据;根据所述待检物品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出所述待检物品的吸收谱数据;
其中,所述参考太赫兹频域波形数据是预先计算得到的:在所述太赫兹光路中,测量没有所述待检物品的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;将所述参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据。
较佳地,所述违禁物品指纹谱库中的各种违禁物品的标准品的吸收谱数据是预先存储的:
针对每种违禁物品,将该种违禁物品的标准品放置于所述预设的太赫兹光路中,测量该标准品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据;根据该标准品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出该标准品的吸收谱数据;
将提取出的各种违禁物品的标准品的吸收谱数据存储到所述违禁物品指纹谱库中;
其中,所述参考太赫兹频域波形数据是预先计算得到的:在所述太赫兹光路中,测量没有所述标准品的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;将所述参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据。
较佳地,所述将所述待检物品的修复后吸收谱数据与违禁物品指纹谱库中的各种违禁物品的标准品的吸收谱数据进行比对,根据比对结果判断所述待检物品是否为其中一种违禁物品,具体包括:
确定出所述待检物品的修复后吸收谱数据中的各波峰;
针对确定出的每个波峰,判断所述修复后吸收谱数据中在该波峰所对应的频点为中心的设定窗口频段内是否存在波峰,若存在,则将该频点记为一个相似点,并计算两个波峰的峰值差,作为该相似点处的峰值差;
针对违禁物品指纹谱库中的每种违禁物品的标准品,根据得到的所述待检物品的修复后吸收谱数据中的波峰总数、相似点的个数、每个相似点处的峰值差、以及该标准品的吸收谱数据的波峰总数,计算所述待检物品与该标准品之间的相似度;
若计算得到的相似度大于设定的相似度阈值,则判断所述待检物品是该种违禁物品;否则,判断所述待检物品不是该种违禁物品。
较佳地,所述根据得到的所述待检物品的修复后吸收谱数据中的波峰总数、相似点的个数、每个相似点处的峰值差、以及该标准品的吸收谱数据的波峰总数,计算所述待检物品与该标准品之间的相似度,具体包括:
根据如下公式2计算得到所述待测物品与该标准品之间的相似度S:
其中,nsim、nsam、nstd分别表示所述相似点的个数、所述待测物品的吸收谱数据的波峰总数、该标准品的吸收谱数据的波峰总数,i为不小于1且不大于nsim的自然数,pi表示第i个相似点处的峰值差,ω1、ω2、ω3分别表示三个预设的权重系数。
较佳地,所述针对确定出的每个波峰,判断所述修复后吸收谱数据中在该波峰所对应的频点为中心的设定窗口频段内是否存在波峰之前,还包括:
从所述待检物品的修复后吸收谱数据中确定出所述待检物品的修复后吸收谱数据中的各波峰;
针对所述违禁物品指纹谱库中的每种违禁物品的标准品,从该标准品的吸收谱数据中确定出该标准品的特征峰数据,并将该标准品的吸收谱数据的特征峰所对应的频点作为中心频点,将该中心频点两边设定范围的频段作为比对窗口;判断所述待检物品的修复后吸收谱数据在所述比对窗口内是否存在波峰;若不存在,则判断所述待检物品不是该种违禁物品;
其中,该标准品的特征峰数据具体为:该标准品的吸收谱数据的各波峰中太赫兹吸收值最高的波峰的数据,或者,针对该标准品的吸收谱数据的每个波峰,以该波峰为中心的设定频段,若该波峰的峰值最高,则将该波峰的数据作为该标准品的特征峰数据。
根据本发明的另一个方面,还提供了一种违禁物品检测装置,包括:
吸收谱数据采集模块,用于将待检物品放置于预设的太赫兹光路中,采集所述待检物品的吸收谱数据;
吸收谱分析模块,用于确定出所述吸收谱数据采集模块采集的所述待检物品的吸收谱数据中的各波峰;针对确定出的每个波峰,根据该波峰对应的太赫兹吸收值、该波峰的相邻点的太赫兹吸收值,计算出该波峰的峰值明显度;并将计算出的波峰明显度与预设的干扰阈值进行比较,根据比较结果判断该峰值点是否为干扰点;
吸收谱去干扰模块,用于利用预设的滑动平均方法对所述吸收谱分析模块判断出的所述待检物品的吸收谱数据中的干扰点进行数据修复,得到所述待检物品的修复后吸收谱数据并输出;
违禁物品检测模块,用于将所述吸收谱去干扰模块输出的所述待检物品的修复后吸收谱数据与违禁物品指纹谱库中的各种违禁物品的标准品的吸收谱数据进行比较,根据比较结果判断所述待检物品是否为其中一种违禁物品。
较佳地,所述吸收谱数据采集模块具体用于将待检物品放置于预设的太赫兹光路中,测量所述待检物品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到所述待检物品的太赫兹频域波形数据;根据所述待检物品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出所述待检物品的吸收谱数据;
其中,所述参考太赫兹频域波形数据是预先计算得到的:在所述太赫兹光路中,测量没有所述待检物品的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;将所述参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据。
较佳地,所述违禁物品检测装置还包括:
违禁物品指纹谱库构建模块,用于针对每种违禁物品,将该种违禁物品的标准品放置于太赫兹光路中,测量该标准品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据;根据该标准品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出该标准品的吸收谱数据;将提取出的各种违禁物品的标准品的吸收谱数据存储到所述违禁物品指纹谱库中;其中,所述参考太赫兹频域波形数据是预先计算得到的:在所述太赫兹光路中,测量没有所述标准品的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;将所述参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据。
较佳地,所述违禁物品检测模块具体用于确定出所述吸收谱去干扰模块输出的所述待检物品的修复后吸收谱数据中的各波峰;针对确定出的每个波峰,判断所述修复后吸收谱数据中在该波峰所对应的频点为中心的设定窗口频段内是否存在波峰,若存在,则将该频点记为一个相似点,并计算两个波峰的峰值差,作为该相似点处的峰值差;针对违禁物品指纹谱库中的每种违禁物品的标准品,根据得到的所述待检物品的修复后吸收谱数据中的波峰总数、相似点的个数、每个相似点处的峰值差、以及该标准品的吸收谱数据的波峰总数,计算所述待检物品与该标准品之间的相似度;若计算得到的相似度大于设定的相似度阈值,则判断所述待检物品是该种违禁物品;否则,判断所述待检物品不是该种违禁物品。
本发明的技术方案中,将待检物品对太赫兹波的吸收谱数据与各种违禁物品的标准品的吸收谱数据进行比对,可以更为准确地检测待检物品是否为违禁物品。进一步地,在将待检物品的吸收谱数据与各种违禁物品的标准品的吸收谱数据进行比对之前,对待检物品的吸收谱数据进行干扰消除,可以提高吸收谱数据的纯净度,有助于提高检测结果的准确度。
附图说明
图1为本发明实施例的参考太赫兹频域波形数据获取方法的流程示意图;
图2为本发明实施例的违禁物品检测方法的流程示意图;
图3为本发明实施例的待检物品的修复后吸收谱数据与一种违禁物品的标准品的吸收谱数据进行比对的实例示意图;
图4为本发明实施例的违禁物品检测装置的内部结构示意图。
具体实施方式
以下将结合附图对本发明的技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施例,都属于本发明所保护的范围。
本发明的发明人发现,在同一太赫兹光路中,不同的违禁物品对太赫兹(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所示,具体可以通过如下方式来构建:
针对每种违禁物品,将该种违禁物品的标准品放置于预设的太赫兹光路中,测量该标准品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据;根据该标准品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出该标准品的吸收谱数据;将提取出的各种违禁物品的标准品的吸收谱数据存储到违禁物品指纹谱库中。
基于上述参考太赫兹频域波形数据和违禁物品指纹谱库,本发明实施例提供的违禁物品检测方法,其流程如图2所示,具体可以包括如下步骤:
S201:将待检物品放置于预设的太赫兹光路中,采集待检物品的吸收谱数据。
具体地,可以将待检物品放置于预设的太赫兹光路中,测量待检物品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到待检物品的太赫兹频域波形数据;根据待检物品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出待检物品的吸收谱数据。
对于待检物品的太赫兹频域波形数据中的每个频点ωj,根据该待检物品的太赫兹频域波形数据中的在该频点的fft值、参考太赫兹频域波形数据在该频点的fft值,以及如下公式1计算该待检物品在该频点的吸收系数αj:
公式1中,j为0到2m-1之间的非负整数;fftj(sam,2m)、fftj(ref,2m)分别表示该待检物品的太赫兹频域波形数据、参考太赫兹频域波形数据在频点ωj的fft值。
之后,将由计算得到的待检物品在各频点的吸收系数组成该待检物品的吸收谱数据。其中,待检物品的吸收谱数据可以记为2m行两列的数组(ω,α),α列中的数据为待检物品在各频点的吸收系数。
S202:确定出采集的待检物品的吸收谱数据中的各波峰。
具体地,可以针对采集的待检物品的吸收谱数据中的每个频点,将该频点的太赫兹吸收值、与该频点左右相邻的两个频点的太赫兹吸收值进行比较,若该频点的太赫兹吸收值都大于上述两个频点的太赫兹吸收值,则说明该频点是一个峰值点,可以将该频点的太赫兹吸收值确定为待检物品的吸收谱数据中的一个波峰所对应的太赫兹吸收值。
S203:针对确定出的每个波峰,根据该波峰对应的太赫兹吸收值、该波峰的相邻点的太赫兹吸收值,计算出该波峰的峰值明显度。
具体地,针对确定出的每个波峰,以该波峰所对应的频点为中心,确定出与该波峰左右相邻的四个频点及其对应的太赫兹吸收值,并计算出该四个频点所对应的太赫兹吸收值的平均值;将该波峰对应的太赫兹吸收值与计算出的平均值的比值作为该波峰的峰值明显度。
S204:针对确定出的每个波峰,将计算出的波峰明显度与预设的干扰阈值进行比较,根据比较结果判断该波峰所对应的频点是否为干扰点。
具体地,针对确定出的每个波峰,将通过步骤S203计算出的该波峰的峰值明显度与预设的干扰阈值进行比较,若比较结果为该波峰的峰值明显度小于预设的干扰阈值,则可以判断该波峰为干扰点。其中,干扰阈值具体由本领域技术人员根据经验进行预先设定的,例如,可以取值为1.2。
S205:利用预设的滑动平均方法对判断出的待检物品的吸收谱数据中的干扰点进行数据修复,得到待检物品的修复后吸收谱数据。
具体地,可以针对判断出的每个干扰点,按照设定的滑动窗口(比如5个频率步长),从采集的待检物品的吸收谱数据中选取出与该干扰点左右相邻的4个频点;计算出该干扰点所对应的太赫兹吸收值与选取出的各频点所对应的太赫兹吸收值的平均值,作为该干扰点的替换太赫兹吸收值。
利用计算出的各干扰点的替换太赫兹吸收值对待检物品的吸收谱数据中的各干扰点所对应的太赫兹吸收值进行替换,得到待检物品的修复后吸收谱数据。这样,仅仅对测量得到的待检物品的吸收谱数据中的干扰点进行数据修复,对其他数据并没有影响,提高了吸收谱的纯净度,进而提高最终比对的可靠性。
S206:将待检物品的修复后吸收谱数据与违禁物品指纹谱库中的各种违禁物品的标准品的吸收谱数据进行比对,根据比对结果判断待检物品是否为其中一种违禁物品。
具体地,可以确定出步骤S205所得到的待检物品的修复后吸收谱数据中的各波峰。继而,针对确定出的每个波峰,判断待检物品的修复后吸收谱数据中在该波峰所对应的频点为中心的设定窗口频段内是否存在波峰,若存在,则可以将该频点记为一个相似点,并计算两个波峰的峰值差,作为该相似点处的峰值差。
例如,如图3所示,横坐标表示频率,纵坐标表示吸收系数;标识为SAMPLE的曲线表示待检物品的修复后吸收谱数据,标识为rdx的曲线表示违禁物品的标准品的吸收谱数据;对于违禁物品的标准品的吸收谱中的一个波峰,以该波峰对应的频点为中心的设定窗口频段内,待检物品的修复后吸收谱数据中存在一个波峰。
继而,针对违禁物品指纹谱库中的每种违禁物品的标准品,根据得到的待检物品的修复后吸收谱数据中的波峰总数、相似点的个数、每个相似点处的峰值差、以及该标准品的吸收谱数据的波峰总数,计算待检物品与该标准品之间的相似度。具体地,根据待检物品的修复后吸收谱数据中波峰总数nsam、相似点的个数nsim、每个相似点处的峰值差、该种违禁物品的标准品的吸收谱数据的波峰总数nstd,以及如下公式2计算出被测物与该标准品之间的相似度S:
公式2中,i为不小于1且不大于nsim的自然数,pi表示第i个相似点处的峰值差,ω1、ω2、ω3分别表示三个预设的权重系数。
若计算得到的相似度大于设定的相似度阈值,则可以判断待检物品是该种违禁物品;否则,判断待检物品不是该种违禁物品。
更优地,为了加快检测速度,可以考虑在针对确定出的每个波峰,判断所述修复后吸收谱数据中在该波峰所对应的频点为中心的设定窗口频段内是否存在波峰之前,判断出待检物品是不是违禁物品。
具体地,可以从待检物品的修复后吸收谱数据中确定出待检物品的修复后吸收谱数据中的各波峰;继而,针对违禁物品指纹谱库中的每种违禁物品的标准品,从该标准品的吸收谱数据中确定出该标准品的特征峰数据,并将该标准品的吸收谱数据的特征峰所对应的频点作为中心频点,将该中心频点两边设定范围的频段作为比对窗口;判断待检物品的修复后吸收谱数据在比对窗口内是否存在波峰;若不存在,则判断待检物品不是该种违禁物品。
其中,该标准品的特征峰数据具体为:该标准品的吸收谱数据的各波峰中太赫兹吸收值最高的波峰的数据,或者,针对该标准品的吸收谱数据的每个波峰,以该波峰为中心的设定频段,若该波峰的峰值最高,则将该波峰的数据作为该标准品的特征峰数据。
基于上述违禁物品检测方法,本发明实施例还提供了一种违禁物品检测装置,其内部结构如图4所示,具体可以包括:吸收谱数据采集模块401、吸收谱分析模块402、吸收谱去干扰模块403、违禁物品检测模块404。
其中,吸收谱数据采集模块401用于将待检物品放置于预设的太赫兹光路中,采集所述待检物品的吸收谱数据。
具体地,吸收谱数据采集模块401具体用于将待检物品放置于预设的太赫兹光路中,测量待检物品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到待检物品的太赫兹频域波形数据;根据待检物品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出待检物品的吸收谱数据。其中,参考太赫兹频域波形数据是预先计算得到的:在太赫兹光路中,测量没有待检物品的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;将参考太赫兹时域波形数据进行快速傅里叶变换,得到参考太赫兹频域波形数据。
吸收谱分析模块402用于确定出吸收谱数据采集模块401采集的待检物品的吸收谱数据中的各波峰;针对确定出的每个波峰,根据该波峰对应的太赫兹吸收值、该波峰的相邻点的太赫兹吸收值,计算出该波峰的峰值明显度;并将计算出的波峰明显度与预设的干扰阈值进行比较,根据比较结果判断该峰值点是否为干扰点。
吸收谱去干扰模块403用于利用预设的滑动平均方法对吸收谱分析模块402判断出的待检物品的吸收谱数据中的干扰点进行数据修复,得到待检物品的修复后吸收谱数据并输出。
违禁物品检测模块404用于将吸收谱去干扰模块403输出的待检物品的修复后吸收谱数据与违禁物品指纹谱库中的各种违禁物品的标准品的吸收谱数据进行比较,根据比较结果判断待检物品是否为其中一种违禁物品。
具体地,违禁物品检测模块404具体用于确定出吸收谱去干扰模块输出的待检物品的修复后吸收谱数据中的各波峰;针对确定出的每个波峰,判断修复后吸收谱数据中在该波峰所对应的频点为中心的设定窗口频段内是否存在波峰,若存在,则将该频点记为一个相似点,并计算两个波峰的峰值差,作为该相似点处的峰值差;针对违禁物品指纹谱库中的每种违禁物品的标准品,根据得到的所述待检物品的修复后吸收谱数据中的波峰总数、相似点的个数、每个相似点处的峰值差、以及该标准品的吸收谱数据的波峰总数,计算所述待检物品与该标准品之间的相似度;若计算得到的相似度大于设定的相似度阈值,则判断待检物品是该种违禁物品;否则,判断待检物品不是该种违禁物品。
更优地,本发明实施例中的违禁物品检测装置还包括:违禁物品指纹谱库构建模块405。
违禁物品指纹谱库构建模块405用于针对每种违禁物品,将该种违禁物品的标准品放置于太赫兹光路中,测量该标准品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据;根据该标准品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出该标准品的吸收谱数据;将提取出的各种违禁物品的标准品的吸收谱数据存储到违禁物品指纹谱库中;其中,所参考太赫兹频域波形数据是预先计算得到的:在太赫兹光路中,测量没有标准品的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;将所参考太赫兹时域波形数据进行快速傅里叶变换,得到参考太赫兹频域波形数据。
本发明实施例中,关于违禁物品检测装置中各模块具体是如何进行待检物品的检测,可以参考上述步骤流程,在此不再详述。
本发明的技术方案中,将待检物品对太赫兹波的吸收谱数据与各种违禁物品的标准品的吸收谱数据进行比对,可以更为准确地检测待检物品是否为违禁物品。进一步地,在将待检物品的吸收谱数据与各种违禁物品的标准品的吸收谱数据进行比对之前,对待检物品的吸收谱数据进行干扰消除,可以提高吸收谱数据的纯净度,有助于提高检测结果的准确度。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (10)
1.一种违禁物品检测方法,其特征在于,包括:
将待检物品放置于预设的太赫兹光路中,采集所述待检物品的吸收谱数据;
确定出采集的所述待检物品的吸收谱数据中的各波峰;
针对确定出的每个波峰,根据该波峰对应的太赫兹吸收值、该波峰的相邻点的太赫兹吸收值,计算出该波峰的峰值明显度;并将计算出的波峰明显度与预设的干扰阈值进行比较,根据比较结果判断该波峰是否为干扰点;
利用预设的滑动平均方法对判断出的所述待检物品的吸收谱数据中的干扰点进行数据修复,得到所述待检物品的修复后吸收谱数据;
将所述待检物品的修复后吸收谱数据与违禁物品指纹谱库中的各种违禁物品的标准品的吸收谱数据进行比对,根据比对结果判断所述待检物品是否为其中一种违禁物品。
2.如权利要求1所述的方法,其特征在于,所述将待检物品放置于预设的太赫兹光路中,采集所述待检物品的吸收谱数据,具体包括:
将待检物品放置于预设的太赫兹光路中,测量所述待检物品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到所述待检物品的太赫兹频域波形数据;根据所述待检物品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出所述待检物品的吸收谱数据;
其中,所述参考太赫兹频域波形数据是预先计算得到的:在所述太赫兹光路中,测量没有所述待检物品的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;将所述参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据。
3.如权利要求1所述的方法,其特征在于,所述违禁物品指纹谱库中的各种违禁物品的标准品的吸收谱数据是预先存储的:
针对每种违禁物品,将该种违禁物品的标准品放置于所述预设的太赫兹光路中,测量该标准品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据;根据该标准品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出该标准品的吸收谱数据;
将提取出的各种违禁物品的标准品的吸收谱数据存储到所述违禁物品指纹谱库中;
其中,所述参考太赫兹频域波形数据是预先计算得到的:在所述太赫兹光路中,测量没有所述标准品的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;将所述参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据。
4.如权利要求1-3任一所述的方法,其特征在于,所述将所述待检物品的修复后吸收谱数据与违禁物品指纹谱库中的各种违禁物品的标准品的吸收谱数据进行比对,根据比对结果判断所述待检物品是否为其中一种违禁物品,具体包括:
确定出所述待检物品的修复后吸收谱数据中的各波峰;
针对确定出的每个波峰,判断所述修复后吸收谱数据中在该波峰所对应的频点为中心的设定窗口频段内是否存在波峰,若存在,则将该频点记为一个相似点,并计算两个波峰的峰值差,作为该相似点处的峰值差;
针对违禁物品指纹谱库中的每种违禁物品的标准品,根据得到的所述待检物品的修复后吸收谱数据中的波峰总数、相似点的个数、每个相似点处的峰值差、以及该标准品的吸收谱数据的波峰总数,计算所述待检物品与该标准品之间的相似度;
若计算得到的相似度大于设定的相似度阈值,则判断所述待检物品是该种违禁物品;否则,判断所述待检物品不是该种违禁物品。
5.如权利要求4所述的方法,其特征在于,所述根据得到的所述待检物品的修复后吸收谱数据中的波峰总数、相似点的个数、每个相似点处的峰值差、以及该标准品的吸收谱数据的波峰总数,计算所述待检物品与该标准品之间的相似度,具体包括:
根据如下公式2计算得到所述待测物品与该标准品之间的相似度S:
其中,nsim、nsam、nstd分别表示所述相似点的个数、所述待测物品的吸收谱数据的波峰总数、该标准品的吸收谱数据的波峰总数,i为不小于1且不大于nsim的自然数,pi表示第i个相似点处的峰值差,ω1、ω2、ω3分别表示三个预设的权重系数。
6.如权利要求4所述的方法,其特征在于,所述针对确定出的每个波峰,判断所述修复后吸收谱数据中在该波峰所对应的频点为中心的设定窗口频段内是否存在波峰之前,还包括:
从所述待检物品的修复后吸收谱数据中确定出所述待检物品的修复后吸收谱数据中的各波峰;
针对所述违禁物品指纹谱库中的每种违禁物品的标准品,从该标准品的吸收谱数据中确定出该标准品的特征峰数据,并将该标准品的吸收谱数据的特征峰所对应的频点作为中心频点,将该中心频点两边设定范围的频段作为比对窗口;判断所述待检物品的修复后吸收谱数据在所述比对窗口内是否存在波峰;若不存在,则判断所述待检物品不是该种违禁物品;
其中,该标准品的特征峰数据具体为:该标准品的吸收谱数据的各波峰中太赫兹吸收值最高的波峰的数据,或者,针对该标准品的吸收谱数据的每个波峰,以该波峰为中心的设定频段,若该波峰的峰值最高,则将该波峰的数据作为该标准品的特征峰数据。
7.一种违禁物品检测装置,其特征在于,包括:
吸收谱数据采集模块,用于将待检物品放置于预设的太赫兹光路中,采集所述待检物品的吸收谱数据;
吸收谱分析模块,用于确定出所述吸收谱数据采集模块采集的所述待检物品的吸收谱数据中的各波峰;针对确定出的每个波峰,根据该波峰对应的太赫兹吸收值、该波峰的相邻点的太赫兹吸收值,计算出该波峰的峰值明显度;并将计算出的波峰明显度与预设的干扰阈值进行比较,根据比较结果判断该峰值点是否为干扰点;
吸收谱去干扰模块,用于利用预设的滑动平均方法对所述吸收谱分析模块判断出的所述待检物品的吸收谱数据中的干扰点进行数据修复,得到所述待检物品的修复后吸收谱数据并输出;
违禁物品检测模块,用于将所述吸收谱去干扰模块输出的所述待检物品的修复后吸收谱数据与违禁物品指纹谱库中的各种违禁物品的标准品的吸收谱数据进行比较,根据比较结果判断所述待检物品是否为其中一种违禁物品。
8.如权利要求7所述的装置,其特征在于,
所述吸收谱数据采集模块具体用于将待检物品放置于预设的太赫兹光路中,测量所述待检物品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到所述待检物品的太赫兹频域波形数据;根据所述待检物品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出所述待检物品的吸收谱数据;
其中,所述参考太赫兹频域波形数据是预先计算得到的:在所述太赫兹光路中,测量没有所述待检物品的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;将所述参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据。
9.如权利要求7所述的装置,其特征在于,还包括:
违禁物品指纹谱库构建模块,用于针对每种违禁物品,将该种违禁物品的标准品放置于太赫兹光路中,测量该标准品的太赫兹时域波形数据;对得到的太赫兹时域波形数据进行快速傅里叶变换,得到该标准品的太赫兹频域波形数据;根据该标准品的太赫兹频域波形数据以及参考太赫兹频域波形数据,计算出该标准品的吸收谱数据;将提取出的各种违禁物品的标准品的吸收谱数据存储到所述违禁物品指纹谱库中;其中,所述参考太赫兹频域波形数据是预先计算得到的:在所述太赫兹光路中,测量没有所述标准品的情况下的太赫兹时域波形数据,作为参考太赫兹时域波形数据;将所述参考太赫兹时域波形数据进行快速傅里叶变换,得到所述参考太赫兹频域波形数据。
10.如权利要求7-9任一所述的装置,其特征在于,
所述违禁物品检测模块具体用于确定出所述吸收谱去干扰模块输出的所述待检物品的修复后吸收谱数据中的各波峰;针对确定出的每个波峰,判断所述修复后吸收谱数据中在该波峰所对应的频点为中心的设定窗口频段内是否存在波峰,若存在,则将该频点记为一个相似点,并计算两个波峰的峰值差,作为该相似点处的峰值差;针对违禁物品指纹谱库中的每种违禁物品的标准品,根据得到的所述待检物品的修复后吸收谱数据中的波峰总数、相似点的个数、每个相似点处的峰值差、以及该标准品的吸收谱数据的波峰总数,计算所述待检物品与该标准品之间的相似度;若计算得到的相似度大于设定的相似度阈值,则判断所述待检物品是该种违禁物品;否则,判断所述待检物品不是该种违禁物品。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410482295.0A CN104316487B (zh) | 2014-09-19 | 2014-09-19 | 违禁物品检测方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410482295.0A CN104316487B (zh) | 2014-09-19 | 2014-09-19 | 违禁物品检测方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104316487A true CN104316487A (zh) | 2015-01-28 |
CN104316487B CN104316487B (zh) | 2017-02-15 |
Family
ID=52371748
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410482295.0A Active CN104316487B (zh) | 2014-09-19 | 2014-09-19 | 违禁物品检测方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104316487B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106198482A (zh) * | 2015-05-04 | 2016-12-07 | 清华大学 | 基于拉曼光谱的检测保健品中是否添加有西药的方法 |
CN107328735A (zh) * | 2017-08-30 | 2017-11-07 | 浙江大学 | 基于太赫兹光谱技术的油菜种类鉴别方法 |
CN108154168A (zh) * | 2016-12-05 | 2018-06-12 | Sap欧洲公司 | 综合货物检验系统及方法 |
CN108226088A (zh) * | 2017-12-30 | 2018-06-29 | 雄安华讯方舟科技有限公司 | 一种药物检测方法及装置 |
CN108732636A (zh) * | 2018-05-21 | 2018-11-02 | 深圳市神飞电子科技有限公司 | 一种对违禁物品建立太赫兹物质指纹谱的方法及系统 |
CN109470647A (zh) * | 2019-01-20 | 2019-03-15 | 南京林业大学 | 一种水蒸气太赫兹吸收谱的测量方法 |
CN109883985A (zh) * | 2019-02-18 | 2019-06-14 | 西安科技大学 | 一种基于分子动力学理论的纳米声学效应分析方法 |
CN109946749A (zh) * | 2019-03-27 | 2019-06-28 | 嘉兴腓特烈太赫科技有限公司 | 检测体内毒品的装置 |
WO2019223364A1 (zh) * | 2018-05-25 | 2019-11-28 | 深圳市太赫兹科技创新研究院有限公司 | 基于太赫兹频段的电磁波检测冬虫夏草质量的方法 |
CN114460161A (zh) * | 2021-12-27 | 2022-05-10 | 中船重工安谱(湖北)仪器有限公司 | 一种基于离子迁移时间的痕量物质检测方法 |
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 | 首都师范大学 | 利用太赫兹频段红外光谱技术检测粮食中农药含量的方法 |
-
2014
- 2014-09-19 CN CN201410482295.0A patent/CN104316487B/zh active Active
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 (1)
Title |
---|
赵晶晶: "基于THz光谱数据库的爆炸物分类识别研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106198482B (zh) * | 2015-05-04 | 2019-07-05 | 清华大学 | 基于拉曼光谱的检测保健品中是否添加有西药的方法 |
CN106198482A (zh) * | 2015-05-04 | 2016-12-07 | 清华大学 | 基于拉曼光谱的检测保健品中是否添加有西药的方法 |
CN108154168A (zh) * | 2016-12-05 | 2018-06-12 | Sap欧洲公司 | 综合货物检验系统及方法 |
CN108154168B (zh) * | 2016-12-05 | 2022-07-08 | Sap欧洲公司 | 综合货物检验系统及方法 |
CN107328735A (zh) * | 2017-08-30 | 2017-11-07 | 浙江大学 | 基于太赫兹光谱技术的油菜种类鉴别方法 |
CN108226088B (zh) * | 2017-12-30 | 2019-09-17 | 雄安华讯方舟科技有限公司 | 一种药物检测方法及装置 |
CN108226088A (zh) * | 2017-12-30 | 2018-06-29 | 雄安华讯方舟科技有限公司 | 一种药物检测方法及装置 |
CN108732636A (zh) * | 2018-05-21 | 2018-11-02 | 深圳市神飞电子科技有限公司 | 一种对违禁物品建立太赫兹物质指纹谱的方法及系统 |
WO2019223364A1 (zh) * | 2018-05-25 | 2019-11-28 | 深圳市太赫兹科技创新研究院有限公司 | 基于太赫兹频段的电磁波检测冬虫夏草质量的方法 |
CN109470647A (zh) * | 2019-01-20 | 2019-03-15 | 南京林业大学 | 一种水蒸气太赫兹吸收谱的测量方法 |
CN109883985A (zh) * | 2019-02-18 | 2019-06-14 | 西安科技大学 | 一种基于分子动力学理论的纳米声学效应分析方法 |
CN109883985B (zh) * | 2019-02-18 | 2021-03-26 | 西安科技大学 | 一种基于分子动力学理论的纳米声学效应分析方法 |
CN109946749A (zh) * | 2019-03-27 | 2019-06-28 | 嘉兴腓特烈太赫科技有限公司 | 检测体内毒品的装置 |
CN114460161A (zh) * | 2021-12-27 | 2022-05-10 | 中船重工安谱(湖北)仪器有限公司 | 一种基于离子迁移时间的痕量物质检测方法 |
CN114460161B (zh) * | 2021-12-27 | 2022-09-27 | 中船重工安谱(湖北)仪器有限公司 | 一种基于离子迁移时间的痕量物质检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104316487B (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104316487A (zh) | 违禁物品检测方法和装置 | |
CN109581546B (zh) | 基于微波链路衰减及偏振信息的降雨类型识别方法 | |
CN106896306B (zh) | 一种gis振荡冲击耐压试验用信号去噪方法 | |
CN104266993A (zh) | 基于太赫兹频段的物品特征提取方法和装置 | |
ITPR20070059A1 (it) | Procedimento per rilevare, identificare e localizzare scariche parziali aventi luogo in un sito di scarica lungo un apparato elettrico | |
CN104736986A (zh) | 使用同位素比测量在人口稠密地区中进行气体泄漏检测和定位的方法 | |
CN104198898A (zh) | 基于脉冲序列分析的局部放电发展过程诊断方法 | |
CN108875706B (zh) | 基于滑动平均与能量归集的海洋结构时频分析方法 | |
FI107081B (fi) | Menetelmä ja järjestely osittaispurkauslähteiden lukumäärän selvittämiseksi | |
CN103941254A (zh) | 一种基于地质雷达的土壤物性类别识别方法和装置 | |
CN104749156A (zh) | 拉曼光谱检测方法 | |
CN105333841A (zh) | 基于反射型太赫兹时域光谱的金属表面粗糙度检测装置及方法 | |
CN102928713B (zh) | 一种磁场天线的本底噪声测量方法 | |
CN103115893A (zh) | 一种检测托卡马克钨第一壁灰尘沉积层成分及厚度的装置 | |
CN104076345B (zh) | 一种针对纯转动拉曼激光雷达测量温度的饱和修正方法 | |
CN108538011B (zh) | 一种激光雷达火灾检测方法 | |
CN104392084A (zh) | 基于太赫兹频段的标准品特征数据库建立方法和装置 | |
CN103512877B (zh) | 一种拉曼光谱物质检测快速样本查找方法 | |
CN104977602A (zh) | 一种地震数据采集施工的控制方法及装置 | |
CN103091259B (zh) | 在线检测托卡马克钨第一壁灰尘成分及厚度的太赫兹方法 | |
CN104864894B (zh) | 一种botdr布里渊散射谱寻峰方法 | |
CN104316488A (zh) | 违禁物品检测方法和装置 | |
CN104819382A (zh) | 一种用于光纤预警系统的自适应恒虚警率振源检测方法 | |
Kohlmann et al. | Estimation of Charge Transfer During Long Continuing Currents in Natural Downward Flashes Using Single‐Station E‐Field Measurements | |
CN105866552A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |