CN115298689A - 双能量中的分辨率改进 - Google Patents

双能量中的分辨率改进 Download PDF

Info

Publication number
CN115298689A
CN115298689A CN202180018359.XA CN202180018359A CN115298689A CN 115298689 A CN115298689 A CN 115298689A CN 202180018359 A CN202180018359 A CN 202180018359A CN 115298689 A CN115298689 A CN 115298689A
Authority
CN
China
Prior art keywords
inspection
pixel intensity
data
intensity values
energy mode
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.)
Pending
Application number
CN202180018359.XA
Other languages
English (en)
Inventor
G·杰古
S·格拉博伊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Smith Detection France Ag
Original Assignee
Smith Detection France Ag
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Smith Detection France Ag filed Critical Smith Detection France Ag
Publication of CN115298689A publication Critical patent/CN115298689A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
    • G01V5/20Detecting prohibited goods, e.g. weapons, explosives, hazardous substances, contraband or smuggled objects
    • G01V5/22Active interrogation, i.e. by irradiating objects or goods using external radiation sources, e.g. using gamma rays or cosmic rays
    • G01V5/224Multiple energy techniques using one type of radiation, e.g. X-rays of different energies
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/40Image enhancement or restoration using histogram techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • G06T5/92Dynamic range modification of images or parts thereof based on global image properties
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/50Extraction of image or video features by performing operations within image blocks; by using histograms, e.g. histogram of oriented gradients [HoG]; by summing image-intensity values; Projection analysis
    • G06V10/507Summing image-intensity values; Histogram projection analysis
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/40Extracting pixel data from image sensors by controlling scanning circuits, e.g. by modifying the number of pixels sampled or to be sampled
    • H04N25/46Extracting pixel data from image sensors by controlling scanning circuits, e.g. by modifying the number of pixels sampled or to be sampled by combining or binning pixels
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20216Image averaging
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • G06T2207/30112Baggage; Luggage; Suitcase

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Signal Processing (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Spectroscopy & Molecular Physics (AREA)

Abstract

在一些实施例中,提供了一种用于处理与由检查辐射的多个N个脉冲照射的货物相关联的检查数据的方法,该方法包括:获得所述检查数据,该检查数据表示所述货物的检查图像的像素强度值,所述货物的检查图像是作为多个检测器检测被发射穿过所述货物的所述多个N个脉冲的结果而生成的,以及所述检查数据包括与较高能量模式相关联的数据和与较低能量模式相关联的数据;生成直方图,所述直方图具有作为第一轴的对应于与所述较高能量模式相关联的像素强度值HM的区间和作为第二轴的对应于与所述较低能量模式相关联的像素强度值LM的区间,其中,生成所述直方图包括针对所述检查图像的与所述多个检测器中的每个脉冲i和每个检测器j相对应的每个像素(i,j),对所获得的检查数据中的出现次数进行合并,使得:
Figure DDA0003827509130000011
在所生成的直方图上,针对所述像素强度值LM的每个区间选择对应于所述像素强度值HM的最频繁区间的区间;以及通过将所述像素强度值LM的每个区间与所述像素强度值HM的所选择的区间进行映射来生成转换表。

Description

双能量中的分辨率改进
技术领域
本发明涉及但不限于用于处理与由多个检查辐射脉冲照射的货物相关联的检查数据的方法。本发明还涉及相关联的设备和计算机程序产品。
背景技术
一些加速器具有双能量模式,例如用于材料辨别。加速器用于以低能模式和高能模式辐射货物。对于给定的货物,模式的传输级别是不同的,并且货物的材料可以被确定。
发明内容
本发明的方面和实施例在所附权利要求中阐述。本发明的这些和其它方面和实施例也在此描述。
附图说明
现在将参考附图仅通过示例的方式描述本公开的实施例,在附图中:
图1示意性地表示了根据本公开的示例方法的流程图;
图2示意性地表示了根据本公开的示例直方图;
图3示意性地表示了根据本公开的示例转换表;以及
图4示意性地示出了包括分析器的检查系统,该分析器被配置成至少部分地执行根据本公开的任何方面的方法。
在附图中,相同的附图标记用于表示相同的元件。
具体实施方式
概述
在本公开的一些实施例中,提供了一种用于处理使用较高能量模式数据和较低能量模式数据从货物检查获得的检查数据的方法。该方法可以包括生成直方图,其示出(i)具有与较高能量模式数据相关联的给定强度和(ii)具有与较低能量模式数据相关联的给定强度的像素在检查数据中的出现次数(occurrence)。所生成的直方图可以用于将与较低能量模式数据相关联的每个给定强度关联到与较高能量模式数据相关联的对应的最频繁强度。可以生成映射相关联的强度的转换表。所述转换表可用于通过转换与较低能量模式数据相关联的强度来确定对应于较高能量模式数据的转换强度。
本公开的实施例使得能够结合与较高能量模式数据相关联的可用强度来使用经转换的强度。替代地或另外,本发明的实施例能够在检查期间有效地使获取频率加倍。替代地或另外,本发明的实施例能够在检查期间在扫描方向上使图像分辨率加倍。替代地或另外,本发明的实施例使得能够通过使用组合的强度来改进检查图像的质量。替代地或另外,本发明的实施例实现检查期间的较快扫描。
替代地或另外,如已经陈述,本发明的实施例使得能够使用从货物本身的检查获得的检查数据产生转换表。替代地或附加地,本公开的实施例使得转换表能够保持准确,即使对于货物,检查辐射例如由于检查系统老化、辐射源能量变化、检测器更换、滤波器的添加等而改变。替代地或另外,本发明的实施例可使得能够避免使用不同材料和/或厚度的校准设备来产生转换表的校准的需要。替代地或另外,本发明的实施例可使得能够改进所产生的转换表的可靠性。
示例性实施例的详细描述
图1示意性地表示了根据本公开的示例方法100的流程图。图1的方法100可以用于处理与由检查辐射的多个(N个)脉冲照射的货物相关联的检查数据。
在一些示例中,检查数据可以表示货物的检查图像的像素强度值。检查图像可以是作为多个检测器检测被发射穿过货物的多个(N个)脉冲的结果而生成的。在使用双能量模式的示例中,检查数据包括与较高能量模式相关联的数据和与较低能量模式相关联的数据。
在一些示例中,与较高能量模式相关联的数据包括使用较高能量模式的脉冲(例如,通过生成检查辐射的脉冲的辐射源)生成的数据,并且与较低能量模式相关联的数据包括使用较低能量模式的脉冲生成的数据。在这些示例中,辐射源可以被配置成生成具有较低能量的脉冲,诸如基本上4MeV的能量,以及具有较高能量的脉冲,诸如基本上6MeV的能量。可以设想用于每种模式的其它能量值。
在一些示例中,多个(N个)脉冲包括检查辐射的二元组,每个检查辐射的二元组包括较高能量模式的脉冲和较低能量模式的脉冲。在一些示例中,辐射源可以被配置为交替地生成具有较低能量和较高能量的脉冲,即较低能量的脉冲跟随较高能量的脉冲,反之亦然。替代地或附加地,辐射源可以被配置为生成具有不同于交替序列的序列的较低能量和较高能量的脉冲,即在较低能量的脉冲之前和之后彼此跟随的两个较高能量的脉冲。可以设想较高能量和较低能量的脉冲的其它序列。
替代地或另外地,在一些示例中,与较高能量模式相关联的数据包括使用包括较高能量模式的检测器生成的数据,并且与较低能量模式相关联的数据包括使用包括较低能量模式的检测器生成的数据。在这些示例中,辐射源可以被配置为生成具有给定能量的脉冲,诸如包括在4MeV-6MeV之间的能量(可以设想针对给定能量的其他能量值),并且具有给定的较高能量的每个脉冲由两个检测器探测,即,检测器包括较高能量模式并且检测器包括较低能量模式两者。
图1的方法100包括在S1获得检查数据。方法100还包括在S3生成直方图。
图2示意性地表示了根据本公开的示例性直方图200。图2的直方图具有:
作为第一轴,对应于与较高能量模式相关联的像素强度值HM的区间(bin),
作为第二轴,对应于与低能量模式相关联的像素强度值LM的区间。
如图2的右手侧的刻度示意性表示的,在S2生成直方图200包括针对检查图像的与多个检测器的每个脉冲i和每个检测器j对应的每个像素(i,j),将获得的检查数据中的出现次数进行合并(binning),使得:
Figure BDA0003827509110000041
被检查货物的厚度相对于脉冲的发射被假定为在脉冲(i-1)、i和(i+1)之间是恒定的。如上所述,HM(i,j)对应于与较高能量模式相关联的像素(i,j)的像素强度值,在一些示例中,较高能量模式诸如较高能量脉冲i。类似地,LM(i-1,j)和LM(i+1,j)分别对应于与较低能量模式相关联的像素(i-1,j)和(i+1,j)的像素强度值,在一些示例中,较低能量模式诸如较低能量脉冲(i-1)或(1+1)。
表达式:
Figure BDA0003827509110000042
对应于与两个较低模式相关联的强度的平均值,并且校正小的厚度转变。
在产生图2的直方图之后,方法100可进行到在S4在所生成的直方图上为像素强度值LM的每个区间选择对应于像素强度值HM的最频繁区间的区间。图2上的虚线表示S4的步骤。在所产生的直方图200上,对于像素强度值LM0的区间,选择对应于像素强度值HM的最频繁区间的区间HM0。
图1的方法100还包括在S5通过将像素强度值LM的每个区间与像素强度值HM的所选区间映射来生成转换表。图3示意性地表示了根据本公开的示例转换表300。图3的转换表300包括具有像素强度值LM的一列和具有映射像素强度值HM(=LM')的对应列。
本公开的实施例使得能够使用从货物本身的检查获得的检查数据来生成转换表。替代地或附加地,如已经陈述的,本公开的实施例使得转换表能够保持准确,即使对于货物,检查辐射例如由于检查系统老化、辐射源能量变化、检测器替换、滤波器的添加等而改变。替代地或另外,本发明的实施例可使得能够避免使用不同材料和/或厚度的校准设备来产生转换表的校准的需要。替代地或另外,本发明的实施例可使得能够改进所产生的转换表的可靠性。
如图3中所说明,所述方法可使得确定转换的像素强度值LM'可包括在所生成的转换表300中内插映射的区间。例如,对于包括在LM1和LM2之间的像素强度值LM(x)301,可以使用HM1和HM2内插转换的像素强度值LM'(x)303。
图1的方法还可包括在S6对于与较低能量模式相关联的检查图像的每个像素,确定强度值是否大于预定阈值。预定阈值可以基本上等于例如0.01。可以设想预定阈值的其它值。预定阈值可以由检查系统的用户设置,货物的检查在该检查系统上执行。在一些示例中,预定阈值使得能够确保与较低能量模式相关联的强度值不会太低而不能被使用,因为例如具有较低能量的脉冲没有充分穿透货物。
如果在S6确定强度值大于预定阈值,则图1的方法100可包括对于具有大于预定阈值的强度值的每个像素,在S7使用如图3所示的生成的转换表300,确定对应于像素强度值HM的转换像素强度值LM'。
如果在S6确定强度值小于预定阈值,则图1的方法100可包括在S8确定对应于像素强度值HM的平均像素强度值HM',使得:
Figure BDA0003827509110000061
图1的方法100还可包括在S9使用与较高能量模式关联的像素强度值HM以及确定的转换像素强度值LM'和/或确定的平均像素强度值HM'中的至少一者来生成增强的检查图像。
本公开的实施例使得能够结合与较高能量模式数据相关联的可用强度来使用经转换的强度。替代地或另外,本发明的实施例能够在检查期间有效地使获取频率加倍。替代地或另外,本发明的实施例能够在检查期间在扫描方向上使图像分辨率加倍。替代地或另外,本发明的实施例使得能够通过使用组合强度来改进检验图像的质量。替代地或另外,本发明的实施例实现检查期间的较快扫描。
在一些示例中,可使用双曲正切函数,并且双曲正切函数在增强检查图像中提供LM'和HM'的预定使用之间的平滑过渡。LM'和HM'的预定使用可以由检查系统的用户确定。也可以设想除双曲正切函数之外的其它函数。
下面解释使用双曲正切函数的例子。
如上所述,如果在S6确定与LM模式相关联的强度值高于预定阈值,则存在足够的与LM模式相关联的强度值来使用转换表。转换表可用于确定LM'。如果在S6确定与LM模式相关联的强度值低于预定阈值,则不存在与LM模式相关联的足够强度值,且可不使用转换表。在对应于较高模式HM的像素中存在更多的强度。可以使用与如上所述的像素强度值HM对应的平均像素强度值HM'。
设T为预定阈值。S可以是与低能模式相关联的像素强度值。W可以是增强检查图像中LM'的预定使用和HM'的预定使用之间的过渡区的预定宽度。因此,宽度W可以由检查系统的用户选择,在该检查系统上执行货物的检查。LM'是从对应于如上所述的S的转换表计算的强度值,HM'是对应于如上所述的S的平均值。
对于所获得的检查图像中的每个像素强度值S,在增强的检查图像中使用的像素值V可以使得:
Figure BDA0003827509110000071
当像素强度值S远低于阈值T时,我们得到
Figure BDA0003827509110000072
使得HM'可以用于增强的图像。
当像素强度值S远高于阈值T时,我们得到:
Figure BDA0003827509110000073
使得LM'可用于增强图像。
相对于脉冲的发射,要检查的货物可以包括不同的材料或不同的厚度。替代地或附加地,在其上执行货物检查的检查系统的用户可能仅对货物的一部分感兴趣。替代地或附加地,由于检查辐射的脉冲的角度特性,货物的图像可以包括不同的特征。在图1的方法100中,在S1获得数据还可包括可选步骤:在S2在检查数据中选择与图像的部分相关联的至少一个区域。
在S2的选择步骤可以包括在检查数据中选择与货物的部分相关联的至少一个区域,该货物的部分对于多个(N个)脉冲具有恒定的材料。在要检查的货物包括不同材料的情况下,可以选择与图像的部分相关联的至少一个区域,使得可以针对具有恒定材料的每个区域生成直方图。
可替换地或附加地,对于多个N脉冲,图像的所述部分可以与具有恒定厚度的货物的部分相关联。在要检查的货物包括不同厚度的情况下,可以选择与图像的部分相关联的至少一个区域,使得可以针对具有恒定厚度的每个区域生成直方图。
替代地或附加地,可以选择与图像的部分相关联的至少一个区域,使得可以针对检查系统的用户感兴趣的每个区域生成直方图。
可替换地或附加地,图像的所述部分可以与以下中的至少一者相关联:图像的上部、图像的中部、图像的下部。在货物的图像由于检查辐射的脉冲的角度特性而包括不同特征的情况下,可以针对图像的不同部分生成不同的直方图。
在一些示例中,例如,因为没有像素具有在0.01和0.02之间的归一化信号,所以可能是检查数据不完整。在本公开的示例方法中,获得检查数据还可以包括获得另外的数据。在一些示例中,所述另外的数据可以包括表示额外货物的另外的检查图像的像素强度值的检查数据,所述另外的检查图像是作为通过发射多个脉冲穿过额外货物而扫描额外货物的结果而生成的。可替换地或附加地,所述另外的数据可以包括校准数据,所述校准数据表示校准图像的像素强度值。所述校准图像作为通过发射N个脉冲穿过校准设备来扫描校准设备的结果而被生成。
在一些实例中,为了从直方图估计给定LM的最可能的HM,可如下文解释而选择直方图的区间。
区间大小可以被选择为足够小以使得能够确定准确的最可能的HM,并且足够大以包含许多像素。在一些实例中,强度范围的区间大小可取决于所述强度范围内的检查图像的像素的数目。换句话说,对于检查图像中具有大量图像像素的强度范围,区间大小可以较小,而在其他地方可以较大。
图4示出了检查系统1,其包括被配置成至少部分地执行根据本公开的任何方面的方法的分析器5。在图4中,要被检查的货物101位于容器100中。图4的检查系统1可以被配置为根据本公开的任何方面生成检查数据。
分析器5可以被配置成例如通过通信网络6(其可以是有线和/或可以是无线的),从系统1接收检查数据。分析器5通常至少包括处理器和存储器,以便执行根据本公开的示例方法。
检查系统1被配置成检查容器100。在一些非限制性示例中,可以通过将检查辐射从检查辐射源穿过容器100发射到检测器来执行检查。
在图4所示的示例中,通信服务器8可以被配置成通过通信网络7(其可以是有线和/或可以是无线的)与系统1和/或分析器5通信。在一些示例中,通信服务器8可以被配置成提供远程数据管理系统。在一些示例中,服务器8可以包括数据库。数据库可以被配置为存储本公开的任何方面的检查数据和/或另外的数据。
类似地,分析器5可以被配置成存储本公开的任何方面的检查数据和/或另外的数据。
变型和修改
应当理解,检查辐射源可以包括其他辐射源,作为非限制性示例,例如电离辐射源(例如伽马射线或中子)。检查辐射源还可以包括不适于由电源激活的源,例如放射源(例如使用Co60或Cs137)。
作为一种可能性,提供了一种计算机程序、计算机程序产品或计算机可读介质,其包括计算机程序指令以使得可编程计算机执行本文描述的方法中的任何一者或多者。在示例性实施方式中,与这里的分析器5相关的活动的至少一些部分可以以软件来实现。应当理解,如果需要,本公开的软件组件可以以ROM(只读存储器)形式实现。如果需要,软件组件通常可以使用传统技术以硬件实现。
在一些示例中,分析器5和/或通信网络6和/或7的组件可以使用专用应用和硬件。
如本领域技术人员将清楚的,服务器8和/或分析器5不应被理解为单个实体,而是指包括至少处理器和存储器的物理和/或虚拟设备,存储器可以被包括在一个或多个服务器中,所述一个或多个服务器可以位于单个位置或者可以彼此远离以形成分布式网络(诸如“服务器场”,例如使用有线或无线技术)。
在一些示例中,一个或多个存储器元件(例如,数据库和/或处理器的存储器)可以存储用于本文描述的操作的数据。这包括能够存储软件、逻辑、代码或处理器指令的存储器元件,所述软件、逻辑、代码或处理器指令被执行以执行本公开中描述的活动。
处理器可执行与数据相关联的任何类型的指令,以实现本公开中在此详述的操作。在一个示例中,处理器可以将元件或物品(例如,数据)从一个状态或事物转换到另一状态或事物。在另一示例中,本文概述的活动可以用固定逻辑或可编程逻辑(例如,由处理器执行的软件/计算机指令)来实现,并且本文标识的元件可以是某种类型的可编程处理器、可编程数字逻辑(例如,现场可编程门阵列(FPGA)、可擦除可编程只读存储器(EPROM)、电可擦除可编程只读存储器(EEPROM))、包括数字逻辑的ASIC、软件、代码、电子指令、闪存、光盘、CD-ROM、DVD ROM、磁卡或光卡、适于存储电子指令的其他类型的机器可读介质、或其任何合适的组合。
通信网络6和通信网络7可以仅形成一个网络。分析器5接收的数据通常可以在可能的通信网络6和/或7的范围上接收,通信网络6和/或7至少诸如:基于卫星的通信网络;基于电缆的通信网络;基于电话的通信网络;基于移动电话的通信网络;因特网协议(IP)通信网络;和/或基于计算机的通信网络。
在一些示例中,通信网络6和/或7和/或分析器5可以包括一个或多个网络。可以以任何形式提供网络,包括但不限于局域网(LAN)、无线局域网(WLAN)、虚拟局域网(VLAN)、城域网(MAN)、广域网(WAN)、虚拟专用网(VPN)、内联网、外联网、任何其它适当的体系结构或系统、或促进网络中的通信的其任何组合。
容器100可以是任何类型的容器,例如支架、器皿或盒子等。因此,作为非限制性示例,容器100可以是拖车和/或托板(例如欧洲标准、美国标准或任何其他标准的托板)和/或火车货车和/或车辆(诸如卡车、货车和/或汽车和/或火车)的罐(tank)和/或行李箱,和/或容器100可以是“运输容器”(例如罐或ISO容器或非ISO容器或单元装载设备(ULD)容器)。因此,应理解,容器100可以是任何类型的容器,并且因此在一些示例中可以是手提箱。辐射源被配置成引起穿过容器100的壁的材料(通常为钢)对货物101的检查,例如用于货物101的检测和/或识别。
系统1被配置成在检查模式中引起容器100的整体(即,检查整个容器100)或部分(即,仅检查容器的选定部分,例如,通常,当检查车辆时,可能不检查车辆的车厢,而检查车辆的后部部分)的检查。
检查系统1可以是可移动的,并且可以从一个位置运输到另一个位置(系统1可以包括机动车辆)。替代地或附加地,检查系统1可以相对于地面静止并且不能移位。
检查辐射源可以包括X射线发生器。X射线的能量可以在100keV和15MeV之间,剂量可以在2mGy和20Gy(戈瑞)之间。对于移动检查系统,对于例如在40mm到400mm之间(通常例如300mm(12英寸))的钢穿透能力,X射线源的功率可以例如在100keV和9.0MeV之间(通常例如2MeV、3.5MeV、4MeV或6MeV)。对于移动检查系统,剂量可以是例如20mGy至120mGy。对于静态检查系统,例如在300mm到450mm之间(通常例如410mm(16.1英寸))的钢穿透能力,X射线源的功率可以例如在4MeV和10MeV之间(通常例如9MeV)。对于静态检查系统,剂量可以是17Gy。
除了其它常规的电气元件之外,检测器可以包括辐射检测线,例如X射线检测线。检测器还可以包括其它类型的检测器,例如可选的伽马和/或中子检测器(例如适于在X射线检查的同时探测容器100内放射性伽马和/或中子发射材料的存在)。对于移动检查系统,检测器还可包括电动液压吊杆,其可在运输模式中在缩回位置中操作且在检查位置中操作。吊杆可由液压致动器(例如液压缸)操作。对于静态检查系统,检测器还可以包括结构和/或台架。检测线可以安装在吊杆或结构和/或台架上,面向容器100的另一侧上的源。
为了检查容器100,系统1可以包括运动生成设备,使得系统1可以被移位,容器100是静止的(该模式有时被称为“扫描”模式)。替代地或附加地,运动生成设备可以使容器100移位,系统1相对于地面静止。在一些实施例中,吞吐量(throughput)(即,单位时间内容器100和/或图像10的数量)可以是20到30个图像/小时。替代地或附加地,在“通过”模式中,系统1不包括运动生成设备,并且容器相对于系统1运动,系统1相对于地面是静止的。在实施例中,通过模式中的吞吐量可以高于扫描模式中的吞吐量,并且在检查通过列车的情况下,通过模式中的吞吐量可以是例如50至200个图像/小时,或者甚至50至几千个图像/小时(例如,超过1000个图像/小时的吞吐量)。
上述实施例应理解为说明性示例,并且设想了另外的实施例。应当理解,关于任何一个实施例描述的任何特征可以单独使用,或者与描述的其它特征组合使用,并且还可以与任何其它实施例的一个或多个特征组合使用,或者与任何其它实施例的任何组合使用。此外,在不背离所附权利要求中限定的本发明的范围的情况下,也可以采用上面未描述的等同物和修改。

Claims (18)

1.一种用于处理与由检查辐射的多个N个脉冲照射的货物相关联的检查数据的方法,所述方法包括:
获得所述检查数据,
所述检查数据表示所述货物的检查图像的像素强度值,所述货物的检查图像是作为多个检测器检测被发射穿过所述货物的所述多个N个脉冲的结果而生成的,以及
所述检查数据包括:与较高能量模式相关联的数据和与较低能量模式相关联的数据;
生成直方图,所述直方图具有作为第一轴的对应于与所述较高能量模式相关联的像素强度值HM的区间和作为第二轴的对应于与所述较低能量模式相关联的像素强度值LM的区间,
其中,生成所述直方图包括:针对所述检查图像的与所述多个检测器中的每个脉冲i和每个检测器j相对应的每个像素(i,j),对所获得的检查数据中的出现次数进行合并,使得:
Figure FDA0003827509100000011
在所生成的直方图上,针对所述像素强度值LM的每个区间,选择对应于所述像素强度值HM的最频繁区间的区间;以及
通过将所述像素强度值LM的每个区间与所述像素强度值HM的所选择的区间进行映射来生成转换表。
2.根据前述权利要求所述的方法,还包括:
针对与所述较低能量模式相关联的所述检查图像的每个像素,确定所述强度值是否大于预定阈值;以及
针对具有大于所述预定阈值的强度值的每个像素,使用所生成的转换表,确定对应于像素强度值HM的转换的像素强度值LM'。
3.根据前述权利要求所述的方法,还包括:
针对具有小于所述预定阈值的强度值的每个像素,确定对应于像素强度值HM的平均像素强度值HM',使得:
Figure FDA0003827509100000021
4.根据前两项权利要求中任一项所述的方法,其中确定所述转换的像素强度值LM'包括:
在所生成的转换表中内插映射的区间。
5.根据前述权利要求中任一项所述的方法,其中获得所述检查数据还包括:
在所述检查数据中选择与所述图像的部分相关联的至少一个区域。
6.根据前述权利要求所述的方法,其中所述图像的所述部分与以下各项中的至少一者相关联:所述图像的上部、所述图像的中部、所述图像的下部。
7.根据前两项权利要求中任一项所述的方法,其中所述图像的所述部分与对于所述多个N个脉冲具有恒定厚度的所述货物的部分相关联。
8.根据前述三个权利要求中的任一项所述的方法,其中获得所述检查数据还包括:
在所述检查数据中选择与针对所述多个N个脉冲具有恒定材料的所述货物的部分相关联的至少一个区域。
9.根据前述权利要求中任一项所述的方法,还包括:
基于强度范围内的所述检查图像的像素的数量来选择针对所述强度范围的区间大小。
10.根据前述权利要求中任一项所述的方法,其中获得所述检查数据还包括获得另外的数据,所述另外的数据包括以下各项中的至少一者:
表示额外货物的另外的检查图像的像素强度值的检查数据,所述另外的检查图像是作为通过发射多个脉冲穿过所述额外货物而扫描所述额外货物的结果而生成的,和/或
表示校准图像的像素强度值的校准数据,所述校准图像是作为通过发射多个N个脉冲穿过校准设备而扫描所述校准设备的结果而生成的。
11.根据前述权利要求中任一项所述的方法,其中:
与所述较高能量模式相关联的所述数据包括使用较高能量模式的脉冲生成的数据,以及
与所述较低能量模式相关联的所述数据包括使用较低能量模式的脉冲生成的数据。
12.根据前述权利要求所述的方法,其中所述多个N个脉冲包括检查辐射的二元组,每个检查辐射的二元组包括较高能量模式的脉冲和较低能量模式的脉冲。
13.权利要求1-11中任一项的方法,其中:
与所述较高能量模式相关联的所述数据包括:使用包括较高能量模式的检测器生成的数据,以及
与所述较低能量模式相关联的所述数据包括:使用包括较低能量模式的检测器生成的数据。
14.根据从属于权利要求2或3的前述权利要求中任一项所述的方法,还包括使用以下各项来生成增强的检查图像:
与所述较高能量模式相关联的所述像素强度值HM;以及
所确定的转换的像素强度值LM'和/或所确定的平均像素强度值HM'中的至少一者。
15.根据从属于权利要求3时的前述权利要求所述的方法,其中生成所述增强的检查图像包括:针对与所获得的检查图像中的所述较低能量模式LM相关联的每个像素强度值S,使用所述增强的检查图像中的强度值V,使得:
Figure FDA0003827509100000041
其中:T为所述预定阈值,
LM'为从所述转换表计算的对应于S的所转换的像素强度值,
HM'是对应于S的所述平均像素强度值,以及
W是LM'的预定使用和HM'的预定使用之间的过渡区域的预定宽度。
16.一种分析器,包括:
处理器;以及
存储指令的存储器,所述指令在由所述处理器执行时使得所述处理器能够执行用于处理与由检查辐射的多个N个脉冲照射的货物相关联的检查数据的方法,所述方法包括:
获得所述检查数据,
所述检查数据表示所述货物的检查图像的像素强度值,所述货物的检查图像是作为多个检测器检测被发射穿过所述货物的所述多个N个脉冲的结果而生成的,以及
所述检查数据包括:与较高能量模式相关联的数据和与较低能量模式相关联的数据;
生成直方图,所述直方图具有作为第一轴的对应于与所述较高能量模式相关联的像素强度值HM的区间和作为第二轴的对应于与所述较低能量模式相关联的像素强度值LM的区间,
其中,生成所述直方图包括:针对所述检查图像的与所述多个检测器中的每个脉冲i和每个检测器j相对应的每个像素(i,j),对所获得的检查数据中的出现次数进行合并,使得:
Figure FDA0003827509100000051
在所生成的直方图上,针对所述像素强度值LM的每个区间,选择对应于所述像素强度值HM的最频繁区间的区间;以及
通过将所述像素强度值LM的每个区间与所述像素强度值HM的所选择的区间进行映射来生成转换表。
17.根据权利要求16所述的分析器,其中所述方法是根据权利要求2至15中任意项所述的方法。
18.一种计算机程序或计算机程序产品,其包括指令,所述指令在由处理器执行时使得所述处理器能够执行根据权利要求1至15中任意项所述的方法。
CN202180018359.XA 2020-03-03 2021-03-02 双能量中的分辨率改进 Pending CN115298689A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
GB2003046.6A GB2592607B (en) 2020-03-03 2020-03-03 Resolution improvement in dual energy
GB2003046.6 2020-03-03
PCT/GB2021/050521 WO2021176207A1 (en) 2020-03-03 2021-03-02 Resolution improvement in dual energy

Publications (1)

Publication Number Publication Date
CN115298689A true CN115298689A (zh) 2022-11-04

Family

ID=70278653

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202180018359.XA Pending CN115298689A (zh) 2020-03-03 2021-03-02 双能量中的分辨率改进

Country Status (5)

Country Link
US (1) US20230089927A1 (zh)
EP (1) EP4115384A1 (zh)
CN (1) CN115298689A (zh)
GB (1) GB2592607B (zh)
WO (1) WO2021176207A1 (zh)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1730502B1 (en) * 2004-03-01 2014-04-09 Varian Medical Systems, Inc. Dual energy radiation scanning of objects
US7983382B2 (en) * 2008-11-26 2011-07-19 General Electric Company System and method for material segmentation utilizing computed tomography scans
US9128584B2 (en) * 2013-02-15 2015-09-08 Carl Zeiss X-ray Microscopy, Inc. Multi energy X-ray microscope data acquisition and image reconstruction system and method

Also Published As

Publication number Publication date
US20230089927A1 (en) 2023-03-23
GB202003046D0 (en) 2020-04-15
WO2021176207A1 (en) 2021-09-10
GB2592607A (en) 2021-09-08
EP4115384A1 (en) 2023-01-11
GB2592607B (en) 2024-06-19

Similar Documents

Publication Publication Date Title
US20090052622A1 (en) Nuclear material detection system
US7492855B2 (en) System and method for detecting an object
US7116751B2 (en) System and method for resolving threats in automated explosives detection in baggage and other parcels
US10713759B2 (en) Denoising and/or zooming of inspection images
US10379252B2 (en) Methods and devices for estimating weight of an object to be inspected in an inspection system
US8090150B2 (en) Method and system for identifying a containment vessel
Jaccard et al. Automated detection of cars in transmission X-ray images of freight containers
US7839971B2 (en) System and method for inspecting containers for target material
CN113508388A (zh) 使用数据生成的分类器
US20230221266A1 (en) Pulse frequency adjustment
CN115298689A (zh) 双能量中的分辨率改进
CN110023990B (zh) 利用配准检测违法物品
Rogers et al. Reduction of wobble artefacts in images from mobile transmission x-ray vehicle scanners
CN107077728B (zh) 确定图像中的均匀度
US20230334639A1 (en) Detection of liquid
US20230147681A1 (en) Methods and Systems for Performing On-The-Fly Automatic Calibration Adjustments of X-Ray Inspection Systems
US8254676B2 (en) Methods and systems for identifying a thin object
CN110199209B (zh) 散射成像
US20090087012A1 (en) Systems and methods for identifying similarities among alarms
CN118235158A (zh) 去噪和超分辨率
Fu et al. A novel algorithm for material discrimination using a dual energy imaging system

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