CN101109820A - 确定辐射检测器的灵敏度的方法 - Google Patents

确定辐射检测器的灵敏度的方法 Download PDF

Info

Publication number
CN101109820A
CN101109820A CN200710110200.2A CN200710110200A CN101109820A CN 101109820 A CN101109820 A CN 101109820A CN 200710110200 A CN200710110200 A CN 200710110200A CN 101109820 A CN101109820 A CN 101109820A
Authority
CN
China
Prior art keywords
detecting device
signal
inventory
detector
value
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
CN200710110200.2A
Other languages
English (en)
Inventor
M·克雷森斯
W·埃克塞尔曼斯
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.)
Agfa HealthCare NV
Original Assignee
Agfa HealthCare NV
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 Agfa HealthCare NV filed Critical Agfa HealthCare NV
Publication of CN101109820A publication Critical patent/CN101109820A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T7/00Details of radiation-measuring instruments
    • G01T7/005Details of radiation-measuring instruments calibration techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/20Measuring radiation intensity with scintillation detectors
    • G01T1/2012Measuring radiation intensity with scintillation detectors using stimulable phosphors, e.g. stimulable phosphor sheets

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Measurement Of Radiation (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

使得检测器清单的每一个检测器经受完全相同的平场曝光,以便在每一个所述检测器中产生辐射图像,根据所述辐射图像确定总体场分布并且在所述图像中抵消所述场分布,接下来确定所述灵敏度。

Description

确定辐射检测器的灵敏度的方法
技术领域
本发明涉及数字放射线摄影系统的质量保证。
更具体来说,本发明涉及与用在数字放射线摄影系统中的一组检测器相关的质量保证,该组检测器被称作检测器清单(detectorinventory)。
背景技术
所设想的检测器可以是粉末磷光体屏幕或针图像板(needle IP)、直接放射线摄影检测器(被设置成用于直接放射线摄影的无定形硅、无定形硒、Cmos、磷光体检测器等等)或类似检测器。
通过把所述检测器(也被称作“平板”)曝光于X射线场而在该检测器上记录辐射图像。在所谓的读出系统(也被称作“数字化器”)中读出由该检测器临时存储的所述辐射图像,在该读出系统中利用具有适当波长的光对所曝光的检测器进行扫描,其中检测到由该检测器在受到激励时发射的经过图像方面调制的光并且将其转换成表示所述辐射图像的数字图像信号。
磷光体屏幕或针图像板通常被放置在盒带中,并且不是所述读出系统的一部分。这意味着需要对于每一个检测器评估该检测器的均匀性。
均匀性的概念指的是所述检测器在该检测器表面的每一点处相对于平均灵敏度的灵敏度。这一概念被称作固有均匀性。其对于该检测器表面的每一点处的恒定剂量被表示为由所述屏幕发射的平均信号的一个百分比。
对于每一个检测器还可以评估所述灵敏度是否在规范之内,即对于给定的剂量评估由该检测器发射的信号是否在规范之内。
美国医学物理学家协会(简称为AAPM)计算机放射线摄影第10任务组讨论了接收器可再现性、密度均匀性和伪像分析。该组织声明,预期所述固有的接收器之间的均匀性是均质且连贯的。该任务组提出一种用于检测器测试的程序,并且还定义了验收水平。
所规定的程序涉及对检测器清单的所有检测器部件进行曝光,从而把所有检测器相对于入射x射线束的轴集中。在整个平板表面上对各检测器进行均匀地照射。随后相对于预定义的验收水平评估由校准后的读出系统从各检测器读出的信号。
IEC公布了一种用来评估平板间灵敏度变化的测试程序。该测试是为了确认存储磷光体系统的平板间灵敏度变化处于所指定的范围内,并且满足制造商的规范或者规章或合同要求。
所述AAPM组已经规定,在不均匀性方面,10%的峰-峰变化是可以接受的,并且对于灵敏度规范,10%的变化是可以接受的。
本发明的一个目的是提供一种确定参与在检测器清单测试中的各辐射检测器的灵敏度的方法。
发明内容
上述目的是通过一种具有如权利要求1所述的特征的方法实现的。在各从属权利要求中阐述了本发明优选实施例的特定特征。
参照下面的描述和附图,本发明的其他优点和实施例将变得显而易见。
附图说明
图1是数字放射线摄影中的一般设置;
图2示出了本发明的方法的相继步骤;
图3示出了数据信号发生器的结果;
图4以图形的方式示出了仿真的每个检测器的ROI信号;
图5对于多个检测器示出了对应于每个感兴趣区域的信号值的位置;
图6对于多个检测器尺寸示出了不同感兴趣区域的位置的分布;
图7A和7B示出了关于补偿x射线拍摄之间(shot-to-shot)的曝光剂量变化的方法步骤以及用于自动检查不均匀性和/或灵敏度要求的程序;
图8是示出了相对于被取作参考的感兴趣中心区域表示的检测器信号比值的曲线图;
图9是示出了对应于检测器清单中的不同检测器的不均匀性值的曲线图;
图10是示出了对应于检测器清单中的不同检测器的相对灵敏度的曲线图。
具体实施方式
为了能够执行根据本发明的方法,要求所述检测器清单(也被称作平板清单)的所有检测器部件都经受基本上均匀并且基本上完全相同的照射场,该照射场是由对应于给定辐射剂量水平的诸如X射线源之类的辐射源发射的。
在图1中示出了用于照射各检测器以及用于读取存储在检测器中的信号的一般设置。
基于通过在读出设备(也被称作数字化器)中从检测器读出信号所获得的数字图像,对所述检测器清单的各检测器进行分析。在下面假设所述检测器是光激励磷光体屏幕,但是也可以设想其他类型的检测器。
由所述数字化器检测到的信号源自撞击到所述检测器上的x射线的照射水平以及该检测器自身的均匀性和/或灵敏度。由于这些因素具有乘法效果,因此要求所述照射水平在所述磷光体屏幕上的每一个位置处是相等的,以便能够对一组检测器的均匀性和/或灵敏度做出正确的评估。
然而,在实践中,在所述照射场中会发生变化。
这些变化源自这一事实:不是所述检测器表面的每一点都与所述辐射源具有相同的距离。在照射轴上的一点最接近所述辐射源的焦点,并且比起位于该检测器边界处的一点将会自动接收到更多照射。
导致不均匀照射的第二种效应是所谓的足跟效应,其导致在所有方向上评估的所发射x射线的频谱和强度并不均匀。这种效应源自这一事实:由于所述焦点不在所述x射线源的阳极的表面上而是在所述阳极之下,因此在该x射线源的阳极处通过距离差而进行滤波,从而产生强度变化以及频谱变化。一般来说,检测器不仅对于照射量是敏感的,而且对于撞击辐射的频谱也是敏感的。因此,所述足跟效应导致对所述检测器表面的照射的不均匀性。
所述照射场的变化的第三种原因是由当光束落在所添加的滤波器上时的光束距离差而造成的,从而导致垂直落在所述滤波器上的光束相对于有角度地穿过这些滤波器的光束受到不同的衰减。
由于上面描述的效应,实际的输入照射场具有非常复杂的分布。实现一个均匀的进入辐射场以用于对检测器进行照射看起来是不可能的。
根据本发明的一个方面,首先确定所述实际的照射场不均匀性,随后在计算所述检测器的灵敏度和/或均匀性(不均匀性)时抵消所述照射场不均匀性,从而最终结果只包括所述检测器自身的效应,而不包括所述照射场的不均匀性的效应。
在现有技术中,已经尝试通过用一半的剂量对所述检测器执行两次照射来解决上述问题,其中在第一次照射期间把该检测器定位在第一位置处,随后在第二次照射之前将其旋转180度。
通过把从受到两次半剂量照射的该检测器读出的信号相加,至少可以部分地消除所述不均匀照射场的效应。
通过仿真产生关于检测器清单测试的数据
为了获得足够多数量的数据以便执行对灵敏度和/或均匀性的测量,可以通过所谓的感兴趣区域(ROI)信号发生器执行如下所述的数据仿真程序。
对于多个测试检测器,针对该检测器表面上的多个感兴趣区域产生一个信号值。在所描述的示例性实施例中,存在五个感兴趣区域。
可以对应于由上面提到的AAPM组织定义的位置来选择所述区域的位置。
第一区域在所述照射场的光束轴上位于所述检测器的中心处。在检测器的特定尺寸的表面区域内定义四个另外的感兴趣区域的位置,这些感兴趣区域当中的每一个形成该检测器的所述表面区域上的一个象限的中心。
首先产生每个感兴趣区域内的平均图像信号。
在实践中,所述数据不是被仿真,而是从所曝光的各检测器的读出中收集的。利用所述辐射源的给定的电压和mAs设置执行所述曝光。
然而还采用一个剂量计,这是由于从实践中已经知道仅仅依赖于从照射管的控制台读取的电压和mAs值是不够的,因为所述x射线发生器和照射设备存在x射线拍摄间的差异。
所述剂量计被定位在所述照射场中的接近所述检测器但是处在该检测器区域外部的所指定的固定位置处,从而利用该剂量计可以确定所述照射设备的x射线拍摄间的不稳定性。
通过执行这一测量,将有可能对于在检测器的所述五个感兴趣区域当中的每个感兴趣区域内测量的信号执行校正。
在所述数据仿真步骤中,对所述x射线拍摄间的差异进行仿真。
在图3的左上列中示出了由信号发生器产生的数据。所述数据是根据一个仿真模型产生的,其中对于多个检测器,对所述照射设备的给定的x射线拍摄间的变化进行仿真。
该仿真示例的检测器部件的数目等于27,所述各检测器被标记为A到Z、AA、AB和AC。
此外,在第二列中对每一个检测器上的不均匀性的量进行仿真。这在该列中被表示为“局部灵敏度变化,峰-峰增益”。在该示例中,所述局部灵敏度变化是“5”,这意味着假设对于一个图像检测器可以测量到5%的峰-峰信号的变化(灵敏度)。
这被用在所述数据生成模型中,以便自动产生ROI信号。
图3中的第三列包括与不均匀的2D曝光场的仿真相关的数据。所述不均匀性是前面所描述的足跟效应、局部的源到检测器距离以及所使用滤波器的组合效应的结果。其代表被分配给每个ROI的遮蔽(shading)(相对于中心感兴趣区域表示)。
将所述遮蔽与对应于所述检测器的峰-峰增益变化相组合。
这样,获得了适用于根据本发明方法的图像分析的数据。这些数据在图3中的标题为“仿真的图像分析”一列中给出。
图3中示出的表的该部分对于所述中心感兴趣区域和四个另外的感兴趣区域包括一个仿真的图像信号值,该仿真的图像信号值是从对上述效应的仿真得到的。所述信号值被表示为“百分比全尺度线性值”。
所述数字化器提供经过量化的信号,所述经过量化的信号的值在零到最大信号之间,其值取决于量化范围(10比特,12比特,对应于饱和水平的最大值)。
在图4中以图形的方式显示了通过所述仿真得到的结果。
图4的横轴示出了不同的检测器,所述检测器由字母A到Z以及AA、AB和AC表示。
该图的纵轴代表所产生的信号值。
图中的各点代表对应于每个感兴趣区域的信号值的位置。
在该图的右侧示出了平均值。
图5对于多个检测器(H、F、N、V)示出了对应于在横轴上表示的每一个ROI(Q1、Q2、Q3、Q4、中心ROI)的信号值的位置。
此外还示出了对应于所述四个检测器的平均信号分布(见点划线)。
上面描述的仿真得到多个数据信号,所述数据信号可以被用来说明和测试本发明的方法。这些数据充当对本发明的算法的虚拟输入。
应当明白,在本发明的真实世界应用中,使用在检测器清单的各检测器部件上测量的信号而不是仿真的数据。
在图1中示出了所述真实世界设置。
在图2中示出了说明本发明的方法的相继步骤的方框图。
确定及提取不均匀照射场
下面对本发明的关于确定不均匀照射影响的一方面进行描述。
一旦提取了该不均匀照射,将在计算所述各检测器的不均匀性和灵敏度的过程中抵消所测量的分布。
所述测试开始于检测器清单的需要被评估的多个(例如20个、30个)检测器部件。
所述各检测器被曝光于完全相同的几乎均匀的照射场。
所曝光的检测器被输入到一个数字化设备中,在该设备中利用具有适当激励波长的激光对它们进行扫描。在受到激励之后,所述各检测器发光,利用一个光电换能器检测到所述光并且将其转换成信号表示,随后将其转换成数字信号表示。
由上述读出设备(也被称作数字化器)在五个感兴趣区域的位置上读出的信号随后被用于如下所述的评估。在该特定实施例中,所述读出信号被变换成线性信号,并且校正数字化器的偏移。
对于每个检测器定义所述五个感兴趣区域,在这些感兴趣区域(在特定应用中,这些感兴趣区域是位于该检测器的中心处以及位于所述四个象限当中的每一个象限中心处的约为2mm的圆)内计算所述平均信号,并且在这五个数据值的基础上对所述检测器清单测试的每一个检测器部件继续进行评估。
图6示出了不同的感兴趣区域的位置分布。该图中的各点示出了对应于可能的不同检测器格式的感兴趣区域的位置。
在一个特定实施例中,不同的检测器格式相对于所述照射轴被集中。
优选地,各感兴趣区域被定位成使得在所述测试中使用的最大检测器格式的各象限的中心位于作为该测试一部分的最小检测器格式的边界内,从而独立于所述格式,即使在最小格式的检测器中也可以找到五个感兴趣区域。
假设作为所述检测器测试的一部分的所有检测器都具有相同的尺寸,从而在图像中的各感兴趣区域被完全相同地定位。
然而,这对于本发明并不是十分重要。
还有可能在单一检测器清单测试中混合具有不同尺寸的平板。在这种情况下,将对于所有检测器收集关于每个感兴趣区域内的信号值的数据,其中所述感兴趣区域是对于所有所述格式定义的。因此,将对于每一个检测器收集对应于经受所述检测器清单测试的检测器的每一种检测器格式的图像数据(在所描述的实施例中是对应于所述检测器清单测试的每个检测器格式部分的5个信号值)。所有这些数据的组合有助于抵消所述不均匀照射场的效应。
为了能够定义及评估所述不均匀照射场,把每个检测器具有一个平均灵敏度这一事实纳入考虑。
检测器的局部灵敏度由以下各项的组合导致:在粉末磷光体平板的情况下是磷光体层的厚度,或者在针图像板的情况下是针晶体层的组成。
由于源自这种所曝光的检测器的对应于给定照射水平的信号被认为是该检测器的局部灵敏度的函数,并且该局部灵敏度具有正态分布,因此在考虑到大量检测器的情况下,我们应用统计量假设平均而言任何检测器的局部灵敏度都是相同的。
通过把所述检测器清单测试的不同检测器部件上的相应感兴趣区域(其在每个检测器中具有相同的位置)的信号求平均,可以对撞击在所述检测器表面上的所述照射场的不均匀性做出估计。
为了不使得描述复杂化,假设所述检测器清单测试仅仅包括相同格式的检测器。
因此,所描述的实施例中的每个检测器具有五个感兴趣区域,所述感兴趣区域位于每个检测器内的相同地方。
这意味着当抵消所述照射场不均匀性时,仅有源自这些感兴趣区域的信息被抵消。
在实际情况中,并不是每个医院都具有足够多的相同格式的检测器来执行对应于该格式的测试,以便能够对所述不均匀照射的效应执行统计抵消。
如上所述,在一个优选实施例中,所述感兴趣区域的位置使得在特定检测器清单中的最大可用检测器格式的各感兴趣区域仍然被定位在所述最小可用格式的检测器场的外缘尺寸内,从而对于每个检测器的每个图像中的所有格式,可以对于所有检测器格式的所有感兴趣区域找到一个平均信号。
因此,对于经受所述测试的每个(不同格式的)检测器具有一个数字信号,其中可以在对于所有可用检测器格式(从而不仅是对于所评估的特定检测器的格式,而且是对于所述检测器清单中的其他格式)定义的每一个感兴趣区域的位置上从该数字信号中提取所述信号值。
例如,在一个示例性测试中(其中有七个格式可用,并且对于每个检测器取5个感兴趣区域的信号),35项数据将可用于通过如上所述的统计求平均来确定所述不均匀照射分布。
剂量补偿
通过图7中示出的示例性数据说明关于x射线拍摄间的曝光剂量变化的方法步骤。
图7B的左列包括表示对于多个检测器测量的辐射剂量的输入数据,所述检测器被标记为A-Z以及AA、AB和AC。
执行所述剂量测量以便允许消除所述辐射发生器和辐射源的x射线拍摄间的不稳定性。
如果一个检测器接收了过多的辐射,则关于该检测器的结果与接收了较少辐射的检测器的结果不可比较。
为了能够消除该效应,使用通过剂量计测量的反馈。重新计算数据,以便适合实际上恒定的辐射水平。
在第一列下面是包括相同数据但是以大小递增的顺序组织的一列。
在一个实施例中,如下计算所述虚拟平均照射水平:首先确定所述值的中值,随后关于围绕该中值的多个检测器值计算所述各值的平均值。
该平均值被用作该检测器清单测试的虚拟参考剂量。
在该值的基础上,通过一个剂量补偿因数来补偿对应于所述测试中的不同检测器的各单独剂量值。
该剂量补偿因数使得对于每个检测器获得的对应于五个感兴趣区域的各信号被校正,从而补偿由曝光的x射线拍摄间的变化导致的不均匀照射。
当检测器的所测量剂量大于所计算的平均剂量时,所述剂量补偿因数小于1;当对于某一检测器检测到的剂量小于所述平均计算剂量时,该补偿因数大于1。
在该图中的图像分析列具有5列,其中包括被表示为百分比全尺度线性的所测量的经过线性化的平均信号(因此是对于所述数字化器的偏移以及对于该数字化器的平方根或对数特性被校正的数据)。所述五列与位于所述四个象限Q1、Q2、Q3、Q4和所述中心感兴趣区域当中的每一个中心处的感兴趣区域相关。通过前面解释的分析获得这些数据,所述分析是对于每个ROI执行的。该表是二维的,这是因为其涉及到与对应于单一特定格式的各检测器的ROI相关的数据。该表可以具有附加的维度,其包括对应于在所述测试中可用的其他检测器格式的ROI的ROI数据(如前面所解释的那样)。
下面的解释将基于与单一特定盒带格式相关的数据。为了解释方便,假设该测试中的所有检测器具有相同的格式。
为了测量所述剂量,把一个剂量计放置在处于所述照射场中的接近所述检测器但是却在该检测器的图像场外部的一个固定位置处的检测器固定器上。该位置是与检测器格式无关的固定位置。
在一个实施例中,不同格式的所有检测器都被放置在一个固定器上,所述固定器处在相对于所述照射光束轴集中的位置处。这是借助于一组孔洞以及与之协作的参考引脚执行的。
由所述剂量计测量的实际剂量并不精确地对应于在所述中心(所述照射轴的位置)位置处的剂量。
对于该效应,计算一个考虑到距离差的校正。在所述测量中考虑该校正。
由于所述检测器清单测试通常是“相对测试”,这意味着在所述测试中的各检测器的灵敏度是相对于彼此被确定的,因此检测器所接收到的有效实际剂量并不很重要。
然而,在诸如图2(右上)中所示出的图的聚类图中,平均剂量值被表示成一个绝对值。当计算所述各检测器的灵敏度的绝对值时,重要的是考虑到由所评估的检测器接收到的有效实际剂量。
为了测量绝对灵敏度,重要的是知道在所述进入辐射的轴上的剂量。对于所述剂量计的轴外位置及其与所述RX轴上的所述检测器上位置的高度差,将确定一个校正因数(基于距离差)。
该校正因数在图7B中被称作“距离和轴外校正因数”,并且其考虑到在所述检测器上的轴上的剂量水平。
在图7B的表的“图像分析”列中的数据对于特定的检测器格式给出在五个感兴趣区域Q1、Q2、Q3、Q4和所述中心位置上的信号值。这些信号值是原始数据(作为全尺度百分比的线性化的平均信号),其尚未对于以下因素被校正:所述辐射发生器的x射线拍摄间的不稳定性,不同检测器上的局部灵敏度对不同感兴趣区域的影响,作为足跟效应的结果的共模辐射因数,以及场分布中的距离差。
最重要的是在每一个所述原始图像数据信号中提取所述共模辐射变化。
为此目的,在每个检测器的感兴趣区域Q1、Q2、Q3和Q4中的信号值将相对于在所述中心感兴趣区域内找到的信号值被表示。在该中心感兴趣区域内的信号被取作参考并且被设置到“1”,其他信号则相对于该中心感兴趣区域内的该信号而被转换。这导致每个感兴趣区域的信号比值。
参见图7A中的表(相对于中心感兴趣区域的局部信号比值)。
这也在图8的曲线图中示出。该检测器信号比值曲线图对于每个检测器示出了每个感兴趣区域相对于所述中心感兴趣区域的信号比值。
该曲线图对于被标记为H、F、N和V的四个检测器表示每个象限相对于中心ROI(被称作参考)的检测器信号比值。每个检测器的作为感兴趣区域的函数的信号分布相对于在所述中心感兴趣区域中检测到的信号被归一化。
在作为参考位置的中心位置处组合所有曲线图。
对于每个检测器获得一种场信号分布。从而消除所述x射线拍摄间的不稳定性。
获得经过剂量校正的遮蔽分布。该信号遮蔽分布还包括由照射场不均匀性产生的照射遮蔽。
通过对这些信号遮蔽分布求平均,可以确定底层共同照射遮蔽分布是什么。
这是通过对于每个感兴趣区域对与所述检测器清单测试中的每个检测器相关的信号(见图7A的表:相对于中心ROI的局部信号比值)求平均来执行的。
在该表的底部示出了每个感兴趣区域的所述平均值。(平均值或者中值)。该值统计地表示每个感兴趣区域的共模信号比值因数。通过所述求平均,可能存在于每个检测器中的不均匀灵敏度的效应被统计地消除,这是由于对大量检测器取了平均值。
在图7A的表(顶部)的下方对于每个感兴趣区域表示出一个中值平均信号比值,其中所述中心ROI被取作参考。
可以通过连接这些中值平均信号比值而形成的曲线图表示二维辐射遮蔽,这是由于通过上述求平均抵销了所述检测器的不均匀灵敏度的效应。
在图8中示出的检测器信号比值的曲线图对于四个检测器示出了在对所述中心感兴趣区域进行了参考校准之后的相应的信号比值。
如果已经对于所述测试中的所有检测器执行了所述操作,则将已经找到之前计算的所述中值平均值。
在图8中通过点划线示出了表示所述中值遮蔽分布的这些值。
如果在所述测试中有足够多的检测器可用,则所述值将对应于所有检测器的照射的变化。
通过提取该曲线图,获得对应于所述不均匀检测器照射的校正因数,这些值可以在实际计算灵敏度和不均匀性时被用于补偿的目的。
直到现在为止,已经根据以下假设提取了所述不均匀照射遮蔽:如果对于大量检测器计算了所述信号比值的中值(其中所述中心感兴趣区域内的信号被取作参考),则得到所述不均匀照射场的签名。
确定不均匀性
现在必须确定每个检测器的不均匀性。
此外还确定特定检测器的信号对于其每个感兴趣区域与所述不均匀照射遮蔽模式的上面确定的值的差异程度。
这在图2中示出。对于每个检测器,以百分比的方式(+或-)确定针对特定位置找到的ROI信号比值相对于相同位置处的中值的偏差。该偏差对于每个检测器的每个ROI表示在抵销了不均匀照射和x射线拍摄间的差异之后的局部的固有不均匀性。
接下来,对于每个检测器,针对每个ROI局部地确定最大偏差。
在图7B的表中,在标题“不均匀性-与平均曝光遮蔽的偏差,%平均信号”下示出了该最大偏差。
这些值被表示为所述检测器平均信号的百分比,并且代表对于不均匀照射在局部ROI信号比值与参考信号比值之间找到的最大偏差。
在计算所述清单平均不均匀性时,下一步是对于所述检测器清单测试的所有图像部分计算所有所找到的检测器固有最大值的平均值。该值被表示为清单平均不均匀性值(在图7B的表的底部)。
在图7B的表中的所述不均匀性列中找到的值是每个检测器的最大值,其以绝对值的形式被表示为所述检测器平均信号的一个百分比。
在图9中以图形的方式对于每个检测器示出了这些值。在横轴上示出了检测器标识,在纵轴上示出了作为对于每个检测器确定的最大值的百分比绝对偏差。
在该图的右侧示出的条表示所述清单平均不均匀性。其是所述不均匀性的平均最大绝对偏差。
确定灵敏虚
当确定所述灵敏度时,考虑到下面的校正值:对所述辐射源和场分布(遮蔽)的x射线拍摄间的不稳定性进行校正的距离和轴外校正因数。前面已经描述了对于这些值的计算。
在图7B的表中,中间列表示经过剂量和遮蔽校正的线性化的检测器平均信号(以全尺度线性值表示)。
这些值源自每个感兴趣区域的所述原始图像信号。该原始图像信号对于照射的局部不均匀性被校正(图7A的表中示出的因数-不均匀性的中值计算)。对于每个ROI执行该操作。
计算这些经过校正的原始信号值的平均值(由于与5个感兴趣区域相关的数据可用性,因此和被除以5)。随后对于x射线拍摄间的变化校正该平均值,这是通过将该值乘以所述剂量补偿因数。
该动作得到一个对于照射异质性以及x射线拍摄间的变化被校正的平均图像信号。
该信号是一个经过剂量归一化的信号,其表示所评估的检测器的平均灵敏度。
该信号被表示为百分比全尺度线性,并且对于所述检测器清单的所有检测器部件被计算。
这样做对于所述测试中的每个检测器产生经过剂量和曝光遮蔽归一化的信号值。清单平均信号被计算为对于每个检测器计算的所述信号的平均值。该值在图7B的表的底部被表示为清单平均信号。
在图10中示出了标题为“灵敏度”的图。(灵敏度被表示为一个百分比,其参考是零百分比,该参考表示经过剂量校正的各单独的检测器平均信号的平均值)。该图对于被标记为A-Z、AA、AB、AC的所有检测器以百分比的形式表示与(被设置在零百分比的)平均灵敏度信号的偏差。
已经发现,平均来说检测器比所述参考更加灵敏或者更不灵敏。对应于整个检测器清单测试的所有成员的这些检测器平均值的平均值被任意地设置到零百分比。
所述零百分比还指代在图2中的聚类图中示出的零百分比,其中,在单个图中示出了所有检测器关于其在不均匀性以及灵敏度方面的特性的位置。
所述清单平均测试信号是一个百分比线性化的信号,其被用于各检测器平均灵敏度的相对比较。
除了所述聚类图之外,示出所述信号所代表的内容也是令人感兴趣的。
因此,相对于利用所述照射轴上的剂量值计算为每单位剂量的信号量的值示出所述清单平均信号。其代表绝对灵敏度。
在图7B的字段“清单平均灵敏度”中示出了所述值,其中所述清单平均被示为一个绝对值,该绝对值是通过把所述清单平均值除以对于所述各检测器上的中心位置(轴上)计算的剂量而获得的。该值表示每单位剂量的信号量。
到现在为止,已经对于整个清单确定了平均不均匀性和平均灵敏度。
所述平均灵敏度信号已经被转化成一个相对于所述平均剂量(在各检测器的中心处计算的平均剂量)的绝对值。(如图2的聚类图中所示)。
在详细描述了本发明的优选实施例之后,本领域技术人员将会明白,在不脱离由所附权利要求书所限定的本发明范围的情况下可以对本发明进行许多修改。

Claims (7)

1.一种检测参与在检测器清单测试中的各辐射检测器的灵敏度的方法,包括以下步骤:
-使得每一个所述检测器经受基本上完全相同的平场辐射曝光,以便在该检测器中产生辐射图像;
-利用剂量计对于每一个所述辐射图像测量该检测器被曝光的相应的辐射剂量;
-对于每一个所述检测器产生表示在预定义数目的感兴趣区域(ROI)处的所述辐射图像的多个信号(ROI信号),针对所述平场曝光在检测器表面上的不均匀特性校正所述ROI信号;
-利用所述校正后的ROI信号对于每个检测器计算检测器信号;
-还针对由所测量的辐射剂量的差异所反映的所述曝光设备的x射线拍摄间的辐射输出变化率校正所述检测器信号;
-通过对所述校正后的检测器信号求平均对于所述清单计算清单平均信号;
-把所述相对检测器灵敏度确定为所述校正后的检测器信号与所述清单平均信号的偏差。
2.根据权利要求1的方法,其中,定义对于所述相对检测器灵敏度的验收标准,并且评估所述相对检测器灵敏度是否满足所述验收标准。
3.根据权利要求1的方法,修改之处在于,替代所述相对检测器灵敏度,通过以下步骤确定检测器清单平均灵敏度:
-基于所述剂量测量数据计算在所述检测器平面处的轴上清单平均剂量,其中所述剂量测量数据是在参与所述清单测试的各检测器的辐射曝光期间收集的;
-把所述轴上清单平均剂量除以所述清单平均信号。
4.根据权利要求1的方法,其中通过以下步骤针对所述x射线拍摄间的曝光变化率校正所述ROI信号:
-通过计算多个所测量的剂量值的平均值来计算对应于所述清单的平均剂量值;
-通过考虑到对应于每个检测器的实际测量的剂量与所述清单平均剂量值的偏差,对于该检测器计算剂量补偿因数;
-把该剂量补偿因数乘以检测器信号,所述检测器信号是通过对每个检测器的经过不均匀曝光校正的ROI信号求平均而获得的,所述每个检测器的经过不均匀曝光校正的ROI信号是通过把所述ROI信号除以其相应的清单中值ROI信号比值来计算的。
5.根据权利要求4的方法,其中通过以下步骤获得所述平均剂量值:
-对于为每个所述检测器所检测到的剂量值进行排序;
-考虑在排序后的剂量值的中值周围的多个剂量值;
-计算所考虑的各值的平均值,所述平均值表示所述清单平均剂量值;
-对于每个检测器计算对应于该检测器的照射的实际测量的剂量与所述清单平均剂量的偏差;
-对于每个检测器在该偏差的基础上计算剂量补偿因数。
6.根据权利要求1的方法,其中,通过以下步骤产生附加的ROI信号数据:把所述各检测器旋转预定义的角度,把所述各检测器曝光于所述基本上完全相同的平场曝光,以及在所述各感兴趣区域的相同位置处检测ROI信号。
7.根据权利要求1的方法,其中,所述各感兴趣区域位于辐射检测器的中心处以及位于每一个所述象限的中心处。
CN200710110200.2A 2006-06-19 2007-06-18 确定辐射检测器的灵敏度的方法 Pending CN101109820A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP06115649A EP1870732A1 (en) 2006-06-19 2006-06-19 Method of determining the sensitivity of a radiation detector.
EP06115649.3 2006-06-19

Publications (1)

Publication Number Publication Date
CN101109820A true CN101109820A (zh) 2008-01-23

Family

ID=37180052

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200710110200.2A Pending CN101109820A (zh) 2006-06-19 2007-06-18 确定辐射检测器的灵敏度的方法

Country Status (3)

Country Link
US (1) US7629574B2 (zh)
EP (1) EP1870732A1 (zh)
CN (1) CN101109820A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110823937A (zh) * 2018-08-08 2020-02-21 株式会社岛津制作所 电子射线显微分析仪、数据处理方法以及存储介质

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2407106B1 (en) * 2010-07-15 2018-06-27 Agfa HealthCare NV Method of determining the spatial response signature of a detector in computed radiography
US9867581B2 (en) * 2014-01-14 2018-01-16 Siemens Medical Solutions Usa, Inc. Cross-calibration for quantitative functional imaging
US11020083B2 (en) * 2016-09-19 2021-06-01 Siemens Medical Solutions Usa, Inc. Flexible application of cross-calibration for quantitative functional imaging

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6118846A (en) * 1999-02-23 2000-09-12 Direct Radiography Corp. Bad pixel column processing in a radiation detection panel
US6350985B1 (en) * 1999-04-26 2002-02-26 Direct Radiography Corp. Method for calculating gain correction factors in a digital imaging system
DE10327036A1 (de) 2003-06-16 2004-09-09 Siemens Ag Verfahren zur Korrektur der Empfindlichkeit einer Röntgenspeicherleuchtstoffschicht und Röntgenvorrichtung
US20040257352A1 (en) * 2003-06-18 2004-12-23 Nuelight Corporation Method and apparatus for controlling

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110823937A (zh) * 2018-08-08 2020-02-21 株式会社岛津制作所 电子射线显微分析仪、数据处理方法以及存储介质
CN110823937B (zh) * 2018-08-08 2022-07-29 株式会社岛津制作所 电子射线显微分析仪、数据处理方法以及存储介质

Also Published As

Publication number Publication date
US7629574B2 (en) 2009-12-08
US20080251708A1 (en) 2008-10-16
EP1870732A1 (en) 2007-12-26

Similar Documents

Publication Publication Date Title
CN105854191B (zh) 一种放射治疗中三维剂量验证系统及验证方法
CN101297221B (zh) 用于谱计算机断层摄影的方法和设备
EP1081947B1 (en) Correction methods and apparatus for digital x-ray imaging
CN102124372B (zh) 用于确定计数结果的处理电子器件和方法以及用于x射线成像设备的探测器
CN104412127B (zh) 数字正电子发射断层摄影(dpet)能量校准系统及方法
WO2017069286A1 (ja) X線装置、データ処理装置及びデータ処理方法
England et al. Proximal sensing for soil carbon accounting
JPH0410034B2 (zh)
CN100358471C (zh) 放射计算机断层摄像装置及断层图像数据生成方法
CN104729718B (zh) 一种用于红外成像系统netd的处理系统及方法
CN108056785A (zh) 一种平板探测器增益校正模板的更新方法
US8077827B2 (en) Method for thickness calibration and measuring thickness of material
CN106796298B (zh) 数据处理装置、求取各像素的特性的方法以及数据处理的方法及程序
CN101109820A (zh) 确定辐射检测器的灵敏度的方法
US5345082A (en) Scintillation camera utilizing energy dependent linearity correction
JP4994225B2 (ja) 2パラメータスペクトルの処理方法及び装置
CN101271162B (zh) 确定检测器清单的不均匀性的方法
CA1134066A (en) Device for determining internal body structures by means of scattered radiation
CN115266779A (zh) 一种ct扫描用平板探测器像元响应非线性校正方法
CN103728323A (zh) 校准设备和方法
CN103445798A (zh) 确定平板x射线图像探测器中传感器的几何偏移的方法
Youn et al. Energy calibration of energy-resolved photon-counting pixel detectors using laboratory polychromatic x-ray beams
EP1935340B1 (en) Method for neutralizing image artifacts prior to the determination of the Signal-to-noise ratio in CR/DR radiography systems
Thomas et al. GafChromic® RTQA film for routine quality assurance of high-energy photon beams
CN107961028A (zh) 一种归一化校正因子获取、确定方法及医学成像方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Open date: 20080123