CN112666186A - 一种双能量x射线安全检查设备的物质分类方法 - Google Patents
一种双能量x射线安全检查设备的物质分类方法 Download PDFInfo
- Publication number
- CN112666186A CN112666186A CN202011635381.2A CN202011635381A CN112666186A CN 112666186 A CN112666186 A CN 112666186A CN 202011635381 A CN202011635381 A CN 202011635381A CN 112666186 A CN112666186 A CN 112666186A
- Authority
- CN
- China
- Prior art keywords
- curve
- mixture
- energy
- confidence
- classification
- 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
- 239000000463 material Substances 0.000 title claims abstract description 37
- 238000007689 inspection Methods 0.000 title claims abstract description 23
- 238000000034 method Methods 0.000 title claims abstract description 19
- 239000000203 mixture Substances 0.000 claims abstract description 43
- 239000005416 organic matter Substances 0.000 claims abstract description 18
- 230000007704 transition Effects 0.000 claims abstract description 15
- 229910052782 aluminium Inorganic materials 0.000 claims abstract description 11
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims abstract description 10
- 239000011521 glass Substances 0.000 claims abstract description 8
- 238000001514 detection method Methods 0.000 claims abstract description 7
- 229910001209 Low-carbon steel Inorganic materials 0.000 claims abstract description 6
- 239000000126 substance Substances 0.000 claims description 31
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 3
- 239000011707 mineral Substances 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims 1
- 230000009977 dual effect Effects 0.000 abstract 1
- 238000004040 coloring Methods 0.000 description 7
- XEEYBQQBJWHFJM-UHFFFAOYSA-N iron Substances [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 7
- 238000005070 sampling Methods 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 229920003229 poly(methyl methacrylate) Polymers 0.000 description 2
- 239000004926 polymethyl methacrylate Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011368 organic material Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Abstract
一种双能量X射线安全检查设备的物质分类方法，包括：1）利用双能量X射线安全检查设备分别对不同厚度的有机玻璃、铝、低碳钢标定块进行扫描，采集探测板所接收每种标定块的像素点数据；2）以像素点的低能值为横坐标，高能值为纵坐标，分别对3组数据进行曲线拟合；3）确定有机物和混合物的边界曲线；4）建立边界曲线的过渡区域；5）对过渡区域设置分类置信度；6）置信度函数归一化；7）提前计算物质分类点存成字典模式，生成材料查找字典。通过本发明方法所得的物质分类方式更合理，物质属性判断的准确，同时预生成字典模式，加速了物质判断速度。
Description
技术领域
本发明涉及图像处理技术领域，尤其涉及一种双能量X射线安全检查设备的物质分类方法。
背景技术
传统安检机的运行速度是0.22米/秒，高速安检机一般会在2.5米/秒以上，是传统安检机的10倍以上。所以，高速安检机在单位时间需要处理的图像数据也是传统安检机的10倍以上，使用传统的图像伪彩色上色方法需要进行计算每一点的物质分界点，大量使用计算资源，使得安检机伪彩色上色较慢，同时传统的伪彩色分界点计算不够准确，造成安检机图像伪彩色上色错误，影响包裹图像的展示和人员的准确判断。
发明内容
有鉴于此，须提供一种判读效率高、效果好的物质分类方法，以便设备能准确判断物质类别进行上色，提高物品识别准确率。
本发明的目的是通过以下技术方案来实现的，一种双能量X射线安全检查设备的物质分类方法，包括：
1）利用双能量X射线安全检查设备分别对不同厚度的有机玻璃、铝、低碳钢标定块进行扫描，采集探测板所接收每种标定块的像素点数据。
2）以像素点的低能值为横坐标，高能值为纵坐标，分别对3组数据进行曲线拟合。
3）确定有机物和混合物的边界曲线为：
f(x)1 = k1 * f(x)pmma + k2 * f(x)Al，其中k1 * 6.5 + k2 * 13 = 10，k1 + k2 =1；
确定无机物和混合物的边界曲线为：
f(x)2 = k3 * f(x)Fe + k4 * f(x)Al，其中k3 * 26 + k4 * 13 = 18，k3 + k4 = 1。
4）建立边界曲线的过渡区域。
确定有机物上界曲线为μo＝(1-s1)* f(x)1；
混合物下界曲线为μm1＝(1+s1)* f(x)1；
混合物上界曲线为μm2＝(1-s2)* f(x)2；
无机物下界曲线为μino＝(1+s1)* f(x)2。
在有机物上界曲线为μorg和混合物下界曲线μmix1之间为f(x)1过渡区域；在混合物上界曲线μmix2和无机物下界曲线μino之间f(x)2过渡区域。
5）对过渡区域设置分类置信度，选取置信度变化函数为：
可得，三种物质的置信度函数为：
其中，
，
；
y为安检机获取X射线穿过物质后的高能值，
、
、
、
分别为有机物上界点，混合物下界点，混合物上界点，无机物下界点，
、
、
分别为有机物、混合物、无机物的置信度函数。
6）置信度函数归一化。
通过计算归一化后的置信度
、
、
，取置信度值较高的作为物质判断的依据。通过上述方法进行判断后可以更加准确地区分边界区域物质的区分，提高识别的准确性。
7）提前计算物质分类点存成字典模式，生成材料查找字典。
通过本发明方法所得的物质分类方式更合理，物质属性判断的准确，同时预生成字典模式，加速了物质判断速度，加快伪彩上色速度，进而提高了物品识别准确率。
附图说明
图1为3种物质不同厚度标定块所得的像素点数据示意图。
图2为拟合曲线示意图。
图3为置信度曲线示意图。
具体实施方式
为了更好地理解本说明书一个或多个技术方案，下面将结合本说明书一个或多个实施例中的附图，对本说明书一个或多个实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本说明书一个或多个一部分实施例，而不是全部的实施例。
当今国内外X射线安检系统根据射线源和探测器的种类可以分为单能和双能系统。由于单能射线系统不具备物质区分的能力，因此仅仅通过外形来判断危险品对安检员的经验有很高的要求。双能系统根据使用的射线源个数又可以分为真双能系统和伪双能系统。所谓真双能系统是指安检机具有两个单能射线源，分别发射高能射线和低能射线，同时接收端有两组探测器，分别接收衰减后的高能射线和低能射线。所谓伪双能系统是指安检机只具有一个射线源，该射线源是一个具有一定能谱宽度的多能射线源，X射线穿过物体到达探测器，先被低能探测器接收得到低能数据，接着通过一个铜片低能滤波器，滤除低能部分，留下的高能部分被下一层的高能探测器接收，得到高能数据。双能系统将得到的高低能信号输入PC机，经过一系列数据处理和物质等效原子序数相关的属性值计算，得到的灰度图像通常要转换成彩色图像显示最终在计算机屏幕上呈现有一定物质类别区分的伪彩色图像，大大减少了安检员的判断难度。
双能系安检物质的分类主要包括有机物、混合物、无机物三种，通常的代表物质为有机玻璃、铝、铁，等效原子序数分别为6.5、13和26。国家标准规定等效原子序数Z≤10时，认为是有机物；10≤Z≤18时，为混合物；Z＞18为无机物。而用有机玻璃、铝、铁进行采样拟合出来的曲线无法直接用于物质分类，需进行一定的比例转换才可用于物质的分类和伪彩上色。
为了更加准确、快速地对通过安检的物品进行上色，本发明实施例提供了一种双能量X射线安全检查设备的物质分类方法，该方法包括以下步骤：
1）选择有机玻璃代表有机物，铝代表混合物，低碳钢代表无机物，分别取上述三种物质的不同厚度的长方形样品作为标定块，每个标定块的厚度统一。在这个实施例中，选用了厚度范围为1mm～120mm的若干聚甲基丙烯酸甲脂标定块表示有机物材料；厚度范围为1mm～60mm的铝标定块表示混合物；厚度范围在0.2mm～14mm的钢标定块表示无机物。标定块的大小及形状不作要求，但需厚度统一。
利用双能量X射线安全检查设备分别对3种物质的标定块进行扫描，采集每种物质不同厚度时的数据绘制曲线如图1所示。图1横轴是双能量X射线低能值，纵轴是双能量X射线高能值。其中图2中最上方的“×”是无机物钢不同厚度标定块所得的采样结果，中间的“*”是混合物铝不同厚度标定块所得的采样结果，最下方的“·”是有机物聚甲基丙烯酸甲脂不同厚度标定块所得的采样结果，且厚度均是从左至右依次增加。
2）对图1图中的三种标定块物质进行曲线拟合。根据统计点的分布情况，本实施例中无机物和混合物选用二次多项式拟合，有机物选用线性拟合，分别得到对应的拟合公式f(x)Fe、f(x)Al和f(x)pmma和最终的拟合曲线图，如图2所示。同样的，横坐标表示的是像素点的低能值，纵坐标表示的是高能值。
3）确定有机物和混合物的边界曲线为：
f(x)1 = k1 * f(x)pmma + k2 * f(x)Al，其中k1 * 6.5 + k2 * 13 = 10，k1 + k2 =1；
确定无机物和混合物的边界曲线为：
f(x)2 = k3 * f(x)Fe + k4 * f(x)Al，其中k3 * 26 + k4 * 13 = 18，k3 + k4 = 1。
4）建立边界曲线的过渡区域。
在实际区分有机物、混合物、无机物时因为噪声的存在，物质种类的区分不会以上述曲线为绝对的标准进行分隔判断。实际应用场景中测试发现，探测板某点在上述边界曲线附近，但落入f(x)1曲线下方时，不一定是有机物，有可能是混合物；探测板某点在曲线f(x)2附近，但落入f(x)1曲线和f(x)2曲线之间时，不一定是混合物，有可能是无机物。因此需要对探测板物质分类曲线设置模糊区域进行物质过度：
确定有机物上界曲线为μo＝(1-s1)* f(x)1；
混合物下界曲线为μm1＝(1+s1)* f(x)1；
混合物上界曲线为μm2＝(1-s2)* f(x)2；
无机物下界曲线为μino＝(1+s1)* f(x)2；
其中，s1，s2为大于0小于1的系数，优选的均为0.05。
即在有机物上界曲线μorg和混合物下界曲线μmix1之间为f(x)1过渡区域；在混合物上界曲线μmix2和无机物下界曲线μino之间为f(x)2过渡区域。
5）对过渡区域设置分类置信度，选取置信度变化函数为：
可得，三种物质的置信度函数为：
其中，
，
；
y为安检机获取X射线穿过物质后的高能值，
、
、
、
分别为有机物上界点，混合物下界点，混合物上界点，无机物下界点，
、
、
分别为有机物、混合物、无机物的置信度函数。对上述置信度曲线作图，如图3所示。
6）置信度函数归一化。
这样，通过计算归一化后的置信度
、
、
，取置信度值较高的作为物质判断的依据。通过上述方法进行判断后可以更加准确地区分边界区域物质的区分，提高识别的准确性。
7）计算物质分类点存成字典模式，生成材料查找字典。
因为高速安检机需要处理的图像时长不宜过长，同时在伪彩上色时每一点都需要重复以上的步骤进行物质分类的判断，导致图像伪彩上色时间过长，所以为了节省时间，提前计算物质分类点，将步骤4）中的每个点都进行置信度值计算，存成字典模式，即可在实际应用中，通过物质高低能数值，得到其准确的物质分类。具体的，字典包含四个部分，有机物上界字典O_Top(x),混合物下界字典Mix_Bottom(x), 混合物上界字典Mix_Top(x),无机物下界字典Ino_Bottom(x)，这些字典可通过代入x，x从0到探测板所能接收的最大值进行按1-4的步骤遍历求解即可获得。经应用验证，该方法在物质判断上准确度有极大提升。
以上公开的本申请优选实施例只是用于帮助阐述本申请。可选实施例并没有详尽叙述所有的细节，也不限制该发明仅为所述的具体实施方式。显然，根据本说明书的内容，可作很多的修改和变化。本说明书选取并具体描述这些实施例，是为了更好地解释本申请的原理和实际应用，从而使所属技术领域技术人员能很好地理解和利用本申请。本申请仅受权利要求书及其全部范围和等效物的限制。
Claims (6)
1.一种双能量X射线安全检查设备的物质分类方法，其特征在于，包括以下步骤：1）利用双能量X射线安全检查设备分别对不同厚度的多个有机玻璃、铝、低碳钢标定块进行扫描，采集探测板所接收每种标定块的像素点数据；
2）以像素点的低能值为横坐标，高能值为纵坐标，分别对所述每种标定块的像素点数据进行曲线拟合；
3）确定有机物和混合物的边界曲线为：
f(x)1 = k1 * f(x)pmma + k2 * f(x)Al，其中k1 * 6.5 + k2 * 13 = 10，k1 + k2 = 1；
确定无机物和混合物的边界曲线为：
f(x)2 = k3 * f(x)Fe + k4 * f(x)Al，其中k3 * 26 + k4 * 13 = 18，k3 + k4 = 1；
4）建立边界曲线的过渡区域；确定有机物上界曲线为μo＝(1-s1)* f(x)1；
混合物下界曲线为μm1＝(1+s1)* f(x)1；
混合物上界曲线为μm2＝(1-s2)* f(x)2；
无机物下界曲线为μino＝(1+s1)* f(x)2；
其中，s1和s2为大于0小于1的系数；
得到在有机物上界曲线μorg和混合物下界曲线μmix1之间为f(x)1过渡区域；在混合物上界曲线μmix2和无机物下界曲线μino之间为f(x)2过渡区域；
5）对过渡区域设置分类置信度，选取置信度变化函数为：
可得，三种物质的置信度函数为：
其中，
，
；
y为安检机获取X射线穿过物质后的高能值，
、
、
、
分别为有机物上界点，混合物下界点，混合物上界点，无机物下界点，
、
、
分别为有机物、混合物、无机物的置信度函数；
6）置信度函数归一化，得到
、
、
即分别在有机物、混合物、无机物的置信度值，取置信度最高的作为分类依据。
2.根据权利要求1所述的一种物质分类方法，其特征在于，还包括步骤7）计算物质分类点存成字典模式，生成材料查找字典。
3.根据权利要求1或2任一所述的一种物质分类方法，其特征在于，所述有机玻璃代表有机物，所述铝代表混合物，所述低碳钢代表无机物。
4.根据权利要求1或2任一所述的一种物质分类方法，其特征在于，所述有机玻璃的厚度范围为1mm～120mm；所述铝的厚度范围为1mm～60mm；所述低碳钢的厚度范围为0.2mm～14mm。
5.根据权利要求1或2任一所述的一种物质分类方法，其特征在于，每个标定块厚度均一。
6.根据权利要求1或2任一所述的一种物质分类方法，其特征在于，所述归一化方法为
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011635381.2A CN112666186B (zh) | 2020-12-31 | 2020-12-31 | 一种双能量x射线安全检查设备的物质分类方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011635381.2A CN112666186B (zh) | 2020-12-31 | 2020-12-31 | 一种双能量x射线安全检查设备的物质分类方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112666186A true CN112666186A (zh) | 2021-04-16 |
| CN112666186B CN112666186B (zh) | 2024-03-15 |
Family
ID=75413340
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011635381.2A Active CN112666186B (zh) | 2020-12-31 | 2020-12-31 | 一种双能量x射线安全检查设备的物质分类方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112666186B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110018524A (zh) * | 2019-01-28 | 2019-07-16 | 同济大学 | 一种基于视觉-属性的x射线安检违禁品识别方法 |
| CN111077171A (zh) * | 2020-02-24 | 2020-04-28 | 南京全设智能科技有限公司 | 一种可调脉冲式平板便携式x射线检查装置及其双能量材料分辨方法 |
-
2020
- 2020-12-31 CN CN202011635381.2A patent/CN112666186B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110018524A (zh) * | 2019-01-28 | 2019-07-16 | 同济大学 | 一种基于视觉-属性的x射线安检违禁品识别方法 |
| CN111077171A (zh) * | 2020-02-24 | 2020-04-28 | 南京全设智能科技有限公司 | 一种可调脉冲式平板便携式x射线检查装置及其双能量材料分辨方法 |
Non-Patent Citations (1)
| Title |
|---|
| 马?;丁军航;谭虹;: "基于神经网络的违禁品检测研究", 青岛大学学报(工程技术版), no. 02, 15 May 2020 (2020-05-15) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112666186B (zh) | 2024-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8861675B2 (en) | Method and apparatus for sorting materials according to relative composition | |
| AU2006252235B2 (en) | Method for Inspecting Object Using Multi-Energy Radiations and Apparatus Thereof | |
| EP1711800B1 (de) | Verfahren und vorrichtung zur bestimmung eines objektmaterials | |
| CN105649615B (zh) | Ct定量、三维可视化测试储层致密油赋存状态的方法 | |
| RU2005141128A (ru) | Способ и устройство для определения параметров частиц и рабочих характеристик процессора в системе обработки каменного угля и минералов | |
| US20170184756A1 (en) | Inspection systems for quarantine and methods thereof | |
| CN113926737A (zh) | 一种利用双能双机x射线的煤矸装置及其识别方法 | |
| NL8902196A (nl) | Geautomatiseerde werkwijze voor het identificeren van mineralen en het karakteriseren van gesteenten. | |
| US20060115037A1 (en) | System and a method of automatically sorting objects | |
| CN112666186B (zh) | 一种双能量x射线安全检查设备的物质分类方法 | |
| CN111539251A (zh) | 一种基于深度学习的安检物品识别方法和系统 | |
| JP2011038939A (ja) | エネルギー分散型ｘ線分析装置のスペクトルの分類方法及び装置 | |
| Hassler et al. | Non-destructive provenance study of chert using infrared reflectance microspectroscopy | |
| CN113406130A (zh) | 一种块状岩石样品矿物含量测量的方法 | |
| Kolkoori et al. | High energy X-ray imaging technology for the detection of dangerous materials in air freight containers | |
| CN115326845A (zh) | 矿物分选装置及矿物分选方法 | |
| CN110084786A (zh) | 一种渐变背景的数字x射线图像缺陷自动识别方法 | |
| CN113781426A (zh) | 一种识别液体成分的智能安检方法 | |
| RU2623692C2 (ru) | Система и способ для обнаружения алмазов в кимберлите и способ предварительного обогащения алмазов с их применением | |
| CN219799278U (zh) | 矿物分选装置及矿物分选系统 | |
| Zhang | Assessment of dual-energy X-ray transmission image analysis process on sulphide ore and coal cases | |
| Leisner et al. | Circular economy: Material sorting in waste streams using dual energy X-ray transmission | |
| CN200947087Y (zh) | 能谱调制装置和识别材料设备 | |
| Modise et al. | A Fast and Noise Rejecting Kolmogorov-Smirnvo Sorting Algorithm in X-ray Diamond sorting | |
| CN113075086B (zh) | 一种基于x射线的木片检测及筛选装置 |
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 |