CN113552640A - 射线检查系统及散射校正方法 - Google Patents
射线检查系统及散射校正方法 Download PDFInfo
- Publication number
- CN113552640A CN113552640A CN202010253098.7A CN202010253098A CN113552640A CN 113552640 A CN113552640 A CN 113552640A CN 202010253098 A CN202010253098 A CN 202010253098A CN 113552640 A CN113552640 A CN 113552640A
- Authority
- CN
- China
- Prior art keywords
- detector array
- detector
- signal
- scatter
- modules
- 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
Links
- 238000007689 inspection Methods 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000005855 radiation Effects 0.000 claims abstract description 58
- 230000005540 biological transmission Effects 0.000 claims abstract description 4
- 238000001514 detection method Methods 0.000 claims description 76
- 230000035515 penetration Effects 0.000 claims description 28
- 238000003708 edge detection Methods 0.000 claims description 8
- 230000014509 gene expression Effects 0.000 claims description 5
- 229910000831 Steel Inorganic materials 0.000 description 10
- 239000010959 steel Substances 0.000 description 10
- 230000000694 effects Effects 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 239000000203 mixture Substances 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 230000005461 Bremsstrahlung Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000010894 electron beam technology Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 229910001385 heavy metal Inorganic materials 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011368 organic material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
- G01V5/20—Detecting prohibited goods, e.g. weapons, explosives, hazardous substances, contraband or smuggled objects
- G01V5/22—Active interrogation, i.e. by irradiating objects or goods using external radiation sources, e.g. using gamma rays or cosmic rays
-
- 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
-
- 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
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
-
- 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
-
- 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/50—Detectors
- G01N2223/501—Detectors array
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Molecular Biology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本公开涉及一种射线检查系统及散射校正方法。射线检查系统包括:射线源(1),被配置为产生射线束流(10);第一探测器阵列(2),至少部分位于所述射线束流(10)的覆盖范围内;第二探测器阵列(3),与所述第一探测器阵列(2)位于所述射线检查系统的检查对象(4)的同侧,且位于所述射线束流(10)的覆盖范围外,所述第二探测器阵列(3)被配置为接收所述射线束流(10)在透过所述检查对象(4)的过程中的散射信号;和处理器(5),与所述第一探测器阵列(2)和所述第二探测器阵列(3)信号连接,被配置为根据所述散射信号对所述第一探测器阵列(2)的接收信号进行散射校正。
Description
技术领域
本公开涉及射线检查领域,尤其涉及一种射线检查系统及散射校正方法。
背景技术
在X射线检查领域,散射通常是造成图像质量下降和产生错误信号的主要原因。相关技术对散射校正的研究主要是针对于CT型的检查系统的研究,这种检查系统所使用的X射线能量一般在几百keV的级别,且X射线为锥形束。散射校正多采用复杂的算法分析实现。
在其他应用场景中,例如海关中使用的集装箱检查系统,往往使用能量级别更高的X射线,例如兆伏级的X射线。为了减少散射影响,相关技术中通过增加屏蔽和准直器来减少散射信号。
发明内容
在本公开的一个方面,提供一种射线检查系统,包括:
射线源,被配置为产生射线束流;
第一探测器阵列,至少部分位于所述射线束流的覆盖范围内;
第二探测器阵列,与所述第一探测器阵列位于所述射线检查系统的检查对象的同侧,且位于所述射线束流的覆盖范围外,所述第二探测器阵列被配置为接收所述射线束流在透过所述检查对象的过程中的散射信号;和
处理器,与所述第一探测器阵列和所述第二探测器阵列信号连接,被配置为根据所述散射信号对所述第一探测器阵列的接收信号进行散射校正。
在一些实施例中,所述第一探测器阵列包括至少一组第一探测器模块,每组第一探测器模块包括沿第一方向排列的多个第一探测器模块,所述第二探测器阵列包括至少一组第二探测器模块,每组第二探测器模块包括沿第一方向排列的多个第二探测器模块;所述第一方向与所述射线束流的束流平面平行。
在一些实施例中,每组第一探测器模块中的第一探测器模块的数量与每组第二探测器模块的数量相同,且每组第一探测器模块中各个第一探测器模块与每组第二探测器模块中各个第二探测器模块一一对应,且与对应的第二探测器模块在所述第一方向上的位置相同。
在一些实施例中,所述第二探测器阵列包括一组第二探测器模块,该组第二探测器模块位于所述第一探测器阵列在第二方向的一侧;或者所述第二探测器阵列包括两组第二探测器模块,所述两组第二探测器模块分别位于所述第一探测器阵列在第二方向的两侧;其中,所述第二方向与所述射线束流的束流平面垂直。
在一些实施例中,所述第一探测器阵列包括的第一探测器模块与所述第二探测器阵列包括的第二探测器模块在规格上均相同。
在一些实施例中,所述第一探测器阵列包括的第一探测器模块与所述第二探测器阵列包括的第二探测器模块在规格、数量和排列位置中的至少一种上不同。
在一些实施例中,所述第一探测器阵列与所述第二探测器阵列在第二方向上的间距大于所述第一探测器阵列的灵敏区和所述第二探测器阵列的灵敏区中的最小像素尺寸,所述第二方向与所述射线束流的束流平面垂直。
在一些实施例中,所述处理器被配置为对所述第一探测器阵列与所述第二探测器阵列的散射信号检测能力之间的预设关系进行标定。
在一些实施例中,所述射线源产生的射线束流为X射线束流,且X射线束流的电子束能量大于等于1.0MeV,和/或所述射线源产生的射线束流的宽度不超过100mm。
根据本公开的一个方面,提供一种基于前述的射线检查系统的散射校正方法,包括:
接收第一探测器阵列探测的探测信号,并接收第二探测器阵列探测的散射信号;
根据所述散射信号对所述探测信号进行散射校正,以获得校正后的探测信号作为射线束流穿过被检对象时的穿透信号。
在一些实施例中,所述第一探测器阵列包括的第一探测器模块与所述第二探测器阵列包括的第二探测器模块的规格、数量和排列位置上均相同,散射校正的步骤包括:
在所述探测信号中减去所述散射信号,以获得所述穿透信号。
在一些实施例中,散射校正的步骤包括:
根据所述第一探测器阵列与所述第二探测器阵列的散射信号检测能力之间的预设关系对所述散射信号进行换算,以获得换算后的散射信号;
在所述探测信号中减去所述换算后的散射信号,以获得所述穿透信号。
在一些实施例中,还包括:
对所述第一探测器阵列与所述第二探测器阵列的散射信号检测能力之间的预设关系进行标定。
在一些实施例中,所述预设关系的表达式为:
RSSIC=a*RSSIO;
其中,RSSIO为所述第二探测器阵列探测到的散射信号的信号强度,RSSIC为换算后的散射信号的信号强度,a为比例系数。
在一些实施例中,所述比例系数a大于0,且小于1。
在一些实施例中,所述第一探测器阵列包括多个第一探测器模块,所述多个第一探测器模块中的至少一个第一探测器模块包括多个第一探测单元,所述第二探测器阵列包括多个第二探测器模块,所述多个第二探测器模块中的至少一个第二探测器模块包括多个第二探测单元;
在对所述散射信号进行换算时,所述多个第二探测单元中的边缘探测单元对应的比例系数大于所述多个第二探测单元中的非边缘探测单元对应的比例系数。
在一些实施例中,在对所述散射信号进行换算之后,还包括:
对所述换算后的散射信号进行平滑去噪,并以平滑去噪后的散射信号作为所述换算后的散射信号。
因此,根据本公开实施例,将第一探测器阵列至少部分地设在射线束流的覆盖范围内,并在第一探测器阵列的侧方且位于射线束流的覆盖范围外设置第二探测器阵列,以利用第二探测器阵列检测到的散射信号对第一探测器阵列接收到的检测信号进行散射校正,将该检测信号中的散射信号尽量移除,从而尽量减少或消除散射信号对射线检查结果的干扰。
附图说明
构成说明书的一部分的附图描述了本公开的实施例,并且连同说明书一起用于解释本公开的原理。
参照附图,根据下面的详细描述,可以更加清楚地理解本公开,其中:
图1是兆伏级X射线检查系统的电子加速器的韧致辐射能谱图;
图2是不同光子能量的X射线作用于钢铁材料时的质量衰减系数图;
图3是根据本公开射线检查系统的一些实施例的俯视角度的剖面结构示意图;
图4是根据本公开射线检查系统的一些实施例的侧视角度的剖面结构示意图;
图5和图6分别是本公开射线检查系统的一些实施例中第一探测器模块和第二探测器模块的结构示意图;
图7-图9分别是本公开散射校正方法的一些实施例的流程示意图;
图10是本公开射线检查系统的一些实施例中第一探测器阵列中部分探测单元检测到的穿透信号、散射信号及其两者总和的信号强度示意图;
图11是本公开射线检查系统的一些实施例中第一探测器阵列中部分探测单元检测到的散射信号与位于不同位置的第二探测器阵列中部分探测单元的散射信号的信号强度示意图;
图12是本公开射线检查系统的一些实施例中第一探测器阵列接收的探测信号、第二探测器阵列探测到并经换算和平滑去噪后的散射信号以及通过散射信号散射校正后得到的穿透信号的信号强度示意图;
图13的(a)和(b)分别是本公开射线检查系统的一些实施例中第一探测器阵列在未经散射校正时生成的穿透力图像和经过散射校正后的穿透力图像。
应当明白,附图中所示出的各个部分的尺寸并不是按照实际的比例关系绘制的。此外,相同或类似的参考标号表示相同或类似的构件。
具体实施方式
现在将参照附图来详细描述本公开的各种示例性实施例。对示例性实施例的描述仅仅是说明性的,决不作为对本公开及其应用或使用的任何限制。本公开可以以许多不同的形式实现,不限于这里所述的实施例。提供这些实施例是为了使本公开透彻且完整,并且向本领域技术人员充分表达本公开的范围。应注意到:除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对布置、材料的组分、数字表达式和数值应被解释为仅仅是示例性的,而不是作为限制。
本公开中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性,而只是用来区分不同的部分。“包括”或者“包含”等类似的词语意指在该词前的要素涵盖在该词后列举的要素,并不排除也涵盖其他要素的可能。“上”、“下”、“左”、“右”等仅用于表示相对位置关系,当被描述对象的绝对位置改变后,则该相对位置关系也可能相应地改变。
在本公开中,当描述到特定器件位于第一器件和第二器件之间时,在该特定器件与第一器件或第二器件之间可以存在居间器件,也可以不存在居间器件。当描述到特定器件连接其它器件时,该特定器件可以与所述其它器件直接连接而不具有居间器件,也可以不与所述其它器件直接连接而具有居间器件。
本公开使用的所有术语(包括技术术语或者科学术语)与本公开所属领域的普通技术人员理解的含义相同,除非另外特别定义。还应当理解,在诸如通用字典中定义的术语应当被解释为具有与它们在相关技术的上下文中的含义相一致的含义,而不应用理想化或极度形式化的意义来解释,除非这里明确地这样定义。
对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论,但在适当情况下,所述技术、方法和设备应当被视为说明书的一部分。
图1是兆伏级X射线检查系统的电子加速器的韧致辐射能谱图。图2是不同光子能量的X射线作用于钢铁材料时的质量衰减系数图。
在图1和图2中可以看到,兆伏级的X射线能量较高,该电子加速器所能实现的最高的X射线光子能量为6MeV,而平均X射线光子能量为1.5MeV。而当这种兆伏级X射线检查系统的电子加速器发出的X射线作用在钢铁材料时,能量在0.9MeV以上(尤其是1MeV以上)的X射线光子与钢铁材料的主要作用以康普顿效应为主,强于光电效应和电子对效应,因此散射信号更加显著。
由于兆伏级的X射线能量较高,在作用于物体时还能够形成多次散射,使得散射干扰更为严重。对于被检对象的组成、形状非常不固定的情形,例如被检对象可能是重金属或有机材料,可能存在多种形状,都会使得散射算法校正存在一定的困难。
另外,散射的强度与X射线束流宽度相关,束流宽度越大,散射干扰越大,理想的束流宽度应基本等于对应的像素尺寸。但对于兆伏级的X射线检查系统来说,其一般用于检查集装箱、车辆等尺寸在数米的物体,射线源与探测器的距离一般也在数米到十多米,因此在实际使用中,难以实现高精度的准直来准确地约束束流宽度,导致在探测器位置的束流宽度一般要大于像素尺寸。这样进一步增加了散射干扰。对于兆伏级的X射线检查系统来说,通过增加屏蔽和准直器也难以完全避免散射影响。
有鉴于此,本公开实施例提供一种射线检查系统及散射校正方法,能够尽量减少或消除散射信号对检查结果的影响。
如图3所示,是根据本公开射线检查系统的一些实施例的俯视角度的剖面结构示意图。参考图3及图4-图6,在一些实施例中,射线检查系统包括:射线源1、第一探测器阵列2、第二探测器阵列3和处理器5。射线源1被配置为产生射线束流10。在一些实施例中,射线源1可以为X射线源,其产生的射线束流10为X射线束流。该X射线束流的电子束能量可大于等于1.0MeV,即该射线源1为兆伏级的射线源。在另一些实施例中,射线源1也可以为其他射线源,例如γ射线。或者,X射线源产生的X射线束流的电子束能量低于1.0MeV。
在图3和图4中,射线源1可包括电子加速器11和准直器12。射线束流10从电子加速器11的靶点T发出,经过准直器12,形成具有一定宽度且呈扇形的射线束流区域。
参考图3,第一探测器阵列2至少部分位于所述射线束流10的覆盖范围内。在图4中,第一探测器阵列2包括至少一组第一探测器模块20,每组第一探测器模块20包括沿第一方向x排列的多个第一探测器模块20。在一些实施例中,射线源1产生的射线束流10的宽度不超过100mm,例如30mm。对于较大尺寸的被检对象来说,射线束流10所形成的射线束流区域可近似为射线束流10的束流平面,而第一方向x与所述射线束流10的束流平面平行。
参考图3和图4,在一些实施例中,第一方向x与所述射线束流10的束流平面平行,第二方向y与所述射线束流10的束流平面垂直,第三方向z与第一方向x和第二方向y均垂直。
在一组第一探测器模块20中,多个第一探测器模块20可被分别布置成不同的摆角,以适应射线束流10的入射角度。在一些实施例中,多个第一探测器模块20的间距可以相同,也可以不同。参考图4,第一方向x可以与射线检查系统的安装表面垂直,而在另一些实施例中,也可以与射线检查系统的安装表面呈倾斜角度。在一些实施例中,第一探测器阵列可包括两组以上第一探测器模块,各组第一探测器模块中多个第一探测器模块所排列的第一方向x相同或不同,但均与所述射线束流10的束流平面平行。
第一探测器模块20可包括多个探测单元(例如16、32、48个),对于一组第一探测器模块20来说,各个第一探测器模块20内的多个探测单元可以依次进行编号,不同编号的探测单元所接收的检测信号对应了被检对象的不同位置。除此之外,第一探测器模块20还包括与多个探测单元对应的数据采集电路和结构件等,从而与多个探测单元组成具有统一对外接口的结构体。探测单元可包括灵敏体和数据读取电子线路。
在图3和图4中,射线束流10能够穿过被检对象4后到达第一探测器阵列2。第一探测器阵列2由于至少部分位于射线束流10的覆盖范围,因此在第一探测器阵列2接收到的检测信号中既包括穿透信号,也包括散射信号。
这里的穿透信号是指由穿过被检物质、以及他非灵敏介质的物质,并且未与其发生作用的射线粒子或光子(例如X射线光子),在探测单元中的灵敏介质产生的信号。该穿透信号的强度反应了被检物质的组成信息。
散射信号是射线粒子或光子(例如X射线光子)与被检物质、空气等发生了相互作用(主要是康普顿效应),所产生的散射粒子或光子到达探测单元并产生的信号,这种信号可能导致检测结果变劣。
图10是本公开射线检查系统的一些实施例中第一探测器阵列中部分探测单元检测到的穿透信号、散射信号及其两者总和的信号强度示意图。该图基于300mm穿透力的蒙特卡洛模拟计算,其中,穿透力的测试方法参考国家标准《辐射型货物和(或)车辆检查系统(GB/T19211-2015)》中第8.1节。
图10内探测单元编号从278到378对应的位置为穿透力钢板的信号。在图10中,编号278到378的散射信号能够基本占到信号总和的2/3,高于穿透信号,因此对检测结果影响比较显著。并且,钢板的上下边缘受散射较为严重;同时,相邻第一探测器模块交界处也更容易收到散射干扰。
为了减少或消除散射信号的影响,第二探测器阵列3可被设置为与所述第一探测器阵列2位于所述射线检查系统的检查对象4的同侧,且位于所述射线束流10的覆盖范围外。第一探测器阵列2和第二探测器阵列3可均设置在臂架6上,也可分别设置在不同的臂架上。
由于第二探测器阵列3位于射线束流10的覆盖范围之外,因此其不能接收到射线束流10在透过所述检查对象4的过程中的穿透信号,而接收的是射线束流10在透过所述检查对象4的过程中的散射信号。
图11是本公开射线检查系统的一些实施例中第一探测器阵列中部分探测单元检测到的散射信号与位于不同位置的第二探测器阵列中部分探测单元的散射信号的信号强度示意图。经过实验验证,在图11中,0mm对应的曲线是第一探测器阵列中编号256到384的探测单元检测到的散射信号的信号强度曲线,而30mm和60mm对应的曲线分别是相对于第一探测器阵列向第二方向y偏移30mm和60mm的第二探测器阵列中编号256到384的探测单元检测到的散射信号的信号强度曲线。
在图11中,穿透力钢板对应于编号278到378的探测单元,可以看到这三个位置的散射信号强度相差不大,随着第二探测器阵列相对于第一探测器阵列的偏移距离的增加,第二探测器阵列中的探测单元所接收的散射信号有所减小。并且考虑到散射一般为前向散射和多次散射形成的背景散射,其信号强度在空间中的分布一般不会突变,并随着与束流间距的增加,散射强度会平滑降低。
由于第一探测器阵列接收到的散射信号与第二探测器阵列接收到的散射信号存在对应关系,本公开实施例则通过与第一探测器阵列2和第二探测器阵列3信号连接的处理器5进行第一探测器阵列2的接收信号的校正工作,即根据所述散射信号对所述第一探测器阵列2的接收信号进行散射校正。通过第二探测器阵列检测到的散射信号对第一探测器阵列接收到的检测信号进行散射校正,能够将该检测信号中的散射信号尽量移除,从而尽量减少或消除散射信号对射线检查结果的干扰。
所述第二探测器阵列3包括至少一组第二探测器模块30。每组第二探测器模块30包括沿第一方向x排列的多个第二探测器模块30,其排列形式可参考图4中的第一探测器阵列2中的多个第一探测器模块20的排列形式。
为了方便探测器阵列的布置以及散射校正计算,在一些实施例中,每组第一探测器模块20中的第一探测器模块20的数量与每组第二探测器模块30的数量相同。每组第一探测器模块20中各个第一探测器模块20与每组第二探测器模块30中各个第二探测器模块30一一对应,且与对应的第二探测器模块30在所述第一方向x上的位置相同。这样每个第一探测器模块20都有一个与其位置对应且有少量偏移的第二探测器模块30,相应地,第二探测器模块30能够接收到的散射信号能够与对应的第一探测器模块20所接收到的散射信号基本一致。
在另一些实施例中,每组第一探测器模块20中的第一探测器模块20的数量与每组第二探测器模块30的数量不同,或者在所述第一方向x上的对应位置不同,在计算时可通过曲线拟合的方式分别拟合出第一探测器阵列和第二探测器阵列各自接收到的信号曲线,再进行运算。
在图3中,所述第二探测器阵列3包括一组第二探测器模块30,该组第二探测器模块30位于所述第一探测器阵列2在第二方向y的一侧。在另一些实施例中,所述第二探测器阵列3包括两组第二探测器模块30,所述两组第二探测器模块30分别位于所述第一探测器阵列2在第二方向y的两侧。两组第二探测器模块30接收的散射信号可以取均值后再参与校正运算。
参考图3,在一些实施例中,第一探测器阵列2与所述第二探测器阵列3在第二方向y上的间距d(即偏移距离)大于所述第一探测器阵列2的灵敏区和所述第二探测器阵列3的灵敏区中的最小像素尺寸Ps,例如10mm。该间距d可根据射线束流10的覆盖范围的宽度确定,例如取值为2~10倍的Ps,以便既能够避开射线束流10的覆盖,又能尽量少的减小接收到的散射信号的强度。
在一些实施例中,所述第一探测器阵列2包括的第一探测器模块20与所述第二探测器阵列3包括的第二探测器模块30在规格上均相同。这里的规格可以包括模块内探测单元的数量、排列位置、性能等。当第一探测器阵列2包括的第一探测器模块20与所述第二探测器阵列3包括的第二探测器模块30在规格上均相同时,第二探测器阵列3所接收到的散射信号能够与第一探测器阵列3所接收到的散射信号更加接近,从而方便散射校正的计算。例如直接在第一探测器阵列接收的探测信号中减去第二探测器阵列接收的散射信号,以此来近似地消除第一探测器阵列所接收的散射信号,从而获得所述穿透信号。
在另一些实施例中,所述第一探测器阵列2包括的第一探测器模块20与所述第二探测器阵列3包括的第二探测器模块30在规格、数量和排列位置中的至少一种上不同。为了在散射校正时尽量消除探测信号中的散射信号,处理器5可根据所述第一探测器阵列2与所述第二探测器阵列3的散射信号检测能力之间的预设关系对所述散射信号进行换算,以获得换算后的散射信号。然后,在所述探测信号中减去所述换算后的散射信号,以获得所述穿透信号。
这里的预设关系可由处理器进行标定。处理器可在第一探测器阵列2与所述第二探测器阵列3对被测对象4进行探测之前,通过多次实验来确定所述第一探测器阵列2与所述第二探测器阵列3的散射信号检测能力之间的预设关系,并对其进行标定。举例来说,对于探测器模块为不同规格或者距离第一探测器阵列较远的第二探测器阵列来说,通过多次实验可确定第一探测器阵列中各个探测单元接收到的散射信号强度的数值RSSI1,并确定第二探测器阵列中各个探测单元接收到的散射信号强度的数值RSSI1,通过统计来确定两者的信号强度的比值a=RSSI1/RSSI2作为比例系数,从而对基于该比值a的表达式进行标定。
在正常探测时,则可根据该预设关系的表达式RSSIC=a*RSSIO进行换算。RSSIO为第二探测器阵列3探测到的散射信号的信号强度,RSSIC为换算后的散射信号的信号强度。
基于本公开射线检查系统的上述实施例,本公开提供了对应的散射校正方法。图7-图9分别是本公开散射校正方法的一些实施例的流程示意图。参考图7,在一些实施例中,散射校正方法包括步骤100-步骤200。在步骤100中,接收第一探测器阵列2探测的探测信号,并接收第二探测器阵列3探测的散射信号。在步骤200中,根据所述散射信号对所述探测信号进行散射校正,以获得校正后的探测信号作为射线束流10穿过被检对象时的穿透信号。
在本实施例中,步骤100和步骤200的执行主体均为与第一探测器阵列2和第二探测器阵列3信号连接的处理器5。探测信号经过校正后,得到的穿透信号由于去除了散射信号的影响,因此能够获得更加准确的检测结果。
对于第一探测器阵列2包括的第一探测器模块20与所述第二探测器阵列3包括的第二探测器模块30的规格、数量和排列位置上均相同的情形,步骤200中散射校正的步骤可包括:在所述探测信号中减去所述散射信号,以获得所述穿透信号。由于第一探测器模块20和第二探测器模块30的位置接近且性能相同,因此接收到的散射信号也非常接近,因此通过在探测信号中减去散射信号的方式能够近似获得穿透信号。
对于第一探测器阵列2与第二探测器阵列3的探测器模块规格不同、数量不同或排列位置不同等情况,或者考虑到第二探测器阵列3相对于第一探测器阵列2的偏移量对散射信号的影响,在另一些实施例中,参考图8,步骤200中散射校正的步骤可包括:步骤210和步骤230。在步骤210中,根据所述第一探测器阵列2与所述第二探测器阵列3的散射信号检测能力之间的预设关系对所述散射信号进行换算,以获得换算后的散射信号。在步骤230中,在所述探测信号中减去所述换算后的散射信号,以获得所述穿透信号。
为了确定上述预设关系,在一些实施例中,散射校正方法还包括:对所述第一探测器阵列2与所述第二探测器阵列3的散射信号检测能力之间的预设关系进行标定。该预设关系的表达式可为:RSSIC=a*RSSIO;其中,RSSIO为所述第二探测器阵列3探测到的散射信号的信号强度,RSSIC为换算后的散射信号的信号强度,a为比例系数。为防止各个探测单元的一致性、以及信号涨落造成的过度校正,在一些实施例中,比例系数a大于0,且小于1。
参考图5和图6,第一探测器模块20可包括多个第一探测单元,第二探测器模块30可包括多个第二探测单元。考虑到第一探测器模块20中的边缘探测单元21(即多个探测单元中沿多个探测单元排列方向最外端的探测单元)比非边缘探测单元22受到散射信号的影响更加显著,因此在对所述散射信号进行换算时,使所述多个第二探测单元中的边缘探测单元31对应的比例系数大于所述多个第二探测单元中的非边缘探测单元32对应的比例系数。
参考图9,为了降低第二探测器阵列中某一个或几个探测单元的散射信号的涨落(即散射噪声),在一些实施例中,可在图8所示的基础上,增加步骤220,即在步骤210之后,对所述换算后的散射信号进行平滑去噪,并以平滑去噪后的散射信号作为所述换算后的散射信号。并在步骤230中,采用平滑去噪后的散射信号作为所述换算后的散射信号从所述探测信号内去除。
图12是本公开射线检查系统的一些实施例中第一探测器阵列接收的探测信号、第二探测器阵列探测到并经换算和平滑去噪后的散射信号以及通过散射信号散射校正后得到的穿透信号的信号强度示意图。
在图12中分别示出了第一探测器阵列(作为主探测器阵列)的探测信号的信号强度曲线和第二探测器阵列(作为散射探测器阵列)的散射信号的信号强度曲线。在散射校正过程中,使散射信号的信号强度乘以比例系数0.9,再对结果进行平滑处理,然后使探测信号的信号强度减去平滑处理后的结果,从而获得校正后的第一探测器阵列的探测信号。
通过本公开射线检查系统实施例的散射校正,图13的(a)和(b)分别示出了第一探测器阵列在未经散射校正时生成的穿透力图像和经过散射校正后的穿透力图像,以便进行对比。从图13的(a)可看到,图像受散射干扰,钢板上下边缘较亮,同时,钢板背侧的铅块的图像有缺失。而图13的(b)的图像经过散射校正后,钢板的亮暗较为均匀,且铅块的图像较为完整。
至此,已经详细描述了本公开的各实施例。为了避免遮蔽本公开的构思,没有描述本领域所公知的一些细节。本领域技术人员根据上面的描述,完全可以明白如何实施这里公开的技术方案。
虽然已经通过示例对本公开的一些特定实施例进行了详细说明,但是本领域的技术人员应该理解,以上示例仅是为了进行说明,而不是为了限制本公开的范围。本领域的技术人员应该理解,可在不脱离本公开的范围和精神的情况下,对以上实施例进行修改或者对部分技术特征进行等同替换。本公开的范围由所附权利要求来限定。
Claims (17)
1.一种射线检查系统,其特征在于,包括:
射线源(1),被配置为产生射线束流(10);
第一探测器阵列(2),至少部分位于所述射线束流(10)的覆盖范围内;
第二探测器阵列(3),与所述第一探测器阵列(2)位于所述射线检查系统的检查对象(4)的同侧,且位于所述射线束流(10)的覆盖范围外,所述第二探测器阵列(3)被配置为接收所述射线束流(10)在透过所述检查对象(4)的过程中的散射信号;和
处理器(5),与所述第一探测器阵列(2)和所述第二探测器阵列(3)信号连接,被配置为根据所述散射信号对所述第一探测器阵列(2)的接收信号进行散射校正。
2.根据权利要求1所述的射线检查系统,其特征在于,所述第一探测器阵列(2)包括至少一组第一探测器模块(20),每组第一探测器模块(20)包括沿第一方向(x)排列的多个第一探测器模块(20),所述第二探测器阵列(3)包括至少一组第二探测器模块(30),每组第二探测器模块(30)包括沿第一方向(x)排列的多个第二探测器模块(30);所述第一方向(x)与所述射线束流(10)的束流平面平行。
3.根据权利要求2所述的射线检查系统,其特征在于,每组第一探测器模块(20)中的第一探测器模块(20)的数量与每组第二探测器模块(30)的数量相同,且每组第一探测器模块(20)中各个第一探测器模块(20)与每组第二探测器模块(30)中各个第二探测器模块(30)一一对应,且与对应的第二探测器模块(30)在所述第一方向(x)上的位置相同。
4.根据权利要求2所述的射线检查系统,其特征在于,所述第二探测器阵列(3)包括一组第二探测器模块(30),该组第二探测器模块(30)位于所述第一探测器阵列(2)在第二方向(y)的一侧;或者所述第二探测器阵列(3)包括两组第二探测器模块(30),所述两组第二探测器模块(30)分别位于所述第一探测器阵列(2)在第二方向(y)的两侧;其中,所述第二方向(y)与所述射线束流(10)的束流平面垂直。
5.根据权利要求2所述的射线检查系统,其特征在于,所述第一探测器阵列(2)包括的第一探测器模块(20)与所述第二探测器阵列(3)包括的第二探测器模块(30)在规格上均相同。
6.根据权利要求2所述的射线检查系统,其特征在于,所述第一探测器阵列(2)包括的第一探测器模块(20)与所述第二探测器阵列(3)包括的第二探测器模块(30)在规格、数量和排列位置中的至少一种上不同。
7.根据权利要求1所述的射线检查系统,其特征在于,所述第一探测器阵列(2)与所述第二探测器阵列(3)在第二方向(y)上的间距(d)大于所述第一探测器阵列(2)的灵敏区和所述第二探测器阵列(3)的灵敏区中的最小像素尺寸,所述第二方向(y)与所述射线束流(10)的束流平面垂直。
8.根据权利要求1所述的射线检查系统,其特征在于,所述处理器(5)被配置为对所述第一探测器阵列(2)与所述第二探测器阵列(3)的散射信号检测能力之间的预设关系进行标定。
9.根据权利要求1所述的射线检查系统,其特征在于,所述射线源(1)产生的射线束流(10)为X射线束流,且X射线束流的电子束能量大于等于1.0MeV,和/或所述射线源(1)产生的射线束流(10)的宽度不超过100mm。
10.一种基于权利要求1~9任一所述的射线检查系统的散射校正方法,其特征在于,包括:
接收第一探测器阵列(2)探测的探测信号,并接收第二探测器阵列(3)探测的散射信号;
根据所述散射信号对所述探测信号进行散射校正,以获得校正后的探测信号作为射线束流(10)穿过被检对象时的穿透信号。
11.根据权利要求10所述的散射校正方法,其特征在于,所述第一探测器阵列(2)包括的第一探测器模块(20)与所述第二探测器阵列(3)包括的第二探测器模块(30)的规格、数量和排列位置上均相同,散射校正的步骤包括:
在所述探测信号中减去所述散射信号,以获得所述穿透信号。
12.根据权利要求10所述的散射校正方法,其特征在于,散射校正的步骤包括:
根据所述第一探测器阵列(2)与所述第二探测器阵列(3)的散射信号检测能力之间的预设关系对所述散射信号进行换算,以获得换算后的散射信号;
在所述探测信号中减去所述换算后的散射信号,以获得所述穿透信号。
13.根据权利要求12所述的散射校正方法,其特征在于,还包括:
对所述第一探测器阵列(2)与所述第二探测器阵列(3)的散射信号检测能力之间的预设关系进行标定。
14.根据权利要求12所述的散射校正方法,其特征在于,所述预设关系的表达式为:
RSSIC=a*RSSIO;
其中,RSSIO为所述第二探测器阵列(3)探测到的散射信号的信号强度,RSSIC为换算后的散射信号的信号强度,a为比例系数。
15.根据权利要求14所述的散射校正方法,其特征在于,所述比例系数a大于0,且小于1。
16.根据权利要求14所述的散射校正方法,其特征在于,所述第一探测器阵列(2)包括多个第一探测器模块(20),所述多个第一探测器模块(20)中的至少一个第一探测器模块(20)包括多个第一探测单元,所述第二探测器阵列(3)包括多个第二探测器模块(30),所述多个第二探测器模块(30)中的至少一个第二探测器模块(30)包括多个第二探测单元;
在对所述散射信号进行换算时,所述多个第二探测单元中的边缘探测单元(31)对应的比例系数大于所述多个第二探测单元中的非边缘探测单元(32)对应的比例系数。
17.根据权利要求12所述的散射校正方法,其特征在于,在对所述散射信号进行换算之后,还包括:
对所述换算后的散射信号进行平滑去噪,并以平滑去噪后的散射信号作为所述换算后的散射信号。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010253098.7A CN113552640A (zh) | 2020-04-02 | 2020-04-02 | 射线检查系统及散射校正方法 |
PCT/CN2021/081348 WO2021197075A1 (zh) | 2020-04-02 | 2021-03-17 | 射线检查系统及散射校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010253098.7A CN113552640A (zh) | 2020-04-02 | 2020-04-02 | 射线检查系统及散射校正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113552640A true CN113552640A (zh) | 2021-10-26 |
Family
ID=77927804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010253098.7A Pending CN113552640A (zh) | 2020-04-02 | 2020-04-02 | 射线检查系统及散射校正方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113552640A (zh) |
WO (1) | WO2021197075A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114113173A (zh) * | 2021-11-18 | 2022-03-01 | 上海联影医疗科技股份有限公司 | 一种x射线设备、应用于x射线设备中的散射校正方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1192900A2 (de) * | 2000-09-27 | 2002-04-03 | Philips Patentverwaltung GmbH | Vorrichtung und Verfahren zur Erzeugung eines Röntgen-Computertomogramms mit einer Streustrahlungskorrektur |
CN1343883A (zh) * | 2001-09-21 | 2002-04-10 | 清华大学 | 一种X或γ辐射成像检测方法与装置 |
US20030076924A1 (en) * | 2001-10-19 | 2003-04-24 | Mario Arthur W. | Tomographic scanning X-ray inspection system using transmitted and compton scattered radiation |
US20050025278A1 (en) * | 2003-07-29 | 2005-02-03 | Akira Hagiwara | X-ray CT system |
CN101065685A (zh) * | 2004-11-23 | 2007-10-31 | 皇家飞利浦电子股份有限公司 | 在x射线系统中的散射补偿 |
US7396162B1 (en) * | 2007-03-30 | 2008-07-08 | General Electric Company | Scatter correction for CT method and apparatus |
CN101262819A (zh) * | 2005-09-13 | 2008-09-10 | 皇家飞利浦电子股份有限公司 | Ct散射的直接测量及校正 |
WO2011036968A1 (ja) * | 2009-09-28 | 2011-03-31 | 株式会社日立メディコ | X線ct装置 |
CN102313753A (zh) * | 2010-06-30 | 2012-01-11 | 清华大学 | 物品检测设备及其检测方法 |
US20140294147A1 (en) * | 2013-03-15 | 2014-10-02 | Varian Medical Systems, Inc. | Systems and methods for multi-view imaging and tomography |
CN110907481A (zh) * | 2018-09-18 | 2020-03-24 | 同方威视技术股份有限公司 | 一种x射线的检测系统和检测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4118535B2 (ja) * | 2001-07-03 | 2008-07-16 | 株式会社日立メディコ | X線検査装置 |
US7072440B2 (en) * | 2001-10-19 | 2006-07-04 | Control Screening, Llc | Tomographic scanning X-ray inspection system using transmitted and Compton scattered radiation |
DE102004035943B4 (de) * | 2004-07-23 | 2007-11-08 | GE Homeland Protection, Inc., , Newark | Röntgencomputertomograph sowie Verfahren zur Untersuchung eines Prüfteils mit einem Röntgencomputertomographen |
-
2020
- 2020-04-02 CN CN202010253098.7A patent/CN113552640A/zh active Pending
-
2021
- 2021-03-17 WO PCT/CN2021/081348 patent/WO2021197075A1/zh active Application Filing
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1192900A2 (de) * | 2000-09-27 | 2002-04-03 | Philips Patentverwaltung GmbH | Vorrichtung und Verfahren zur Erzeugung eines Röntgen-Computertomogramms mit einer Streustrahlungskorrektur |
CN1343883A (zh) * | 2001-09-21 | 2002-04-10 | 清华大学 | 一种X或γ辐射成像检测方法与装置 |
US20030076924A1 (en) * | 2001-10-19 | 2003-04-24 | Mario Arthur W. | Tomographic scanning X-ray inspection system using transmitted and compton scattered radiation |
US20050025278A1 (en) * | 2003-07-29 | 2005-02-03 | Akira Hagiwara | X-ray CT system |
CN101065685A (zh) * | 2004-11-23 | 2007-10-31 | 皇家飞利浦电子股份有限公司 | 在x射线系统中的散射补偿 |
CN101262819A (zh) * | 2005-09-13 | 2008-09-10 | 皇家飞利浦电子股份有限公司 | Ct散射的直接测量及校正 |
US7396162B1 (en) * | 2007-03-30 | 2008-07-08 | General Electric Company | Scatter correction for CT method and apparatus |
WO2011036968A1 (ja) * | 2009-09-28 | 2011-03-31 | 株式会社日立メディコ | X線ct装置 |
CN102313753A (zh) * | 2010-06-30 | 2012-01-11 | 清华大学 | 物品检测设备及其检测方法 |
US20140294147A1 (en) * | 2013-03-15 | 2014-10-02 | Varian Medical Systems, Inc. | Systems and methods for multi-view imaging and tomography |
CN110907481A (zh) * | 2018-09-18 | 2020-03-24 | 同方威视技术股份有限公司 | 一种x射线的检测系统和检测方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114113173A (zh) * | 2021-11-18 | 2022-03-01 | 上海联影医疗科技股份有限公司 | 一种x射线设备、应用于x射线设备中的散射校正方法 |
CN114113173B (zh) * | 2021-11-18 | 2024-06-18 | 上海联影医疗科技股份有限公司 | 一种x射线设备、应用于x射线设备中的散射校正方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2021197075A1 (zh) | 2021-10-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110072459B (zh) | 用于自校准的自校准ct检测器、系统和方法 | |
EP3242126A1 (en) | Dual-energy ray imaging method and system | |
US10521936B2 (en) | Device and method for image reconstruction at different X-ray energies, and device and method for X-ray three-dimensional measurement | |
US10386532B2 (en) | Radiation signal processing system | |
US7583783B2 (en) | X-ray computer tomograph and method for examining a test piece using an x-ray computer tomograph | |
JP2008510132A (ja) | 放射線検出器用の抗散乱グリッド | |
JP6596414B2 (ja) | X線検出器におけるx線線量パラメータを測定するための装置、及び放射線学的情報を取り扱うためのコンピュータネットワーク | |
CN101960297B (zh) | X射线计算机断层成像装置和采用x射线计算机断层成像检测对象的方法 | |
JP2015021913A (ja) | 放射線エネルギー分布の検出方法及びその検出装置 | |
Osipov et al. | Numerical modeling of radiographic images as the basis for correctly designing digital radiography systems of large-sized objects | |
Udod et al. | State-of-the art and development prospects of digital radiography systems for nondestructive testing, evaluation, and inspection of objects: a review | |
CN113552640A (zh) | 射线检查系统及散射校正方法 | |
Stritt et al. | Performance quantification of a flat-panel imager in industrial mega-voltage X-ray imaging systems | |
JP5946012B2 (ja) | 放射線検査装置 | |
CN112033983A (zh) | 一种x射线发射系统及检测系统 | |
US11872071B2 (en) | Method for correcting a spectral image | |
JP4037590B2 (ja) | 陽電子消滅γ線測定方法および装置 | |
CN212780578U (zh) | X射线发射系统及检测系统 | |
JP7520772B2 (ja) | 粒子線到達位置モニタ、粒子線治療システムおよび粒子線到達位置モニタ方法 | |
CN114113173B (zh) | 一种x射线设备、应用于x射线设备中的散射校正方法 | |
JP7449821B2 (ja) | 内部状態検査システム及び内部状態検査方法 | |
JPH09108207A (ja) | X線検出方法 | |
Chen et al. | Validation of HADES-based simulations of radiographic experiments | |
JP2017142956A (ja) | 電子増幅器および放射線検出器 |
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 |