CN106353343B - 一种小像素透射多信息融合高速安检识别方法 - Google Patents
一种小像素透射多信息融合高速安检识别方法 Download PDFInfo
- Publication number
- CN106353343B CN106353343B CN201610843896.9A CN201610843896A CN106353343B CN 106353343 B CN106353343 B CN 106353343B CN 201610843896 A CN201610843896 A CN 201610843896A CN 106353343 B CN106353343 B CN 106353343B
- Authority
- CN
- China
- Prior art keywords
- ray
- signal
- energy
- formula
- 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
- 230000005540 biological transmission Effects 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000004927 fusion Effects 0.000 title claims abstract description 17
- 238000007689 inspection Methods 0.000 claims abstract description 48
- 239000000463 material Substances 0.000 claims abstract description 34
- 238000001514 detection method Methods 0.000 claims abstract description 32
- 239000000203 mixture Substances 0.000 claims abstract description 30
- 239000005416 organic matter Substances 0.000 claims abstract description 24
- 230000004888 barrier function Effects 0.000 claims abstract description 16
- 230000011218 segmentation Effects 0.000 claims abstract description 7
- 230000009977 dual effect Effects 0.000 claims description 34
- 239000000126 substance Substances 0.000 claims description 32
- 239000010936 titanium Substances 0.000 claims description 15
- 238000010521 absorption reaction Methods 0.000 claims description 11
- 238000013178 mathematical model Methods 0.000 claims description 11
- 230000035943 smell Effects 0.000 claims description 11
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 claims description 9
- 239000004411 aluminium Substances 0.000 claims description 9
- 229910052782 aluminium Inorganic materials 0.000 claims description 9
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000010354 integration Effects 0.000 claims description 9
- HBMJWWWQQXIZIP-UHFFFAOYSA-N silicon carbide Chemical compound [Si+]#[C-] HBMJWWWQQXIZIP-UHFFFAOYSA-N 0.000 claims description 9
- 229910010271 silicon carbide Inorganic materials 0.000 claims description 9
- 229910052719 titanium Inorganic materials 0.000 claims description 9
- 229910052720 vanadium Inorganic materials 0.000 claims description 9
- LEONUFNNVUYDNQ-UHFFFAOYSA-N vanadium atom Chemical compound [V] LEONUFNNVUYDNQ-UHFFFAOYSA-N 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 238000012544 monitoring process Methods 0.000 claims description 6
- 230000002285 radioactive effect Effects 0.000 claims description 6
- 229910000831 Steel Inorganic materials 0.000 claims description 5
- 239000010959 steel Substances 0.000 claims description 5
- 238000010586 diagram Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- PLXMOAALOJOTIY-FPTXNFDTSA-N Aesculin Natural products OC[C@@H]1[C@@H](O)[C@H](O)[C@@H](O)[C@H](O)[C@H]1Oc2cc3C=CC(=O)Oc3cc2O PLXMOAALOJOTIY-FPTXNFDTSA-N 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 230000000149 penetrating effect Effects 0.000 claims description 3
- 230000035515 penetration Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 208000033748 Device issues Diseases 0.000 claims description 2
- 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 claims 1
- -1 inorganic matter Substances 0.000 abstract description 5
- 230000006872 improvement Effects 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 4
- 239000002360 explosive Substances 0.000 description 3
- 229910052751 metal Inorganic materials 0.000 description 3
- 239000002184 metal Substances 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000000763 evoking effect Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 230000008092 positive effect Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/005—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using neutrons
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/06—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption
- G01N23/083—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption the radiation being X-rays
- G01N23/087—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption the radiation being X-rays using polyenergetic X-rays
-
- G01V5/20—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06K—GRAPHICAL DATA READING; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
- G06K17/00—Methods or arrangements for effecting co-operative working between equipments covered by two or more of main groups G06K1/00 - G06K15/00, e.g. automatic card files incorporating conveying and reading operations
-
- G06Q50/60—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/03—Investigating materials by wave or particle radiation by transmission
- G01N2223/04—Investigating materials by wave or particle radiation by transmission and measuring absorption
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/045—Investigating materials by wave or particle radiation combination of at least 2 measurements (transmission and scatter)
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/10—Different kinds of radiation or particles
- G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
- G01N2223/1016—X-ray
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/10—Different kinds of radiation or particles
- G01N2223/106—Different kinds of radiation or particles neutrons
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/30—Accessories, mechanical or electrical features
- G01N2223/33—Accessories, mechanical or electrical features scanning, i.e. relative motion for measurement of successive object-parts
- G01N2223/3307—Accessories, mechanical or electrical features scanning, i.e. relative motion for measurement of successive object-parts source and detector fixed; object moves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/421—Imaging digitised image, analysed in real time (recognition algorithms)
Abstract
本发明公开了一种小像素透射多信息融合高速安检系统及其识别方法；由交叉带式高速分拣机和控制交叉带式高速分拣机运行的控制装置组成；交叉带式高速分拣机包括分类装置、输送装置及分拣道口；控制装置包括全方位条形码扫描仪、X光检查装置、光障开关检测子系统、外部图像采集子系统、物理参数采集子系统和主控子系统通信联系，本发明提供了利用上述小像素透射多信息融合高速安检系统对邮件进行识别的方法，本系统实现了危险品检查仪设备状态的自检功能，具有将有机物、无机物、混合物进行科学区域分割的功能。
Description
技术领域
本发明涉及到邮政分拣设备系统技术领域，具体地说涉及一种小像素透射多信息融合高速安检系统及其识别方法。
背景技术
邮政分拣是国家公共安全的重点行业，百姓有意无意邮寄受其知识局限违禁品包裹投入到邮寄流程将对本人及他人生命和国家财产构成一种潜在的威胁，极易造成严重事故。更值得注意的是一些犯罪分子和恐怖分子为达到破坏社会安定的目的，有意邮寄违禁品来“考核”国家公共安全的建设水平更是一种严重的威胁，发生在国内外的恐怖事件已证明了这一点。随着国际反恐形势的不断深化，恐怖分子的破坏手段变得越来越高科技、越来越复杂化，使得国家职能机构应对不暇。采用高科技手段来提高安全检测的准确性、查堵破坏手段的复杂性变得越来越重要。
多年来广泛使用的安全检查设备，如磁力针、金属武器检测门、X射线检测仪等，能发现武器和普通炸药等危险品，在安全检查工作中发挥了重要作用。但是，这些设备受原技术条件所限，不尽如人意，漏检失误较多。再加之当今随着科学技术的发展，犯罪分子和恐怖分子也利用高新技术，制造新的武器、爆炸物等。如利用集成电路技术制造的高精密炸弹、塑性炸药及毒品等，上述传统检测手段就无能为力。
发明内容
本发明的目的在于针对上述现有技术的缺陷，提供一种小像素透射多信息融合高速安检系统及其识别方法。
为了实现上述目的，本发明的技术方案是：
一种小像素透射多信息融合高速安检系统，该小像素透射多信息融合高速安检系统由交叉带式高速分拣机和控制交叉带式高速分拣机运行的控制装置组成；
所述交叉带式高速分拣机包括输送装置、分类装置及分拣道口；所述控制装置按时间序列分别读取邮件的条码、邮件面单信号，接收和处理分拣信号；输送装置采用交叉带小车在固定轨道上运行的输送方式，运载着待分拣的邮件鱼贯通过控制装置；控制装置识别邮件条形码的信息，向分类装置发送邮件的分拣指令，并对邮件进行分类；分类装置根据控制装置发出的分拣指令，改变邮件在输送装置上运行的方向；交叉带小车按照分类装置的指令将邮件送达到指定的分拣道口处，小车托盘倾斜邮件滑下轨道台车进入分拣道口；
所述控制装置包括用于处理接收到的信息并发出指令信息的主控子系统、接收小车托盘的“车到”信号以及“高包裹”判别信号的光障开关检测子系统、采集邮件外部形状图像信息的外部图像采集子系统、扫描邮件并获取邮件地藉信息的全方位条形码扫描仪、识别邮件包裹的内部物质属性的X光检查装置、采集邮件放射性指标数据、挥发性指标数据的物理参数采集子系统；所述X光检查装置为可以发出和接收两种能级的双能量的X光检查装置，所述全方位条形码扫描仪、X光检查装置、光障开关检测子系统、外部图像采集子系统、物理参数采集子系统与主控子系统通信联系；
所述光障开关检测子系统监测“车到”和“高包裹”信号，并发送给主控子系统作为整个系统的周期控制启动信号；邮件鱼贯通过全方位条形码扫描仪，全方位条形码扫描仪按时间序列分别读取邮件的条形码、邮件面单信号，对载物小车身份识别；所述外部图像采集子系统监测采集邮件外部形状图像信息；物理参数采集子系统采用有效原子序数物质分类计算对被检目标实现定性识别，获取放射性指标数据和挥发性指标数据，实现ppM级气味识别；所述X光检查装置识别邮件包裹的内部物质属性，对被检目标进行线扫描透照图像获取；所述主控子系统对外部形状图像信号、原始X射线阵列信号进行数字化图形图像处理，并由物质分类实验建立材料分类区域及分类阈值，界定有机物、无机物、混合物，分辨出被检包裹内部物质的属性。
作为对上述技术方案的改进，所述主控子系统包括由若干个电脑组成的用于处理和显示X光图像的X光检查工作站，一台电脑组成的用于实现系统总体控制以及外部图像和其它物理参数采集的主控工作站，一台电脑组成的负载X 光透照数据采集、以及与主控工作站通信实现控制流程、并且控制X光检查工作站工作的X光透照控制工作站，所述主控工作站、X光检查工作站、X光透照控制工作站通过局域网相连接。
作为对上述技术方案的改进，所述X光检查装置包括X射线发生器和L形的X射线接收器，该X光检查装置安装在条形码扫描仪后部，X射线发生器位于输送装置的底部，L形的X射线接收器位于输送装置的顶部；所述X光检查装置与X光透照控制工作站连接。
作为对上述技术方案的改进，所述物理参数采集子系统包括γ中子检测仪、识别简单和复杂气味的、毫秒级响应时间的电子鼻气味传感器，该电子鼻气味传感器为PID光离子传感器，对待检测邮件进行气体检测；该γ中子检测仪检测邮件中是否含有射线。
作为对上述技术方案的改进，所述光障开关检测子系统包括光障开关；所述外部图像采集子系统包括工业数字相机，所述工业数字相机内置CCD传感器。
本发明并提供了一种利用上述小像素透射多信息融合高速安检系统对邮件进行识别的方法，其特征在于，该识别方法的步骤是：
S1、采用双能量的X光检查装置对邮件检测，得到被扫描邮件的平均原子数；
S2、将物质分类，设Z为物质的原子序数，根据定义，当1＜Z≤10时，为有机物，当Z＞20时，为无机物，当10＜Z≤20时为混合物；Z＝10和Z＝20 分别为有机物与混合物、无机物与混合物的边界；
S3、根据有效原子序数计算公式找到边界物质，该有效原子序数计算公式：
式中:wi为化合物质种类当量个数；为物质的原子序数。
根据上式，计算得到碳化硅的有效原子序数大约等于10，钛六铝四钒的有效原子序数大约等于20，这样我们就拿碳化硅和钛六铝四钒代替有机物与无机物的两个边界；
不同物质的质量吸收系数不同，X射线所能穿透的最大厚度也不同，有如下关系：
当双能量的X光检查装置穿透率为22mm钢板时，在100KeV时，
(μ/ρ)Fe＝0.3713(μ/ρ)Al＝0.1704(μ/ρ)Si＝0.1835
(μ/ρ)C＝0.1514(μ/ρ)V＝0.2877(μ/ρ)Ti＝0.21
ρFe＝7.874ρAl＝2.699ρsi＝2.330ρC＝1.7ρV＝6.11ρTi＝4.450
可以计算出碳化硅和钛六铝四钒的最大穿透厚度分别为108mm和36mm；
S4、将实验数据用曲线专家分别进行拟合得到两条边界曲线及其数学模型；
混合物与无机物的边界曲线数学模型为：
T(x2)＝-50587+26.8056·x2-0.00463·x2 2+3.7233·10-7·x2 3- 1.1318·10-11·x2 4 (3)
混合物与有机物的边界曲线数学模型为：
S(x4)＝-16403+9.6105·x4-0.00143·x4 2+1.06446·10-7·x4 3- 3.055·10-12·x4 4 (4)
S5、计算物质的有效原子序数和密度；
S6、双能量的X光检查装置的分类探测信号分析；
A1、透射信号，对于透射图像，探测信号与光电效应和散射都有关，有如下方程：
其中T(x，y)是在投影图中对应于坐标(x，y)探测到的透射信号，N(E) 是能量位于E和E+dE之间的X射线光子数，μt(x′，y′，z′，E)是沿着X射线方向的总的线吸收系数，Pd(E)是探测器的光子能量为E时的探测概率， E是光子能量；
A2、背散射信号；对于向后散射，探测信号B(x，y)估计为：
其中，与公式(11)中相同的符号其含义相同，其余的，B(x，y)是在投影图中(x，y)点的散射光子信号；Z′是发生散射时，笔状射线束沿射线方向进入物体的深度；σs(E)是散射的横截面积，是能量的函数；ρ(x′，y′，z′)是包裹中的物体沿X射线方向的密度(g/cm3)；E′是散射后的光子能量。射；
A3、向前散射信号的表达式与向后散射信号的表达式相似，此式中第二个积分项的积分限为从Z′到D；
除了F(x，y)表示在投影图中(x，y)位置的散射光子信号外，其它的符号与 (11)、(12)两式中的含义相同。
A4、提取R值，由双能量的X光检查装置透射图像的低、高能灰度值计算出R值。由散射图像的灰度值可以计算出L值；
双能量的X光检查装置使用将TH和TL相结合的方法，得到R值：
R称作双能量，TH0和TL0是在自由空间获得的高能、低能透射信号。
S7、考虑物质厚度对双能量透射图像的影响；
一定输入能量Ein下，总的透射信号是能量范围内透射信号的积分，
将透射信号流程改写为：
对于简单的物体，总的线吸收系数为μt(x′，y′，E)，厚度为t，于是方程(17)变为：
对于在图像中任意位置(x，y)的物体，从方程(18)中去除空间坐标，得到：
从方程(19)可知，总的透射信号强度取决于物体的厚度；
双能量的X光检查装置受到X射线源的多色谱线的影响；对于某一特定的物质，输入能量分别为EH，EL，于是有：
从方程(20)可知，R不仅与物质的原子序数Z，原子重量A，密度ρ有关，也于物质的厚度有关；
S8、克服厚度影响的数字化方法
对于给定的双能量的X光检查装置，在方程(20)中，入射的X射线光子数N(E)，X射线光子能量E，光子的探测概率Pd(E)都是已知的或是通过仿真可以得到的参数；
让K(E)＝N(E)EPd(E),在特定点的透射信号可写成：
被积函数K(E)仅在离散点上已知；
1)复化求积公式
积分限[0，Ein]可以被分成几个小积分段，在函数线下的区域被积分间隔分成若干个细长条的区域；
0＝E0＜E1＜…＜EM＝Ein
分割方法确定，方程(21)变成：
其中
2)辛甫森公式
辛甫森公式是在数字化积分中应用最普遍的简化公式，辛甫森公式又称三点求积公式、抛物线公式，是一个近似积分，等于步长乘以被积函数分别在左边界、中心点、右边界点的加权平均。
其中步长h＝Em+1-Em，是截断误差。
在方程(23)中使用相等的子区间，将辛甫森公式求和，变成复化辛甫森公式：
当M为偶数时，
当M为奇数时，
其中M＝Ein/h是子积分段的数目。总的误差是每一子积分段的误差之和：
通常将方程(25)写成下面的形式
其中
ψm＝μt(Em)t (28)
其中不能用方程(28)来计算总的透射信号T，只能在特定的能级Em下，用测量的T值来估计面积吸收参数ψm。
与现有技术相比，本发明具有的优点和积极效果是：
本发明的小像素透射多信息融合高速安检系统，
1)高分辨率双能量L型探测器，它具有将有机物、无机物、混合物进行科学区域分割的功能，且探测器成像的像素尺寸是目前国际及其国内最小的。它可分辨出0.09mm的金属丝，并能穿透28mm的钢板；
2)建立了具有自主知识产权的常见危险品数据库。本数据库可达到自动识别旅客行包内常见危险品的目的；
3)实现了旅客行包的图形图像处理软件系统。可在双显示器上显示出旅客行包的黑白图像和彩色图像，并包括图形的边缘增强、图形放大、图形局部放大、有机物剔出、无机物剔出、黑白图像反转、伪彩色图像、图像增暗、图像增亮、危险品报警等功能；
4)实现了危险品检查仪设备状态的自检功能。在显示器底部设置一滚动条状态监测检查仪的自身信息，并设有442接口，做到检测数据的远程通讯和设备自身状态的远程诊断；
5)机械系统具有可折叠性，便于运输和机械系统的维护及其机械部件的更换；
6)实现细化分析、区域分析及特殊图像处理等数学模型；
7)推证出一个提取R值的算法并研究减小物体厚度对R值的影响的方法；
R是与有效原子序数直接相关的量，通过R可以将有机物从金属和无机物中分离出来。这是很重要的，因为大多数违禁品包含在有机物中，这就大大缩小了检查范围。进一步推证出减小物体厚度对R值的影响的方法，使R值的变化相对较小，这可提高物质识别的准确性。
附图说明
图1为本发明的交叉带式高速分拣机的结构示意简图；
图2为本发明的控制装置的构成系统图；
图3为本发明的主控子系统的构成系统图。
具体实施方式
下面结合附图和具体实施方式对本发明技术作进一步详细说明。
如图1、2、3所示，本发明的小像素透射多信息融合高速安检系统，该小像素透射多信息融合高速安检系统由交叉带式高速分拣机和控制交叉带式高速分拣机运行的控制装置1组成；
所述交叉带式高速分拣机包括输送装置3、分类装置2及分拣道口4；所述控制装置1按时间序列分别读取邮件的条码、邮件面单信号，接收和处理分拣信号；输送装置3采用交叉带小车在固定轨道上运行的输送方式，运载着待分拣的邮件鱼贯通过控制装置1；控制装置1识别邮件条形码的信息，向分类装置2发送邮件的分拣指令，并对邮件进行分类；分类装置2根据控制装置1 发出的分拣指令，改变邮件在输送装置上3运行的方向；交叉带小车按照分类装置2的指令将邮件送达到指定的分拣道口4处，小车托盘倾斜邮件滑下轨道台车进入分拣道口4；
所述控制装置2包括用于处理接收到的信息并发出指令信息的主控子系统、接收小车托盘的“车到”信号以及“高包裹”判别信号的光障开关检测子系统、采集邮件外部形状图像信息的外部图像采集子系统、扫描邮件并获取邮件地藉信息的全方位条形码扫描仪、识别邮件包裹的内部物质属性的X光检查装置、采集邮件放射性指标数据、挥发性指标数据的物理参数采集子系统；所述X光检查装置为可以发出和接收两种能级的双能量的X光检查装置，所述全方位条形码扫描仪、X光检查装置、光障开关检测子系统、外部图像采集子系统、物理参数采集子系统与主控子系统通信联系；
所述光障开关检测子系统监测“车到”和“高包裹”信号，并发送给主控子系统作为整个系统的周期控制启动信号；邮件鱼贯通过全方位条形码扫描仪，全方位条形码扫描仪按时间序列分别读取邮件的条形码、邮件面单信号，对载物小车身份识别；所述外部图像采集子系统监测采集邮件外部形状图像信息；物理参数采集子系统采用有效原子序数物质分类计算对被检目标实现定性识别，获取放射性指标数据和挥发性指标数据，实现ppM级气味识别；所述X光检查装置识别邮件包裹的内部物质属性，对被检目标进行线扫描透照图像获取；所述主控子系统对外部形状图像信号、原始X射线阵列信号进行数字化图形图像处理，并由物质分类实验建立材料分类区域及分类阈值，界定有机物、无机物、混合物，分辨出被检包裹内部物质的属性。
作为对上述技术方案的改进，所述主控子系统包括由若干个电脑组成的用于处理和显示X光图像的X光检查工作站，一台电脑组成的用于实现系统总体控制以及外部图像和其它物理参数采集的主控工作站，一台电脑组成的负载X 光透照数据采集、以及与主控工作站通信实现控制流程、并且控制X光检查工作站工作的X光透照控制工作站，所述主控工作站、X光检查工作站、X光透照控制工作站通过局域网相连接。
作为对上述技术方案的改进，所述X光检查装置包括X射线发生器和L形的X射线接收器，该X光检查装置安装在条形码扫描仪后部，X射线发生器位于输送装置的底部，L形的X射线接收器位于输送装置的顶部；所述X光检查装置与X光透照控制工作站连接。
作为对上述技术方案的改进，所述物理参数采集子系统包括γ中子检测仪、识别简单和复杂气味的、毫秒级响应时间的电子鼻气味传感器，该电子鼻气味传感器为PID光离子传感器，对待检测邮件进行气体检测；该γ中子检测仪检测邮件中是否含有射线。
作为对上述技术方案的改进，所述光障开关检测子系统包括光障开关；所述外部图像采集子系统包括工业数字相机，所述工业数字相机内置CCD传感器。
本发明并提供了一种利用上述小像素透射多信息融合高速安检系统对邮件进行识别的方法，其特征在于，该识别方法的步骤是：
S1、采用双能量的X光检查装置对邮件检测，得到被扫描邮件的平均原子数；
S2、将物质分类，设Z为物质的原子序数，根据定义，当1＜Z≤10时，为有机物，当Z＞20时，为无机物，当10＜Z≤20时为混合物；Z＝10和Z＝20 分别为有机物与混合物、无机物与混合物的边界；
S3、根据有效原子序数计算公式找到边界物质，该有效原子序数计算公式：
式中：wi为化合物质种类当量个数；为物质的原子序数。
根据上式，计算得到碳化硅的有效原子序数大约等于10，钛六铝四钒的有效原子序数大约等于20，这样我们就拿碳化硅和钛六铝四钒代替有机物与无机物的两个边界；
不同物质的质量吸收系数不同，X射线所能穿透的最大厚度也不同，有如下关系：
当双能量的X光检查装置穿透率为22mm钢板时，在100KeV时，
(μ/ρ)Fe＝0.3713(μ/ρ)Al＝0.1704(μ/ρ)Si＝0.1835
(μ/ρ)C＝0.1514(μ/ρ)V＝0.2877(μ/ρ)Ti＝0.21
ρFe＝7.874ρAl＝2.699ρSi＝2.330ρC＝1.7ρV＝6.11ρTi＝4.450
可以计算出碳化硅和钛六铝四钒的最大穿透厚度分别为108mm和36mm；
S4、将实验数据用曲线专家分别进行拟合得到两条边界曲线及其数学模型；
混合物与无机物的边界曲线数学模型为：
T(x2)＝-50587+26.8056·x2-0.00463·x2 2+3.7233·10-7·x2 3- 1.1318·10-11·x2 4 (3)
混合物与有机物的边界曲线数学模型为：
S(x4)＝-16403+9.6105·x4-0.00143·x4 2+1.06446·10-7·x4 3- 3.055·10-12·x4 4 (4)
S5、计算物质的有效原子序数和密度；
S6、双能量的X光检查装置的分类探测信号分析；
A1、透射信号，对于透射图像，探测信号与光电效应和散射都有关，有如下方程：
其中T(x，y)是在投影图中对应于坐标(x，y)探测到的透射信号，N(E) 是能量位于E和E+dE之间的X射线光子数，μt(x′，y′，z′，E)是沿着X射线方向的总的线吸收系数(cm-1)，Pd(E)是探测器的光子能量为E时的探测概率，E是光子能量；
A2、背散射信号；对于向后散射，探测信号B(x，y)估计为：
其中，与公式(11)中相同的符号其含义相同，其余的，B(x，y)是在投影图中(x，y)点的散射光子信号；Z′是发生散射时，笔状射线束沿射线方向进入物体的深度；σs(E)是散射的横截面积，是能量的函数；ρ(x′，y′，z′)是包裹中的物体沿X射线方向的密度(g/cm3)；E′是散射后的光子能量。射；
A3、向前散射信号的表达式与向后散射信号的表达式相似，此式中第二个积分项的积分限为从Z′到D；
除了F(x，y)表示在投影图中(x，y)位置的散射光子信号外，其它的符号与 (11)、(12)两式中的含义相同。
A4、提取R值，由双能量的X光检查装置透射图像的低、高能灰度值计算出R值。由散射图像的灰度值可以计算出L值；
双能量的X光检查装置使用将TH和TL相结合的方法，得到R值：
R称作双能量，TH0和TL0是在自由空间获得的高能、低能透射信号。
S7、考虑物质厚度对双能量透射图像的影响；
一定输入能量Ein下，总的透射信号是能量范围内透射信号的积分，
将透射信号流程改写为：
对于简单的物体，总的线吸收系数为μt(x′，y′，E)，厚度为t，于是方程(17)变为：
对于在图像中任意位置(x，y)的物体，从方程(18)中去除空间坐标，得到：
从方程(19)可知，总的透射信号强度取决于物体的厚度；
双能量的X光检查装置受到X射线源的多色谱线的影响；对于某一特定的物质，输入能量分别为EH，EL，于是有：
从方程(20)可知，R不仅与物质的原子序数Z，原子重量A，密度ρ有关，也于物质的厚度有关；
S8、克服厚度影响的数字化方法
对于给定的双能量的X光检查装置，在方程(20)中，入射的X射线光子数N(E)，X射线光子能量E，光子的探测概率Pd(E)都是已知的或是通过仿真可以得到的参数；
让K(E)＝N(E)EPd(E),在特定点的透射信号可写成：
被积函数K(E)仅在离散点上已知；
1)复化求积公式
积分限[0，Ein]可以被分成几个小积分段，在函数线下的区域被积分间隔分成若干个细长条的区域；
0＝E0＜E1＜…＜EM＝Ein
分割方法确定，方程(21)变成：
其中
2)辛甫森公式
辛甫森公式是在数字化积分中应用最普遍的简化公式，辛甫森公式又称三点求积公式、抛物线公式，是一个近似积分，等于步长乘以被积函数分别在左边界、中心点、右边界点的加权平均。
其中步长h＝Em+1-Em，是截断误差。
在方程(23)中使用相等的子区间，将辛甫森公式求和，变成复化辛甫森公式：
当M为偶数时，
当M为奇数时，
其中M＝Ein/h是子积分段的数目。总的误差是每一子积分段的误差之和：
通常将方程(25)写成下面的形式
其中
ψm＝μt(Em)t (28)
其中不能用方程(28)来计算总的透射信号T，只能在特定的能级Em下，用测量的T值来估计面积吸收参数ψm。
与现有技术相比，本发明具有的优点和积极效果是：
本发明的小像素透射多信息融合高速安检系统，
1)高分辨率双能量L型探测器，它具有将有机物、无机物、混合物进行科学区域分割的功能，且探测器成像的像素尺寸是目前国际及其国内最小的。它可分辨出0.09mm的金属丝，并能穿透28mm的钢板；
2)建立了具有自主知识产权的常见危险品数据库。本数据库可达到自动识别旅客行包内常见危险品的目的；
3)实现了旅客行包的图形图像处理软件系统。可在双显示器上显示出旅客行包的黑白图像和彩色图像，并包括图形的边缘增强、图形放大、图形局部放大、有机物剔出、无机物剔出、黑白图像反转、伪彩色图像、图像增暗、图像增亮、危险品报警等功能；
4)实现了危险品检查仪设备状态的自检功能。在显示器底部设置一滚动条状态监测检查仪的自身信息，并设有442接口，做到检测数据的远程通讯和设备自身状态的远程诊断；
5)机械系统具有可折叠性，便于运输和机械系统的维护及其机械部件的更换；
6)实现细化分析、区域分析及特殊图像处理等数学模型；
7)推证出一个提取R值的算法并研究减小物体厚度对R值的影响的方法；
R是与有效原子序数直接相关的量，通过R可以将有机物从金属和无机物中分离出来。这是很重要的，因为大多数违禁品包含在有机物中，这就大大缩小了检查范围。进一步推证出减小物体厚度对R值的影响的方法，使R值的变化相对较小，这可提高物质识别的准确性。
以上所述仅为本发明的较佳实施案例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (5)
1.一种小像素透射多信息融合高速安检系统对邮件进行识别的方法，该小像素透射多信息融合高速安检系统由交叉带式高速分拣机和控制交叉带式高速分拣机运行的控制装置组成；
所述交叉带式高速分拣机包括输送装置、分类装置及分拣道口；所述控制装置按时间序列分别读取邮件的条码、邮件面单信号，接收和处理分拣信号；输送装置采用交叉带小车在固定轨道上运行的输送方式，运载着待分拣的邮件鱼贯通过控制装置；控制装置识别邮件条形码的信息，向分类装置发送邮件的分拣指令，并对邮件进行分类；分类装置根据控制装置发出的分拣指令，改变邮件在输送装置上运行的方向；交叉带小车按照分类装置的指令将邮件送达到指定的分拣道口处，小车托盘倾斜邮件滑下轨道台车进入分拣道口；
所述控制装置包括用于处理接收到的信息并发出指令信息的主控子系统、接收小车托盘的“车到”信号以及“高包裹”判别信号的光障开关检测子系统、采集邮件外部形状图像信息的外部图像采集子系统、扫描邮件并获取邮件地藉信息的全方位条形码扫描仪、识别邮件包裹的内部物质属性的X光检查装置、采集邮件放射性指标数据、挥发性指标数据的物理参数采集子系统；所述X光检查装置为可以发出和接收两种能级的双能量的X光检查装置，所述全方位条形码扫描仪、X光检查装置、光障开关检测子系统、外部图像采集子系统、物理参数采集子系统与主控子系统通信联系；
所述光障开关检测子系统监测“车到”和“高包裹”信号，并发送给主控子系统作为整个系统的周期控制启动信号；邮件鱼贯通过全方位条形码扫描仪，全方位条形码扫描仪按时间序列分别读取邮件的条形码、邮件面单信号，对载物小车身份识别；所述外部图像采集子系统监测采集邮件外部形状图像信息；物理参数采集子系统采用有效原子序数物质分类计算对被检目标实现定性识别，获取放射性指标数据和挥发性指标数据，实现ppM级气味识别；所述X光检查装置识别邮件包裹的内部物质属性，对被检目标进行线扫描透照图像获取；所述主控子系统对外部形状图像信号、原始X射线阵列信号进行数字化图形图像处理，并由物质分类实验建立材料分类区域及分类阈值，界定有机物、无机物、混合物，分辨出被检包裹内部物质的属性；
其特征在于，该识别方法的步骤是：
S1、采用双能量的X光检查装置对邮件检测，得到被扫描邮件的平均原子数；
S2、将物质分类，设Z为物质的原子序数，根据定义，当1＜Z≤10时，为有机物，当Z＞20时，为无机物，当10＜Z≤20时为混合物；Z＝10和Z＝20分别为有机物与混合物、无机物与混合物的边界；
S3、根据有效原子序数计算公式找到边界物质，该有效原子序数计算公式：
式中:wi为化合物质种类当量个数；为物质的原子序数；
根据上式，计算得到碳化硅的有效原子序数大约等于10，钛六铝四钒的有效原子序数大约等于20，这样我们就拿碳化硅和钛六铝四钒代替有机物与无机物的两个边界；
不同物质的质量吸收系数不同，X射线所能穿透的最大厚度也不同，有如下关系：
当双能量的X光检查装置穿透率为22mm钢板时，在100KeV时，
(μ/ρ)Fe＝0·3713(μ/ρ)Al＝0.1704(μ/ρ)Si＝0.1835
(μ/ρ)C＝0.1514(μ/ρ)V＝0.2877(μ/ρ)Ti＝0.21
ρFe＝7.874ρAl＝2.699ρsi＝2.330ρC＝1.7ρV＝6.11ρTi＝4.450
可以计算出碳化硅和钛六铝四钒的最大穿透厚度分别为108mm和36mm；
S4、将实验数据用曲线专家分别进行拟合得到两条边界曲线及其数学模型；
混合物与无机物的边界曲线数学模型为：
T(x2)＝-50587+26.8056·x2-0.00463·x2 2+3.723333·10-7·x2 3-1.1318·10-11·x2 4(3)
混合物与有机物的边界曲线数学模型为：
S(x4)＝-16403+9.6105·x4-0.00143·x4 2+1.06446·10-7·x4 3-3.055.10-12·x4 4(4)
S5、计算物质的有效原子序数和密度；
S6、双能量的X光检查装置的分类探测信号分析；
A1、透射信号，对于透射图像，探测信号与光电效应和散射都有关，有如下方程：
其中T(x，y)是在投影图中对应于坐标(x，y)探测到的透射信号，N(E
是能量位于E和E+dE之间的X射线光子数，μt(x′，y′，z′，E)是沿着X射线方向的总的线吸收系数，Pd(E)是探测器的光子能量为E时的探测概率，E是光子能量；
A2、背散射信号；对于向后散射，探测信号B(x，y)估计为：
其中，与公式(11)中相同的符号其含义相同，其余的，B(x，y)是在投影图中(x，y)点的散射光子信号；Z′是发生散射时，笔状射线束沿射线方向进入物体的深度；σs(E)是散射的横截面积，是能量的函数；ρ(x′，y′，z′)是包裹中的物体沿X射线方向的密度,g/cm3；E′是散射后的光子能量；
A3、向前散射信号的表达式与向后散射信号的表达式相似，此式中第二个积分项的积分限为从Z′到D；
除了F(x，y)表示在投影图中(x，y)位置的散射光子信号外，其它的符号与(11)、(12)两式中的含义相同；
A4、提取R值，由双能量的X光检查装置透射图像的低、高能灰度值计算出R值；由散射图像的灰度值可以计算出L值；
双能量的X光检查装置使用将TH和TL相结合的方法，得到R值：
R称作双能量，TH0和TL0是在自由空间获得的高能、低能透射信号；
S7、考虑物质厚度对双能量透射图像的影响；
一定输入能量Ein下，总的透射信号是能量范围内透射信号的积分，
将透射信号流程改写为：
对于简单的物体，总的线吸收系数为c(x′，y′，E)，厚度为t，于是方程(17)变为：
对于在图像中任意位置(x，y)的物体，从方程(18)中去除空间坐标，得到：
从方程(19)可知，总的透射信号强度取决于物体的厚度；
双能量的X光检查装置受到X射线源的多色谱线的影响；对于某一特定的物质，输入能量分别为EH，EL，于是有：
从方程(20)可知，R不仅与物质的原子序数Z，原子重量A，密度ρ有关，也于物质的厚度有关；
S8、克服厚度影响的数字化方法
对于给定的双能量的X光检查装置，在方程(20)中，入射的X射线光子数N(E)，X射线光子能量E，光子的探测概率Pd(E)都是已知的或是通过仿真可以得到的参数；
让K(E)＝N(E)EPd(E),在特定点的透射信号可写成：
被积函数K(E)仅在离散点上已知；
1)复化求积公式
积分限[0，Ein]可以被分成几个小积分段，在函数线下的区域被积分间隔分成若干个细长条的区域；
0＝E0＜E1＜…＜EM＝Ein
分割方法确定，方程(21)变成：
其中
2)辛甫森公式
辛甫森公式是在数字化积分中应用最普遍的简化公式，辛甫森公式又称三点求积公式、抛物线公式，是一个近似积分，等于步长乘以被积函数分别在左边界、中心点、右边界点的加权平均；
其中步长h＝Em+1-Em，是截断误差；
在方程(23)中使用相等的子区间，将辛甫森公式求和，变成复化辛甫森公式：
当M为偶数时，
当M为奇数时，
其中M＝Ein/h是子积分段的数目；总的误差是每一子积分段的误差之和：
通常将方程(25)写成下面的形式
其中
ψm＝μt(Em)t (28)
其中不能用方程(28)来计算总的透射信号T，只能在特定的能级Em下，用测量的T值来估计面积吸收参数ψm。
2.根据权利要求1所述的小像素透射多信息融合高速安检系统对邮件进行识别的方法，其特征在于：所述主控子系统包括由若干个电脑组成的用于处理和显示X光图像的X光检查工作站，一台电脑组成的用于实现系统总体控制以及外部图像和其它物理参数采集的主控工作站，一台电脑组成的负载X光透照数据采集、以及与主控工作站通信实现控制流程、并且控制X光检查工作站工作的X光透照控制工作站，所述主控工作站、X光检查工作站、X光透照控制工作站通过局域网相连接。
3.根据权利要求1所述的小像素透射多信息融合高速安检系统对邮件进行识别的方法，其特征在于：所述X光检查装置包括X射线发生器和L形的X射线接收器，该X光检查装置安装在条形码扫描仪后部，X射线发生器位于输送装置的底部，L形的X射线接收器位于输送装置的顶部；所述X光检查装置与X光透照控制工作站连接。
4.根据权利要求1所述的小像素透射多信息融合高速安检系统对邮件进行识别的方法，其特征在于：所述物理参数采集子系统包括γ中子检测仪、识别简单和复杂气味的、毫秒级响应时间的电子鼻气味传感器，该电子鼻气味传感器为PID光离子传感器，对待检测邮件进行气体检测；该γ中子检测仪检测邮件中是否含有射线。
5.根据权利要求1所述的小像素透射多信息融合高速安检系统对邮件进行识别的方法，其特征在于：所述光障开关检测子系统包括光障开关；所述外部图像采集子系统包括工业数字相机，所述工业数字相机内置CCD传感器。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610843896.9A CN106353343B (zh) | 2016-09-22 | 2016-09-22 | 一种小像素透射多信息融合高速安检识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610843896.9A CN106353343B (zh) | 2016-09-22 | 2016-09-22 | 一种小像素透射多信息融合高速安检识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106353343A CN106353343A (zh) | 2017-01-25 |
| CN106353343B true CN106353343B (zh) | 2019-05-07 |
Family
ID=57859387
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610843896.9A Active CN106353343B (zh) | 2016-09-22 | 2016-09-22 | 一种小像素透射多信息融合高速安检识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106353343B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108254394B (zh) * | 2017-12-28 | 2020-09-01 | 合肥美亚光电技术股份有限公司 | X射线双能检测方法及系统 |
| CN112950877B (zh) * | 2021-02-18 | 2021-10-15 | 黑龙江省机场管理集团有限公司 | 行李安检x光机防入侵系统 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104444185A (zh) * | 2014-11-10 | 2015-03-25 | 镇江威速自动化系统技术有限公司 | 交叉带分拣机 |
| CN104502368B (zh) * | 2014-12-29 | 2018-11-23 | 同方威视技术股份有限公司 | 一体化安全检查系统 |
| CN105523372B (zh) * | 2015-10-26 | 2018-05-18 | 浙江朗奥物流科技有限公司 | 一种交叉带分拣机的上包装置 |
| CN205236453U (zh) * | 2015-12-15 | 2016-05-18 | 田华 | 一种自动化分拣装置 |
-
2016
- 2016-09-22 CN CN201610843896.9A patent/CN106353343B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106353343A (zh) | 2017-01-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9915752B2 (en) | Inspection systems with two X-ray scanners in a first stage inspection system | |
| US7366282B2 (en) | Methods and systems for rapid detection of concealed objects using fluorescence | |
| JP5173405B2 (ja) | ２エネルギー放射線走査および遅発中性子検出による物体調査 | |
| US9239404B2 (en) | Device and method for continuous monitoring of persons, vehicles, containers or packets | |
| GB2424065A (en) | Radiation detection apparatus | |
| CN105074440A (zh) | μ介子断层扫描检查中的初级和次级扫描 | |
| US7877340B2 (en) | System and method for resolving gamma-ray spectra | |
| US10001571B2 (en) | Fission meter and neutron detection using poisson distribution comparison | |
| KR102187231B1 (ko) | 검사 설비 및 검사 방법 | |
| US8401270B2 (en) | Examination of a region using dual-energy radiation | |
| CN106353343B (zh) | 一种小像素透射多信息融合高速安检识别方法 | |
| CN110415409A (zh) | 一种安检与检票一体化的闸机系统 | |
| US8194814B2 (en) | Detecting fission from special nuclear material sources | |
| US20160282481A1 (en) | System and method for detecting neutron, gamma and muon radiations with contiguous plastics scintillators | |
| EP3008498A1 (en) | A method for detecting and identifying nuclear material in a container and radiation monitoring arrangement for conducting the method | |
| US20110193711A1 (en) | Method and device for detecting the presence, in a load, of objects suspected of containing at least one material having a given atomic weight | |
| US20140240134A1 (en) | System and method for the detection and control of illicit trafficking of special nuclear materials | |
| CN219143113U (zh) | 人物一体查验系统 | |
| CN108398727B (zh) | 一种基于图像反馈的x射线成像能量选择装置和选择方法 | |
| van der Ende et al. | An integrated system for non-destructive detection of shielded nuclear material using multiple modes of passive detection | |
| Yifan | Research on a material discrimination method by cosmic ray muon tomography | |
| Stadnikia et al. | Situational Information Guidance for Revised Radiological Detection Thresholds | |
| RU2339021C2 (ru) | Тест-образец | |
| CN101539557B (zh) | 用于放射性物质检测和x光辐射成像的集成系统 | |
| ROBERTO et al. | SYSTEM AND METHOD FOR THE |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |