Quantitative SPECT reconstruction using CT-derived corrections

Single-centre experience with peptide receptor radionuclide therapy for neuroendocrine tumours (NETs): results using a theranostic molecular imaging-guided approach

AIM

To summarise our centre's experience managing patients with neuroendocrine tumours (NETs) in the first 5 years after the introduction of peptide receptor radionuclide therapy (PRRT) with [177Lu]Lu-DOTA-octreotate (LUTATE). The report emphasises aspects of the patient management related to functional imaging and use of radionuclide therapy.

METHODS

We describe the criteria for treatment with LUTATE at our centre, the methodology for patient selection, and the results of an audit of clinical measures, imaging results and patient-reported outcomes. Subjects are treated initially with four cycles of ~ 8 GBq of LUTATE administered as an outpatient every 8 weeks.

RESULTS

In the first 5 years offering LUTATE, we treated 143 individuals with a variety of NETs of which approx. 70% were gastroentero-pancreatic in origin (small bowel: 42%, pancreas: 28%). Males and females were equally represented. Mean age at first treatment with LUTATE was 61 ± 13 years with range 28-87 years. The radiation dose to the organs considered most at risk, the kidneys, averaged 10.6 ± 4.0 Gy in total. Median overall survival (OS) from first receiving LUTATE was 72.5 months with a median progression-free survival (PFS) of 32.3 months. No evidence of renal toxicity was seen. The major long-term complication seen was myelodysplastic syndrome (MDS) with a 5% incidence.

CONCLUSIONS

LUTATE treatment for NETs is a safe and effective treatment. Our approach relies heavily on functional and morphological imaging informing the multidisciplinary team of NET specialists to guide appropriate therapy, which we suggest has contributed to the favourable outcomes seen.

Standardized Radiomics Analysis of Clinical Myocardial Perfusion Stress SPECT Images to Identify Coronary Artery Calcification

PURPOSE

Myocardial perfusion (MP) stress single-photon emission computed tomography (SPECT) is an established diagnostic test for patients suspected of coronary artery disease (CAD). Meanwhile, coronary artery calcification (CAC) scoring obtained from diagnostic CT is a highly sensitive test, offering incremental diagnostic information in identifying patients with significant CAD yet normal MP stress SPECT (MPSS) scans. However, after decades of wide utilization of MPSS, CAC is not commonly reimbursed (e.g. by the CMS), nor widely deployed in community settings. We studied the potential of complementary information deduced from the radiomics analysis of normal MPSS scans in predicting the CAC score.

METHODS

We collected data from 428 patients with normal (non-ischemic) MPSS (99mTc-sestamibi; consensus reading). A nuclear medicine physician verified iteratively reconstructed images (attenuation-corrected) to be free from fixed perfusion defects and artifactual attenuation. Three-dimensional images were automatically segmented into four regions of interest (ROIs), including myocardium and three vascular segments (left anterior descending [LAD]-left circumference [LCX]-right coronary artery [RCA]). We used our software package, standardized environment for radiomics analysis (SERA), to extract 487 radiomic features in compliance with the image biomarker standardization initiative (IBSI). Isotropic cubic voxels were discretized using fixed bin-number discretization (eight schemes). We first performed blind-to-outcome feature selection focusing on a priori usefulness, dynamic range, and redundancy of features. Subsequently, we performed univariate and multivariate machine learning analyses to predict CAC scores from i) selected radiomic features, ii) 10 clinical features, and iii) combined radiomics + clinical features. Univariate analysis invoked Spearman correlation with Benjamini-Hotchberg false-discovery correction. The multivariate analysis incorporated stepwise linear regression, where we randomly selected a 15% test set and divided the other 85% of data into 70% training and 30% validation sets. Training started from a constant (intercept) model, iteratively adding/removing features (stepwise regression), invoking the Akaike information criterion (AIC) to discourage overfitting. Validation was run similarly, except that the training output model was used as the initial model. We randomized training/validation sets 20 times, selecting the best model using log-likelihood for evaluation in the test set. Assessment in the test set was performed thoroughly by running the entire operation 50 times, subsequently employing Fisher's method to verify the significance of independent tests.

RESULTS

Unsupervised feature selection significantly reduced 8×487 features to 56. In univariate analysis, no feature survived the false-discovery rate (FDR) to directly correlate with CAC scores. Applying Fisher's method to the multivariate regression results demonstrated combining radiomics with the clinical features to enhance the significance of the prediction model across all cardiac segments.  Conclusions: Our standardized and statistically robust multivariate analysis demonstrated significant prediction of the CAC score for all cardiac segments when combining MPSS radiomic features with clinical features, suggesting radiomics analysis can add diagnostic or prognostic value to standard MPSS for wide clinical usage.

Validation of [Formula: see text]Tc and [Formula: see text]Lu quantification parameters for a Monte Carlo modelled gamma camera

PURPOSE

Monte Carlo (MC) simulation in Nuclear Medicine is a powerful tool for modeling many physical phenomena which are difficult to track or measure directly. MC simulation in SPECT/CT imaging is particularly suitable for optimizing the quantification of activity in a patient, and, consequently, the absorbed dose to each organ. To do so, validating MC results with real data acquired with gamma camera is mandatory. The aim of this study was the validation of the calibration factor (CF) and the recovery coefficient (RC) obtained with SIMIND Monte Carlo code for modeling a Siemens Symbia Intevo Excel SPECT-CT gamma camera to ensure optimal [Formula: see text]Tc and [Formula: see text]Lu SPECT quantification.

METHODS

Phantom experiments using [Formula: see text]Tc and [Formula: see text]Lu have been performed to measure spatial resolution and sensitivity, as well as to evaluate the CF and RC from acquired data. The geometries used for 2D planar imaging were (1) Petri dish and (2) capillary source while for 3D volumetric imaging were (3) a uniform filled cylinder phantom and (4) a Jaszczack phantom with spheres of different volumes. The experimental results have been compared with the results obtained from Monte Carlo simulations performed in the same geometries.

RESULTS

Comparison shows good accordance between simulated and experimental data. The measured planar spatial resolution was 8.3[Formula: see text] mm for [Formula: see text]Tc and 11.8±0.6 mm for [Formula: see text]Lu. The corresponding data obtained by SIMIND for [Formula: see text]Tc was 7.8±0.1 mm, while for [Formula: see text]Lu was 12.4±0.4 mm. The CF was 110.1±5.5 cps/MBq for Technetium and 18.3±1.0 cps/MBq for Lutetium. The corresponding CF obtained by SIMIND for [Formula: see text]Tc was 107.3±0.3 cps/MBq, while for [Formula: see text]Lu 20.4±0.7 cps/MBq. Moreover, a complete curve RCs vs Volume (ml) both for Technetium and Lutetium was determined to correct the PVE for all volumes of clinical interest. In none of the cases, a RC coefficient equal to 100 was found.

CONCLUSIONS

The validation of quantification parameters shows that SIMIND can be used for simulating both gamma camera planar and SPECT images of Siemens Symbia Intevo using [Formula: see text]Tc and [Formula: see text]Lu radionuclides for different medical purposes and treatments.

Quantitative assessment of ventilation-perfusion relationships with gallium-68 positron emission tomography/computed tomography imaging in lung cancer patients

•

Voxel-wise correlations between gallium-68 (⁶⁸Ga) positron emission tomography/computed tomography (PET/CT)-measured ventilation and perfusion varied widely among 19 patients with lung cancer (range: 0.26–0.88).

•

⁶⁸Ga PET/CT-measured percent gas exchanging lung volume was moderately correlated with diffusing capacity of the lung for carbon monoxide (D_LCO) as the reference standard, with the highest correlation coefficient of 0.59 (P < 0.01).

•

⁶⁸Ga PET/CT ventilation/perfusion imaging may provide a reasonable surrogate for regional gas exchange.

Pulmonary functional imaging has demonstrated potential to improve thoracic radiotherapy. The purpose of this study was twofold: 1) to quantify ventilation/perfusion relationships in lung cancer patients using a new functional imaging approach, gallium-68 (⁶⁸Ga)-positron emission tomography/computed tomography (PET/CT); and 2) to compare ventilation/perfusion matching with diffusing capacity of the lung for carbon monoxide (D_LCO). Voxel-wise correlations between ventilation and perfusion varied widely among 19 patients (range: 0.26–0.88). ⁶⁸Ga-PET/CT-measured percent gas exchanging lung volume was moderately correlated with D_LCO (≤0.59). Our findings suggested that ⁶⁸Ga-PET/CT ventilation/perfusion imaging provided complementary information and a reasonable surrogate for gas exchange in lung cancer patients.

Quantitative SPECT/CT for dosimetry of peptide receptor radionuclide therapy

Neuroendocrine tumors (NETs) are uncommon malignancies of increasing incidence and prevalence. As these slow growing tumors usually overexpress somatostatin receptors (SSTRs), the use of ⁶⁸Ga-DOTA-peptides (gallium-68 chelated with dodecane tetra-acetic acid to somatostatin), which bind to the SSTRs, allows for PET based imaging and selection of patients for peptide receptor radionuclide therapy (PRRT). PRRT with radiolabeled somatostatin analogues such as ¹⁷⁷Lu-DOTATATE (lutetium-177-[DOTA,Tyr³]-octreotate), is mainly used for the treatment of metastatic or inoperable NETs. However, PRRT is generally administered at a fixed injected activity in order not to exceed dose limits in critical organs, which is suboptimal given the variability in radiopharmaceutical uptake among patients. Advances in SPECT (single photon emission computed tomography) imaging enable the absolute quantitative measure of the true radiopharmaceutical distribution providing for PRRT dosimetry in each patient. Personalized PRRT based on patient-specific dosimetry could improve therapeutic efficacy by optimizing effective tumor absorbed dose while limiting treatment related radiotoxicity.

Prediction of Absorbed Dose to Normal Organs with Endocrine Tumors for I-131 by use of 99mTC Single Photon Emission Computed Tomography/Computed Tomography and Geant4 Application for Tomographic Emission Simulation

Introduction:

This study aimed to predict the dose absorbed by normal organs with neuroendocrine tumors for ¹³¹I using single photon emission computed tomography/computed tomography (SPECT/CT) images and Geant4 application for tomographic emission (GATE) simulation.

Materials and Methods:

Four to 5 whole-body planar scan series, along with one SPECT/CT image, were taken from four patients following ^99mTc-hynic-Tyr³-octreotide radiotracer injection. After image quantification, the residence time of each organ was calculated using the image analysis and the activity time curves. The energy deposit and dose conversion (S-value) were extracted from the GATE simulation for the target organs of each patient. Using the residence times and S-values, the mean absorbed dose for the target organs of each patient was calculated and compared with the data obtained from the standard method.

Results:

Very close agreement was obtained between the S-value of the self–organ irradiation. The mean percentage difference between the two methods (i.e. GATE and Medical Internal Radiation Dose [MIRD]) was 1.8%, while a weak agreement was observed for cross-organ irradiation. The percentage difference between the total absorbed doses by the organs was 2%. The percentage difference between the absorbed doses obtained for tumors and three considered normal organs estimated by the GATE method was slightly higher than the MIRD method (about 11% on average for tumors).

Conclusion:

Regardless of the small difference between the obtained results for the organs and absorbed doses of the tumors in the present study, patient-specific dosimetry by the GATE methods is useful and essential for therapeutic radionuclides such as ¹³¹I due to high cross-dose effects, especially for young adult patients, to ensure the radiation safety and increase the effectiveness of the treatment.

Quantitative evaluation of single-photon emission computed tomography findings in lower extremity possible without computed tomography-based attenuation correction

OBJECTIVE

For performing accurate quantitative analysis of single-photon emission computed tomography (SPECT)/computed tomography (CT) images, CT-based attenuation correction (CTAC) is considered to be necessary. However, the effect on quantitative values for an examined area close to the body surface, such as in the lower extremity, has yet to be elucidated. We performed the present investigation to determine the possibility of quantitative evaluation using a SPECT standalone device without CT.

METHODS

Validation was performed using clinical data of patients who underwent a lower extremity SPECT/CT examination, with grouping based on presence or absence of CTAC, scatter correction and resolution recovery. Using a reference group in which all types of correction were applied, standardized uptake values (SUVs), including maximum (SUVmax) and peak (SUVpeak), were examined in each group and compared.

RESULTS

As compared to the reference group, the difference in quantitative values became smaller in the order of the applied scatter correction and resolution recovery, applied resolution recovery, applied scatter correction, and neither scatter correction or resolution recovery applied groups, with no significant difference between the reference group and that with neither scatter correction or resolution recovery applied. A similar tendency was seen for both SUVmax and SUVpeak.

CONCLUSIONS

In bone SPECT quantitative examinations of the lower extremity, quantitative evaluation without CTAC is possible without the use of scatter correction or resolution recovery. Thus, quantitative evaluation can be performed with use of a standalone SPECT device without CT.

Evaluation of Iodine-123 and Iodine-131 SPECT activity quantification: a Monte Carlo study

Purpose

The quantitative accuracy of Nuclear Medicine images, acquired for both planar and SPECT studies, is influenced by the isotope-collimator combination as well as image corrections incorporated in the iterative reconstruction process. These factors can be investigated and optimised using Monte Carlo simulations. This study aimed to evaluate SPECT quantification accuracy for ¹²³I with both the low-energy high resolution (LEHR) and medium-energy (ME) collimators and ¹³¹I with the high-energy (HE) collimator.

Methods

Simulated SPECT projection images were reconstructed using the OS-EM iterative algorithm, which was optimised for the number of updates, with appropriate corrections for scatter, attenuation and collimator detector response (CDR), including septal scatter and penetration compensation. An appropriate calibration factor (CF) was determined from four different source geometries (activity-filled: water-filled cylindrical phantom, sphere in water-filled (cold) cylindrical phantom, sphere in air and point-like source), investigated with different volume of interest (VOI) diameters. Recovery curves were constructed from recovery coefficients to correct for partial volume effects (PVEs). The quantitative method was evaluated for spheres in voxel-based digital cylindrical and patient phantoms.

Results

The optimal number of OS-EM updates was 60 for all isotope-collimator combinations. The CF_point with a VOI diameter equal to the physical size plus a 3.0-cm margin was selected, for all isotope-collimator geometries. The spheres’ quantification errors in the voxel-based digital cylindrical and patient phantoms were less than 3.2% and 5.4%, respectively, for all isotope-collimator combinations.

Conclusion

The study showed that quantification errors of less than 6.0% could be attained, for all isotope-collimator combinations, if corrections for; scatter, attenuation, CDR (including septal scatter and penetration) and PVEs are performed. ¹²³I LEHR and ¹²³I ME quantification accuracies compared well when appropriate corrections for septal scatter and penetration were applied. This can be useful in departments that perform ¹²³I studies and may not have access to ME collimators.

Optimization of SPECT/CT imaging protocols for quantitative and qualitative 99mTc SPECT

Background

The introduction of hybrid SPECT/CT devices enables quantitative imaging in SPECT, providing a methodological setup for quantitation using SPECT tracers comparable to PET/CT. We evaluated a specific quantitative reconstruction algorithm for SPECT data using a ^99mTc-filled NEMA phantom. Quantitative and qualitative image parameters were evaluated for different parametrizations of the acquisition and reconstruction protocol to identify an optimized quantitative protocol.

Results

The reconstructed activity concentration (AC_rec) and the signal-to-noise ratio (SNR) of all examined protocols (n = 16) were significantly affected by the parametrization of the weighting factor k used in scatter correction, the total number of iterations and the sphere volume (all, p < 0.0001). The two examined SPECT acquisition protocols (with 60 or 120 projections) had a minor impact on the AC_rec and no significant impact on the SNR.

In comparison to the known AC, the use of default scatter correction (k = 0.47) or object-specific scatter correction (k = 0.18) resulted in an underestimation of AC_rec in the largest sphere volume (26.5 ml) by − 13.9 kBq/ml (− 16.3%) and − 7.1 kBq/ml (− 8.4%), respectively. An increase in total iterations leads to an increase in estimated AC and a decrease in SNR. The mean difference between AC_rec and known AC decreased with an increasing number of total iterations (e.g., for 20 iterations (2 iterations/10 subsets) = − 14.6 kBq/ml (− 17.1%), 240 iterations (24i/10s) = − 8.0 kBq/ml (− 9.4%), p < 0.0001). In parallel, the mean SNR decreased significantly from 2i/10s to 24i/10s by 76% (p < 0.0001).

Conclusion

Quantitative SPECT imaging is feasible with the used reconstruction algorithm and hybrid SPECT/CT, and its consistent implementation in diagnostics may provide perspectives for quantification in routine clinical practice (e.g., assessment of bone metabolism). When combining quantitative analysis and diagnostic imaging, we recommend using two different reconstruction protocols with task-specific optimized setups (quantitative vs. qualitative reconstruction). Furthermore, individual scatter correction significantly improves both quantitative and qualitative results.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-021-00405-3.

[Modern methods in diagnostics of lacrimal system vertical part obliteration]

A review of the literature is devoted to the description of modern diagnostic methods that are directed at detecting lacrimal pathway obstruction. The presented analysis is based on the study of data from Russian and foreign literary sources on both routine and high-tech lacrimal system examination methods, as well as on authors' experience in studying the diagnostic efficiency of imaging methods such as lacrimal scintigraphy, computed tomography with lacrimal pathway contrasting, and also single photon emission computed tomography. The article gives information about the advantages and disadvantages of the presented methods, highlights the questions of the prospects of their usage in practical dacryology.

The effect of calibration factors and recovery coefficients on 177Lu SPECT activity quantification accuracy: a Monte Carlo study

BACKGROUND

Different gamma camera calibration factor (CF) geometries have been proposed to convert SPECT data into units of activity concentration. However, no consensus has been reached on a standardised geometry. The CF is dependent on the selected geometry and is further affected by partial volume effects. This study investigated the effect of two CF geometries and their corresponding recovery coefficients (RCs) on the quantification accuracy of 177Lu SPECT images using Monte Carlo simulations.

METHODS

The CF geometries investigated were (i) a radioactive-sphere surrounded by non-radioactive water (sphere-CF) and (ii) a cylindrical phantom uniformly filled with radioactive water (cylinder-CF). Recovery coefficients were obtained using the sphere-CF and cylinder-CF, yielding the sphere-RC and cylinder-RC values, respectively, for partial volume correction (PVC). The quantification accuracy was evaluated using four different-sized spheres (15.6-65.4 ml) and a kidney model with known activity concentrations inside a cylindrical, torso and patient phantom. Images were reconstructed with the 3D OS-EM algorithm incorporating attenuation, scatter and detector-response corrections. Segmentation was performed using the physical size and a small cylindrical volume inside the cylinder for the sphere-CF and cylinder-CF, respectively.

RESULTS

The sphere quantification error (without PVC) was better for the sphere-CF (≤ - 5.54%) compared to the cylinder-CF (≤ - 20.90%), attributed to the similar geometry of the quantified and CF spheres. Partial volume correction yielded comparable results for the sphere-CF-RC (≤ 3.47%) and cylinder-CF-RC (≤ 3.53%). The accuracy of the kidney quantification was poorer (≤ 22.34%) for the sphere-CF without PVC compared to the cylinder-CF (≤ 2.44%). With PVC, the kidney quantification results improved and compared well for the sphere-CF-RC (≤ 3.50%) and the cylinder-CF-RC (≤ 3.45%).

CONCLUSION

The study demonstrated that upon careful selection of CF-RC combinations, comparable quantification errors (≤ 3.53%) were obtained between the sphere-CF-RC and cylinder-CF-RC, when all corrections were applied.

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

Background

SPECT-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.

Methods

We 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 ⁶⁸Ga/¹⁷⁷Lu IEC phantom study and a clinical example using ⁶⁴Cu/⁶⁷Cu.

Results

SPECTRE is shown to produce images comparable in accuracy and recovery to PET with minimal introduction of artefacts and amplification of noise.

Conclusion

The 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.

Phantom and clinical evaluation of bone SPECT/CT image reconstruction with xSPECT algorithm

BACKGROUND

Two novel methods of image reconstruction, xSPECT Quant (xQ) and xSPECT Bone (xB), that use an ordered subset conjugate gradient minimizer (OSCGM) for SPECT/CT reconstruction have been proposed. The present study compares the performance characteristics of xQ, xB, and conventional Flash3D (F3D) reconstruction using images derived from phantoms and patients.

METHODS

A custom-designed body phantom for bone SPECT was scanned using a Symbia Intevo (Siemens Healthineers), and reconstructed xSPECT images were evaluated. The phantom experiments proceeded twice with different activity concentrations and sphere sizes. A phantom with 28-mm spheres containing a 99mTc-background and tumor-to-normal bone ratios (TBR) of 1, 2, 4, and 10 were generated, and convergence property against various TBR was evaluated across 96 iterations. A phantom with four spheres (13-, 17-, 22-, and 28-mm diameters), containing a 99mTc-background at TBR4, was also generated. The full width at half maximum of an imaged spinous process (10 mm), coefficients of variance (CV), contrast-to-noise ratio (CNR), and recovery coefficients (RC) were evaluated after reconstructing images of a spine using Flash 3D (F3D), xQ, and xB. We retrospectively analyzed images from 20 patients with suspected bone metastases (male, n = 13) which were acquired using [99mTc]Tc-(H)MDP SPECT/CT, then CV and standardized uptake values (SUV) at the 4th vertebral body (L4) were compared after xQ and xB reconstruction in a clinical setup.

RESULTS

Mean activity concentrations with various TBR converged according to increasing numbers of iterations. The spatial resolution of xB was considerably superior to xQ and F3D, and it approached almost the actual size regardless of the iteration numbers during reconstruction. The CV and RC were better for xQ and xB than for F3D. The CNR peaked at 24 iterations for xQ and 48 iterations for F3D and xB, respectively. The RC between xQ and xB significantly differed at lower numbers of iterations but were almost equivalent at higher numbers of iterations. The reconstructed xQ and xB images of the clinical patients showed a significant difference in the SUVmax and SUVpeak.

CONCLUSIONS

The reconstructed xQ and xB images were more accurate than those reconstructed conventionally using F3D. The xB for bone SPECT imaging offered essentially unchanged spatial resolution even when the numbers of iterations did not converge. The xB reconstruction further enhanced SPECT image quality using CT data. Our findings provide important information for understanding the performance characteristics of the novel xQ and xB algorithms.

Validation of a Monte Carlo modelled gamma camera for Lutetium-177 imaging

This study validated a model of the Siemens Symbia T16 dual-head SPECT/CT gamma camera created using the Monte Carlo program SIMIND for ¹⁷⁷Lu. The validation was done by comparing experimental and simulated gamma camera performance criteria tests for the ¹⁷⁷Lu 208 keV photopeak with a medium-energy collimator. Results showed good agreement between the experimental and simulated values. These results illustrated that SIMIND could emulate the Symbia T16 successfully and therefore, can be used with confidence to model ¹⁷⁷Lu images.

Head-to-head comparison of lung perfusion with dual-energy CT and SPECT-CT

PURPOSE

To compare the quantitative and qualitative lung perfusion data acquired with dual energy CT (DECT) to that acquired with a large field-of-view cadmium-zinc-telluride camera single-photon emission CT coupled to a CT system (SPECT-CT).

MATERIALS AND METHODS

A total of 53 patients who underwent both dual-layer DECT angiography and perfusion SPECT-CT for pulmonary hypertension or pre-operative lobar resection surgery were retrospectively included. There were 30 men and 23 women with a mean age of 65.4±17.5 (SD)years (range: 18-88years). Relative lobar perfusion was calculated by dividing the amount (of radiotracer or iodinated contrast agent) per lobe by the total amount in both lungs. Linear regression, Bland-Altman analysis, and Pearson's correlation coefficient were also calculated. Kappa test was used to test agreements in morphology and severity of perfusion defects assessed on SPECT-CT and on DECT iodine maps with a one-month interval. Wilcoxon rank sum test was used to compare the sharpness of perfusion defects and radiation dose among modalities.

RESULTS

Strong correlations for relative lobar perfusion using linear regression analysis and Pearson's correlation coefficient (r=0.93) were found. Bland-Altman analysis revealed a -0.10 bias, with limits of agreement between [-6.01; 5.81]. With respect to SPECT- CT as standard of reference, the sensitivity, specificity, PPV, NPV, accuracy for lobar perfusion defects were 89.4% (95% 
CI: 82.6-93.4%), 96.5% (95% CI: 92.1-98.5%), 95.6% (95% CI: 
90.9-97.8%), 91.4% (95% CI: 85.6-94.9%) and 93.0% (95% CI: 
87.6-96.1%) respectively. High level of agreement was found for morphology and severity of perfusion defects between modalities (Kappa=0.84 and 0.86 respectively) and on DECT images among readers (Kappa=0.94 and 0.89 respectively). A significantly sharper delineation of perfusion defects was found on DECT images (P<0.0001) using a significantly lower equivalent dose of 4.1±2.3 (SD) mSv (range: 1.9-11.85mSv) compared to an equivalent dose of 5.3±1.1 (SD) mSv (range: 2.8-7.3mSv) for SPECT-CT, corresponding to a 21.2% dose reduction (P=0.0004).

CONCLUSION

DECT imaging shows strong quantitative correlations and qualitative agreements with SPECT-CT for the evaluation of lung perfusion.

Comprehensive SPECT/CT system characterization and calibration for 177Lu quantitative SPECT (QSPECT) with dead-time correction

Background

Personalization of ¹⁷⁷Lu-based radionuclide therapy requires implementation of dosimetry methods that are both accurate and practical enough for routine clinical use. Quantitative single-photon emission computed tomography/computed tomography (QSPECT/CT) is the preferred scanning modality to achieve this and necessitates characterizing the response of the camera, and calibrating it, over the full range of therapeutic activities and system capacity. Various methods to determine the camera calibration factor (CF) and the deadtime constant (τ) were investigated, with the aim to design a simple and robust protocol for quantitative ¹⁷⁷Lu imaging.

Methods

The SPECT/CT camera was equipped with a medium energy collimator. Multiple phantoms were used to reproduce various attenuation conditions: rod sources in air or water-equivalent media, as well as a Jaszczak phantom with inserts. Planar and tomographic images of a wide range of activities were acquired, with multiple energy windows for scatter correction (double or triple energy window technique) as well as count rate monitoring over a large spectrum of energy. Dead time was modelled using the paralysable model. CF and τ were deduced by curve fitting either separately in two steps (CF determined first using a subset of low-activity acquisitions, then τ determined using the full range of activity) or at once (both CF and τ determined using the full range of activity). Total or segmented activity in the SPECT field of view was computed. Finally, these methods were compared in terms of accuracy to recover the known activity, in particular when planar-derived parameters were applied to the SPECT data.

Results

The SPECT camera was shown to operate as expected on a finite count rate range (up to ~ 350 kcps over the entire energy spectrum). CF and τ from planar (sources in air) and SPECT segmented Jaszczak data yielded a very good agreement (CF < 1% and τ < 3%). Determining CF and τ from a single curve fit made dead-time-corrected images less prone to overestimating recovered activity. Using triple-energy window scatter correction while acquiring one or more additional energy window(s) to enable wide-spectrum count rate monitoring (i.e. ranging 55–250 or 18–680 keV) yielded the most consistent results across the various geometries. The final, planar-derived calibration parameters for our system were a CF of 9.36 ± 0.01 cps/MBq and a τ of 0.550 ± 0.003 μs. Using the latter, the activity in a Jaszczak phantom could be quantified by QSPECT with an accuracy of 0.02 ± 1.10%.

Conclusions

Serial planar acquisitions of sources in air using an activity range covering the full operational capacity of the SPECT/CT system, with multiple energy windows for wide-spectrum count rate monitoring, and followed by simultaneous determination of CF and τ using a single equation derived from the paralysable model, constitutes a practical method to enable accurate dead-time-corrected QSPECT imaging in a post-¹⁷⁷Lu radionuclide therapy setting.

Quantitative SPECT/CT-Technique and Clinical Applications

The continuous development of SPECT over the past 50 years has led to improved image quality and increased diagnostic confidence. The most influential developments include the realization of hybrid SPECT/CT devices, as well as the implementation of attenuation correction and iterative image reconstruction techniques. These developments have led to a preference for SPECT/CT devices over SPECT-only systems and to the widespread adoption of the former, strengthening the role of SPECT/CT as the workhorse of Nuclear Medicine imaging. New trends in the ongoing development of SPECT/CT are diverse. For example, whole-body SPECT/CT images, consisting of acquisitions from multiple consecutive bed positions in the manner of PET/CT, are increasingly performed. Additionally, in recent years, some interesting approaches in detector technology have found their way into commercial products. For example, some SPECT cameras dedicated to specific organs employ semiconductor detectors made of cadmium telluride or cadmium zinc telluride, which have been shown to increase the obtainable image quality by offering a higher sensitivity and energy resolution. However, the advent of quantitative SPECT/CT which, like PET, can quantify the amount of tracer in terms of Bq/mL or as a standardized uptake value could be regarded as most important development. It is a major innovation that will lead to increased diagnostic accuracy and confidence, especially in longitudinal studies and in the monitoring of treatment response. The current work comprises two main aspects. At first, physical and technical fundamentals of SPECT image formation are described and necessary prerequisites of quantitative SPECT/CT are reviewed. Additionally, the typically achievable quantitative accuracy based on reports from the literature is given. Second, an extensive list of studies reporting on clinical applications of quantitative SPECT/CT is provided and reviewed.

Analysis of quantitative [I-123] mIBG SPECT/CT in a phantom and in patients with neuroblastoma

Purpose

To determine the accuracy of quantitative SPECT, intersystem and interpatient standardized uptake value (SUV) calculation consistency for a manufacturer-independent quantitative SPECT/CT reconstruction algorithm, and the range of SUVs of normal and neoplastic tissue.

Methods

A NEMA body phantom with 6 spheres (ranging 10–37 mm) was filled with a known activity-to-volume ratio and used to determine the contrast recovery coefficient (CRC) for each visible sphere, and the measured SUV accuracy of those spheres and background water solution. One hundred eleven 123I-metaiodobenzylguanidine ([I-123] mIBG) SPECT/CT examinations from 43 patients were reconstructed using SUV SPECT® (HERMES Medical Solutions Inc.); 42 examinations were acquired using a GE Infinia Hawkeye 4 SPECT/CT, and 69 were acquired on a Siemens Symbia Intevo SPECT/CT. Inter scanner SUV analysis of 9 regions of normal [I-123] mIBG tissue uptake was conducted. Intrapatient mean SUV variability was calculated by measuring normal liver uptake within patients scanned on both cameras. The intensity of uptake by neoplastic tissue in the images was quantified using maximum SUV and, if present, compared over time.

Results

The phantom results of the visible spheres and background resulted in accuracy calculations better than 5–10% with CRC correction. Interscanner SUV variability showed no statistical difference (average p value 0.559; range 0.066–1.0) among the 9 normal tissues analyzed. Intrapatient liver mean SUV varied ≤ 16% as calculated for 28 patients (87 examinations) studied on both scanners. In one patient, a thoracic tumor evaluated over 10 time points (18 months) underwent a 74% (3.1/12.0) reduction in maximum SUV with treatment.

Conclusion

The results demonstrate quantitative accuracy to better than 10%, and both consistent SUV calculation between 2 different SPECT/CT scanners for 9 tissues, and low intrapatient measurement variability for quantitative SPECT/CT analysis in a pediatric population with neuroblastoma. Quantitative SPECT/CT offers the opportunity for objective analysis of tumor response using [I-123] mIBG by normalizing the uptake to injected dose and patient weight, as is done for PET.

What You See Is Not What You Get: On the Accuracy of Voxel-Based Dosimetry in Molecular Radiotherapy

Improvements in quantitative SPECT/CT have aroused growing interest in voxel-based dosimetry for radionuclide therapies, because it promises visualization of absorbed doses at a voxel level. In this work, SPECT/CT-based voxel-level dosimetry of a 3-dimensional (3D) printed 2-compartment kidney phantom was performed, and the resulting absorbed dose distributions were examined. Additionally, the potential of the PETPVC partial-volume correction tool was investigated. Methods: Both kidney compartments (70% cortex, 30% medulla) were filled with different activity concentrations, and SPECT/CT imaging was performed. The images were reconstructed using varying settings (iterations, subsets, and postfiltering). On the basis of these activity concentration maps, absorbed dose distributions were calculated with precalculated ¹⁷⁷Lu voxel S values and an empiric kidney half-life. An additional set of absorbed doses was calculated after applying PETPVC for partial-volume correction of the SPECT reconstructions. Results: SPECT/CT imaging blurs the 2 discrete suborgan absorbed dose values into a continuous distribution. Although this effect is slightly improved by applying more iterations, it is enhanced by additional postfiltering. By applying PETPVC, the absorbed dose values are separated into 2 peaks. Although this leads to a better agreement between SPECT/CT-based and nominal values, considerable discrepancies remain. In contrast to the calculated nominal absorbed doses of 7.8 and 1.6 Gy (in the cortex and medulla, respectively), SPECT/CT-based voxel-level dosimetry resulted in mean absorbed doses of 3.0-6.6 Gy (cortex) and 2.7-5.1 Gy (medulla). PETPVC led to improved ranges of 6.1-8.9 Gy (cortex) and 2.1-5.4 Gy (medulla). Conclusion: Our study showed that ¹⁷⁷Lu quantitative SPECT/CT imaging leads to voxel-based dose distributions largely differing from the real organ distribution. SPECT/CT imaging and reconstruction deficiencies might directly translate into unrealistic absorbed dose distributions, thus questioning the reliability of SPECT-based voxel-level dosimetry. Therefore, SPECT/CT reconstructions should be adapted to ensure an accurate quantification of the underlying activity and, therefore, absorbed dose in a volume of interest of the expected object size (e.g., organs, organ substructures, lesions, or voxels). As an example, PETPVC largely improves the match between SPECT/CT-based and nominal dose distributions. In conclusion, the concept of voxel-based dosimetry should be treated with caution. Specifically, one should remember that the absorbed dose distribution is mainly a convolved version of the underlying SPECT reconstruction.

Protocol-driven multidetector SPECT/CT: integration of hybrid imaging into the routine workflow of whole-body bone scintigraphy in oncology patients

AIM

To analyse the additional clinical value of protocol-driven and selective use of multidetector single-photon-emission tomography/computed tomography (SPECT/CT) in oncology patients undergoing whole-body bone scintigraphy (BS) and to analyse reporter confidence in diagnosis with and without SPECT/CT.

MATERIALS AND METHODS

During a 2-year period, 2,954 whole-body BS examinations were performed in oncology patients, with 444 (15%) undergoing additional protocol-driven SPECT/CT. Retrospective evaluation of planar BS and SPECT/CT images was performed by two experienced dual-trained nuclear medicine radiologists. The BS and SPECT/CT images were graded blindly using a five-point scale designed to evaluate the likelihood of a lesion being benign or malignant. Interpretation was applied on a per-patient basis.

RESULTS

There was a 74.5% increase in definitive diagnostic classification and a 26.6% reduction in equivocal findings with SPECT/CT when compared to BS alone (p<0001). Of cases initially classified as "probably benign" on BS, 5.1% (10/193) were reclassified to "probably malignant" (1%) or "malignant" (4.1%) using the SPECT/CT data. The highest impact in reporter confidence was seen with SPECT/CT in the interpretation of lesions within the pelvis (34%), ribs (23%), lumbar spine (22%), and thoracic spine (21%).

CONCLUSION

Protocol-driven, selective use of SPECT/CT imaging to augment planar BS reduces equivocal findings and improves reporter confidence whilst minimising the impact on patient and reporting workflows.

A Brief History of Lung Ventilation and Perfusion Imaging Over the 50-Year Tenure of the Editors of Seminars in Nuclear Medicine

The ventilation/perfusion lung scan has been in continuous use for approximately half a century, the same lifetime as Seminars in Nuclear Medicine. Remarkably, the founding Editors-in-Chief have continued to guide the journal over this entire period. In this Feschrift issue celebrating their enormous contribution, we review the history of the lung scan, its highs and lows, the transition from planar to SPECT/CT V/Q scans, and the future that is in store in this age of multimodality functional imaging. We concur with the published view of one of the retiring editors (LMF) that V/Q scintigraphy is indeed alive and well and has a definite future in clinical medicine.

Feasibility and accuracy of single time point imaging for renal dosimetry following 177Lu-DOTATATE ('Lutate') therapy

Background

This study aims to assess both feasibility and accuracy of renal dosimetry imaging protocols in patients receiving Lutate therapy for neuroendocrine tumours (NETs), when data acquisition over multiple days is not possible on all cycles.

Method

Patients who had received a full 4 cycles of Lutate therapy with complete imaging at each cycle were included. Imaging consisted of quantitative SPECT/CT of the kidneys at 4, 24 and 96–120 h post injection. Renal absorbed dose was calculated for each data set, and in addition, five alternative methods were explored for comparison. Method 1: a patient average clearance time (t_1/2 average) derived from the first half of contributing patient data was used to estimate absorbed dose for subsequent patients based on 4 h imaging alone; method 2: t_1/2 average was applied to subsequent patients on 24 h imaging alone; method 3: a patient-specific clearance rate (t_1/2 patient) was determined from complete image data of cycle 1 and applied subsequently to remaining cycles using 4 h image data alone; method 4: t_1/2 patient was applied to 24 h imaging alone in subsequent cycles; method 5: the 120 h data was estimated on subsequent cycles based on the cycle 1 fraction of injected activity (%IA) at 24 and 120 h.

Results

Twenty treatments from 18 patients, resulting in 80 cycles of therapy, were analysed. The measured average renal absorbed dose per cycle of treatment was 0.38 ± 0.19 Gy/GBq when derived from full imaging data. The use of t_1/2 average applied to a single time point led to large deviations of dose estimates from true values (on average 59% and 30%, when using 4 h data and 24 h data, respectively). The use of complete image data on cycle 1 and the derivation of t_1/2 patient led to improved dose estimates, with an average deviation from true values of 13% and 2% when using 4 h data only and 24 h data only, respectively. The use of a 120 h %IA derived from cycle 1 led to an average deviation from true dose estimates of 14%.

Conclusion

In instances where demands on both patients and facilities make multiple time point data acquisition impractical, renal dosimetry is best derived through complete imaging at cycle 1 only followed by a single 24 h imaging time point on subsequent cycles, assuming no significant changes in renal function during the time course of therapy.

A Comparison of 2D and 3D Kidney Absorbed Dose Measures in Patients Receiving 177Lu-DOTATATE

Objective(s):

To investigate and compare quantitative accuracy of kidney absorbed dose measures made from both 2D and 3D imaging in patients receiving ¹⁷⁷Lu-DOTATATE (Lutate) for treatment of neuroendocrine tumours (NETs).

Methods:

Patients receiving Lutate therapy underwent both whole body planar imaging and SPECT/CT imaging over the kidneys at time points 0.5, 4, 24, and 96-120 hours after injection. Planar data were corrected for attenuation using transmission data, and were converted to units of absolute activity via two methods, using either a calibration standard in the field of view or relative to pre-voiding image total counts. Hand drawn regions of interest were used to generate time activity curves and kidney absorbed dose estimates in OLINDA-EXM. Fully quantitative SPECT data were generated using CT-derived corrections for both scatter and attenuation, before correction for dead time and application of a camera specific sensitivity factor to convert data to units of absolute activity. Volumes of interest were defined for kidney using the co-registered x-ray CT, before time activity curves and absorbed dose measures were generated in OLINDA-EXM, both with and without corrections made to the model for patient specific kidney volumes. Quantitative SPECT data were also used to derive dose maps through dose kernel convolution (DKC), which was treated as the gold standard.

Results:

A total of 50 studies were analysed, corresponding to various cycles of treatment from 21 patients. Planar absorbed dose estimates were consistently higher than SPECT derived estimates by, on average, a factor of 3.

Conclusion:

Quantitative SPECT is considered the gold standard approach for organ specific dosimetry however often relies on in house software. As such planar methods for estimating absorbed dose are much more widely available, and in particular, are often the only source of reference in previously published data. For the case of Lutate dosimetry, planar measures may lead to a three-fold increase in measures of kidney absorbed dose.

Image-based SPECT calibration based on the evaluation of the Fraction of Activity in the Field of View

Background

SPECT quantification is important for dosimetry in targeted radionuclide therapy (TRT) and the calibration of SPECT images is a crucial stage for image quantification. The current standardized calibration protocol (MIRD 23) uses phantom acquisitions to derive a global calibration factor in specific conditions. It thus requires specific acquisitions for every clinical protocols. We proposed an alternative and complementary image-based calibration method that allows to determine a calibration factor adapted to each patient, radionuclide, and acquisition protocol and that may also be used as an additional independent calibration.

Results

The proposed method relies on a SPECT/CT acquisition of a given region of interest and an initial whole-body (WB) planar image. First, the conjugate view of WB planar images is computed after scatter and attenuation correction. 3D SPECT images are reconstructed with scatter, attenuation, and collimator-detector response (CDR) corrections and corrected from apparent dead-time. The field of view (FOV) of the SPECT image is then projected on the corrected WB planar image. The fraction of activity located in the area corresponding to the SPECT FOV is then calculated based on the counts on the corrected WB planar image. The Fraction of Activity in Field Of View (FAF) is then proposed to compute the calibration factor as the total number of counts in the SPECT image divided by this activity. Quantification accuracy was compared with the standard calibration method both with phantom experiments and on patient data.

Both standard and image-based calibrations give good accuracy on large region of interest on phantom experiments (less than 7% of relative difference compared to ground truth). Apparent dead-time correction allows to reduce the uncertainty associated with standard calibration from 2.5 to 1%. The differences found between both methods were lower than the uncertainty range of the standard calibration (<3%). In patient data, although no ground truth was available, both methods give similar calibration factor (average difference 3.64%).

Conclusions

A calibration factor may be computed directly from the acquired SPECT image providing that a WB planar image is also available and if both acquisitions are performed before biological elimination. This method does not require to perform phantom acquisition for every different acquisition conditions and may serve to double check the calibration with an independent factor.

Determination of gamma camera calibration factors for quantitation of therapeutic radioisotopes

Background

Camera calibration, which translates reconstructed count map into absolute activity map, is a prerequisite procedure for quantitative SPECT imaging. Both planar and tomographic scans using different phantom geometries have been proposed for the determination of the camera calibration factor (CF). However, there is no consensus on which approach is the best. The aim of this study is to evaluate all these calibration methods, compare their performance, and propose a practical and accurate calibration method for SPECT quantitation of therapeutic radioisotopes. Twenty-one phantom experiments (Siemens Symbia SPECT/CT) and 12 Monte Carlo simulations (GATE v6.1) using three therapy isotopes (¹³¹I, ¹⁷⁷Lu, and ¹⁸⁸Re) have been performed. The following phantom geometries were used: (1) planar scans of point source in air (PS), (2) tomographic scans of insert(s) filled with activity placed in non-radioactive water (HS + CB), (3) tomographic scans of hot insert(s) in radioactive water (HS + WB), and (4) tomographic scans of cylinders uniformly filled with activity (HC). Tomographic data were reconstructed using OSEM with CT-based attenuation correction and triple energy window (TEW) scatter correction, and CF was determined using total counts in the reconstructed image, while for planar scans, the photopeak counts, corrected for scatter and background with TEW, were used. Additionally, for simulated data, CF obtained from primary photons only was analyzed.

Results

For phantom experiments, CF obtained from PS and HS + WB agreed to within 6% (below 3% if experiments performed on the same day are considered). However, CF from HS + CB exceeded those from PS by 4–12%. Similar trend was found in simulation studies. Analysis of CFs from primary photons helped us to understand this discrepancy. It was due to underestimation of scatter by the TEW method, further enhanced by attenuation correction. This effect becomes less important when the source is distributed over the entire phantom volume (HS + WB and HC).

Conclusions

Camera CF could be determined using planar scans of a point source, provided that the scatter and background contributions are removed, for example using the clinically available TEW method. This approach is simple and yet provides CF with sufficient accuracy (~ 5%) to be used in clinics for radiotracer quantification.

Variations of clinical SPECT/CT operations

Aim: This survey gathers information about clinical SPECT/CT operations worldwide to help guide standardization of clinical SPECT/ CT imaging. Methods: An international, webbased survey of SPECT/CT users was initiated in 12/2010 through an e-mail distribution. Users were asked 71 questions related to (A) demographics, (B) SPECT/CT operations/utilization and (C) variations in imaging protocols. Results: Collected responses originated from 117 imaging centers in the Americas (66%), Europe (20%), Asia-Pacific (11%) and the Middle-East (3%), with the majority of responding sites representing public health care institutions (69%). Most sites operate 1–2 SPECT/CT-systems (74%), typically installed in Nuclear Medicine departments (84%) with extensive prior SPECT-only experience (82%). Only 14% of SPECT/CTs are installed in Radiology departments. Clinical SPECT/CT imaging is performed either as routine (51%) or ad-hoc “add-on” procedure (49%) with a high inter-site and inter-examination variability. The main application of the integrated CT is to provide anatomical localization of the tracer uptake rather than to produce contrast enhanced or other high-quality CT images. Consequently, in only 22% of the sites a CT contrast injector is installed. Only 6% of centers use SPECT/CT devices for stand-alone CT procedures. Conclusion: An international survey among clinical SPECT/CT users revealed that SPECT/CT is a not a routine component of nuclear medicine procedures. The majority of the centers responding do not fully utilize the diagnostic potential of the CT components. Significant variations in standard imaging protocols were observed. These findings illustrate the need for training and standardization and underscore the need for revisiting the role of SPECT/CT in diagnostic imaging.

Optimizing Image Quantification for 177Lu SPECT/CT Based on a 3D Printed 2-Compartment Kidney Phantom

The aim of this work was to find an optimal setup for activity determination of ¹⁷⁷Lu-based SPECT/CT imaging reconstructed with 2 commercially available methods (xSPECT Quant and Flash3D). For this purpose, 3-dimensional (3D)-printed phantoms of different geometries were manufactured, different partial-volume correction (PVC) methods were applied, and the accuracy of the activity determination was evaluated. Methods: A 2-compartment kidney phantom (70% cortical and 30% medullary compartment), a sphere, and an ellipsoid of equal volumes were 3D printed, filled with ¹⁷⁷Lu, and scanned with a SPECT/CT system. Reconstructions were performed with xSPECT and Flash3D. Different PVC methods were applied to find an optimal quantification setup: method 1 was a geometry-specific recovery coefficient based on the 3D printing model, method 2 was a geometry-specific recovery coefficient based on the low-dose CT scan, method 3 was an enlarged volume of interest including spilled-out counts, method 4 was activity concentration in the peak milliliter applied to the entire CT-based volume, and method 5 was a fixed threshold of 42% of the maximum in a large volume containing the object of interest. Additionally, the influence of postreconstruction gaussian filtering was investigated. Results: Although the recovery coefficients of sphere and ellipsoid differed by only 0.7%, a difference of 31.7% was observed between the sphere and the renal cortex phantoms. Without postfiltering, the model-based recovery coefficients (methods 1 and 2) resulted in the best accuracies (xSPECT, 1.5%; Flash3D, 10.3%), followed by the enlarged volume (method 3) (xSPECT, 8.5%; Flash3D, 13.0%). The peak-milliliter method (method 4) showed large errors only for sphere and ellipsoid (xSPECT, 23.4%; Flash3D, 21.6%). Applying a 42% threshold (method 5) led to the largest quantification errors (xSPECT, 32.3%; Flash3D, 46.7%). After postfiltering, a general increase in the errors was observed. Conclusion: In this work, 3D printing was used as a prototyping technique for a geometry-specific investigation of SPECT/CT reconstruction parameters and PVC methods. The optimal setup for activity determination was found to be an unsmoothed SPECT/CT reconstruction in combination with a recovery coefficient based on the low-dose CT. The difference between spheric and renal recovery coefficients suggests that the typically applied volume-dependent but only sphere-based recovery coefficient lookup tables should be replaced by a more geometry-specific alternative.

Validation of a calibration method using the cross-calibration factor and system planar sensitivity in quantitative single-photon emission computed tomography imaging

The present study aimed to validate the absolute quantitative accuracy of a calibration method for single-photon emission computed tomography (SPECT) using cross-calibration factor (CCF)- and system sensitivity-based calibration methods. The CCF obtained with different reconstruction parameters was evaluated using a cylindrical phantom (diameter 20 cm, height 20 cm). SPECT images were acquired with a positron emission tomography/computed tomography (CT) phantom. Subsequently, they were reconstructed by using ordered subset expectation maximization with resolution recovery, scatter, and CT-based attenuation correction. All reconstructed SPECT counts were converted to activity concentrations based on the CCF and system planar sensitivity. We placed 12 circular regions of interest, 37 mm in diameter, on the phantom background, and the converted activity concentration and relative measurement error were assessed. The CCF obtained using a cylindrical phantom was affected by the iterative update number and post-smoothing filter function. The activity concentration calibrated using the CCF showed over- and underestimation. However, the activity concentration obtained from the system planar sensitivity was similar to that gained using the phantom. The values obtained using the system planar sensitivity were within 10% of the activity concentrations obtained with the phantom. These findings demonstrated that the calibration method using system planar sensitivity provides accurate quantification within 10% of the true activity concentration. Further clinical examination is required to validate the present results.

A Prospective Study Comparing 99mTc-Hydroxyethylene-Diphosphonate Planar Bone Scintigraphy and Whole-Body SPECT/CT with 18F-Fluoride PET/CT and 18F-Fluoride PET/MRI for Diagnosing Bone Metastases

We prospectively evaluated and compared the diagnostic performance of ^99mTc-hydroxyethylene-diphosphonate (^99mTc-HDP) planar bone scintigraphy (pBS), ^99mTc-HDP SPECT/CT, ¹⁸F-NaF PET/CT, and ¹⁸F-NaF PET/MRI for the detection of bone metastases. Methods: One hundred seventeen patients with histologically proven malignancy referred for clinical pBS were prospectively enrolled. pBS and whole-body SPECT/CT were performed followed by ¹⁸F-NaF PET/CT within 9 d. ¹⁸F-NaF PET/MRI was also performed in 46 patients. Results: Bone metastases were confirmed in 16 patients and excluded in 101, which was lower than expected. The number of equivocal scans was significantly higher for pBS than for SPECT/CT and PET/CT (18 vs. 5 and 6, respectively; P = 0.004 and 0.01, respectively). When equivocal readings were excluded, no statistically significant difference in sensitivity, specificity, positive predictive value, negative predictive value, or overall accuracy were found when comparing the different imaging techniques. In the per-patient analysis, equivocal scans were either assumed positive for metastases ("pessimistic analysis") or assumed negative for metastases ("optimistic analysis"). The percentages of misdiagnosed patients for the pessimistic analysis were 21%, 15%, 9%, and 7% for pBS, SPECT/CT, PET/CT, and PET/MRI, respectively. Corresponding figures for the optimistic analysis were 9%, 12%, 5%, and 7%. In those patients identified as having bone metastases according to the reference standard, SPECT/CT, ¹⁸F-NaF PET/CT, and PET/MRI detected additional lesions compared with pBS in 31%, 63%, and 71%, respectively. Conclusion:¹⁸F-NaF PET/CT and whole-body SPECT/CT resulted in a significant reduction of equivocal readings compared with pBS, which implies an improved diagnostic confidence. However, the clinical benefit of using, for example, ¹⁸F-NaF PET/CT or PET/MRI as compared with SPECT/CT and pBS in this patient population with a relatively low prevalence of bone metastases (14%) is likely limited. This conclusion is influenced by the low prevalence of patients with osseous metastases. There may well be significant differences in the sensitivity of SPECT/CT, PET/CT, and PET/MRI compared with pBS, but a larger patient population or a patient population with a higher prevalence of bone metastases would have to be studied to demonstrate this.

Fast analytical approach of application specific dose efficient spectrum selection for diagnostic CT imaging and PET attenuation correction

Computed tomography (CT) has been used for a variety of applications, two of which include diagnostic imaging and attenuation correction for PET or SPECT imaging. Ideally, the x-ray tube spectrum should be optimized for the specific application to minimize the patient radiation dose while still providing the necessary information. In this study, we proposed a projection-based analytic approach for the analysis of contrast, noise, and bias. Dose normalized contrast to noise ratio (CNRD), inverse noise normalized by dose (IND) and bias are used as evaluation metrics to determine the optimal x-ray spectrum. Our simulation investigated the dose efficiency of the x-ray spectrum ranging from 40 kVp to 200 kVp. Water cylinders with diameters of 15 cm, 24 cm, and 35 cm were used in the simulation to cover a variety of patient sizes. The effects of electronic noise and pre-patient copper filtration were also evaluated. A customized 24 cm CTDI-like phantom with 13 mm diameter inserts filled with iodine (10 mg ml^-1), tantalum (10 mg ml^-1), water, and PMMA was measured with both standard (1.5 mGy) and ultra-low (0.2 mGy) dose to verify the simulation results at tube voltages of 80, 100, 120, and 140 kVp. For contrast-enhanced diagnostic imaging, the simulation results indicated that for high dose without filtration, the optimal kVp for water contrast is approximately 100 kVp for a 15 cm water cylinder. However, the 60 kVp spectrum produces the highest CNRD for bone and iodine. The optimal kVp for tantalum has two selections: approximately 50 and 100 kVp. The kVp that maximizes CNRD increases when the object size increases. The trend in the CTDI phantom measurements agrees with the simulation results, which also agrees with previous studies. Copper filtration improved the dose efficiency for water and tantalum, but reduced the iodine and bone dose efficiency in a clinically-relevant range (70-140 kVp). Our study also shows that for CT-based attenuation correction applications for PET or SPECT, a higher-kVp spectrum with copper filtration is preferable. This method is developed based on filter back projection and does not require image reconstruction or Monte Carlo dose estimates; thus, it could potentially be used for patient-specific and task-based on-the-fly protocol optimization.

Technical Pitfalls and Limitations of SPECT/CT

The synergy of functional and anatomic information in hybrid systems has undoubtedly enhanced the diagnostic potential of radionuclide imaging in recent years, contributing to the advancement of SPECT/CT in clinical practice. Since the introduction of commercial SPECT/CT in the late 1990 s, the field has seen rapid expansion and development toward multidetector CT subsystems, establishing the role of SPECT/CT as a routine imaging tool. It is, however, important to discuss possible challenges and technical limitations of such systems and how these influence imaging outcomes. In particular, the issues of patient motion and spatial misalignment of the SPECT and CT modalities, data corrections such as those for photon attenuation, and the choice of CT acquisition protocols in relation to radiation exposure are discussed in the article.

Quantitative SPECT/CT Imaging of (177)Lu with In Vivo Validation in Patients Undergoing Peptide Receptor Radionuclide Therapy

PURPOSE

The purpose of this study is to extend an established SPECT/CT quantitation protocol to (177)Lu and validate it in vivo using urine samples, thus providing a basis for 3D dosimetry of (177)Lu radiotherapy and improvement over current planar methods which improperly account for anatomical variations, attenuation, and overlapping organs.

PROCEDURES

In our quantitation protocol, counts in images reconstructed using an ordered subset-expectation maximization algorithm are converted to kilobecquerels per milliliter using a calibration factor derived from a phantom experiment. While varying reconstruction parameters, we tracked the ratio of image to true activity concentration (recovery coefficient, RC) in hot spheres and a noise measure in a homogeneous region. The optimal parameter set was selected as the point where recovery in the largest three spheres (16, 8, and 4 ml) stagnated, while the noise continued to increase. Urine samples were collected following 12 SPECT/CT acquisitions of patients undergoing [(177)Lu]DOTATATE therapy, and activity concentrations were measured in a well counter. Data was reconstructed using parameters chosen in the phantom experiment, and estimated activity concentration from the images was compared to the urine values to derive RCs.

RESULTS

In phantom data, our chosen parameter set yielded RCs in 16, 8, and 4 ml spheres of 80.0, 74.1, and 64.5 %, respectively. For patients, the mean bladder RC was 96.1 ± 13.2% (range, 80.6-122.4 %), with a 95 % confidence interval between 88.6 and 103.6 %. The mean error of SPECT/CT concentrations was 10.1 ± 8.3% (range, -19.4-22.4 %).

CONCLUSIONS

Our results show that quantitative (177)Lu SPECT/CT in vivo is feasible but could benefit from improved reconstruction methods. Quantifying bladder activity is analogous to determining the amount of activity in the kidneys, an important task in dosimetry, and our results provide a useful benchmark for future efforts.

[Development of a Novel Body Phantom with Bone Equivalent Density for Evaluation of Bone SPECT]

We developed a custom-designed phantom for bone single photon emission computed tomography (SPECT)-specific radioactivity distribution and linear attenuation coefficient. The aim of this study was to evaluate the accuracy of the phantom. The lumbar phantom consisted of the trunk of a body phantom (background) containing a cylinder (vertebral body), a sphere (tumor), and a T-shaped container (processus). The vertebral body, tumor, and processus phantoms contained a K(2)HPO(4) solution of bone equivalent density and 50, 300 and 50 kBq/mL of (99m)Tc, respectively. The body phantom contained 8 kBq/mL of (99m)Tc solution. SPECT images were acquired using low-energy high-resolution collimation, a 128 × 128 matrix and 120 projections over 360° with a dwell time of 15 sec/view × 4 times. Thereafter, CT images were acquired at 130 kV and 70 ref mAs using adaptive dose modulation. The SPECT data were reconstructed with ordered subset expectation maximization with three-dimensional, scatter, and CT-based attenuation correction. Count ratio, linear attenuation coefficient (LAC), and full-width at half-maximum (FWHM) were measured. Count ratios between the background, the vertebral body, and the tumor in SPECT images were 463.8: 2888.0: 15150.3 (1: 6.23: 32.7). The LAC of the background and vertebral body in the CT-derived attenuation map were 0.155 cm⁻¹ and 0.284 cm⁻¹, respectively, and the FWHM measured from the processus was 15.27 mm. The precise counts and LAC indicated that the phantom was accurate and could serve as a tool for evaluating acquisition, reconstruction parameters, and quantitation in bone SPECT images.

Intravenous contrast-enhanced CT can be used for CT-based attenuation correction in clinical (111)In-octreotide SPECT/CT

BACKGROUND

CT-based attenuation correction (CT-AC) using contrast-enhancement CT impacts (111)In-SPECT image quality and quantification. In this study we assessed and evaluated the effect.

METHODS

A phantom (5.15 L) was filled with an aqueous solution of In-111. Three SPECT/CT scans were performed: (A) no IV contrast, (B) with 100-mL IV contrast, and (C) with 200-mL IV contrast added. Scan protocol included a localization CT, a low-dose CT (LD), and a full-dose CT (FD). Phantom, LD and FD scan series were performed at 90, 120, and 140 kVp. Phantom data were evaluated looking at mean counts in a central volume. Ten patients referred for (111)In-octreotide scintigraphy were scanned according to our clinical (111)In-SPECT/CT protocol including a topogram, a LD (140 kVp), and a FD (120 kVp). The FD/contrast-enhanced CT was acquired in both arterial (FDAP) and venous phase (FDVP) following a mono-phasic IV injection of 125-mL Optiray (4.5 mL/s). For patient data, we report image quality, Krenning scores, and mean/max values for liver and tumor regions.

RESULTS

Phantoms: in uncorrected emission data, mean counts (average ± SD) decreased with increasing IV concentration: (A) 119 ± 9, (B) 113 ± 8, and (C) 110 ± 9. For all attenuation correction (AC) scans, the mean values increased with increasing iodine concentration.

PATIENTS

there were no visible artifacts in single photon emission computed tomography (SPECT) following CT-AC with contrast-enhanced CT. The average score of image quality was 4.1 ± 0.3, 3.8 ± 0.4, and 4.2 ± 0.4 for LD, arterial phase, and venous phase, respectively. A total of 16 lesions were detected. The Krenning scores of 13/16 lesions were identical across all scan series. The max pixel values for the 16 lesions showed generally lower values for LD than for contrast-enhanced CT.

CONCLUSIONS

In (111)In-SPECT/CT imaging of phantoms and patients, the use of IV CT contrast did neither degrade the SPECT image quality nor affect the clinical Krenning score. Reconstructed counts in healthy liver tissues were unaffected, and there was a generally lower count value in lesions following CT-AC based on the LD non-enhanced images. Overall, for clinical interpretation, no separate low-dose CT is required for CT-AC in (111)In-SPECT/CT.

In vivo quantification of (177)Lu with planar whole-body and SPECT/CT gamma camera imaging

Background

Advances in gamma camera technology and the emergence of a number of new theranostic radiopharmaceutical pairings have re-awakened interest in in vivo quantification with single-photon-emitting radionuclides. We have implemented and validated methodology to provide quantitative imaging of ¹⁷⁷Lu for 2D whole-body planar studies and for 3D tomographic imaging with single-photon emission computed tomography (SPECT)/CT.

Methods

Whole-body planar scans were performed on subjects to whom a known amount of [¹⁷⁷Lu]-DOTA-octreotate had been administered for therapy. The total radioactivity estimated from the images was compared with the known amount of the radionuclide therapy administered. In separate studies, venous blood samples were withdrawn from subjects after administration of [¹⁷⁷Lu]-DOTA-octreotate while a SPECT acquisition was in progress and the concentration of the radionuclide in the venous blood sample compared with that estimated from large blood pool structures in the SPECT reconstruction. The total radioactivity contained within an internal SPECT calibration standard was also assessed.

Results

In the whole-body planar scans (n = 28), the estimated total body radioactivity was accurate to within +4.6 ± 5.9 % (range −17.1 to +11.2 %) of the correct value. In the SPECT reconstructions (n = 12), the radioactivity concentration in the cardiac blood pool was accurate to within −4.0 ± 7.8 % (range −16.1 to +7.5 %) of the true value and the internal standard measurements (n = 89) were within 2.0 ± 8.5 % (range −16.3 to +24.2 %) of the known amount of radioactivity contained.

Conclusions

In our hands, state-of-the-art hybrid SPECT/CT gamma cameras were able to provide accurate estimates of in vivo radioactivity to better than, on average, ±10 % for use in biodistribution and radionuclide dosimetry calculations.

Repeatability of Radiotracer Uptake in Normal Abdominal Organs with ¹¹¹In-Pentetreotide Quantitative SPECT/CT

With an increasing emphasis on quantitation of SPECT imaging and its use in dosimetry to guide therapies, it is desirable to understand the repeatability in normal-organ SPECT uptake values (SPECT-UVs). We investigated the variability of normal abdominal organ uptake in repeated (111)In-pentetreotide SPECT studies.

METHODS

Nine patients with multiple (111)In-pentetreotide SPECT/CT studies for clinical purposes were evaluated. Volumes of interest were drawn for the abdominal organs and applied to SPECT-UVs. The variability of those values was assessed.

RESULTS

The average SPECT-UV for the liver (1.7 ± 0.6) was much lower than for the kidneys (right, 8.0 ± 2.4; left, 7.5 ± 1.7). Interpatient and intrapatient variability was similar (intraclass correlation coefficients, 0.40-0.59) for all organs. The average coefficients of variation for each organ for each patient were obtained and averaged across all patients (0.26 for liver, 0.22 for right kidney, and 0.20 for left kidney). The coefficients of variation for the organs across all scans were 0.33 (liver), 0.30 (right kidney), and 0.22 (left kidney).

CONCLUSION

Variability across all patients and all scans for the liver was higher than reported with (18)F-FDG PET, though left kidney variability was similar to PET liver variability and left kidney uptake may be able to serve as an internal metric for determining the quantifiability of an (111)In-pentetreotide SPECT study.

Estimation of radioactivity in single-photon emission computed tomography for sentinel lymph node biopsy in a torso phantom study

OBJECTIVE

The number of lymph nodes to be removed is determined from residual counts. Advance estimation of residual radioactivity in lymphatic nodes before a biopsy is useful for reducing surgical operation time. The purpose of this study was to estimate the total radioactivity of a small hotspot in single-photon emission computed tomography (SPECT) of a torso phantom.

METHODS

A cross-calibration study was performed to convert counts in SPECT images to radioactivity. A simulation study was performed to estimate the size of the volume of interest (VOI) covering a hotspot corrupted with full-width at half-maximum between 8 and 16 mm. The estimation of total radioactivity was validated in a torso phantom study using small sources.

RESULTS

True radioactivity was approximately equal to integrated values of hotspots using the VOI with a diameter of 40 mm in our simulation study. The difference was less than 18% in cases of more than 9.4 kBq.

CONCLUSION

The total radioactivity in small sources simulating a typical sentinel node was estimated from SPECT images using a VOI of 40 mm in a torso phantom study. Because the difference from actual values was less than 10% on average when radioactivities were more than 9.4 kBq, the total radioactivity of a lymph node can be estimated in a clinical examination.

Factors affecting the repeatability of gamma camera calibration for quantitative imaging applications using a sealed source

Several applications in nuclear medicine require absolute activity quantification of single photon emission computed tomography images. Obtaining a repeatable calibration factor that converts voxel values to activity units is essential for these applications. Because source preparation and measurement of the source activity using a radionuclide activity meter are potential sources of variability, this work investigated instrumentation and acquisition factors affecting repeatability using planar acquisition of sealed sources. The calibration factor was calculated for different acquisition and geometry conditions to evaluate the effect of the source size, lateral position of the source in the camera field-of-view (FOV), source-to-camera distance (SCD), and variability over time using sealed Ba-133 sources. A small region of interest (ROI) based on the source dimensions and collimator resolution was investigated to decrease the background effect. A statistical analysis with a mixed-effects model was used to evaluate quantitatively the effect of each variable on the global calibration factor variability. A variation of 1 cm in the measurement of the SCD from the assumed distance of 17 cm led to a variation of 1-2% in the calibration factor measurement using a small disc source (0.4 cm diameter) and less than 1% with a larger rod source (2.9 cm diameter). The lateral position of the source in the FOV and the variability over time had small impacts on calibration factor variability. The residual error component was well estimated by Poisson noise. Repeatability of better than 1% in a calibration factor measurement using a planar acquisition of a sealed source can be reasonably achieved. The best reproducibility was obtained with the largest source with a count rate much higher than the average background in the ROI, and when the SCD was positioned within 5 mm of the desired position. In this case, calibration source variability was limited by the quantum noise.

A method for estimating radiation interaction coefficients for tissues from single energy CT

A parametric model for the x-ray linear attenuation coefficient is used to describe the compositional dependence of Hounsfield numbers measured by medical CT scanners. Measurements with materials of known density and composition, that span and evenly sample the compositional range of tissues, are written as linear simultaneous equations and solved for model coefficients. An algorithm is identified for this purpose. Results are expressed as atomic cross-sections in units of barn per electron divided by the attenuation coefficient for water. With the CT scanner characterised, a virtual CT scan can be simulated to predict HN for tissues based upon their known density and composition. Similar calculations using the tabulations and mixture rule deliver attenuation coefficients and mass energy absorption coefficients for mono-energetic radiation 10 keV to 20 MeV. Results are presented for measurements with a radiotherapy CT simulator, the RMI-467 phantom with tissue substitute materials, plus common polymer materials and silicon. Published measurements with earlier generations of the phantom and tissue substitutes using different CT scanners are also considered. Measured atomic cross-sections differ from expectations for mono-energetic radiation due to the use of a filtered spectrum and energy integrating detection system. The cross-sections for different CT scanners are similar, without large variations with kVp. Results are presented showing the relationship between predicted HN for tissues, electron density and photon interaction coefficients for healthy tissues and mono-energetic radiation. A strategy is suggested for accommodating strongly attenuating materials such as calculi and metallic implants.

Nonlinear dual reconstruction of SPECT activity and attenuation images

In single photon emission computed tomography (SPECT), accurate attenuation maps are needed to perform essential attenuation compensation for high quality radioactivity estimation. Formulating the SPECT activity and attenuation reconstruction tasks as coupled signal estimation and system parameter identification problems, where the activity distribution and the attenuation parameter are treated as random variables with known prior statistics, we present a nonlinear dual reconstruction scheme based on the unscented Kalman filtering (UKF) principles. In this effort, the dynamic changes of the organ radioactivity distribution are described through state space evolution equations, while the photon-counting SPECT projection data are measured through the observation equations. Activity distribution is then estimated with sub-optimal fixed attenuation parameters, followed by attenuation map reconstruction given these activity estimates. Such coupled estimation processes are iteratively repeated as necessary until convergence. The results obtained from Monte Carlo simulated data, physical phantom, and real SPECT scans demonstrate the improved performance of the proposed method both from visual inspection of the images and a quantitative evaluation, compared to the widely used EM-ML algorithms. The dual estimation framework has the potential to be useful for estimating the attenuation map from emission data only and thus benefit the radioactivity reconstruction.

A free-response evaluation determining value in the computed tomography attenuation correction image for revealing pulmonary incidental findings: a phantom study

RATIONALE AND OBJECTIVES

The purpose of this study was to compare lesion-detection performance when interpreting computed tomography (CT) images that are acquired for attenuation correction when performing single photon emission computed tomography/computed tomography (SPECT/CT) myocardial perfusion studies. In the United Kingdom, there is a requirement that these images be interpreted; thus, it is necessary to understand observer performance on these images.

MATERIALS AND METHODS

An anthropomorphic chest phantom with inserted spherical lesions of different sizes and contrasts was scanned on five different SPECT/CT systems using site-specific CT protocols for SPECT/CT myocardial perfusion imaging. Twenty-one observers (0-4 years of CT experience) searched 26 image slices (17 abnormal, containing 1-3 lesions, and 9 normal, containing no lesions) for each CT acquisition. The observers marked and rated perceived lesions under the free-response paradigm. Four analyses were conducted using jackknife alternative free-response receiver operating characteristic (JAFROC) analysis: (1) 20-pixel acceptance radius (AR) with all 21 readers, abbreviated to 20/ALL analysis, (2) 40-pixel AR with 21 readers (40/ALL), (3) 20-pixel AR with 14 readers experienced in CT (20/EXP), and (4) 20-pixel AR with 7 readers with no CT experience (20/NOT). The significance level of the test was set so as to conservatively control the overall probability of a type I error to <0.05.

RESULTS

The mean JAFROC figure of merit (FOM) for the five CT acquisitions for the 20/ALL study were 0.602, 0.639, 0.372, 0.475, and 0.719 with a significant difference in lesion-detection performance evident between all individual treatment pairs (P < .0001) with the exception of the 1-2 pairing, which was not significant (these differed only in milliamp seconds). System 5, which had the highest performance, had the smallest slice thickness and the largest matrix size. For the other analyses, the system orderings remained unchanged, and the significance of FOM difference findings remained identical to those for 20/ALL, with one exception: for 20/EXP analysis the 1-2 difference became significant with the higher milliamp seconds superior. Improved detection performance was associated with a smaller slice thickness, increased matrix size, and, to a lesser extent, increased tube charge.

CONCLUSIONS

Protocol variations for CT-based attenuation correction (AC) in SPECT/CT imaging have a measurable impact on lesion-detection performance. The results imply that z-axis resolution and matrix size had the greatest impact on lesion detection, with a weaker but detectable dependence on the product of milliamp and seconds.