Heemskerk Jan W T, Defrise Michel
Department of Nuclear Medicine, Universitair Ziekenhuis Brussel, Brussel, Belgium.
EJNMMI Phys. 2020 May 11;7(1):27. doi: 10.1186/s40658-020-00296-w.
For therapeutic applications of several isotopes (e.g., I, Sm, Lu) in nuclear medicine, the high activities typically applied require accurate dead time correction in early time point imaging. We present a novel, straightforward dead time correction method using the Lambert W function, which is in principle exact for the paralyzable detector model with a single parameter τ (i.e., dead time).
As a proof of concept, the method is validated with a simple model: a commonly used isotope, Tc, with a single photopeak. We measured count rates of a gamma camera both intrinsically and extrinsically (i.e., with collimators) with point sources in air and in a scatter phantom (extrinsic only). τ was estimated for both open window (τ) and a Tc photopeak window (τ), using a "graphical" method for fitting the count rate of decaying sources. These values for τ were subsequently used for dead time correction. τ varied significantly between the different geometries for both energy windows, but τ was more reproducible than τ, particularly for the scatter phantom measurements. τ measured from the phantom measurements was approximately 30% lower than τ from the intrinsic measurement but corresponded within 15% with the extrinsic point source measurements. Accordingly, using the intrinsic τ led to an overcorrection of 8% at high count rates; τ from the extrinsic point source measurements corrected the phantom measurement to within 2%. However, significant differences were observed between τ values. All measured τ values underestimated dead time losses in a second independent phantom measurement, with even τ from the first phantom measurement underestimating activity with 5-10% at the highest count rates. Based on measurements of the effect of energy window settings and geometry, we tentatively attribute the added dead time losses to pulse pile-up.
Analytic dead time correction based on the Lambert W function is accurate for the range in which gamma detectors behave as paralyzable systems. However, further investigation indicated measured τ values to be variable with geometry as well as window fraction. We propose that dead time correction should be based on the open window value, τ, corrected for window fraction.
对于核医学中几种同位素(如碘、钐、镥)的治疗应用,早期时间点成像中通常应用的高活度需要精确的死时间校正。我们提出了一种使用兰伯特W函数的新颖、直接的死时间校正方法,该方法原则上对于具有单个参数τ(即死时间)的可麻痹探测器模型是精确的。
作为概念验证,该方法用一个简单模型进行了验证:一种常用的同位素锝,具有单个光电峰。我们在空气中和散射体模(仅外部)中使用点源,对γ相机进行了本征和外部(即使用准直器)计数率测量。使用“图形”方法拟合衰变源的计数率,对开窗(τ)和锝光电峰窗(τ)估计τ。这些τ值随后用于死时间校正。对于两个能量窗,τ在不同几何结构之间有显著差异,但τ比τ更具可重复性,特别是对于散射体模测量。从体模测量得到的τ比本征测量得到的τ低约30%,但与外部点源测量结果相差15%以内。因此,使用本征τ在高计数率下会导致高估8%;外部点源测量得到的τ将体模测量校正到2%以内。然而,观察到τ值之间存在显著差异。在第二次独立的体模测量中,所有测量的τ值都低估了死时间损失,甚至第一次体模测量得到的τ在最高计数率下也低估了活度5 - 10%。基于能量窗设置和几何结构影响的测量,我们初步将额外的死时间损失归因于脉冲堆积。
基于兰伯特W函数的解析死时间校正在γ探测器表现为可麻痹系统的范围内是准确的。然而,进一步研究表明,测量的τ值随几何结构以及窗分数而变化。我们建议死时间校正应基于开窗值τ,并对窗分数进行校正。