一种通道式危险液体检测装置与方法
【专利摘要】本发明涉及一种通道式危险液体检测装置与方法。所述装置包括透视扫描探测装置、直线断层扫描探测装置和处理与控制单元。由透视扫描探测装置获取传送带上的容器内液体的透视图，由直线断层扫描探测装置获取容器内液体的断层扫描数据，由处理与控制单元对所述直线断层扫描数据进行数据重排和数据补偿获得容器内液体的投影数据，采用双能重建算法重建液体某一位置的断层图像，并计算高、低能衰减系数、有效原子序数和电子密度，通过将这四个指标与阈值进行比较判断容器内液体是否为危险液体。本发明检测装置能够进行连续探测，探测器成本低，通过对缺失的投影数据进行补偿，并根据四个指标判断液体是否为危险液体，提高了危险液体的识别精度。
【专利说明】
一种通道式危险液体检测装置与方法
技术领域
[0001] 本发明属于危险品安全检查技术领域，涉及一种通道式危险液体检测装置与方 法。
【背景技术】
[0002] 危险液体探测技术包括X射线探测技术、中子探测技术、电磁探测技术和蒸汽痕迹 探测技术等。但是，由于液体的独有特性使其探测具有更高的难度和复杂性，如许多危险液 体的物理属性和日常安全液体很接近，这必然增加液体探测的难度。
[0003] 常见液体通常包装在一个密封容器中，液体分子难以挥发至包装之外，采用蒸汽 痕迹探测技术探测液体时，通常需要开瓶检测，增加了检测的操作难度。因此蒸汽痕迹探测 技术并不适用于密封液体的检测。而中子探测技术又由于设备造价、中子防护等原因目前 很少在公共场所实际应用。电磁探测技术是目前市场上的液体探测设备常用的技术，如通 过介电断层技术探测液体的电子特性(介电常数和传导性）；通过微波探测液体的传导性、 分子极性;通过近红外技术来探测常用液体炸药原料过氧化氢的含量等。不过通常上述设 备只能探测塑料容器包装的液体，对于高密度容器包装的液体探测往往无能为力。
[0004] 利用X射线CT技术进行液体探测可以对被探测材料的几何形状进行判断，结合X射 线衰减理论可以对被探测物质的线性衰减系数、等效原子序数、电子密度等参数进行探测。 同方威视公司和公安部第一研究所都已经有了产品样机并分别申报了各自的专利。但是这 类设备均为基于工业CT的三代CT扫描模式，系统操作需要人工将液体放入检测仓，检测完 毕后还需要人工从检测仓取出液体，操作复杂，不能连续检测，效率低下。
【发明内容】
[0005] 为了解决现有技术中存在的上述问题，本发明提出一种用于识别容器中危险液体 的通道式危险液体检测装置与方法，利用容器的对称信息对液体的断层扫描数据进行补 偿，然后重建断层图像，从而实现危险液体的识别。
[0006] 为达到上述目的，本发明采用如下技术方案：
[0007] -种通道式危险液体检测装置，包括透视扫描探测装置、直线断层扫描探测装置 和处理与控制单元。所述透视扫描探测装置包括透视扫描X射线源和透视扫描探测器，用于 获取传送带上容器内液体的透视图；所述直线断层扫描探测装置包括直线断层扫描X射线 源和直线断层扫描探测器，用于获取所述容器内液体的直线断层扫描数据。所述处理与控 制单元对所述直线断层扫描数据进行数据重排和数据补偿获得容器内液体的投影数据，重 建液体某一位置的断层图像从而实现危险液体的识别。所述处理与控制单元还输出控制信 号控制所述透视扫描探测装置、直线断层扫描探测装置和传送带的运行。
[0008] 所述透视扫描探测器是由多个探测单元无间隙排列而成的线阵;所述直线断层扫 描探测器是由多个探测单元排列成的间距大于1厘米的线阵。
[0009] 进一步地，所述透视扫描探测器为直线型、L型、U型或弧形线阵探测器。
[0010] 进一步地，所述透视扫描探测器的探测单元的排列方向与所述直线断层扫描探测 器的探测单元的排列方向在空间上成90度角。
[0011] 进一步地，所述透视扫描探测器为单能探测器或双能夹层探测器；
[0012] 进一步地，所述直线断层扫描探测器为双能夹层探测器或光子计数探测器。
[0013] 进一步地，所述检测装置还包括与所述处理与控制单元相连的用于显示透视图像 和断层图像的显示单元。
[0014] -种识别危险液体的方法，包括以下步骤：
[0015] 步骤1，由透视扫描探测装置获取传送带上容器内液体的透视图像，从而获得容器 内液体的整体轮廓结构信息。
[0016] 步骤2,由直线断层扫描探测装置获取容器内液体的直线断层扫描数据，采用双能 重建算法，利用经数据重排和数据补偿获得的容器内液体的投影数据，重建液体某一位置 的断层图像，包括高、低能衰减系数图像、有效原子序数图像和电子密度图像。并根据所述 图像的数值计算高、低能衰减系数、有效原子序数和电子密度。
[0017] 步骤3,将高、低能衰减系数、有效原子序数和电子密度与设定的阈值分别进行比 较，根据比较结果判断容器内液体是否为危险液体。
[0018] 进一步地，步骤2采用双能重建算法重建液体某一位置的断层图像的方法包括以 下步骤：
[0019] 步骤2.1，获取直线断层高、低能投影数据Pro jH、ProjL，通过数据重排获得平行束 投影数据SinoH、SinoL，对SinoH、SinoL进行数据补偿得到SinoH-C、SinoL_C;
[0020] 步骤2.2，根据基材料模型，利用5丨11〇11-(：、5111〇1-(：进行双能投影分解获得基材料 Sin 〇A、Sin〇B;或根据基效应模型，利用Sin〇H-C、Sin〇L-C进行双能投影分解获得基效应投 影SinoA、SinoB;
[0021] 步骤2.3，对SinoA、SinoB进行重建，获得断层图像B1和B2;对SinoH、SinoL进行重 建，获得SI iceH和SI iceL;如果所述容器为对称容器，采取滤波反投影重建算法；如果所述 容器为非对称容器，采取迭代重建算法；
[0022]步骤2.4,利用断层图像B1和B2,计算液体的有效原子序数断层图像和电子密度断 层图像。
[0023] 进一步地，步骤2.1所述数据重排的方法包括以下步骤：
[0024] (1)根据由直线断层扫描探测装置获取的容器内液体的直线断层扫描数据 1 inogram的起始位置和结束位置，确定数据重排时选择的旋转中心的横坐标xo;
[0025] (2)根据linogram中各个视角的投影数据的重心及投影方向确定数据重排时选择 的旋转中心的纵坐标y〇;
[0026] (3)以（XQ，yo)为旋转中心对linogram进行投影空间变换，得到平行束投影数据正 弦图 sinogram。
[0027]进一步地，步骤2.1所述数据补偿是对sinogram中缺失的数据（即在容器投影死角 缺失的投影数据)进行数据补偿，具体方法包括以下步骤：
[0028] (1)对sinogram进行二值化处理；
[0029] (2)计算sinogram中行数据的宽度；
[0030] (3)计算行数据宽度的标准差以及最大宽度和最小宽度的比值；
[0031] (4)将所述标准差和所述比值分别与设定的阈值进行比较，如果所述标准差和所 述比值均小于设定的阈值，所述容器为对称容器，用缺失数据点的对称点的数据对缺失数 据进行补偿;否则，所述容器为非对称容器，维持原数据。
[0032] 与现有技术相比，本发明具有以下有益效果：
[0033] 本发明提出的通道式危险液体检测装置包括透视扫描探测装置、直线断层扫描探 测装置和处理与控制单元。由透视扫描探测装置获取传送带上的容器内液体的透视图，由 直线断层扫描探测装置获取容器内液体的断层扫描数据，由处理与控制单元对所述直线断 层扫描数据进行数据重排和数据补偿获得容器内液体的投影数据，采用双能重建算法重建 液体某一位置的断层图像，根据图像的数值计算高低能衰减系数、有效原子序数和电子密 度，通过将这四个指标与设定的阈值进行比较判断容器内液体是否为危险液体。该检测装 置能够进行高效率的连续探测；由于直线断层扫描探测装置的探测器采用稀疏排列的线阵 探测器组，因此本发明能够有效降低探测器成本;本发明采用透视扫描探测装置，便于观察 容器内液体的内部结构信息；另外，本发明通过对对称容器内液体缺失的投影数据进行补 偿，而且根据高、低能衰减系数、有效原子序数和电子密度四个指标判断液体是否为危险液 体，大大提高了危险液体的识别精度。
【附图说明】
[0034] 图1为本发明实施例涉及的所述检测装置的结构示意图；
[0035] 图2为直线断层扫描探测器中探测单元的排列示意图。
[0036] 图中：1-透视扫描探测装置，11-透视扫描X射线源，12-透视扫描探测器，2-直线断 层扫描探测装置，21-直线断层扫描X射线源，22-直线断层扫描探测器，3-处理与控制单元， 4-显示单元，5-容器，6-传送带。
【具体实施方式】
[0037] 下面结合附图和实施例对本发明做进一步说明。
[0038] -种通道式危险液体检测装置，如图1、2所示，包括透视扫描探测装置1、直线断层 扫描探测装置2和处理与控制单元3。所述透视扫描探测装置1包括透视扫描X射线源11和透 视扫描探测器12,用于获取传送带6上容器5内液体的透视图；所述直线断层扫描探测装置2 包括直线断层扫描X射线源21和直线断层扫描探测器22,用于获取所述容器5内液体的直线 断层扫描数据。所述处理与控制单元3对所述直线断层扫描数据进行数据重排和数据补偿 获得容器内液体的投影数据，重建液体某一位置的断层图像从而实现危险液体的识别。所 述处理与控制单元3还输出控制信号控制所述透视扫描探测装置1、直线断层扫描探测装置 2和传送带6的运行。所述透视扫描探测器12是由多个探测单元无间隙排列而成的线阵;所 述直线断层扫描探测器22是由多个探测单元排列成的间距大于1厘米的线稀疏阵，如图2所 示。采用稀疏排列的线阵，可以减少探测器的数量，从而有效降低探测器的能耗及装置的成 本。
[0039]所述透视扫描探测器12为直线型、L型、U型或弧形线阵探测器。图2中的透视扫描 探测器12采用L型。
[0040]所述透视扫描探测器12的探测单元的排列方向与所述直线断层扫描探测器22的 探测单元的行排列方向在空间上成90度角。如图2所示。
[0041] 所述透视扫描探测器12为单能探测器或双能夹层探测器。双能夹层探测器包括高 能探测器和低能探测器，在高能探测器和低能探测器之间设有一层金属片，一般为铜片，通 过金属片的过滤使高低能探测器接收到不同能量的射线。
[0042] 所述直线断层扫描探测器22为双能夹层探测器或光子计数探测器。光子计数探测 器通过对不同能量的光子计数实现成像，可以分辨不同能量段的光子并进行计数，具有很 强的能量分辨能力。其不足是受计数率等参数制约信噪比难以保证。双能夹层探测器可以 通过一次扫描实现双能成像，但能谱区分度不及光子技术探测器。
[0043] 所述检测装置还包括与所述处理与控制单元3相连的显示单元4。显示单元4用于 显示透视图像和断层图像。
[0044] -种识别危险液体的方法，包括以下步骤：
[0045] 步骤1，由透视扫描探测装置获取传送带上容器内液体的透视图像，从而获得容器 内液体的整体轮廓结构信息；
[0046] 步骤2,由直线断层扫描探测装置获取容器内液体的直线断层扫描数据，采用双能 重建算法，利用经数据重排和数据补偿获得的容器内液体的投影数据，重建液体某一位置 的断层图像，包括高、低能衰减系数图像、有效原子序数图像和电子密度图像。并根据所述 图像的数值计算高、低能衰减系数、有效原子序数和电子密度。
[0047] 步骤3,将高低能衰减系数、有效原子序数和电子密度与设定的阈值分别进行比 较，根据比较结果判断容器内液体是否为危险液体。
[0048] 步骤2采用双能重建算法重建液体某一位置的断层图像的方法包括以下步骤： [0049]步骤2.1，获取直线断层高、低能投影数据Pro jH、ProjL，通过数据重排获得平行束 投影数据SinoH、SinoL，对SinoH、SinoL进行数据补偿得到SinoH-C、SinoL_C;
[0050] 步骤2.2，根据基材料模型，利用5111〇11-(：、5111〇1-(：进行双能投影分解获得基材料 Sin 〇A、Sin〇B;或根据基效应模型，利用Sin〇H-C、Sin〇L-C进行双能投影分解获得基效应投 影SinoA、SinoB;
[0051] 基于双能成像原理由SinoH、SinoL获得SinoA、SinoB，方法如下：
[0052]在200keV以内的射线能量范围内，射线与物质的相互作用遵循康普顿散射和光电 效应。物质的线性衰减系数y(E)满足如下模型：
[0053] y(E) =acfKN(E)+aPfp(E)
[0054] 其中，fP(E)、f?(E)为只与能量E有关而与材质无关的分解系数，aP、ac是独立于能 量只与材质有关的物理量，a P表示光电效应系数，a。为康普顿散射效应系数，且：
L0058J 其中，〇 = ￡/51〇.9751?^，11、12为常数，0为物质密度，2为原子序数4为原子量。该 模型表示在一定的射线能量范围内，物质的衰减是光电效应和康普顿散射共同作用的结 果。此模型通常称之为基效应模型(也可称双效应模型）。
[0059] 与基效应模型对应的还有一种关于物质衰减系数的物理模型，即基材料模型：
[0060] y(E) =biiii(E)+b2li2(E)
[0061 ]其中，m (E)、y2 (E)分别为两种基材料的线性衰减系数，h、b2分别为对应两种基材 料的分解系数，对于某一固定的物质，是两个常数。基材料模型表示任何一种物质的 线性衰减系数都可由两种基材料的线性衰减系数线性叠加而成。
[0062] 令：
[0063] Ac = /acdl ,AP = /aPdl ,Bi = /bidl ,B2 = /b2dl
[0064] Ac、AP、Bi、B2为a c、aP、h、b2的线积分投影值。根据宽能谱射线条件下的BEER定律：
[0065] p//^..\n[js/.(E)cxp\~AJK^E)^ ⑴ ",=-||1吖5,(/:.>邛卜/心，.（/_0-/又"，(/：)16/￡1+1"(^,(/:>/￡
[0066] [/^-in|j5//(￡)expM!，|(￡)-? 2/,(￡)l^HlnJ ⑵
[0067]其中，Sh(E)、Sl(E)分别为高低能系统能谱，Ph、Pl分别为高低能投影（即SinoH、 3111〇〇。基于投影分解的双能(：1'预处理重建算法的核心为根据式（1)、式(2)求解^1)、8 1、 B2(即基材料投影Sin〇A、Sin〇B和基效应投影Sin 〇A、Sin〇B)，称这个求解过程为投影分解过 程。然后根据CT重建的原理，利用滤波反投影图像重建算法计算~、^、^、^，并由此计算材 质的等效原子序数Zm和电子密度&，公式如下：
[0069] pe = K2ac = blPel+b2Pe2
[0070] 其中，心上为常数，n取3或HZ2分别为两种基材料的原子序数，分别为两 种基材料的电子密度。
[0071] 步骤2.3，对SinoA、SinoB进行重建，获得断层图像B1和B2;对SinoH、SinoL进行重 建，获得SI iceH和SI iceL;如果所述容器为对称容器，采取滤波反投影重建算法；如果所述 容器为非对称容器，采取迭代重建算法；
[0072]步骤2.4,利用断层图像B1和B2,计算液体的有效原子序数断层图像和电子密度断 层图像。
[0073] 步骤2.1所述数据重排的方法包括以下步骤：
[0074] (1)根据由直线断层扫描探测装置获取的容器内液体的直线断层扫描数据 1 inogram的起始位置和结束位置，确定数据重排时选择的旋转中心的横坐标xo;
[0075] (2)根据linogram中各个视角的投影数据的重心及投影方向确定数据重排时选择 的旋转中心的纵坐标y〇;
[0076] (3)以（XQ，yo)为旋转中心对linogram进行投影空间变换，得到平行束投影数据正 弦图 sinogram。
[0077] 步骤2.1所述数据补偿是对sinogram中缺失的数据进行数据补偿，具体方法包括 以下步骤：
[0078] (1)对sinogram进行二值化处理；
[0079]将每个像素点的sinogram值与设定的阈值进行比较，将小于所述阈值的sinogram 值设为〇，大于或等于所述阈值的s inogram值设为1。所述阈值稍大于空气背景的阈值。 [0080] (2)计算sinogram中行数据的宽度；
[0081]某一视角下，sinogram值为1的最左边的点与最右边的点之间的距离（以像素点的 个数表示)为行数据的宽度。
[0082] (3)计算行数据宽度的标准差以及最大宽度和最小宽度的比值；
[0083] (4)将所述标准差和所述比值分别与设定的阈值进行比较，如果所述标准差和所 述比值均小于设定的阈值，所述容器为对称容器，用缺失数据点的对称点的数据对缺失数 据进行补偿;否则，所述容器为非对称容器，不进行数据补充定，维持原数据。
[0084] 设待重建图像为f(x，y)或/(/%(柱面坐标），f(x，y)在不同视角巾下的投影为 (xr)，则平行束重建公式为：
[0085] f(r, S'[/' cos(6^ - (p), (p\d(p 0
[0086] 迭代重建算法为基于全变分最小约束代数重建算法，是在代数重建算法的基础上 增加一个全变分最小的约束，将重建问题转化为一个约束优化问题。
[0087] 可以采用多种方法解约束优化问题，本实施例采用最为常用的梯度下降法，首先 计算全变差的梯度，然后沿着负梯度方向寻找优化结果。梯度的计算公式如下：
[0089]其中，f表示断层图像，fs,t表示坐标为（s，t)的像素点的图像值，| |fs,t| |TV表示断 层图像的全变分，其表达式为：
[0091] (4)利用断层图像B1和B2,计算液体的有效原子序数断层图像和电子密度断层图 像。
[0092] 本发明不限于上述实施方式，本领域技术人员所做出的对上述实施方式任何显而 易见的改进或变更，都不会超出本发明的构思和所附权利要求的保护范围。
【主权项】
1. 一种通道式危险液体检测装置，其特征在于，包括透视扫描探测装置、直线断层扫描 探测装置和处理与控制单元;所述透视扫描探测装置包括透视扫描X射线源和透视扫描探 测器，用于获取传送带上容器内液体的透视图；所述直线断层扫描探测装置包括直线断层 扫描X射线源和直线断层扫描探测器，用于获取所述容器内液体的直线断层扫描数据;所述 处理与控制单元对所述直线断层扫描数据进行数据重排和数据补偿获得容器内液体的投 影数据，重建液体某一位置的断层图像从而实现危险液体的识别;所述处理与控制单元还 输出控制信号控制所述透视扫描探测装置、直线断层扫描探测装置和传送带的运行； 所述透视扫描探测器是由多个探测单元无间隙排列而成的线阵;所述直线断层扫描探 测器是由多个探测单元排列成的间距大于1厘米的线阵。2. 根据权利要求1所述的通道式危险液体检测装置，其特征在于，所述透视扫描探测器 为直线型、L型、U型或弧形线阵探测器。3. 根据权利要求1所述的通道式危险液体检测装置，其特征在于，所述透视扫描探测器 的探测单元的排列方向与所述直线断层扫描探测器的探测单元的排列方向在空间上成90 度角。4. 根据权利要求1所述的通道式危险液体检测装置，其特征在于，所述透视扫描探测器 为单能探测器或双能夹层探测器。5. 根据权利要求1所述的通道式危险液体检测装置，其特征在于，所述直线断层扫描探 测器为双能夹层探测器或光子计数探测器。6. 根据权利要求1~5任意一项所述的通道式危险液体检测装置，其特征在于，所述检 测装置还包括与所述处理与控制单元相连的用于显示透视图像和断层图像的显示单元。7. -种识别危险液体的方法，其特征在于，包括以下步骤： 步骤1，由透视扫描探测装置获取传送带上容器内液体的透视图像，从而获得容器内液 体的整体轮廓结构信息； 步骤2,由直线断层扫描探测装置获取容器内液体的直线断层扫描数据，采用双能重建 算法，利用经数据重排和数据补偿获得的容器内液体的投影数据，重建液体某一位置的断 层图像，包括高、低能衰减系数图像、有效原子序数图像和电子密度图像;并根据所述图像 的数值计算高、低能衰减系数、有效原子序数和电子密度； 步骤3,将高、低能衰减系数、有效原子序数和电子密度与设定的阈值分别进行比较，根 据比较结果判断容器内液体是否为危险液体。8. 根据权利要求7所述的识别危险液体的方法，其特征在于，所述步骤2采用双能重建 算法重建液体某一位置的断层图像的方法包括以下步骤： 步骤2.1，获取直线断层高、低能投影数据?仰犯、？仰儿，通过数据重排获得平行束投影 数据 SinoH、SinoL，对 SinoH、SinoL 进行数据补偿得到SinoH-C、SinoL-C; 步骤2.2，根据基材料模型，利用SinoH-C、SinoL-C进行双能投影分解获得基材料 Sin〇A、Sin〇B;或根据基效应模型，利用Sin〇H-C、Sin 〇L-C进行双能投影分解获得基效应投 影SinoA、SinoB; 步骤2.3，对3111(^、3111〇8进行重建，获得断层图像81和82;对3111〇!1、3111〇1^进行重建，获 得SI i ceH和SI i ceL;如果所述容器为对称容器，采取滤波反投影重建算法;如果所述容器为 非对称容器，采取迭代重建算法； 步骤2.4,利用断层图像BI和B2,计算液体的有效原子序数断层图像和电子密度断层图 像。9. 根据权利要求8所述的识别危险液体的方法，其特征在于，步骤2.1所述数据重排的 方法包括以下步骤： (1) 根据由直线断层扫描探测装置获取的容器内液体的直线断层扫描数据Iinogram的 起始位置和结束位置，确定数据重排时选择的旋转中心的横坐标XO; (2) 根据Iinogram中各个视角的投影数据的重心及投影方向确定数据重排时选择的旋 转中心的纵坐标yo; (3) 以（XQ，yQ)为旋转中心对Iinogram进行投影空间变换，得到平行束投影数据正弦图 sinogram〇10. 根据权利要求8所述的识别危险液体的方法，其特征在于，步骤2.1所述数据补偿是 对s inogram中缺失的数据进行数据补偿，具体方法包括以下步骤： (1) 对sinogram进行二值化处理； (2) 计算sinogram中行数据的宽度； (3) 计算行数据宽度的标准差以及最大宽度和最小宽度的比值； (4) 将所述标准差和所述比值分别与设定的阈值进行比较，如果所述标准差和所述比 值均小于设定的阈值，所述容器为对称容器，用缺失数据点的对称点的数据对缺失数据进 行补偿;否则，所述容器为非对称容器，维持原数据。
【文档编号】G01V5/00GK105911604SQ201610371858
【公开日】2016年8月31日
【申请日】2016年5月30日
【发明人】李保磊, 李斌, 陈学亮, 李蒙, 张耀军, 郭双茂, 莫阳, 陈力, 李永清, 张萍宇, 王海鹏, 刘斌, 孟博, 张福, 曹琴琴
【申请人】公安部第研究所, 公安部第一研究所, 北京中盾安民分析技术有限公司