EJNMMI Physics
Positron range in combination with point-spread-function correction: an evaluation of different implementations for [124I]-PET imaging
Hunor Kertész ORCID: orcid.org/0000-0002-4626-930X1,
Maurizio Conti2,
Vladimir Panin2,
Jorge Cabello2,
Deepak Bharkhada2,
Thomas Beyer1,
Laszlo Papp1,
Walter Jentzen3,
Jacobo Cal-Gonzalez1,4,
Joaquín L. Herraiz5,6,
Alejandro López-Montes5 &
Ivo Rausch1
EJNMMI Physics volume 9, Article number: 56 (2022) Cite this article
To evaluate the effect of combining positron range correction (PRC) with point-spread-function (PSF) correction and to compare different methods of implementation into iterative image reconstruction for 124I-PET imaging.
Uniform PR blurring kernels of 124I were generated using the GATE (GEANT4) framework in various material environments (lung, water, and bone) and matched to a 3D matrix. The kernels size was set to 11 × 11 × 11 based on the maximum PR in water and the voxel size of the PET system. PET image reconstruction was performed using the standard OSEM algorithm, OSEM with PRC implemented before the forward projection (OSEM+PRC simplified) and OSEM with PRC implemented in both forward- and back-projection steps (full implementation) (OSEM+PRC). Reconstructions were repeated with resolution recovery, point-spread function (PSF) included. The methods and kernel variation were validated using different phantoms filled with 124I acquired on a Siemens mCT PET/CT system. The data was evaluated for contrast recovery and image noise.
Contrast recovery improved by 2–10% and 4–37% with OSEM+PRC simplified and OSEM+PRC, respectively, depending on the sphere size of the NEMA IQ phantom. Including PSF in the reconstructions further improved contrast by 4–19% and 3–16% with the PSF+PRC simplified and PSF+PRC, respectively. The benefit of PRC was more pronounced within low-density material. OSEM-PRC and OSEM-PSF as well as OSEM-PSF+PRC in its full- and simplified implementation showed comparable noise and convergence. OSEM-PRC simplified showed comparably faster convergence but at the cost of increased image noise.
The combination of the PSF and PRC leads to increased contrast recovery with reduced image noise compared to stand-alone PSF or PRC reconstruction. For OSEM-PRC reconstructions, a full implementation in the reconstruction is necessary to handle image noise. For the combination of PRC with PSF, a simplified PRC implementation can be used to reduce reconstruction times.
Positron emission tomography (PET) is based on the detection of photon radiation arising from the annihilation of positrons originating from the PET tracers with electrons within the patient [1]. To obtain a 3D representation of the tracer distribution, image reconstruction is necessary. This is done today using iterative reconstruction algorithms where the detector signals are estimated from a guess of the tracer distribution using a mathematical model of the system (system matrix) and compared to the actual measured signal in projection space. Standard reconstruction algorithms, such as OSEM, typically use simple geometric assumptions to model the system matrix [2]. However, to optimize spatial resolution a more realistic model of the system accounting for all physical processes is needed. A realistic system response is, for example, derived from point source measurements [3]. The implementation of such a realistic detector response function into the projectors of the iterative image reconstruction process is called resolution recovery or point-spread-function (PSF) correction. The benefit of including PSF in the image reconstruction by means of reduced image noise and increased contrast was shown by various studies [4,5,6,7].
However, besides the imaging systems properties, spatial resolution in PET is further limited by positron range (PR), the distance traveled by the positron from the emission position to the annihilation position [8, 9]. For the most commonly used PET isotope (18F), median PR in soft tissue or water is of the order of 0.4 mm and thus [10] can be practically ignored in clinical PET studies. However, positron emitters with high positron energies such as 68Ga and 124I, which are increasingly applied in the diagnosis of prostate cancer [11] or used for therapy planning in patients with thyroid cancer [12], respectively, present with average PRs of ~ 10 mm in soft tissues which lead to substantial degradation of spatial resolution and image quality [13].
Positron range effects can be corrected with a similar implementation to the one used to incorporate PSF into the image reconstruction algorithm. This is done by adding the PR-induced blurring to the forward projector, for example, by convolving the current image guess with the PR-based blurring kernels before projecting the current image guess into projection space. In addition, the PR blurring has to be taken into account in the back-projection step within an iterative reconstruction since the back-projection is a transpose operation of the forward projection [14]. The main challenges of PRC are deriving adequate PR kernels and extensive computational demand due to the use of non-isotropic and non-shift-invariant PR kernels. Monte Carlo simulations are usually used for defining the ground truth of tissue-dependent PR kernels [10].
However, due to the computational challenges, PR kernels are generally based on simply using isotropic PR kernels for different tissue types [15,16,17]. Only a limited number of studies used spatially variant and tissue-dependent estimations of PR kernels including effects on tissue borders [15, 18]. These studies suggest implementing the PR blurring kernels only prior to the forward projection step in standard iterative reconstructions as an image blur to reduce the computational demand [15, 19]. Only one recent study showed the value of using a full implementation of the PRC into an OSEM algorithm [14].
Therefore, our study aimed to evaluate the effect of combined PSF correction and PRC within the reconstruction on image contrast and convergence. Further, the influence of using simplified PR correction implementation methods limiting computational demand in combination with PSF was assessed.
Positron range kernel calculation
The uniform positron range distribution of 124I was simulated using GATE 9.0 (GEANT4 10.06.p02). The simulation setup consisted of a point source with a radius of 5 nm centered in a uniform phantom with a radius of 30 cm. For the phantom material, three settings were simulated bone material (mass density 1.92 g/cm3), water (1.00 g/cm3), and lung equivalent material (0.26 g/cm3). The predefined "empenelope" physics list was used. The initial activity of 124I was 10 MBq and the emission and annihilation coordinates were recorded for 20 M annihilation events. From the recorded data, the 3D PR distribution was generated by mapping the distribution to a 3D matrix. The size of the kernel (11 × 11 × 11 voxels) was chosen based on the voxel size of the Siemens mCT PET/CT system (2 × 2 × 2 mm3) (Siemens Medical Solutions Inc., Knoxville, TN, USA) [20] and the maximum PR in water for 124I (~ 10 mm), since the use of larger kernels is leading to artifacts [14].
For each voxel of the investigated object, spatially variant and tissue-dependent PR kernels based on the underlying material distribution were created as described before [14]. In short, the attenuation map was segmented into a material map containing lung, soft tissue, and bone. Regions containing air were treated as lung material. Then, for every voxel, PR kernels were generated by combining pre-calculated uniform PR kernels for lung, water, and bone from GATE Monte Carlo simulations.
Positron range correction implementations
The PR correction was implemented into the vendor-based software in combination with OSEM and PSF algorithms as an additional PR-dependent PSF applied in image space. This was done as full implementation and, in addition, in a simplified version applying the PR-dependent blurring only before the forward projector of the reconstruction algorithm, as suggested by [19] to reduce computational demand.
The standard OSEM algorithm taken as the basis can be written as:
$$f_{j}^{n + 1} = \frac{{f_{j}^{n} }}{{\sum\nolimits_{{i^{\prime} \in S^{n} }} {a_{i^{\prime}j} } }}\sum_{{i \in S^{n} }} a_{ij} \frac{{m_{i} }}{{\sum_{k} a_{ik} f_{k}^{n} }}$$
where fjn+1 is the next image estimate of voxel j based on the current image estimate fkn. The measured projection data is mi. The system matrix describing the probability that the emission from voxel j will be detected along the line of response (LOR) i is given by Aij = (aij)IxJ. Only a subset Sn of the data was used in each update. Of note: for simplification, the random and scatter are not included in this description. The system matrix can be factorized to account for the finite resolution effects, in our case the positron range effect by matrix H = (hj'j)JxJ, a matrix X = (xij)IxJ expressing the intersection length and a matrix W = (wii)IxI describing the photon attenuation and geometrical sensitivity variations. Thus, the system matrix can be written as:
$$A =WXH$$
Combining Eq. (1) and Eq. (2), the OSEM algorithm with PRC correction can be written as [14]:
$${f}_{j}^{n+1}=\frac{{f}_{j}^{n}}{{\sum }_{b}{h}_{bj}{\sum }_{i\in {S}^{n}}{w}_{ii}{x}_{ib}}{\sum }_{b}{h}_{bj}{\sum }_{i\in {S}^{n}}{x}_{ib}\frac{{m}_{i}}{{\sum }_{p}{x}_{ip}{\sum }_{v}{h}_{pv}{f}_{v}^{n}}$$
This method is referred to as "PRC" in this work. Since PRC is computationally highly expensive, a simplified version of the PRC implementation was suggested in previous studies [18, 21], where the PR kernels are applied to the current image estimate as a convolution in image space before the forward projection. In this way, Eq. (1) can be written as:
$$f_{j}^{n + 1} = \frac{{f_{j}^{n} }}{{\sum\nolimits_{{i^{\prime} \in S^{n} }} {a_{i^{\prime}j} } }}\sum_{{i \in S^{n} }} a_{ij} \frac{{m_{i} }}{{\sum_{k} a_{ik} f_{k}^{{n^{ \sim } }} }}$$
where \(f_{k}^{{n^{ \sim } }}\) is the current image estimate blurred with the spatially variant and tissue-dependent PR kernel (\(\rho\)) calculated for the given voxel [19]:
$$f_{k}^{{n^{ \sim } }} = f \otimes \rho = \frac{{\sum_{h} f_{k} - \rho_{h} }}{{\mathop \sum \nolimits_{h}^{{}} \rho_{h} }}$$
This method is referred to as "PRC simplified" in this study. Both PRC methods (PRC and PRC simplified) together with the generated PR kernels were implemented into the Siemens e7tools (Siemens Medical Solutions USA, Inc., Knoxville, TN, USA) image reconstruction framework using MATLAB R2019a (MathWorks Inc, USA). The image reconstructions with both PRC implementations were performed with both OSEM and PSF corrected OSEM (PSF): OSEM algorithm with PR implemented in the forward projection only (OSEM+PRC simplified), OSEM with PRC implemented in both forward- and back-projection (OSEM+PRC). The same naming convention was used for PSF reconstructions (PSF+PRC simplified and PSF+PRC).
All reconstructions were done with 1–10 iterations and 12 subsets, except the reconstruction for which PSF was combined with PRC in the back-projection (PSF+PRC). In this case, image reconstructions were performed with 1–20 iterations and 12 subsets. The image size was 400 × 400 x 109 voxels, with a voxel size of 2 × 2 × 2 mm3. All emission data was corrected for attenuation, normalization, scatter, and randoms. The time-of-flight (TOF) information was included and no post-reconstruction filters were applied.
Comparison of the PRC implementations
The PRC implementations were tested and compared by means of contrast, noise, and convergence properties using the standard "NEMA Image Quality (IQ)" phantom, a modified NEMQ IQ phantom with changed hot spheres referred to as "Small-tumor" phantom and a "Bone–lung" phantom. Of note, the NEMA IQ phantom simulates hot lesions embedded in soft tissue (water) only (Additional file 1: Fig. S1a). The Small-tumor phantom is a specially-modified phantom using the housing of the NEMA IQ phantom, with fillable spheres with diameters of 3.7, 4.8, 6.5, 7.7, 8.9, and 9.7 mm (Additional file 1: Fig. S1b). The dedicated Bone–lung phantom simulates hot lesions embedded in lung and bone mimicking material, and thus, allows a comparison of the PRCs in a tissue-dependent environment (Additional file 1: Fig. S1c). All phantom acquisitions were performed at the University Clinic Essen, Germany, using the Siemens mCT PET/CT system [20].
NEMA IQ
In the NEMA IQ phantom, all the spheres (diameter of 10, 13, 17, 22, 28, and 37 mm) were filled with a 124I activity concentration of 30 kBq/ml, and the background region was filled with an activity concentration of 6 kBq/ml. The emission acquisition time was 60 min. The phantom was positioned with the centers of the spheres in the center of FOV. For the evaluation, 12 volume-of-interests (VOIs) were used in the background region (each with a diameter of 37 mm) and 6 individual VOIs covering every sphere (with the diameter corresponding to the sphere size) (Additional file 1: Fig. S1d). For placing the background VOIs, the center of the FOV slice was selected, where all the spheres are centered as well. The reconstructed images were evaluated for contrast recovery calculated as:
$${\text{Contrast}}\; {\text{recovery}} = \frac{{\frac{{{\text{MEAN}}_{{{\text{signal}}}} }}{{{\text{MEAN}}_{{{\text{background}}}} }}}}{{{\text{Activity }}\;{\text{ratio}}}}$$
where the mean value in the background was calculated as the mean overall VOIs.
Image noise was defined in percentage following the guidelines of EFOMP [22]:
$${\text{Noise}} = \frac{{{\text{STDEV}}_{{{\text{background}}}} }}{{{\text{MEAN}}_{{{\text{background}}}} }}*100$$
The convergence of the investigated reconstruction algorithms in combination with the various PRC implementations was assessed by plotting the contrast recovery vs. noise and how fast contrast recovery and noise is changing with each iteration [14, 23]. For a direct comparison of the reconstructed images, the reconstruction settings were selected to match a 10% background noise as defined as a clinical acceptable noise level as gained from phantom scans [22].
Small-tumor
The Small-tumor phantom was filled with a 124I activity concentration of 25 kBq/ml in the hot lesions and a background activity of 1.2 kBq/ml leading to a signal-to-noise ratio of 20:1. The acquisition time was 30 min [7]. The reconstructed images were evaluated for contrast recovery (Eq. 5) and image noise (Eq. 6) by defining 10 VOIs in the background region (diameter of 20 mm) and individual VOIs for every sphere (with the diameter corresponding to the sphere size) (Additional file 1: Fig. S1e).
Bone–lung
The Bone–lung phantom was composed of 3 different cylinders (each with a diameter of 50 mm) made of different materials (lung with -800 HU, bone with 500 HU, and 1000 HU) that are housed within the standard NEMA IQ phantom [24]. Within every cold cylinder, two fillable spheres were inserted, one with a diameter of 8.5 mm and one with 19.4 mm. The spheres were filled with a 124I activity concentration of 30 kBq/ml and the background with 6 kBq/ml. Similar to the NEMA IQ evaluations, 6 background VOIs (diameter of 37 mm) and VOIs covering the hot spheres (diameter of 8.5 and 19.4 mm) were defined (Additional file 1: Fig. S1f). The acquisition time was 60 min, and the reconstructed images were evaluated for contrast recovery (Eq. 5) and noise (Eq. 6).
NEMA IQ phantom
Substantial differences in convergence, image contrast, and noise properties were noticed between the different PRC implementations.
In general, three different convergence patterns were observed that can be described as follows: OSEM+PRC simplified had similar convergence as OSEM with higher achievable contrast but also noise levels (Fig. 1). In both cases, 2 iterations and 12 subsets resulted in the background noise of approximately 10%.
Noise versus contrast recovery for the NEMA IQ phantom: (a) 10 mm, (b) 13 mm, (c) 17 mm, (d) 22 mm, (e) 28 mm and (f) 37 mm spheres filled with 124I and reconstructed with OSEM (1–5 iterations), OSEM+PRC simplified (1–3 iterations), OSEM+PRC (1–10 iterations), PSF (1–10 iterations), PSF+PRC simplified (1–10 iterations), and PSF+PRC (1–20 iterations). Every point corresponds to one iteration (all reconstructions done with 12 subsets). The behavior of the OSEM+PRC is similar to the stand-alone PSF reconstructions. The selected number of iterations for the comparison of the images is marked with red
OSEM+PRC and PSF showed similar convergence, contrast, and noise level with improved contrast and noise properties, albeit with a slower enhancement of noise (the 10% noise level was reached with 8 iterations and 12 subsets) compared to OSEM and OSEM+PRC simplified (Fig. 1). PRC in combination with PSF resulted in similar convergence and contract recovery for the full and simplified implementation; however, convergence was notably slower by means of gain in contrast per iteration for PSF+PRC compared to PSF+PRC simplified (Fig. 1).
The visual inspection of the reconstructed NEMA IQ phantoms is shown in Fig. 2. With OSEM+PRC and PSF-PRC simplified, Gibbs artifacts were introduced, which became more pronounced for PSF+PRC (Fig. 3). Regardless of the PRC implementation, the effect was more pronounced for the smaller spheres (10, 13, and 17 mm).
Comparison of the central axial slice of the reconstructed (top) NEMA IQ phantom using: (a) OSEM; (b) OSEM+PRC simplified; (c) OSEM+PRC; (d) PSF; (e) PSF+PRC simplified; and (f) PSF+PRC. The images with similar noise levels measured in the background are marked with dashed boxes. These noise levels were achieved using 2 iterations for OSEM and OSEM+PRC simplified, 8 iterations from OSEM+PRC and PSF, 7 iterations for PSF+PSF simplified, and 15 iterations for PSF+PRC. All reconstructions were with 12 subsets. The OSEM+PRC is producing similar images then the stand-alone PSF; however, more pronounced Gibbs artifacts are produced. (bottom) Small-tumor phantom using: (g) OSEM; (h) OSEM+PRC simplified; (i) OSEM+PRC; (j) PSF; (k) PSF+PRC simplified, and l) PSF+PRC. The images with similar noise levels measured in the background (~ 33%) are marked with dashed boxes. These noise levels were achieved using 2 iterations for OSEM and OSEM+PRC simplified, 9 iterations from OSEM+PRC, 8 iterations for PSF, 7 iterations for PSF+PSF simplified, and 18 iterations for PSF+PRC All reconstructions were with 12 subsets
Comparison of the reconstructed (a) NEMA Image Quality phantom with OSEM (2 iterations and 12 subsets), OSEM+PRC simplified (2 iterations and 12 subsets), OSEM+PRC (8 iterations and 12 subsets), PSF (8 iterations and 12 subsets), PSF+PRC simplified (7 iterations and 12 subsets), and PSF+PRC (15 iterations and 12 subsets) and (b) Small-tumor phantom with OSEM (2 iterations and 12 subsets), OSEM+PRC simplified (2 iterations and 12 subsets), OSEM+PRC (9 iterations and 12 subsets), PSF (8 iterations and 12 subsets), PSF+PRC simplified (7 iterations and 12 subsets) and PSF+PRC (18 iterations and 12 subsets). For the NEMA IQ phantom, the highest differences (%) seen were for the smaller spheres 10, 13, and 17 mm and the edge of the phantom
Table 1 summarizes the lesion contrasts with all PRC implementations for a clinically acceptable image noise of ~ 10%. For subsequent evaluations, those reconstruction settings were selected that yielded a background noise level of ~ 10% following the convergence analyses (Fig. 1).
Table 1 Recovery coefficient for the different image reconstruction and kernel combinations. The image reconstruction settings were defined to match the background noise around ~ 10% for the NEMA IQ phantom and ~ 33% for the Small-tumor phantom. Relative deviations (%) apply relative to the OSEM and for PSF reconstructions
Small-tumor phantom
For this phantom, similar convergence patterns were observed as in the case of the NEMA IQ phantom (Fig. 4). The effect of PRC was noticeable for hot lesions down to 4.8 mm. When comparing images with similar noise levels (~ 33%) (Fig. 2), OSEM+PRC improved contrast recovery by 80% and 86% for the 4.8 and 6.5 mm spheres, respectively (Table 1), which is also clearly seen visually (Fig. 3b). Table 1 summarizes the contrasts improvements with all PRC implementations for a comparable image noise level of ~ 33%.
Noise vs. contrast recovery for the Small-tumor phantom (a) 3.7 mm, (b) 4.8 mm, (c) 6.5 mm, (d) 7.7 mm, (e) 8.9 mm, 9.7 mm spheres filled with 124I and reconstructed with OSEM (1–5 iterations), OSEM+PRC simplified (1–3 iterations), OSEM+PRC (1–10 iterations), PSF (1–10 iterations), PSF+PRC simplified (1–10 iterations), and PSF+PRC (1–20 iterations). Every point corresponds to one iteration (all reconstructions done with 12 subsets). The behavior of the OSEM+PRC is similar to the stand-alone PSF reconstructions. The selected number of iterations for the comparison of the images are marked with red
Bone–lung phantom
The convergence of the reconstructions was highly affected by the surrounding medium. Within the bone inserts, almost no change in contrast improvement was achieved for either the small or big spheres with the OSEM+PRC simplified (Fig. 5). The convergence was similar for OSEM+PRC, PSF, and PSF+PRC reconstructions. Within the lung medium, the effect of PRC was clearly visible with substantially improved contrasts (Table 2). The OSEM+PRC simplified showed similar convergence as OSEM with a substantially higher contrast recovery. The trend was similar with OSEM+PRC and PSF reconstruction with reduced image noise. PSF+PRC simplified and PSF+PRC had almost identical convergence, with the simplified implementation leading to convergence substantially faster. The effect of PRC within the lung medium is also clearly visible on the difference images shown in Fig. 6. Table 2 summarizes the contrast recoveries for all reconstruction and PRC combinations at comparable image noise (~ 10%).
Contrast recovery vs. noise calculated for the small sphere (8.5 mm) in (a) lung, (b) bone (500 HU), (c) bone (1000 HU) and for the large sphere (19.4 mm) in (d) lung, (e) bone (500 HU), (f) bone (1000 HU) reconstructed with PRC methods in combination with OSEM and PSF. The changes are most pronounced within the lung insert
Table 2 Contrast recovery for every reconstruction method and kernel variation using the Bone–lung phantom. Relative deviations (%) in relation to the standard OSEM and for the PSF reconstructions
Comparison of the reconstructed Bone–lung phantom with (a) OSEM and (b) PSF the central slices are shown through the small (8.5 mm) and (b) the central slices are shown through large (19.4 mm) spheres. The recovery is slightly overestimated for the 8.5 mm sphere within the lung medium. The black dashed lines represent the activity concentrations in the hot spheres and background
In this paper, different PRC implementation methods in combination with PSF image reconstruction were evaluated for PET imaging using 124I. The implementation of the PR-dependent blurring into the PET image reconstruction was done for the OSEM algorithm and for a PSF corrected OSEM as full and simplified implementation [14]. Both PRC implementations, simplified and full implementation, lead to increased contrast recovery compared to standard OSEM. However, the full implementations also led to a slower convergence of the algorithm although with higher achievable contrasts compared to with PRC simplified. The simplified implementation demonstrated a faster convergence although with substantially higher noise levels (Fig. 1). The reason for the faster convergence and the increased noise levels are expected due to result from the strong mismatch of the forward and backward projectors [3].
The image reconstruction with OSEM+PRC showed a similar convergence to the PSF reconstruction (Fig. 1). This can be explained by the similarity of the full implementation of the PRC and the vendor-based PSF image reconstructions [3].
When adding the PRC to the PSF reconstruction, the image contrast recovery was almost identical for the full and the simplified PRC implementation. This observation may be explained by the higher similarity of the forward and backward projector in the case of the combination of the PRC with PSF. This observation is also strengthened by the results for the OSEM-PRC simplified, which indicates that the mismatch of forward and backward projector substantially influences the noise propagation within iterative PRC implementations (Figs. 2 and 3).
The effect of PRC was the most pronounced within low-density materials, such as lung, where the reconstructions with PRC outperformed the standard reconstructions (Fig. 4). In this case, also OSEM+PRC was leading to higher contrast recovery compared to PSF (Fig. 5). On the other hand, within the bone medium, relatively small changes were seen for PRC-based reconstructions. This was expected due to the high positron range of 124I within the lung (max. 30 mm) and short PR in bone. Gibbs artifacts well known from PSF reconstructions were also present in the PRC reconstructions and more prominent when PSF and PRC were combined (Figs. 3 and 6). As shown in Fig. 6 in the line profile through the small spheres, a slight overestimation can be caused by the Gibbs artifacts. The effect is more visible in the edges of the bigger spheres. However, such Gibbs artifacts as known from PSF are mainly occurring at sharp activity concentration borders as present in phantom studies and are expected to be less pronounced in patients. Further, as shown for PSF reconstructions, such artifacts can be handled using post-reconstruction filters or by applying regularization methods [25, 26].
The composition of the spatially variant and tissue-dependent kernels does not take into account the positron energy loss, in particular when a positron is emitted from a higher density material. The solution for these limitations would require a more complex PR such as based on direct MC simulations; however, due to an extensive number of kernels that have to be computed during the image reconstruction process, the full PRC application is still a challenge especially with the aim of clinical application.
The data acquisition of the NEMA IQ phantom was 60 min. For obtaining similar count statistic as for a standard 15 min NEMA IQ scan with 18F given the reduced branching ratio for positrons in 124I. In the clinical scenario, the reconstruction settings may have to be adjusted to account for the reduced count statistics due to the shorter acquisition times.
The evaluated PRC methods are based on the same basic methodology of calculating the spatially variant and tissue-dependent kernels, by analyzing the underlying material composition from the AC maps [14]. Thus, the errors in the attenuation map directly translate into artifacts caused by the PRC. This holds also true in the case of PET/MRI, were artifacts in MR-based AC might occur more frequent that for PET/CT. Furthermore, when applying the evaluated methods for the PET/MRI, the effect of the magnetic field on the PR has to be carefully considered [15, 27].
Combining PSF and PRC notably increases achievable contrast compared to PSF or PRC alone in 124I PET and appears to be preferred for formal correctness, contrast recovery, and noise level. Our study showed that the use of a simplified PRC implementation applying the PR-based blurring exclusively to the forward projector within an OSEM algorithm, resulted in a seriously increase in image noise compared to the full implementation of PRC in both, forward- and back-projection. However, a simplified implementation of PRC in the forward projector appears acceptable, when combined with a PSF correction.
The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
Cal-Gonzalez J, Rausch I, Shiyam Sundar LK, Lassen ML, Muzik O, Moser E, et al. Hybrid imaging: instrumentation and data processing. Front Phys. 2018;6:47.
Reader AJ, Ally S, Bakatselos F, Manavaki R, Walledge RJ, Jeavons AP, et al. One-pass list-mode EM algorithm for high-resolution 3-D PET image reconstruction into large arrays. IEEE Trans Nucl Sci. 2002;49:693–9.
Panin VY, Kehren F, Michel C, Casey M. Fully 3-D PET reconstruction with system matrix derived from point source measurements. IEEE Trans Med Imag. 2006;25:907–21.
Jakoby BW, Bercier Y, Watson CC, Bendriem B, Townsend DW. Performance characteristics of a New LSO PET/CT scanner with extended axial field-of-view and PSF reconstruction. IEEE Trans Nucl Sci. 2009;56:633–9.
Tong S, Alessio AM, Kinahan PE. Noise and signal properties in PSF-based fully 3D PET image reconstruction: an experimental evaluation. Phys Med Biol. 2010;55:1453–73.
Ashrafinia S, Mohy-Ud-Din H, Karakatsanis NA, Jha AK, Casey ME, Kadrmas DJ, et al. Generalized PSF modeling for optimized quantitation in PET imaging. Phys Med Biol. 2017;62:5149–79.
Kersting D, Jentzen W, Fragoso Costa P, Sraieb M, Sandach P, Umutlu L, et al. Silicon-photomultiplier-based PET/CT reduces the minimum detectable activity of iodine-124. Sci Rep. 2021;11:17477.
Rahmim A, Qi J, Sossi V. Resolution modeling in PET imaging: theory, practice, benefits, and pitfalls. Med Phys. 2013;40: 064301.
Cal-González J, Herraiz JL, España S, Corzo PMG, Vaquero JJ, Desco M, et al. Positron range estimations with PeneloPET. Phys Med Biol. 2013;58:5127–52.
Carter LM, Kesner AL, Pratt EC, Sanders VA, Massicano AVF, Cutler CS, et al. The impact of positron range on PET resolution, evaluated with phantoms and PHITS Monte Carlo simulations for conventional and non-conventional radionuclides. Mol Imag Biol. 2020;22:73–84.
Hoffmann MA, Wieler HJ, Baues C, Kuntz NJ, Richardsen I, Schreckenberger M. The impact of 68Ga-PSMA PET/CT and PET/MRI on the management of prostate cancer. Urology. 2019;130:1–12.
Freudenberg LS, Jentzen W, Stahl A, Bockisch A, Rosenbaum-Krumme SJ. Clinical applications of 124I-PET/CT in patients with differentiated thyroid cancer. Eur J Nucl Med Mol Imag. 2011;38(Suppl 1):S48-56.
Conti M, Eriksson L. Physics of pure and non-pure positron emitters for PET: a review and a discussion. EJNMMI Phys. 2016;3:8.
Kertész H, Beyer T, Panin V, Jentzen W, Cal-Gonzalez J, Berger A, et al. Implementation of a spatially-variant and tissue-dependent position range correction for PET/CT imaging. Front Physiol. 2022. https://doi.org/10.3389/fphys.2022.818463.
Bertolli O, Eleftheriou A, Cecchetti M, Camarlinghi N, Belcari N, Tsoumpas C. PET iterative reconstruction incorporating an efficient positron range correction method. Phys Med. 2016;32:323–30.
Fraile LM, Herraiz JL, Udías JM, Cal-González J, Corzo PMG, España S, et al. Experimental validation of gallium production and isotope-dependent positron range correction in PET. Nucl Instrum Methods Phys Res A. 2016;814:110–6.
Anton-Rodriguez JM, Krokos G, Kotasidis F, Asselin M-C, Morris O, Julyan P, et al. Experimental validation of estimated spatially variant radioisotope-specific point spread functions using published positron range simulations and fluorine-18 measurements. Phys Med Biol. 2018;63:24NT01.
Cal-González J, Pérez-Liva M, Herraiz JL, Vaquero JJ, Desco M, Udías JM. Tissue-Dependent and spatially-variant positron range correction in 3D PET. IEEE Trans Med Imag. 2015;34:2394–403.
Cal-González J, Herraiz JL, España S, Vicente E, Herranz E, Desco M, et al. Study of CT-based positron range correction in high resolution 3D PET imaging. Nucl Instrum Methods Phys Res A. 2011;648:S172–5.
Rausch I, Cal-González J, Dapra D, Gallowitsch HJ, Lind P, Beyer T, et al. Performance evaluation of the Biograph mCT Flow PET/CT system according to the NEMA NU2-2012 standard. EJNMMI Phys. 2015;2:26.
Cal-Gonzalez J, Vaquero JJ, Herraiz JL, Pérez-Liva M, Soto-Montenegro ML, Peña-Zalbidea S, et al. Improving PET quantification of small animal [68Ga]DOTA-labeled PET/CT studies by using a CT-based positron range correction. Mol Imag Biol. 2018;20:584–93.
Quality controls in PET/CT and PET/MR. EFOMP; 2022.
Lois C, Jakoby BW, Long MJ, Hubner KF, Barker DW, Casey ME, et al. An assessment of the impact of incorporating time-of-flight information into clinical PET/CT imaging. J Nucl Med. 2010;51:237–45.
Fragoso Costa P, Jentzen W, Süßelbeck F, Fendler WP, Rischpler C, Herrmann K, et al. Reduction of emission time for [68Ga]Ga-PSMA PET/CT using the digital biograph vision: a Phantom study. Q J Nucl Med Mol Imag. 2021. https://doi.org/10.23736/S1824-4785.21.03300-8.
Herraiz JL, España S, Vicente E, Vaquero JJ, Desco M, Udias JM. Noise and physical limits to maximum resolution of PET images. Nucl Instrum Methods Phys Res A. 2007;580:934–7.
Knäusl B, Rausch IF, Bergmann H, Dudczak R, Hirtl A, Georg D. Influence of PET reconstruction para meters on the TrueX algorithm. Nuklearmedizin Schattauer GmbH. 2013;52:28–35.
Eleftheriou A, Tsoumpas C, Bertolli O, Stiliaris Ε. Effect of the magnetic field on positron range using GATE for PET-MR. mMR Whole-Body PET/MR System. IEEE Trans Nucl Sci. 2012;59:1900–9.
The financial support of the Austrian Science Funds (FWF) Project I3451-N32 is gratefully acknowledged. The authors acknowledge the support of Siemens Medical Solutions USA, Inc. (Knoxville, TN, USA).
QIMP Team, Center for Medical Physics and Biomedical Engineering, Medical University of Vienna, Währinger Gürtel 18-20, 1090, Vienna, Austria
Hunor Kertész, Thomas Beyer, Laszlo Papp, Jacobo Cal-Gonzalez & Ivo Rausch
Siemens Medical Solutions USA, Inc., Knoxville, TN, USA
Maurizio Conti, Vladimir Panin, Jorge Cabello & Deepak Bharkhada
Clinic for Nuclear Medicine, University Hospital Essen, Essen, Germany
Walter Jentzen
Ion Beam Applications, Protontherapy Center Quironsalud, Madrid, Spain
Jacobo Cal-Gonzalez
Nuclear Physics Group and IPARCOS, Faculty of Physical Sciences, University Complutense of Madrid, Madrid, Spain
Joaquín L. Herraiz & Alejandro López-Montes
Health Research Institute of the Hospital Clínico San Carlos, Madrid, Spain
Joaquín L. Herraiz
Hunor Kertész
Maurizio Conti
Vladimir Panin
Jorge Cabello
Deepak Bharkhada
Thomas Beyer
Laszlo Papp
Alejandro López-Montes
Ivo Rausch
All authors contributed to writing, critically reviewing and approving the manuscript. Specific author contributions are as follows: HK implemented the PRC methods, did the image reconstructions, analyzed the data, and wrote the manuscript draft. MC, VP, and JC provided expertise on the PRC implementations, DB provided the image reconstruction tools, LP optimized the PRC implementation, WJ performed the phantom measurements, JH and ALM provided with help and expertise with the kernel calculations, TB, JCG, and IR designed, guided, and supervised the project. All authors read and approved the final manuscript.
Correspondence to Hunor Kertész.
No patient information/data was used in the manuscript.
Consent to participate
Dr. Maurizio Conti, Dr. Vladimir Panin, Dr. Jorge Cabello, and Dr. Deepak Bharkhada are employees of Siemens Medical Solutions USA, Inc., (Knoxville, TN, USA) and report no competing interests with this study.
Additional file 1
. Employed phantoms: a) NEMA IQ - Image quality phantom; b) Small tumor phantom and c) Bone-lung phantom. d-f) VOI delineations for the image analysis.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Kertész, H., Conti, M., Panin, V. et al. Positron range in combination with point-spread-function correction: an evaluation of different implementations for [124I]-PET imaging. EJNMMI Phys 9, 56 (2022). https://doi.org/10.1186/s40658-022-00482-y
Positron range correction