CN108267465A - 多视角成像数据处理方法和设备 - Google Patents

Abstract

公开了一种用于多视角成像数据处理的方法和设备。在一实施例中，通过X射线从两个或更多个视角对被检物体进行照射得到所述多视角成像数据。该方法包括：从多视角成像数据中提取代表被检物体中同一区域的数据；根据所提取的数据，对该区域进行截面重建；基于梯度下降，对重建的截面进行处理，得到处理后的重建截面；根据处理后的重建截面，确定该区域的等效原子序数和密度；以及基于所确定的等效原子序数和特征密度，确定该区域是否包含违禁品。

Description

多视角成像数据处理方法和设备

技术领域

本发明涉及安检技术领域，尤其涉及一种对多视角成像数据进行处理的方法和设备。

背景技术

对于行李安全检查来说，X射线安检技术是目前应用最为广泛的一类非接触式安检技术。早期使用的X射线检查设备多为采用单视角单能量的X射线技术。随着X射线技术的不断发展，已经出现了可以估计有效原子序数的双能X射线技术。另外，为了满足高探测精度的安检需求，X射线-断层扫描(X-CT)技术已经被引入到安检领域。但是，X-CT技术需要处理的数据量非常巨大，获得图像的时间较长，制造成本较高。

多视角设备是介于普通单视角设备与CT设备之间的设备，例如，视角一般在3个以上、6个以下。虽然多视角设备不像CT那样可以通过360度扫描的完整数据做到准确重建，但可以利用有限几个视角的数据对被检物的3D形状做出部分重建，并可以在部分重建的基础上能估计出被检物的等效原子序数和密度，能部分克服物体之间相互遮挡的问题，从而达到对违禁物品更加准确的检测和识别。

多视角设备除了具有单视角产品的优点以外，还能够对被检物中的危险物质做出更准确的检测和识别，降低对安检员的依赖，检测速度也更高，其成本低于CT，通道尺寸和检测速度不像CT那样受到限制，有很大的市场前景。

发明内容

本公开的目的至少部分地在于提供一种用于多视角成像数据处理的方法和设备，能够更好地从多视角成像数据中进行图像重建。

根据本公开的一个方面，提供了一种对多视角成像数据进行处理的方法，其中，通过X射线从两个或更多个视角对被检物体进行照射得到所述多视角成像数据，该方法包括：从多视角成像数据中提取代表被检物体中同一区域的数据；根据所提取的数据，对该区域进行截面重建；基于梯度下降，对重建的截面进行处理，得到处理后的重建截面；根据处理后的重建截面，确定该区域的等效原子序数和密度；以及基于所确定的等效原子序数和特征密度，确定该区域是否包含违禁品。

根据本公开的实施例，可以通过双能X射线，获得所述多视角成像数据。此时，可以更好地获得材料特性。例如，从多视角成像数据中提取代表被检物体中同一区域的数据可以包括：基于基材料系数表，根据多视角成像数据得到基材料分解系数。

根据本公开的实施例，可以利用代数重建技术(ART)来进行截面重建。例如，可以将待重建的截面离散化为N个单元格像素x_j，其中j＝1，2，...，N，且可以根据如下公式进行重建：

其中，k＝0，1，...，K是ART的迭代次数，x(k)_j是第k次迭代后单元格x_j的值，x(k+1)_j是第k+1次迭代后单元格x_j的值，x(k+1)₂，...，x(k+1)_N}代表第k+1次迭代后重建的截面，λ为松弛参数，p_i为射线i对应的投影值，投影系数w_ij是射线i穿过单元格x_j的长度。

根据本公开的实施例，可以利用全变差(TV)算法来对重建的截面进行处理。例如，可以根据如下公式对重建的截面进行处理：

其中，l＝0，1，...，L是TV的迭代次数，x(0，k+1)是ART中第k+1次迭代后重建的截面x(L，k+1)作为处理后的重建截面并作为ART中第k+2次迭代的初值，α是控制梯度下降速率的变量，d表示使用ART第k+1次迭代之前和之后的重建截面之间的差别，表示为 是归一化的TV梯度，即

根据本公开的实施例，该方法还可以包括：对处理后的重建截面中的像素进行2-均值分类。

根据本公开的另一方面，提供了一种成像数据处理设备，用于对多视角成像数据进行处理，其中，通过X射线从两个或更多个视角对被检物体进行照射得到所述多视角成像数据，该成像数据处理包括：通信接口，用于接收多视角成像数据；存储器，用于存储接收的多视角成像数据；以及处理器，用于：从多视角成像数据中提取代表被检物体中同一区域的数据；根据所提取的数据，对该区域进行截面重建；基于梯度下降，对重建的截面进行处理，得到处理后的重建截面；根据处理后的重建截面，确定该区域的等效原子序数和特征密度；以及基于所确定的等效原子序数和特征密度，确定该区域是否包含违禁品。

根据本公开的实施例，该成像数据处理设备还可以包括：显示器，用于显示重建结果。

根据本公开的实施例，处理器可以如上所述利用代数重建技术(ART)来进行截面重建，并可以利用全变差(TV)算法来对重建的截面进行处理。

根据本公开的又一方面，提供了一种计算机可读记录介质，其上存储有可执行指令，所述指令在由一个或多个处理器执行时，使所述一个或多个处理器执行以下操作：从多视角成像数据中提取代表被检物体中同一区域的数据；根据所提取的数据，对该区域进行截面重建；基于梯度下降，对重建的截面进行处理，得到处理后的重建截面；根据处理后的重建截面，确定该区域的等效原子序数和特征密度；以及基于所确定的等效原子序数和特征密度，确定该区域是否包含违禁品。

根据本公开的实施例，在截面重建时，基于梯度下降，对重建的截面进行处理。这样，可以使得像素之间更为平滑。

附图说明

通过以下参照附图对本公开实施例的描述，本公开的上述以及其他目的、特征和优点将更为清楚，在附图中：

图1是示意性示出了根据本公开实施例的安检系统的透视图；

图2(a)是示意性示出了根据本公开实施例的辐射源和探测器布局的透视图；

图2(b)是示意性示出了根据本公开实施例的辐射源和探测器布局的平面图；

图3A是示意性示出了根据本公开实施例的多视角成像数据处理方法的流程图；

图3B是示意性示出了根据本公开实施例的截面重建方法的流程图；

图4是示出了代数重建技术(ART)原理的示意图；

图5是示意性示出了根据本公开实施例的多视角成像数据处理设备的框图。

贯穿附图，相同或相似的附图标记表示相同或相似的部件。

具体实施方式

以下，将参照附图来描述本公开的实施例。但是应该理解，这些描述只是示例性的，而并非要限制本公开的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本公开的概念。

在此使用的术语仅仅是为了描述具体实施例，而并非意在限制本公开。这里使用的词语“一”、“一个(种)”和“该”等也应包括“多个”、“多种”的意思，除非上下文另外明确指出。此外，在此使用的术语“包括”、“包含”等表明了所述特征、步骤、操作和/或部件的存在，但是并不排除存在或添加一个或多个其他特征、步骤、操作或部件。

在此使用的所有术语(包括技术和科学术语)具有本领域技术人员通常所理解的含义，除非另外定义。应注意，这里使用的术语应解释为具有与本说明书的上下文相一致的含义，而不应以理想化或过于刻板的方式来解释。

附图中示出了一些方框图和/或流程图。应理解，方框图和/或流程图中的一些方框或其组合可以由计算机程序指令来实现。这些计算机程序指令可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器，从而这些指令在由该处理器执行时可以创建用于实现这些方框图和/或流程图中所说明的功能/操作的装置。

因此，本公开的技术可以硬件和/或软件(包括固件、微代码等)的形式来实现。另外，本公开的技术可以采取存储有指令的计算机可读介质上的计算机程序产品的形式，该计算机程序产品可供指令执行系统使用或者结合指令执行系统使用。在本公开的上下文中，计算机可读介质可以是能够包含、存储、传送、传播或传输指令的任意介质。例如，计算机可读介质可以包括但不限于电、磁、光、电磁、红外或半导体系统、装置、器件或传播介质。计算机可读介质的具体示例包括：磁存储装置，如磁带或硬盘(HDD)；光存储装置，如光盘(CD-ROM)；存储器，如随机存取存储器(RAM)或闪存；和/或有线/无线通信链路。

图1是示意性示出了根据本公开实施例的安检系统的透视图。

如图1所示，根据该实施例的安检系统100可以包括安检通道110。被检物体120例如可以通过传送装置(例如，传送带)穿过安检通道110。安检通道110内可以布设有辐射源和探测器。在被检物体120行进通过安检通道110时，辐射源发出的辐射(通常为X射线)照射在被检物体120上，且探测器可以探测透射过被检物体120的辐射(可选地，还可以探测由被检物体120散射的辐射)。

根据本公开的实施例，可以布设多个(例如，两个或更多)辐射源以及相应的多个探测器，以便从多个(例如，两个或更多个)视角对被检物体120进行扫描成像。以下，将进一步详细描述辐射源和探测器的布置。

安检系统100还可以包括控制装置130。控制装置130(例如，计算机)可以连接(有线或无线方式)到安检通道110，以便控制安检通道中的部件如传输装置、辐射源和探测器等，以及对探测器探测到的数据进行处理等。根据本公开的实施例，控制装置130可以对多视角成像数据进行处理，并重建被检物体120(或其一部分)的图像。重建的图像可以显示在输出装置(例如，显示器)上，以供安检人员查看。另外，如果识别到违禁品如爆炸物、武器等，还可以输出警报(例如，在显示器上高亮显示，或者发出提示音)。

这种安检系统100可以适用于车站、机场、海关等需要对物品进行检查以便防止违禁品进入的场所。

图2(a)是示意性示出了根据本公开实施例的辐射源和探测器布局的透视图，且图2(b)是示意性示出了根据本公开实施例的辐射源和探测器布局的平面图。

如图2(a)和2(b)所示，在安检通道110内，布设了辐射源V1、V2、V3和相应的探测器112-1、112-2和112-3。

辐射源V1、V2和V3可以是X射线发生器，并包括相应的驱动电路。另外，与各辐射源V1、V2和V3相对应，可以设置探测器112-1、112-2和112-3。这些探测器例如可以是L型探测器。各辐射源可以发射具有一定张角的扇形束，且相应的探测器可以与该扇形束相对，以便至少部分地接收透射辐射。这样，辐射源与相应的探测器一起限定了辐射源的视角。安检通道110内所布设的所有辐射源的视角一起可以实现对被检物体的大致360°覆盖。

各辐射源和探测器对(例如，(V1，112-1)、(V2，112-2)、(V3，112-3))可以位于平面内(探测视野也在该平面内)。各辐射源和探测器对所在的平面可以彼此平行且相互隔开。例如，这些平面可以是与被检物体120的行进方向(安检通道110的纵向方向，即z轴)实质上正交的平面(x-y平面)，且沿着z轴方向彼此隔开。

如图2(a)和2(b)所示，在该示例中，辐射源V1可以设置于安检通道110的左下方，且相应的探测器112-1可以设置于安检通道110的右侧壁和顶壁上；辐射源V2可以设置于安检通道110的右下方，且相应的探测器112-2可以设置于安检通道110的左侧壁和顶壁上；辐射源V3可以设置于安检通道110的右侧，且相应的探测器112-3可以设置于安检通道110的左侧壁和顶壁上。当然，辐射源和探测器可以按不同方式来布置。

根据本公开的实施例，可以实现双能扫描。例如，辐射源V1、V2和V3中每一个均可以辐射高、低能两种X射线，相应的探测器112-1、112-2和112-3可以包括探测相应X射线的探测阵列。双能扫描有助于识别材料信息。

探测器112-1、112-2和112-3可以探测从被检物体120透射的辐射，并可以将探测到的辐射转换为电信号。这种探测信号可以送到控制装置130进行处理。由于在多个视角进行扫描成像，在此将这种探测信号称为“多视角成像数据”。随着被检物体120行进通过安检通道110，可以获得该被检物体120不同截面处的多视角扫描成像数据。根据这些多视角扫描成像数据，可以重建被检物体120的各截面，并可以由此重建被检物体120的三维轮廓。

在该示例中，示出了三个辐射源V1、V2和V3(以及相应的探测器112-1、112-2和112-3)，从而实现三个视角下的扫描成像。但是，本公开不限于此。例如，也可以双视角或者更多视角(例如，4-6个视角)下的扫描成像。

图3A是示意性示出了根据本公开实施例的多视角成像数据处理方法的流程图。例如，该多视角成像数据处理方法可以由控制装置130来执行。

如图3A所示，根据该实施例的多视角成像数据处理方法300可以包括在操作310处获得多视角成像数据。例如，这些多视角成像数据如上所述是在被检物体120行进通过安检通道110时，由辐射源发出的X射线照射到被检物体120上，并由探测器探测透射的X射线而得到的。控制装置130可以通过有线或无线连接来接收多视角成像数据。

在操作320，可以对多视角成像数据进行图像分割，例如，分割为前景和背景。本领域存在众多技术来进行图像分割。例如，可以通过边缘检测技术来检测被检物体120(例如，行李箱)内各物品的轮廓或者说边缘，从而得到前景图像(并因此可以剔除背景)。针对得到的前景，可以识别其中可能的危险区域。例如，在双能扫描成像的情况下，可以根据材料对高能X射线的吸收和低能X射线的吸收之间的差异，来识别材料特性。例如，可以基于基材料系数表，对多视角成像数据进行基材料分解，得到基材料分解系数，并基于基材料分解系数来识别可能的危险区域。以下的处理可以针对这些识别出的可能危险区域来进行，以便提高处理效率。

在针对这些可能危险区域的数据中，可以提取出不同视角下针对同一区域的数据(例如，经过基材料分解后的数据)。例如，如果一个视角下的某一区域与另一视角下的某一区域存在一定的匹配关系(例如，在几何形状、采集到的高低能数据等方面)，那么可以认为这两个区域事实上是同一区域在不同视角下的表现。

然后，在操作330，可以基于所提取的针对同一区域的不同视角下的数据，对该区域进行截面重建。接下来，将结合图3B进一步详细描述截面重建的处理。

在重建截面之后，在操作340，可以基于重建的截面执行各种处理。例如，可以根据同一区域的多个重建截面，可以构建该区域的三维轮廓(例如，通过串接各截面)。另外，可以根据重建的截面，计算相应区域处的等效原子序数和密度。可以根据所确定的等效原子序数和密度来确定该区域是否包含违禁品。例如，可以对多个截面中的原子序数和密度求平均。如果该平均值落入违禁品数据库中的违禁品区，则可以判定该区域包含违禁品，并可以报警；如果落入违禁品区之外的其他区域，则可以判定该区域为安全品，不予报警。

图3B是示意性示出了根据本公开实施例的截面重建方法的流程图。

如图3B所示，根据该实施例的截面重建方法330可以包括在操作331中根据多视角成像数据(如上所述，针对同一区域的不同视角下的数据)，重建截面中的像素值。在此，该重建可以称作“初步重建”，且重建的像素值可以称作“初步重建(像素)值”。

本领域存在多种截面重建技术。根据本公开的实施例，可以使用代数重建技术(ART)。在ART中，通过多次迭代来重建像素值。

具体地，参见图4，对于待重建的截面421，可以将其离散化为N个单元格像素x_j(j＝1，2，...，N)。对于各单元格，可以通过如下迭代来计算其重建值：

其中，k＝0，1，...，K是ART的迭代次数(可以事先指定，或者可以自适应地确定)，x(k)_j是第k次迭代后单元格x_j的值，x(k+1)_j是第k+1次迭代后单元格x_j的值，代表第k+1次迭代后重建的截面，λ为松弛参数(可以根据经验选取，例如可以根据实际情况提前评估选取，如λ＝0.5)，p_i为射线i对应的投影值(也即，射线i透射过被检物体之后由探测器412探测到的信号，可以是经过了上述的基材料分解之后的数据)，投影系数w_ij是射线i穿过单元格x_j的长度。这些射线i是由辐射源发出的各条X射线，它们(相对于探测器)的角度和位置由物理布局而定。

根据本公开的实施例，在操作333，可以基于梯度下降，对初步重建值进行处理，以在像素之间进行平滑。例如，可以利用全变差(TV)算法，来进行这种处理。

具体地，针对初步重建过程中的每次迭代结果，可以进行TV算法。例如，初步重建的第k+1次迭代后得到的截面可以作为TV算法中的初值，经TV算法处理之后得到的值可以作为初步重建过程中第k+2次迭代的初值，以此类推。

在初步重建如上所述采用ART的情况下，可以如下进行TV算法：

其中，l＝0，1，...，L是TV的迭代次数(可以事先指定，或者可以自适应地确定)，x(0，k+1)是ART中第k次迭代后重建的截面x(L，k+1)作为处理后的重建截面并可用作ART中第k+2次迭代时的初值(即，处理后的重建截面中各单元格的值x(k+1)_j可以代替以上在ART中计算的相应初步重建值x(k+1)_j)，α是控制梯度下降速率的变量(一般而言，α值越小，TV平滑效果越弱，而α值越大，TV平滑效果越强；另一方面，α值与TV的最大迭代次数L之间有相互平衡制约的关系，因此可以根据需求选定α值，例如α＝0.2)，d表示使用ART第k+1次迭代之前和之后的重建截面之间的差别，表示为 是归一化的TV梯度，即(其中，||x(l，k+1)||_TV表示x(l，k+1)的全变差)

需要指出的是，以上截面重建假设物体内部的密度是基本均匀的。即，在该物体的范围内，每个像素具有基本上相同的值；而在该物体范围之外，像素值基本上为零(0)。也就是说，图像上第一物体所在区域的像素值是相同的第一正值，第二物体所在的区域的像素值可以是相同的第二正值(可以与第一正值不同)；物体之间的空白区域(背景)的像素值可以为0。在这种假设下，在截面重建之后，可以对像素值进行2-均值分类(例如，分类为0像素和1像素)。

图5是示意性示出了根据本公开实施例的多视角成像数据处理设备的框图。

如图5所示，根据该实施例的处理设备530(例如，上述控制装置130)可以包括总线531、处理器532、存储器533、输入/输出(I/O)接口534、显示器535和通信接口536。在本公开的各种实施例中，可以省略上述元件中的至少一个，或可以向处理设备530添加其他元件。

总线531可以将上述部件532至536彼此连接，并且在上述部件之间进行通信(例如控制消息和/或数据)。

处理器532例如可以包括央处理单元(CPU)、应用处理器(AP)或者通信处理器(CP)等。处理器532可以执行数据处理或与处理设备530的至少一个其它部件的通信和/或控制相关的操作。

存储器533可以包括易失性存储器和/或非易失性存储器。存储器533可以存储与处理设备530中的至少一个其它部件相关的指令或数据。根据本公开的实施例，存储器533可以存储软件和/或程序。处理器532可以通过执行存储器533中存储的软件和/或程序来执行相应的操作，例如以上结合图3A和3B所述的操作。另外，存储器533中还可以存储有基材料系数表。

I/O接口534可以用于向处理设备530的其他部件传送从用户或其他外部设备输入的指令或数据。此外，I/O接口534可以向用户或其他外部设备输出从处理设备530的其他部件接收到的指令或数据。根据实施例，I/O接口534可以通过输出装置，例如处理音频信号的扬声器或耳机之类的音频设备，进行报警(例如，在存在违禁品的情况下)。

显示器535可以包括例如液晶显示器(LCD)、发光二极管(LED)显示器、有机发光二极管(OLED)显示器等。显示器535可以向例如用户展示各种内容(例如，文本、图像、视频等)。显示器535可以包括触摸屏，并可以接收触摸输入。如上所述，显示器535可以向用户展示重建结果。

通信接口536可以在处理设备530和其他设备(例如，探测器)之间进行通信。例如，通信接口536可以经由无线通信或有线通信与网络如计算机网络(例如局域网(LAN)、广域网(WAN)、互联网)相连，以便与其他设备进行通信。

无线通信可以使用蜂窝通信协议如长期演进(LET)、码分多址(CDMA)、宽带CDMA(WCDMA)或全球移动通信系统(GSM)等。无线通信可以包括例如短距离通信，如WiFi、蓝牙、近场通信(NFC)等。有线通信可以包括通用串行总线(USB)等。

以上对本公开的实施例进行了描述。但是，这些实施例仅仅是为了说明的目的，而并非为了限制本公开的范围。尽管在以上分别描述了各实施例，但是这并不意味着各个实施例中的措施不能有利地结合使用。本公开的范围由所附权利要求及其等价物限定。不脱离本公开的范围，本领域技术人员可以做出多种替代和修改，这些替代和修改都应落在本公开的范围之内。

Claims (15)

1.一种对多视角成像数据进行处理的方法，其中，通过X射线从两个或更多个视角对被检物体进行照射得到所述多视角成像数据，该方法包括：

从多视角成像数据中提取代表被检物体中同一区域的数据；

根据所提取的数据，对该区域进行截面重建；

基于梯度下降，对重建的截面进行处理，得到处理后的重建截面；

根据处理后的重建截面，确定该区域的等效原子序数和密度；以及

基于所确定的等效原子序数和特征密度，确定该区域是否包含违禁品。

2.根据权利要求1所述的方法，其中，通过双能X射线，获得所述多视角成像数据。

3.根据权利要求2所述的方法，其中，从多视角成像数据中提取代表被检物体中同一区域的数据包括：

基于基材料系数表，根据多视角成像数据得到基材料分解系数。

4.根据权利要求1所述的方法，其中，利用代数重建技术ART来进行截面重建。

5.根据权利要求4所述的方法，其中，

将待重建的截面离散化为N个单元格像素x_j，其中j＝1，2，...，N，

且根据如下公式进行重建：

其中，k＝0，1，...，K是ART的迭代次数，x(k)_j是第k次迭代后单元格x_j的值，x(k+1)_j是第k+1次迭代后单元格x_j的值， 代表第k+1次迭代后重建的截面，λ为松弛参数，p_i为射线i对应的投影值，投影系数w_ij是射线i穿过单元格x_j的长度。

6.根据权利要求4或5所述的方法，其中，利用全变差TV算法来对重建的截面进行处理。

7.根据权利要求6所述的方法，其中，根据如下公式对重建的截面进行处理：

其中，l＝0，1，...，L是TV的迭代次数，x(0，k+1)是ART中第k+1次迭代后重建的截面x(L，k+1)作为处理后的重建截面并作为ART中第k+2次迭代的初值，α是控制梯度下降速率的变量，d表示使用ART第k+1次迭代之前和之后的重建截面之间的差别，表示为 是归一化的TV梯度，即

8.根据权利要求1所述的方法，还包括：对处理后的重建截面中的像素进行2-均值分类。

9.一种成像数据处理设备，用于对多视角成像数据进行处理，其中，通过X射线从两个或更多个视角对被检物体进行照射得到所述多视角成像数据，该成像数据处理包括：

通信接口，用于接收多视角成像数据；

存储器，用于存储接收的多视角成像数据；以及

处理器，用于：

从多视角成像数据中提取代表被检物体中同一区域的数据；

根据所提取的数据，对该区域进行截面重建；

基于梯度下降，对重建的截面进行处理，得到处理后的重建截面；

根据处理后的重建截面，确定该区域的等效原子序数和特征密度；以及

基于所确定的等效原子序数和特征密度，确定该区域是否包含违禁品。

10.根据权利要求9所述的成像数据处理设备，还包括：

显示器，用于显示重建结果。

11.根据权利要求9所述的成像数据处理设备，其中，处理器利用代数重建技术ART来进行截面重建。

12.根据权利要求11所述的成像数据处理设备，其中，

处理器将待重建的截面离散化为N个单元格像素x_j，其中j＝1，2，...N，

且根据如下公式进行重建：

其中，k＝0，1，...，K是ART的迭代次数，x(k)_j是第k次迭代后单元格x_j的值，x(k+1)_j是第k+1次迭代后单元格x_j的值， 代表第k+1次迭代后重建的截面，λ为松弛参数，p_i为射线i对应的投影值，投影系数w_ij是射线i穿过单元格x_j的长度。

13.根据权利要求11或12所述的成像数据处理设备，其中，处理器利用全变差TV算法来对重建的截面进行处理。

14.根据权利要求13所述的成像数据处理设备，其中，处理器根据如下公式对重建的截面进行处理：

其中，l＝0，1，...，L是TV的迭代次数，x(0，k+1)是ART中第k+1次迭代后重建的截面x(L，k+1)作为处理后的重建截面并作为ART中第k+2次迭代的初值，α是控制梯度下降速率的变量，d表示使用ART第k+1次迭代之前和之后的重建截面之间的差别，表示为 是归一化的TV梯度，即

15.一种计算机可读记录介质，其上存储有可执行指令，所述指令在由一个或多个处理器执行时，使所述一个或多个处理器执行以下操作：

从多视角成像数据中提取代表被检物体中同一区域的数据；

根据所提取的数据，对该区域进行截面重建；

基于梯度下降，对重建的截面进行处理，得到处理后的重建截面；

根据处理后的重建截面，确定该区域的等效原子序数和特征密度；以及

基于所确定的等效原子序数和特征密度，确定该区域是否包含违禁品。