IEC Phantom Research Articles

Anatomy-based correction of kidney PVE on 177Lu\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{177}\ ext{Lu}$$\\end{document} SPECT images

BackgroundIn peptide receptor radionuclide therapy (PRRT), accurate quantification of kidney activity on post-treatment SPECT images paves the way for patient-specific treatment. Due to the limited spatial resolution of SPECT images, the partial volume effect (PVE) is a significant source of quantitative bias. In this study, we aimed to evaluate the performance and robustness of anatomy-based partial volume correction (PVC) algorithms to recover the accurate activity concentration of realistic kidney geometries on 177\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{177}$$\\end{document}Lu SPECT images recorded under clinical conditions.MethodsBased on the CT scan data from patients, three sets of fillable kidneys with surface-to-volume (S:V) ratios ranging from 1.5 to 2.8 cm−1, were 3D printed and attached in a IEC phantom. Quantitative 177\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{177}$$\\end{document}Lu SPECT/CT acquisitions were performed on a GE Discovery NM CT 870 DR camera for the three modified IEC phantoms and for 6 different Target-To-Background ratios (TBRs: 2, 4, 6, 8, 10, 12). Two region-based (GTM and Labbé) and five voxel-based (GTM + MTC, Labbé + MTC, GTM + RBV, Labbé + RBV and IY) methods were evaluated with this data set. Additionally, the robustness of PVC methods to Point Spread Function (PSF) discrepancies, registration mismatches and background heterogeneity was evaluated.ResultsWithout PVC, the average kidney RCs across all TBRs ranged from 0.66 ± 0.05 (smallest kidney) to 0.80 ± 0.03 (largest kidney). For a TBR of 12, all anatomy-based method were able to recover the kidneys activity concentration with an error < 6%. All methods result in a comparable decline in RC restoration with decreasing TBR. The Labbé method was the most robust against PSF and registration mismatches but was also the most sensitive to background heterogeneity. Among the voxel-based methods, MTC images were less uniform than RBV and IY images at the outer edge of high uptake areas (kidneys and spheres).ConclusionAnatomy-based PVE correction allows for accurate SPECT quantification of the 177\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{177}$$\\end{document}Lu activity concentration with realistic kidney geometries. Combined with recent progress in deep-learning algorithms for automatic anatomic segmentation of whole-body CT, these methods could be of particular interest for a fully automated OAR dosimetry pipeline with PVE correction.

The Role of Lung Density in the Voxel-Based Dosimetry of 90Y-TARE Evaluated with the Voxel S-Value (VSV) Method and Fast Monte Carlo Simulation

(1) Background: In 90Y-TARE treatments, lung-absorbed doses should be calculated according to the manufacturer’s instructions, using the MIRD-scheme. This scheme is derived from the assumption that 90Y-microspheres deliver the dose in a water-equivalent medium. Since the density of the lungs is quite different from that of the liver, the absorbed dose to the lungs could vary considerably, especially at the liver/lungs interface. The aim of this work is to compare the dosimetric results obtained by two dedicated software packages implementing a water-equivalent dose calculation and a Monte Carlo (MC) simulation, respectively. (2) Methods: An anthropomorphic IEC phantom and a retrospective selection of 24 patients with a diagnosis of HCC were taken into account. In the phantom study, starting from a 90Y-PET/CT acquisition, the liver cavity was manually fixed with a uniform activity concentration on PET series, while the lung compartment was manually expanded on a CT series to simulate a realistic situation in which the liver and lungs are adjacent. These steps were performed by using MIM 90Y SurePlan. Then, a first simulation was carried out with only the liver cavity filled, while a second one was carried out, in which the lung compartment was also manually fixed with a uniform activity concentration corresponding to 10% lung shunt fraction. MIM 90Y SurePlan was used to obtain Voxel S-Value (VSV) approach dose values; instead, Torch was used to obtain MC approach dose values for both the phantom and the patients. (3) Results: In the phantom study, the percentage mean dose differences (∆D%) between VSV and MC in the first and second simulation, respectively were found to be 1.2 and 0.5% (absolute dose variation, ∆D, of 0.7 and 0.3 Gy) for the liver, −56 and 70% (∆D of −0.3 and −16.2 Gy) for the lungs, and −48 and −60% (∆D of −4.3 and −16.5 Gy) for the Liver/Lungs Edge region. The patient study reports similar results with ∆D% between VSV and MC of 7.0%, 4.1% and 6.7% for the whole liver, healthy liver, and tumor, respectively, while the result was −61.2% for the left lung and −61.1% for both the right lung and lungs. (4) Conclusion: Both VSV and MC allowed accurate radiation dose estimation with small differences (&lt;7%) in regions of uniform water-equivalent density (i.e., within the liver). Larger differences between the two methods (&gt;50%) were observed for air-equivalent regions in the phantom simulation and the patient study.

The determination of the optimal threshold on measurement of thyroid volume using quantitative SPECT/CT for Graves' hyperthyroidism

PurposeTo investigate the optimal threshold for measuring thyroid volume in patients with Grave's hyperthyroidism (GH) by SPECT/CT.Materials and methodsA 53 mL butterfly-shaped hollow container made of two 45-degree transparent elbows was put into a NEMA IEC phantom tank. The butterfly-shaped container and the tank were then filled with Na99mTcO4 of different radioactive concentrations, respectively, which could simulate thyroid gland with GH by different target-to-background ratios (T/B) (200:1, 600:1, 1000:1). The different T/B of planar imaging and SPECT/CT were acquired by a Discovery NM/CT 670 Pro SPECT/CT. With Thyroid software (Version 4.0) of GE-Xeleris workstation, the region of the thyroid gland in planar imaging was delineated. The thyroid area and average long diameter of both lobes were substituted into the Allen formula to calculate the thyroid volume. The calculation error was compared with the actual volume. Q-Metrix software was used to perform CT-based attenuation correction, scatter correction, resolution recovery. Ordered-subsets expectation maximization was used to reconstruct SPECT data. 20%, 25%, 30%, 40%, 50%, 60% thresholds were selected to automatically delineate the volume of interest and compared with the real volume, which determinated the optimal threshold. We measured the thyroid volume of 40 GH patients using the threshold and compared the volumes obtained by planar imaging and ultrasound three-dimensional. The differences of the volumes with different T/B and thresholds were compared by the ANOVA and least significant difference t test. The volumes delineated by SPECT/CT were evaluated using ANOVA, least significant difference t test, correlation analysis and, linear regression and Bland–Altman concordance test plot. The differences and consistency of thyroid volume were compared among the above three methods.ResultsThere was no significant difference in the results between different T/B models (P > 0.05). The thyroid volume calculated by the planar imaging formula method was higher than the real volume, with an average overestimation of 22.81%. The volumes delineated by SPECT/CT threshold automatically decreased while the threshold increased. There were significant differences between groups with different thresholds (P < 0.001). With an average error of 3.73%, the thyroid volume analyzed by the threshold of 25% was close to the results of ultrasound measurement (P > 0.05). Thyroid volume measured by planar imaging method was significantly higher than ultrasound and SPECT/CT threshold automatic delineation method (P < 0.05). The agreement between the SPECT/CT 25% threshold and ultrasound (r = 0.956, b = 0.961) was better than that between the planar imaging and ultrasound (r = 0.590, b = 0.574). The Bland–Altman plot also showed that the thyroid volume measured by the 25% threshold automatic delineation method was in good agreement with the ultrasound measurement.ConclusionsThe T/B has no effect on the measurement of thyroid volume in GH patients; planar imaging method can significantly overestimate thyroid volume in GH patients, and 25% threshold automatic delineation method can obtain more accurate thyroid volume in GH patients.

Quantitative 177Lu SPECT/CT imaging for personalized dosimetry using a ring-shaped CZT-based camera

BackgroundDosimetry after radiopharmaceutical therapy with 177Lu (177Lu-RPT) relies on quantitative SPECT/CT imaging, for which suitable reconstruction protocols are required. In this study, we characterized for the first time the quantitative performance of a ring-shaped CZT-based camera using two different reconstruction algorithms: an ordered subset expectation maximization (OSEM) and a block sequential regularized expectation maximization (BSREM) combined with noise reduction regularization. This study lays the foundations for the definition of a reconstruction protocol enabling accurate dosimetry for patients treated with 177Lu-RPT.MethodsA series of 177Lu-filled phantoms were acquired on a StarGuide™ (GE HealthCare), with energy and scatter windows centred at 208 (± 6%) keV and 185 (± 5%) keV, respectively. Images were reconstructed with the manufacturer implementations of OSEM (GE-OSEM) and BSREM (Q.Clear) algorithms, and various combinations of iterations and subsets. Additionally, the manufacturer-recommended Q.Clear-based reconstruction protocol was evaluated. Quantification accuracy, measured as the difference between the SPECT-based and the radionuclide calibrator-based activity, and noise were evaluated in a large cylinder. Recovery coefficients (RCs) and spatial resolution were assessed in a NEMA IEC phantom with sphere inserts. The reconstruction protocols considered suitable for clinical applications were tested on a cohort of patients treated with [177Lu]Lu-PSMA-I&T.ResultsThe accuracy of the activity from the cylinder, although affected by septal penetration, was < 10% for all reconstructions. Both algorithms featured improved spatial resolution and higher RCs with increasing updates at the cost of noise build-up, but Q.Clear outperformed GE-OSEM in reducing noise accumulation. When the reconstruction parameters were carefully selected, similar values for noise (~0.15), spatial resolution (~1 cm) and RCs were found, irrespective of the reconstruction algorithm. Analogue results were found in patients.ConclusionsAccurate activity quantification is possible when imaging 177Lu with StarGuide™. However, the impact of septal penetration requires further investigations. GE-OSEM is a valid alternative to the recommended Q.Clear reconstruction algorithm, featuring comparable performances assessed on phantoms and patients.

Impact of the maximum ring difference on image quality and noise characteristics of a total-body PET/CT scanner

The sensitivity of a PET system highly depends on the axial acceptance angle or maximum ring difference (MRD), which can be particularly high for total-body scanners due to their larger axial field of views (aFOVs). This study aims to evaluate the impact on image quality (IQ) and noise performance when MRD85 (18°), the current standard for clinical use, is increased to MRD322 (52°) for the Biograph Vision Quadra (Siemens Healthineers). MethodsStudies with a cylindrical phantom covering the 106 cm aFOV and an IEC phantom filled with 18F, 68Ga and 89Zr were performed for acquisition times from 60 to 1800 s and activity concentrations from 0.4 to 3 kBq/ml to assess uniformity, contrast recovery coefficients (CRCs) and to characterize noise by coefficient of variation (CV). Spatial resolution was compared for both MRDs by sampling a quadrant of the FOV with a point source. Further IQ, CV, liver SUVmean and SUVmax were compared for a cohort of 5 patients scanned with [18F]FDG (3 MBq/kg, 1 h p.i.) from 30 to 300 s. ResultsCV was improved by a factor of up to 1.49 and is highest for short acquisition times, peaks at the center field of view and mitigates parabolic in axial direction with no difference to MRD85 beyond the central 80 cm. No substantial differences between the two evaluated MRDs in regards to uniformity, SUVmean or CRC for the different isotopes were observed. A degradation of the average spatial resolution of 0.9 ± 0.2 mm in the central 40 cm FOV was determined with MRD322. Depending on the acquisition time MRD322 resulted in a decrease of SUVmax between 23.8% (30 s) and 9.0% (300 s). ConclusionPatient and phantom studies revealed that scan time could be lowered by approximately a factor of two with MRD322 while maintaining similar noise performance. The moderate degradation in spatial resolution for MRD322 is worth to exploit the full potential of the Quadra by either shorten scan times or leverage noise performance in particular for low count scenarios such as ultra-late imaging or dynamic studies with high temporal resolution.

Multicenter quantitative 18F-fluorodeoxyglucose positron emission tomography performance harmonization: use of hottest voxels towards more robust quantification.

Harmonization methods reduce variability between different make and models of positron emission tomography (PET) scanners. The study aims to explore harmonization strategies that lead to comparable and robust quantitative metrics in a multicenter setting. NEMA IEC Phantom data acquisition was performed for low and high spheres-to-background ratios (SBR4:1 and 10:1) on six PET/CT (computed tomography) scanners. Different reconstruction sets, including the number of sub-iterations, number of subsets, and full width at half maximum (FWHM) for each scanner, were evaluated towards optimized and harmonized reconstruction settings. Recovery coefficients (RCs) of four quantitative metrics, including standardized uptake value (SUV)max, SUVISO-50 (SUVmean in 50% isocontour), SUVpeak, and mean uptake of 10 highest concentration voxels were evaluated as RCmax, RCISO-50, RCpeak, and RC10V, representing percent difference relative to the static ground truth case as functions of sphere sizes. A set of image reconstruction parameters was proposed for harmonized reconstruction to minimize variability between scanners. The root mean square error (RMSE), curvature, and reproducibility were examined. The proposed reconstruction protocols for harmonization and standard clinical reconstruction settings were compared to each other across all scanners. A significant difference (P value <0.0001) was observed in the aforementioned quantitative metrics between SBR10 and SBR4. Reconstruction parameter sets with the smallest RMSE and RC values within 10% bias were identified as the best candidate for harmonization. The coefficient of variation of the mean value of RCs (CVMRC) shows a remarkable reduction of about 28%, 26%, 32%, and 19% in harmonized reconstruction settings for MRCmax, MRCISO-50, MRCpeak, and MRC10V, respectively. CVMRC for MRC10V in the harmonized reconstruction setting was 5.9% in SBR4, while the smallest value in SBR10 belongs to MRCpeak, with a value of 5.8%. The reproducibility of RC is improved by deriving the value from ten hottest voxels and is equally reproducible with RCpeak. Compared to RCmax and RCISO-50, the variability is reduced by 18% and 22% if ten voxels are pooled. Harmonizing PET/CT systems with and without point spread function/time of flight (PSF/TOF) using various vendor-developed image reconstruction algorithms improves the quantification reproducibility. RC10V, likewise RCpeak, is superior to the rest of the quantitative indices in terms of accuracy and reproducibility and helpful in quantifying lesion volume below 1 mL.

Transformation of PET raw data into images for event classification using convolutional neural networks.

In positron emission tomography (PET) studies, convolutional neural networks (CNNs) may be applied directly to the reconstructed distribution of radioactive tracers injected into the patient's body, as a pattern recognition tool. Nonetheless, unprocessed PET coincidence data exist in tabular format. This paper develops the transformation of tabular data into n-dimensional matrices, as a preparation stage for classification based on CNNs. This method explicitly introduces a nonlinear transformation at the feature engineering stage and then uses principal component analysis to create the images. We apply the proposed methodology to the classification of simulated PET coincidence events originating from NEMA IEC and anthropomorphic XCAT phantom. Comparative studies of neural network architectures, including multilayer perceptron and convolutional networks, were conducted. The developed method increased the initial number of features from 6 to 209 and gave the best precision results (79.8) for all tested neural network architectures; it also showed the smallest decrease when changing the test data to another phantom.

Clinical and phantom validation of a deep learning based denoising algorithm for F-18-FDG PET images from lower detection counting in comparison with the standard acquisition

BackgroundPET/CT image quality is directly influenced by the F-18-FDG injected activity. The higher the injected activity, the less noise in the reconstructed images but the more radioactive staff exposition. A new FDA cleared software has been introduced to obtain clinical PET images, acquired at 25% of the count statistics considering US practices. Our aim is to determine the limits of a deep learning based denoising algorithm (SubtlePET) applied to statistically reduced PET raw data from 3 different last generation PET scanners in comparison to the regular acquisition in phantom and patients, considering the European guidelines for radiotracer injection activities. Images of low and high contrasted (SBR = 2 and 5) spheres of the IEC phantom and high contrast (SBR = 5) of micro-spheres of Jaszczak phantom were acquired on 3 different PET devices. 110 patients with different pathologies were included. The data was acquired in list-mode and retrospectively reconstructed with the regular acquisition count statistic (PET100), 50% reduction in counts (PET50) and 66% reduction in counts (PET33). These count reduced images were post-processed with SubtlePET to obtain PET50 + SP and PET33 + SP images. Patient image quality was scored by 2 senior nuclear physicians. Peak-signal-to-Noise and Structural similarity metrics were computed to compare the low count images to regular acquisition (PET100).ResultsSubtlePET reliably denoised the images and maintained the SUVmax values in PET50 + SP. SubtlePET enhanced images (PET33 + SP) had slightly increased noise compared to PET100 and could lead to a potential loss of information in terms of lesion detectability. Regarding the patient datasets, the PET100 and PET50 + SP were qualitatively comparable. The SubtlePET algorithm was able to correctly recover the SUVmax values of the lesions and maintain a noise level equivalent to full-time images.ConclusionBased on our results, SubtlePET is adapted in clinical practice for half-time or half-dose acquisitions based on European recommended injected dose of 3 MBq/kg without diagnostic confidence loss.

Evaluating the impact of different positron emitters on the performance of a clinical PET/MR system.

The positron range and prompt gamma emission are distinctive with different positron emitters. The performance assessment of an integrated PET/MR scanner with these positron emitters is required for related applications, as the magnetic field interferes with the positron propagation. Such an assessment is to be performed on the United Imaging uPMR 790-integrated PET/MRsystem. The performance measurement methods were modified based on NEMA NU 2-2012, involving 18 F, 64 Cu, 68 Ga, 89 Zr, and 124 I as positron emitters. The NEMA IEC phantom was used for evaluations of image qualities. An agarose cap was wrapped around the point source for tissue-simulating spatial resolution measurement. The count rate performance was assessed with selected positron emitters. Images of a 3D-printed Derenzo phantom and representative patients were alsoacquired. The image quality measurement showed that all five positron emitters were suitable for the PET/MR system studied. However, due to the magnetic field, the image of the point source showed an elongated comet-tail feature, which could be eliminated by a tissue-simulating cap. This effect is more obvious in 124 I and 68 Ga, due to their long positron ranges. The imaging ability with various positron emitters was further validated with the count rate assessment, the Derenzo phantom, and the clinicalimages. Different positron emitters could be effectively imaged by the PET/MR system tested. The resolution measurement strategy proposed could be applied to measure PET spatial resolution in the magneticfield.

Influence of reconstruction parameters of positron emission tomograph scanning on the effect of partial volume of the pathological lesion

In this work, the following tasks were solved: to perform a comparative analysis of data processing methods when calculating recovery factors; to evaluate the influence of time-of-flight technology and PSF function on the recovery factor and the forecast of recovery factor deviation for potential pathological foci with a diameter of 6–8 mm; to evaluate the influence of parameters of iterative reconstruction algorithms, Gaussian filter and axial filters on the recovery factor. The calculation of the recovery factors was carried out on the basis of quantitative characteristics obtained in the analysis of reconstructions of images of the IEC phantom with six spheres installed inside and filled with a radiopharmaceutical. Eight series of experiments with background / sphere activity ratios 1/3, 1/4, 1/6, 1/8, 1/12, 1/14, 1/16, 1/20 were carried out with the same concentration of activity in the spheres during each separate experiment. The forecast of the effect of the partial volume effect on lesions with a diameter of 6 to 8 mm was carried out, taking into account the used reconstruction algorithms. It is advisable to use the results obtained to harmonize diagnostic protocols for scanning with positron emission tomographs using the input parameters of reconstruction algorithms and filters, which will minimize the error in the quantitative assessment of a radiopharmaceutical when analyzing the dynamics of the development of a pathological process, as well as the response of pathology to therapy.

Evaluation of attenuation correction in PET/MRI with synthetic lesion insertion.

Purpose: One major challenge facing simultaneous positron emission tomography (PET)/ magnetic resonance imaging (MRI) is PET attenuation correction (AC) measurement and evaluation of its accuracy. There is a crucial need for the evaluation of current and emergent PET AC methodologies in terms of absolute quantitative accuracy in the reconstructed PET images. Approach: To address this need, we developed and evaluated a lesion insertion tool for PET/MRI that will facilitate this evaluation process. This tool was developed for the Biograph mMR and evaluated using phantom and patient data. Contrast recovery coefficients (CRC) from the NEMA IEC phantom of synthesized lesions were compared to measurements. In addition, SUV biases of lesions inserted in human brain and pelvis images were assessed from PET images reconstructed with MRI-based AC (MRAC) and CT-based AC (CTAC). Results: For cross-comparison PET/MRI scanners AC evaluation, we demonstrated that the developed lesion insertion tool can be harmonized with the GE-SIGNA lesion insertion tool. About CRC curves difference between simulation and measurement was achieved. An average of 1.6% between harmonized simulated CRC curves obtained with mMR and SIGNA lesion insertion tools was achieved. A range of to 12% MRAC to CTAC SUV bias was respectively achieved in the vicinity and inside bone tissues in patient images in two anatomical regions, the brain, and pelvis. Conclusions: A lesion insertion tool was developed for the Biograph mMR PET/MRI scanner and harmonized with the SIGNA PET/MRI lesion insertion tool. These tools will allow for an accurate evaluation of different PET/MRI AC approaches and permit exploration of subtle attenuation correction differences across systems.

Bone Scintigraphy of Vertebral Fractures With a Whole-Body CZT Camera in a PET-Like Utilization.

Objective: An image display with a standardized uptake value (SUV) scale is recommended for analyzing PET exams, thus requiring the reconstruction of accurate images for both SUV measurement and visual analysis. This study aimed to determine whether such images may also be obtained with a high-speed CZT-SPECT/CT system, with a further application for the longitudinal monitoring of vertebral fractures. Materials and Methods: SPECT image reconstruction was optimized with an IEC phantom according to both image quality parameters and accuracy of measured activity. The optimized reconstruction process was applied to ≤15 min 99mTc-HDP SPECT spine recordings previously acquired from 25 patients (74 ± 12 years old) at both early (1.3 ± 1.1 months) and late (5.2 ± 2.3 months) stages after an acute vertebral fracture. Results: A SPECT reconstruction with 32 equivalent iterations was selected based on the association of high detectability for spheres down to 0.6 ml in volume, with accurate measured activity, although the latter was affected by partial volume effect for spheres ≤5.6 ml. Coherent measurements were obtained on these high-quality SPECT images for the SUVmax from the intact vertebrae of patients, which were stable between basal SPECT/CT and follow-up SPECT/CT (for T1 vertebrae: 5.7 ± 1.1 vs. 5.8 ± 1.1, p = 0.76), and from initially fractured vertebrae, which were dramatically higher on the basal compared with the follow-up SPECT (21.0 ± 8.5 vs. 11.2 ± 4.2, p < 0.001), whereas inverse changes in SUVmax were observed for newly compacted fractures identified on follow-up SPECT (74.4 ± 2.0 vs. 21.8 ± 10.3, p = 0.002). Finally, an image display with an SUV scale was shown to be advantageous for highlighting areas with >7.5 SUV, a level reached by 98% of vertebral fractures of <7 months and 4% of reference intact vertebrae. Conclusion: Bone scintigraphy of vertebral fractures may be obtained with this CZT-SPECT/CT system with fast 3D acquisitions and high-quality images displayed with a reliable SUV scale, approaching what is achieved and recommended for PET imaging.

Simulating NEMA characteristics of the modular total-body J-PET scanner—an economic total-body PET from plastic scintillators

The purpose of the presented research is estimation of the performance characteristics of the economic total-body Jagiellonian-PET system (TB-J-PET) constructed from plastic scintillators. The characteristics are estimated according to the NEMA NU-2-2018 standards utilizing the GATE package. The simulated detector consists of 24 modules, each built out of 32 plastic scintillator strips (each with cross section of 6 mm times 30 mm and length of 140 or 200 cm) arranged in two layers in regular 24-sided polygon circumscribing a circle with the diameter of 78.6 cm. For the TB-J-PET with an axial field-of-view (AFOV) of 200 cm, a spatial resolutions (SRs) of 3.7 mm (transversal) and 4.9 mm (axial) are achieved. The noise equivalent count rate (NECR) peak of 630 kcps is expected at 30 kBq cc−1. Activity concentration and the sensitivity at the center amounts to 38 cps kBq−1. The scatter fraction (SF) is estimated to 36.2 %. The values of SF and SR are comparable to those obtained for the state-of-the-art clinical PET scanners and the first total-body tomographs: uExplorer and PennPET. With respect to the standard PET systems with AFOV in the range from 16 to 26 cm, the TB-J-PET is characterized by an increase in NECR approximately by factor of 4 and by the increase of the whole-body sensitivity by factor of 12.6 to 38. The time-of-flight resolution for the TB-J-PET is expected to be at the level of CRT = 240 ps full width at half maximum. For the TB-J-PET with an AFOV of 140 cm, an image quality of the reconstructed images of a NEMA IEC phantom was presented with a contrast recovery coefficient and a background variability parameters. The increase of the whole-body sensitivity and NECR estimated for the TB-J-PET with respect to current commercial PET systems makes the TB-J-PET a promising cost-effective solution for the broad clinical applications of total-body PET scanners. TB-J-PET may constitute an economic alternative for the crystal TB-PET scanners, since plastic scintillators are much cheaper than BGO or LYSO crystals and axial arrangement of the strips significantly reduces the costs of readout electronics and SiPMs.

Optimisation of the event-based TOF filtered back-projection for online imaging in total-body J-PET

We perform a parametric study of the newly developed time-of-flight (TOF) image reconstruction algorithm, proposed for the real-time imaging in total-body Jagiellonian PET (J-PET) scanners. The asymmetric 3D filtering kernel is applied at each most likely position of electron-positron annihilation, estimated from the emissions of back-to-back γ-photons. The optimisation of its parameters is studied using Monte Carlo simulations of a 1-mm spherical source, NEMA IEC and XCAT phantoms inside the ideal J-PET scanner. The combination of high-pass filters which included the TOF filtered back-projection (FBP), resulted in spatial resolution, 1.5 times higher in the axial direction than for the conventional 3D FBP. For realistic 10-minute scans of NEMA IEC and XCAT, which require a trade-off between the noise and spatial resolution, the need for Gaussian TOF kernel components, coupled with median post-filtering, is demonstrated. The best sets of 3D filter parameters were obtained by the Nelder-Mead minimisation of the mean squared error between the resulting and reference images. The approach allows training the reconstruction algorithm for custom scans, using the IEC phantom, when the temporal resolution is below 50 ps. The image quality parameters, estimated for the best outcomes, were systematically better than for the non-TOF FBP.

Optimizing reconstruction parameters for quantitative 124I-PET in the presence of therapeutic doses of 131I

BackgroundThe goal of this work was to determine the quantitative accuracy and optimal reconstruction parameters for 124I-PET imaging in the presence of therapeutic levels of 131I. In this effort, images were acquired on a GE D710 PET/CT scanner using a NEMA IEC phantom with spheres containing 124I and increasing amounts of 131I activity in the background. At each activity level, two scans were acquired, one with the phantom centered in the field of view (FOV) and one 11.2 cm off-center. Reconstructions used an ordered subset expectation maximization algorithm with up to 100 iterations of 16 subsets, with and without time-of-flight (TOF) information. Results were evaluated visually and by comparing the 124I activity relative to the scan performed in the absence of 131I.Results131I within the FOV added to the randoms rate, to dead time, and to pile-up within the detectors. Using our standard clinical reconstruction parameters, the image quality and quantitative accuracy suffered at 131I activities above 1.4 GBq. Convergence rates slowed progressively in the presence of increasing amounts of 131I for both TOF and nonTOF reconstructions. TOF reconstructions converged more quickly than nonTOF but often towards erroneous concentrations. Iterating nonTOF reconstructions to convergence produced quantitatively accurate images except for the off-center phantom at the very highest level of background 131I tested.ConclusionsThis study shows that quantitative PET is feasible in the presence of large amounts of 131I. The high randoms fractions resulted in slow reconstruction convergence and negatively impacted TOF corrections and/or the accuracy of TOF information. Therefore, increased iterations and nonTOF reconstructions are recommended.

Theranostic SPECT reconstruction for improved resolution: application to radionuclide therapy dosimetry

BackgroundSPECT-derived dose estimates in tissues of diameter less than 3× system resolution are subject to significant losses due to the limited spatial resolution of the gamma camera. Incorporating resolution modelling (RM) into the SPECT reconstruction has been proposed as a possible solution; however, the images produced are prone to noise amplification and Gibbs artefacts. We propose a novel approach to SPECT reconstruction in a theranostic setting, which we term SPECTRE (single photon emission computed theranostic reconstruction); using a diagnostic PET image, with its superior resolution, to guide the SPECT reconstruction of the therapeutic equivalent. This report demonstrates a proof in principle of this approach.MethodsWe have employed the hybrid kernelised expectation maximisation (HKEM) algorithm implemented in STIR, with the aim of producing SPECT images with PET-equivalent resolution. We demonstrate its application in both a dual 68Ga/177Lu IEC phantom study and a clinical example using 64Cu/67Cu.ResultsSPECTRE is shown to produce images comparable in accuracy and recovery to PET with minimal introduction of artefacts and amplification of noise.ConclusionThe SPECTRE approach to image reconstruction shows improved quantitative accuracy with a reduction in noise amplification. SPECTRE shows great promise as a method of improving SPECT radioactivity concentrations, directly leading to more accurate dosimetry estimates in small structures and target lesions. Further investigation and optimisation of the algorithm parameters is needed before this reconstruction method can be utilised in a clinical setting.

3D TOF-PET image reconstruction using total variation regularization.

In this paper we introduce a semi-analytic algorithm for 3-dimensional image reconstruction for positron emission tomography (PET). The method consists of the back-projection of the acquired data into the most likely image voxel according to time-of-flight (TOF) information, followed by the filtering step in the image space using an iterative optimization algorithm with a total variation (TV) regularization. TV regularization in image space is more computationally efficient than usual iterative optimization methods for PET reconstruction with full system matrix that use TV regularization. The efficiency comes from the one-time TOF back-projection step that might also be described as a reformatting of the acquired data. An important aspect of our work concerns the evaluation of the filter operator of the linear transform mapping an original radioactive tracer distribution into the TOF back-projected image. We obtain concise, closed-form analytical formula for the filter operator. The proposed method is validated with the Monte Carlo simulations of the NEMA IEC phantom using a one-layer, 50cm-long cylindrical device called Jagiellonian PET scanner. The results show a better image quality compared with the reference TOF maximum likelihood expectation maximization algorithm.

Asphericity of tumor FDG uptake in non-small cell lung cancer: reproducibility and implications for harmonization in multicenter studies

BackgroundAsphericity (ASP) of the primary tumor’s metabolic tumor volume (MTV) in FDG-PET/CT is independently predictive for survival in patients with non-small cell lung cancer (NSCLC). However, comparability between PET systems may be limited. Therefore, reproducibility of ASP was evaluated at varying image reconstruction and acquisition times to assess feasibility of ASP assessment in multicenter studies.MethodsThis is a retrospective study of 50 patients with NSCLC (female 20; median age 69 years) undergoing pretherapeutic FDG-PET/CT (median 3.7 MBq/kg; 180 s/bed position). Reconstruction used OSEM with TOF4/16 (iterations 4; subsets 16; in-plane filter 2.0, 6.4 or 9.5 mm), TOF4/8 (4 it; 8 ss; filter 2.0/6.0/9.5 mm), PSF + TOF2/17 (2 it; 17 ss; filter 2.0/7.0/10.0 mm) or Bayesian-penalized likelihood (Q.Clear; beta, 600/1750/4000). Resulting reconstructed spatial resolution (FWHM) was determined from hot sphere inserts of a NEMA IEC phantom. Data with approx. 5-mm FWHM were retrospectively smoothed to achieve 7-mm FWHM. List mode data were rebinned for acquisition times of 120/90/60 s. Threshold-based delineation of primary tumor MTV was followed by evaluation of relative ASP/SUVmax/MTV differences between datasets and resulting proportions of discordantly classified cases.ResultsReconstructed resolution for narrow/medium/wide in-plane filter (or low/medium/high beta) was approx. 5/7/9 mm FWHM. Comparing different pairs of reconstructed resolution between TOF4/8, PSF + TOF2/17, Q.Clear and the reference algorithm TOF4/16, ASP differences was lowest at FWHM of 7 versus 7 mm. Proportions of discordant cases (ASP > 19.5% vs. ≤ 19.5%) were also lowest at 7 mm (TOF4/8, 2%; PSF + TOF2/17, 4%; Q.Clear, 10%). Smoothing of 5-mm data to 7-mm FWHM significantly reduced discordant cases (TOF4/8, 38% reduced to 2%; PSF + TOF2/17, 12% to 4%; Q.Clear, 10% to 6%), resulting in proportions comparable to original 7-mm data. Shorter acquisition time only increased proportions of discordant cases at < 90 s.ConclusionsASP differences were mainly determined by reconstructed spatial resolution, and multicenter studies should aim at comparable FWHM (e.g., 7 mm; determined by in-plane filter width). This reduces discordant cases (high vs. low ASP) to an acceptable proportion for TOF and PSF + TOF of < 5% (Q.Clear: 10%). Data with better resolution (i.e., lower FWHM) could be retrospectively smoothed to the desired FWHM, resulting in a comparable number of discordant cases.

Implementation of patient dosimetry in the clinical practice after targeted radiotherapy using [177Lu-[DOTA0, Tyr3]-octreotate

BackgroundThis study’s aim was to develop our dosimetric methodology using a commercial workstation for the routine evaluation of the organs at risk during peptide receptor radionuclide therapy (PRRT) with 177Lu.MethodsFirst, planar and SPECT sensitivity factors were determined on phantoms. The reconstruction parameters were optimized by SPECT/CT image acquisition using a NEMA IEC phantom containing a 500 ml bottle of 177Lu, to simulate a kidney. The recovery coefficients were determined on various phantoms. For the red marrow, this was calculated using a NEMA IEC phantom that contained a centrally placed bottle of 80 ml of 177Lu (to model the L2-L4 red marrow) flanked by two 200 ml bottles with 177Lu to simulate the kidneys.Then, SPECT/CT images were acquired at 4, 24, 72, and 192 h after injection in 12 patients with neuroendocrine tumors who underwent PRRT with 177Lu-DOTATATE. SPECT data were reconstructed using the iterative ordered subset expectation maximization (OSEM) method, with six iterations and ten subsets, attenuation, scatter, recovery resolution corrections, and a Gaussian post-filter of 0.11 cm. The liver, spleen, kidneys, and red marrow dose per administered activity (AD/A admin) values were calculated with the Medical Internal Radiation Dose (MIRD) formalism and the residence times (Dosimetry toolkit® application) using standard and CT imaging-based organ masses (OLINDA/EXM® V1.0 software).ResultsSensitivity factors of 6.11 ± 0.01 and 5.67 ± 0.08 counts/s/MBq were obtained with planar and SPECT/CT acquisitions, respectively. A recovery coefficient of 0.78 was obtained for the modeled L2–L4 red marrow. The mean AD/A admin values were 0.43 ± 0.13 mGy/MBq [0.27–0.91] for kidneys, 0.54 ± 0.58 mGy/MBq [0.12–2.26] for liver, 0.61 ± 0.13 mGy/MBq [0.42–0.89] for spleen, and 0.04 ± 0.02 mGy/MBq [0.01–0.09] for red marrow. The AD/A admin values varied when calculated using the personalized and standard organ mass, particularly for kidneys (p = 1 × 10−7), spleen (p = 0.0069), and red marrow (p = 0.0027). Intra-patient differences were observed especially in organs close to or including tumor cells or metastases.ConclusionsThe obtained AD/A admin values were in agreement with the literature data. This study shows the technical feasibility of patient dosimetry in clinical practice and the need to obtain patient-specific information.

Dose volume histogram-based optimization of image reconstruction parameters for quantitative 90 Y-PET imaging.

90 Y-microsphere radioembolization or selective internal radiation therapy is increasingly being used as a treatment option for tumors that are not candidates for surgery and external beam radiation therapy. Recently, volumetric 90 Y-dosimetry techniques have been implemented to explore tumor dose-response on the basis of 3D 90 Y-activity distribution from PET imaging. Despite being a theranostic study, the optimization of quantitative 90 Y-PET image reconstruction still uses the mean activity concentration recovery coefficient (RC) as the objective function, which is more relevant to diagnostic and detection tasks than is to dosimetry. The aim of this study was to optimize 90 Y-PET image reconstruction by minimizing errors in volumetric dosimetry via the dose volume histogram (DVH). We propose a joint optimization of the number of equivalent iterations (the product of the iterations and subsets) and the postreconstruction filtration (FWHM) to improve the accuracy of voxel-level 90 Y dosimetry. A modified NEMA IEC phantom was used to emulate clinically relevant 90 Y-PET imaging conditions through various combinations of acquisition durations, activity concentrations, sphere-to-background ratios, and sphere diameters. PET data were acquired in list mode for 300 min in a single-bed position; we then rebinned the list mode PET data to 60, 45, 30, 15, and 5 min per bed, with 10 different realizations. Errors in the DVH were calculated as root mean square errors (RMSE) of the differences in the image-based DVH and the expected DVH. The new optimization approach was tested in a phantom study, and the results were compared with the more commonly used objective function of the mean activity concentration RC. In a wide range of clinically relevant imaging conditions, using 36 equivalent iterations with a 5.2-mm filtration resulted in decreased systematic errors in volumetric 90 Y dosimetry, quantified as image-based DVH, in 90 Y-PET images reconstructed using the ordered subset expectation maximization (OSEM) iterative reconstruction algorithm with time of flight (TOF) and point spread function (PSF) modeling. Our proposed objective function of minimizing errors in DVH, which allows for joint optimization of 90 Y-PET iterations and filtration for volumetric quantification of the 90 Y dose, was shown to be superior to conventional RC-based optimization approaches for image-based absorbed dose quantification. Our proposed objective function of minimizing errors in DVH, which allows for joint optimization of iterations and filtration to reduce errors in the PET-based volumetric quantification 90 Y dose, is relevant to dosimetry in therapy procedures. The proposed optimization method using DVH as the objective function could be applied to any imaging modality used to assess voxel-level quantitative information.