Implementation of an iterative scatter correction, the influence of attenuation map quality and their effect on absolute quantitation in SPECT

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.

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.

Classification of the Multiple Stages of Parkinson's Disease by a Deep Convolution Neural Network Based on 99mTc-TRODAT-1 SPECT Images

Single photon emission computed tomography (SPECT) has been employed to detect Parkinson's disease (PD). However, analysis of the SPECT PD images was mostly based on the region of interest (ROI) approach. Due to limited size of the ROI, especially in the multi-stage classification of PD, this study utilizes deep learning methods to establish a multiple stages classification model of PD. In the retrospective study, the 99mTc-TRODAT-1 was used for brain SPECT imaging. A total of 202 cases were collected, and five slices were selected for analysis from each subject. The total number of images was thus 1010. According to the Hoehn and Yahr Scale standards, all the cases were divided into healthy, early, middle, late four stages, and HYS I~V six stages. Deep learning is compared with five convolutional neural networks (CNNs). The input images included grayscale and pseudo color of two types. The training and validation sets were 70% and 30%. The accuracy, recall, precision, F-score, and Kappa values were used to evaluate the models' performance. The best accuracy of the models based on grayscale and color images in four and six stages were 0.83 (AlexNet), 0.85 (VGG), 0.78 (DenseNet) and 0.78 (DenseNet).

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.

Towards standardization of absolute SPECT/CT quantification: a multi-center and multi-vendor phantom study

Abstract

Absolute quantification of radiotracer distribution using SPECT/CT imaging is of great importance for dosimetry aimed at personalized radionuclide precision treatment. However, its accuracy depends on many factors. Using phantom measurements, this multi-vendor and multi-center study evaluates the quantitative accuracy and inter-system variability of various SPECT/CT systems as well as the effect of patient size, processing software and reconstruction algorithms on recovery coefficients (RC).

Methods

Five SPECT/CT systems were included: Discovery™ NM/CT 670 Pro (GE Healthcare), Precedence™ 6 (Philips Healthcare), Symbia Intevo™, and Symbia™ T16 (twice) (Siemens Healthineers). Three phantoms were used based on the NEMA IEC body phantom without lung insert simulating body mass indexes (BMI) of 25, 28, and 47 kg/m². Six spheres (0.5–26.5 mL) and background were filled with 0.1 and 0.01 MBq/mL ^99mTc-pertechnetate, respectively. Volumes of interest (VOI) of spheres were obtained by a region growing technique using a 50% threshold of the maximum voxel value corrected for background activity. RC, defined as imaged activity concentration divided by actual activity concentration, were determined for maximum (RC_max) and mean voxel value (RC_mean) in the VOI for each sphere diameter. Inter-system variability was expressed as median absolute deviation (MAD) of RC. Acquisition settings were standardized. Images were reconstructed using vendor-specific 3D iterative reconstruction algorithms with institute-specific settings used in clinical practice and processed using a standardized, in-house developed processing tool based on the SimpleITK framework. Additionally, all data were reconstructed with a vendor-neutral reconstruction algorithm (Hybrid Recon™; Hermes Medical Solutions).

Results

RC decreased with decreasing sphere diameter for each system. Inter-system variability (MAD) was 16 and 17% for RC_mean and RC_max, respectively. Standardized reconstruction decreased this variability to 4 and 5%. High BMI hampers quantification of small lesions (< 10 ml).

Conclusion

Absolute SPECT quantification in a multi-center and multi-vendor setting is feasible, especially when reconstruction protocols are standardized, paving the way for a standard for absolute quantitative SPECT.

Absolute radiotracer concentration measurement using whole-body solid-state SPECT/CT technology: in vivo/in vitro validation

The accuracy of recently approved quantitative clinical software was determined by comparing in vivo/in vitro measurements for a solid-state cadmium-zinc-telluride SPECT/CT (single photon emission computed tomography/x-ray computed tomography) camera. Bone SPECT/CT, including the pelvic region in the field of view, was performed on 16 patients using technetium-99m methylene diphosphonic acid as a radiotracer. After imaging, urine samples from each patient provided for the measurement of in vitro radiopharmaceutical concentrations. From the SPECT/CT images, three users measured in vivo radiotracer concentration and standardized uptake value (SUV) for the bladder using quantitative software (Q.Metrix, GE Healthcare). Linear regression was used to validate any in vivo/in vitro identity relations (ideally slope = 1, intercept = 0), within a 95% confidence interval (CI). Thirteen in vivo/in vitro pairs were available for further analysis, after rejecting two as clinically irrelevant (SUVs > 100 g/mL) and one as an outlier (via Cook's distance calculations). All linear regressions (R² ≥ 0.85, P < 0.0001) provided identity in vivo/in vitro relations (95% CI), with SUV averages from all users giving a slope of 0.99 ± 0.25 and intercept of 0.14 ± 5.15 g/mL. The average in vivo/in vitro residual difference was < 20%. Solid-state SPECT/CT imaging can reliably provide in vivo urinary bladder radiotracer concentrations within approximately 20% accuracy. This practical, non-invasive, in vivo quantitation method can potentially improve diagnosis and assessment of response to treatment. Graphical abstract.

Clinical Applications of Technetium-99m Quantitative Single-Photon Emission Computed Tomography/Computed Tomography

Single-photon emission computed tomography/computed tomography (SPECT/CT) is an already established nuclear imaging modality. Co-registration of functional information (SPECT) with anatomical images (CT) paved the way to the wider application of SPECT. Recent advancements in quantitative SPECT/CT have made it possible to incorporate quantitative parameters, such as standardized uptake value (SUV) or %injected dose (%ID), in gamma camera imaging. This is indeed a paradigm shift in gamma camera imaging from qualitative to quantitative evaluation. In fact, such quantitative approaches of nuclear imaging have only been accomplished for positron emission tomography (PET) technology. Attenuation correction, scatter correction, and resolution recovery are the three main features that enabled quantitative SPECT/CT. Further technical improvements are being achieved for partial-volume correction, motion correction, and dead-time correction. The reported clinical applications for quantitative SPECT/CT are mainly related to Tc-99m-labeled radiopharmaceuticals: Tc-99m diphosphonate for bone/joint diseases, Tc-99m pertechnetate for thyroid function, and Tc-99m diethylenetriaminepentaacetic acid for measurement of glomerular filtration rate. Dosimetry before trans-arterial radio-embolization is also a promising application for Tc-99m macro-aggregated albumin. In this review, clinical applications of Tc-99m quantitative SPECT/CT will be discussed.

Accuracy of 177Lu activity quantification in SPECT imaging: a phantom study

BACKGROUND

The aim of the study is to assess accuracy of activity quantification of 177Lu studies performed according to recommendations provided by the committee on Medical Internal Radiation Dose (MIRD) pamphlets 23 and 26. The performances of two scatter correction and three segmentation methods were compared. Additionally, the accuracy of tomographic and planar methods for determination of the camera normalization factor (CNF) was evaluated. Eight phantoms containing inserts of different sizes and shapes placed in air, water, and radioactive background were scanned using a Siemens SymbiaT SPECT/CT camera. Planar and tomographic scans with 177Lu sources were used to measure CNF. Images were reconstructed with our SPEQToR software using resolution recovery, attenuation, and two scatter correction methods (analytical photon distribution interpolated (APDI) and triple energy window (TEW)). Segmentation was performed using a fixed threshold method for both air and cold water scans. For hot water experiments three segmentation methods were compared as folows: a 40% fixed threshold, segmentation based on CT images, and our iterative adaptive dual thresholding (IADT). Quantification error, defined as the percent difference between experimental and true activities, was evaluated.

RESULTS

Quantification error for scans in air was better for TEW scatter correction (<6%) than for APDI (<11%). This trend was reversed for scans in water (<10% for APDI and <14% for TEW). For hot water, the best results (<18% for small objects and <5% for objects >100 ml) were obtained when APDI and IADT were used for scatter correction and segmentation, respectively. Additionally, we showed that planar acquisitions with scatter correction and tomographic scans provide similar CNF values. This is an important finding because planar acquisitions are easier to perform than tomographic scans. TEW and APDI resulted in similar quantification errors with APDI showing a small advantage for objects placed in medium with non-uniform density.

CONCLUSIONS

Following the MIRD recommendations for data acquisition and reconstruction resulted in accurate activity quantification (errors <5% for large objects). However, techniques for better organ/tumor segmentation must still be developed.

Microdosing, isotopic labeling, radiotracers and metabolomics: relevance in drug discovery, development and safety

This review discusses the use of stable (¹³C, ²D) or radioactive isotopes (¹⁴C, ¹¹C, ¹⁸F, ¹³¹I, ⁶⁴Cu, ⁶⁸Ga) incorporated into the molecular structure of new drug entities for the purpose of pharmacokinetic or -dynamic studies. Metabolite in safety testing requires the administration of pharmacologically active doses. In such studies, radiotracers find application mainly in preclinical animal investigations, whereby LC-MS/MS is used to identify metabolite structure and drug-related effects. In contrast, first-in-human metabolite studies have to be carried out at nonpharmacological doses not exceeding 100 μg (microdose), which is generally too low for metabolite detection by LC-MS/MS. This short-coming can be overcome by specific radio- or isotopic labeling of the drug of interest and measurements using accelerator mass spectroscopy, single-photon emission computed tomography and positron emission tomography. Such combined radioisotope-based approaches permit Phase 0, first-in-human metabolite study.

Analytically based photon scatter modeling for a multipinhole cardiac SPECT camera

PURPOSE

Dedicated cardiac SPECT scanners have improved performance over standard gamma cameras allowing reductions in acquisition times and/or injected activity. One approach to improving performance has been to use pinhole collimators, but this can cause position-dependent variations in attenuation, sensitivity, and spatial resolution. CT attenuation correction (AC) and an accurate system model can compensate for many of these effects; however, scatter correction (SC) remains an outstanding issue. In addition, in cameras using cadmium-zinc-telluride-based detectors, a large portion of unscattered photons is detected with reduced energy (low-energy tail). Consequently, application of energy-based SC approaches in these cameras leads to a higher increase in noise than with standard cameras due to the subtraction of true counts detected in the low-energy tail. Model-based approaches with parallel-hole collimator systems accurately calculate scatter based on the physics of photon interactions in the patient and camera and generate lower-noise estimates of scatter than energy-based SC. In this study, the accuracy of a model-based SC method was assessed using physical phantom studies on the GE-Discovery NM530c and its performance was compared to a dual energy window (DEW)-SC method.

METHODS

The analytical photon distribution (APD) method was used to calculate the distribution of probabilities that emitted photons will scatter in the surrounding scattering medium and be subsequently detected. APD scatter calculations for ^99mTc-SPECT (140 ± 14 keV) were validated with point-source measurements and 15 anthropomorphic cardiac-torso phantom experiments and varying levels of extra-cardiac activity causing scatter inside the heart. The activity inserted into the myocardial compartment of the phantom was first measured using a dose calibrator. CT images were acquired on an Infinia Hawkeye (GE Healthcare) SPECT/CT and coregistered with emission data for AC. For comparison, DEW scatter projections (120 ± 6 keV ) were also extracted from the acquired list-mode SPECT data. Either APD or DEW scatter projections were subtracted from corresponding 140 keV measured projections and then reconstructed with AC (APD-SC and DEW-SC). Quantitative accuracy of the activity measured in the heart for the APD-SC and DEW-SC images was assessed against dose calibrator measurements. The difference between modeled and acquired projections was measured as the root-mean-squared-error (RMSE). APD-modeled projections for a clinical cardiac study were also evaluated.

RESULTS

APD-modeled projections showed good agreement with SPECT measurements and had reduced noise compared to DEW scatter estimates. APD-SC reduced mean error in activity measurement compared to DEW-SC in images and the reduction was statistically significant where the scatter fraction (SF) was large (mean SF = 28.5%, T-test p = 0.007). APD-SC reduced measurement uncertainties as well; however, the difference was not found to be statistically significant (F-test p > 0.5). RMSE comparisons showed that elevated levels of scatter did not significantly contribute to a change in RMSE (p > 0.2).

CONCLUSIONS

Model-based APD scatter estimation is feasible for dedicated cardiac SPECT scanners with pinhole collimators. APD-SC images performed better than DEW-SC images and improved the accuracy of activity measurement in high-scatter scenarios.

JADA: a graphical user interface for comprehensive internal dose assessment in nuclear medicine

PURPOSE

The main objective of this work was to design a comprehensive dosimetry package that would keep all aspects of internal dose calculation within the framework of a single software environment and that would be applicable for a variety of dose calculation approaches.

METHODS

Our MATLAB-based graphical user interface (GUI) can be used for processing data obtained using pure planar, pure SPECT, or hybrid planar/SPECT imaging. Time-activity data for source regions are obtained using a set of tools that allow the user to reconstruct SPECT images, load images, coregister a series of planar images, and to perform two-dimensional and three-dimensional image segmentation. Curve fits are applied to the acquired time-activity data to construct time-activity curves, which are then integrated to obtain time-integrated activity coefficients. Subsequently, dose estimates are made using one of three methods.

RESULTS

The organ level dose calculation subGUI calculates mean organ doses that are equivalent to dose assessment performed by OLINDA/EXM. Voxelized dose calculation options, which include the voxel S value approach and Monte Carlo simulation using the EGSnrc user code DOSXYZnrc, are available within the process 3D image data subGUI.

CONCLUSIONS

The developed internal dosimetry software package provides an assortment of tools for every step in the dose calculation process, eliminating the need for manual data transfer between programs. This saves times and minimizes user errors, while offering a versatility that can be used to efficiently perform patient-specific internal dose calculations in a variety of clinical situations.

An evidence-based review of quantitative SPECT imaging and potential clinical applications

SPECT has traditionally been regarded as nonquantitative. Advances in multimodality γ-cameras (SPECT/CT), algorithms for image reconstruction, and sophisticated compensation techniques to correct for photon attenuation and scattering have, however, now made quantitative SPECT viable in a manner similar to quantitative PET (i.e., kBq cm(-3), standardized uptake value). This review examines the evidence for quantitative SPECT and demonstrates clinical studies in which the accuracy of the reconstructed SPECT data has been assessed in vivo. SPECT reconstructions using CT-based compensation corrections readily achieve accuracy for (99m)Tc to within ± 10% of the known concentration of the radiotracer in vivo. Quantification with other radionuclides is also being introduced. SPECT continues to suffer from poorer photon detection efficiency (sensitivity) and spatial resolution than PET; however, it has the benefit in some situations of longer radionuclide half-lives, which may better suit the biologic process under examination, as well as the ability to perform multitracer studies using pulse height spectroscopy to separate different radiolabels.

Two methods to generate templates for template-based partial volume effect correction: SPECT phantom experiments

In this paper, we explore the applicability of template-based compensation for the partial volume effect (PVE) for situations where (i) the image has multiple uptake sites (tumors and organs) but only one of them is treated as a region of interest (ROI) with the boundaries available from a high-resolution modality and (ii) no information regarding activities inside or outside this ROI is a priori available. We modeled this situation by performing SPECT acquisitions of phantoms containing 21 containers, which had different shapes and sizes and were filled with different levels of activity. In our analysis, each of these containers was treated as an individual ROI. We compared the performance of two methods of template construction. In method 1, the ROI template value was obtained from a conventionally reconstructed (without PVEC) image. In method 2, the ROI template value was directly (bypassing the PVE-affected conventional image) calculated from projections using region-based reconstruction. Our processing shows that method 1 resulted in consistent (activities for all 21 ROIs were improved) but relatively weak PVE compensation (errors of recovered total activities were equal to or lower than 10% for 5 ROIs only). Application of method 2 resulted in a selective (activities for 19 ROIs were improved) but considerably better compensation when compared to method 1 (errors of recovered total activities were equal to or lower than 10% for 10 ROIs).

Quantitative SPECT/CT

Conventional nuclear medical imaging uses radiopharmaceuticals labeled by single-photon emitters such as Tc-99m, I-123, or I-131 in vivo. Classical clinical examples are the study of bone metabolism by bone scintigraphy with the Tc-99m-labeled polyphosphonates or of iodine transport into the thyroid gland using Tc-99m-pertechnetate. With single-photon emission-computed tomography (SPECT), the distribution of these radiopharmaceuticals within the human body is three-dimensionally visualized. Contrary to positron emission tomography (PET), current SPECT technology does not allow the quantification of regional values of radioactivity tissue concentration as SPECT images are grossly compromised by artifacts caused by photon scatter and attenuation. With the advent of hybrid imaging systems combining a SPECT camera with an X-ray computerized (CT) scanner in one gantry, reliable corrections for these artifacts seem possible, allowing truly quantitative SPECT.

A multivendor phantom study comparing the image quality produced from three state-of-the-art SPECT-CT systems

OBJECTIVE

Ongoing advancements in single photon emission computed tomography with on-board X-ray computed tomography (SPECT-CT) hardware and software raise important questions regarding the relative performances of various cameras and their respective image-processing software. This phantom-based study compares images produced from three state-of-the-art cameras using four image quality measurements.

METHODS

A thorax phantom modeling the spine, lungs, a healthy heart, and three tumors (cylindrical bottles) was scanned using the following SPECT-CT systems: Philips' Precedence (PP), GE's Infinia-Hawkeye (GH), and Siemens' Symbia-T6 (SS). For each scan, Tc-99m solutions were injected into the heart, three bottles, and thorax to yield activity concentration ratios of roughly 6:1 for both heart:thorax and tumor:thorax. The data were reconstructed using the most advanced software available on the cameras, namely, Evolution for Bone and Evolution for Cardiac (EVB and EVC, respectively), Astonish (AST), and Flash3D (FLA) for GH, PP, and SS, respectively. In addition, all sets of data were reconstructed using our in-house software. The mean values of activity error, uniformity, signal-to-noise ratio, and contrast error were used as figures of merit (FOM).

RESULTS

No significant differences were observed for all FOM between all in-house reconstructions using PP, GH, and SS acquisition data. The mean activity error for the AST reconstruction (-24.0±1.6%) was significantly closer to the truth relative to EVB (-38.0±1.6%), EVC (-34.5±2.3%), and FLA (-33.8±1.6%). No significant differences were found between EVC and FLA for all FOM.

CONCLUSION

In this phantom-based study, Philips' AST provided the most quantitatively accurate and highest contrast images, whereas Siemens' FLA and GE's EVC provided relatively higher signal-to-noise ratios and more uniform images.

Quantitative image reconstruction for dual-isotope parathyroid SPECT/CT: phantom experiments and sample patient studies

We investigated the quantitative accuracy of the model-based dual-isotope single-photon emission computed tomography (DI-SPECT) reconstructions that use Klein-Nishina expressions to estimate the scattered photon contributions to the projection data. Our objective was to examine the ability of the method to recover the absolute activities pertaining to both radiotracers: Tc-99m and I-123. We validated our method through a series of phantom experiments performed using a clinical hybrid SPECT/CT camera (Infinia Hawkeye, GE Healthcare). Different activity ratios and different attenuating media were used in these experiments to create cross-talk effects of varying severity, which can occur in clinical studies. Accurate model-based corrections for scatter and cross-talk with CT attenuation maps allowed for the recovery of the absolute activities from DI-SPECT/CT scans with errors that ranged 0-10% for both radiotracers. The unfavorable activity ratios increased the computational burden but practically did not affect the resulting accuracy. The visual analysis of parathyroid patient data demonstrated that our model-based processing improved adenoma/background contrast and enhanced localization of small or faint adenomas.

Quantification and reduction of the collimator-detector response effect in SPECT by applying a system model during iterative image reconstruction: a simulation study

BACKGROUND

Detector blurring and non-ideal collimation decrease the spatial resolution of the single-photon emission computed tomography (SPECT) images. Iterative reconstruction algorithms such as ordered subsets expectation maximization (OSEM) can incorporate degrading factors during reconstruction. We investigated the quantitative errors associated with poor SPECT resolution and evaluated the importance of two-dimensional (2D) and three-dimensional (3D) resolution recovery by modelling system response during iterative image reconstruction.

METHODS

Different phantoms consisted of the NURBS-based cardiac-torso (NCAT) liver phantom with small tumors, the Zubal brain phantom and the NCAT heart phantom were used in this study. Monte Carlo simulation was used to create SPECT projections. Gaussian functions were used to model collimator detector response (CDR). Modeled CDRs were applied during OSEM. Both noise-free and noisy projections were created.

RESULTS

Even with noise-free projections, conventional OSEM algorithm provided limited quantitative accuracy compared to both 2D and 3D resolution recovery. The 3D implementation of resolution recovery, however, yielded superior results compared to its 2D implementation. For the liver phantom, the ability to distinguish small tumors in both transverse and axial planes was improved. For the brain phantom, gray to white matter activity ratio was increased from 3.14 ± 0.04 in simple OSEM to 3.84 ± 0.06 in 3D OSEM. For the NCAT heart phantom, 3D resolution recovery, results in images with thinner wall and higher contrast for different noise levels.

CONCLUSION

There are considerable quantitative errors associated with CDR, especially when the size of the target is comparable with the spatial resolution of the system. Between different reconstruction algorithms, 3D OSEM that consider the 3D nature of CDR, improve both the visual quality and the quantitative accuracy of any SPECT studies.

Toward a practical template-based approach to semiquantitative SPECT myocardial perfusion imaging

PURPOSE

Our template-based quantitative perfusion single photon emission computed tomography (SPECT) method (T-QPS) performs semiquantitative analysis for myocardial perfusion imaging (MPI) without the use of normal databases. However, in its current form, T-QPS requires extensive calculations, which limits its clinical application. In the interest of clinical feasibility, the authors examine the trade-off between accuracy and processing time as the method is simplified.

METHODS

The T-QPS method uses the reconstructed SPECT image of the patient to create a 3D digital template of his∕her healthy heart. This template is then projected, reconstructed, and sampled into the bulls-eye map domain. A ratio of the patient and template images produces a final corrected image in which a threshold is applied to identify perfusion defects. In principle, the template should be constructed with the heart and all extracardiac activity, and the projection step should include primary and scatter components; however, this leads to lengthy calculations. In an attempt to shorten the processing time, the authors analyzed the performance of four template (T) generation methods: T(P-HRT), T(PS-HRT), T(P-HRTBKG), and T(PS-HRTBKG), where P and S represent primary and scattered photons included in the projection step, respectively; and HRT and HRTBKG represent template constructed with the heart only and the heart with background activity, respectively. Forty-eight thorax phantoms and 21 randomly selected patient studies were analyzed using each approach. All studies used GE's Infinia Hawkeye SPECT∕CT system and followed a standard cardiac acquisition protocol.

RESULTS

Approximate processing times for the T(P-HRT), T(PS-HRT), T(P-HRTBKG), and T(PS-HRTBKG) methods were less than a minute, 2-3 h, less than a minute and 3-4 h, respectively. In both the simulation and patient studies, a significant reduction in the quality of perfusion defect definition was exhibited by the T(P-HRT) method relative to the other three methods. The optimal method with respect to perfusion defect definition and processing time was T(P-HRTBKG) with a sensitivity, specificity, and accuracy in spatially defining the perfusion defects (simulation study) of 80%, 84%, and 83%, respectively.

CONCLUSIONS

The T-QPS method using T(P-HRTBKG) leads to accurate and fast semiquantitative analysis of SPECT MPI, without the use of normal databases.

Quantitative evaluation of scintillation camera imaging characteristics of isotopes used in liver radioembolization

Background

Scintillation camera imaging is used for treatment planning and post-treatment dosimetry in liver radioembolization (RE). In yttrium-90 (90Y) RE, scintigraphic images of technetium-99m (99mTc) are used for treatment planning, while 90Y Bremsstrahlung images are used for post-treatment dosimetry. In holmium-166 (166Ho) RE, scintigraphic images of 166Ho can be used for both treatment planning and post-treatment dosimetry. The aim of this study is to quantitatively evaluate and compare the imaging characteristics of these three isotopes, in order that imaging protocols can be optimized and RE studies with varying isotopes can be compared.

Methodology/Principal Findings

Phantom experiments were performed in line with NEMA guidelines to assess the spatial resolution, sensitivity, count rate linearity, and contrast recovery of 99mTc, 90Y and 166Ho. In addition, Monte Carlo simulations were performed to obtain detailed information about the history of detected photons. The results showed that the use of a broad energy window and the high-energy collimator gave optimal combination of sensitivity, spatial resolution, and primary photon fraction for 90Y Bremsstrahlung imaging, although differences with the medium-energy collimator were small. For 166Ho, the high-energy collimator also slightly outperformed the medium-energy collimator. In comparison with 99mTc, the image quality of both 90Y and 166Ho is degraded by a lower spatial resolution, a lower sensitivity, and larger scatter and collimator penetration fractions.

Conclusions/Significance

The quantitative evaluation of the scintillation camera characteristics presented in this study helps to optimize acquisition parameters and supports future analysis of clinical comparisons between RE studies.

Quantitative (177)Lu SPECT (QSPECT) imaging using a commercially available SPECT/CT system

Purpose: The combination of single photon emission computed tomography (SPECT) and computer tomography (CT) that incorporates iterative reconstruction algorithms with attenuation and scatter correction should facilitate accurate non-invasive quantitative imaging. Quantitative SPECT (QSPECT) may improve diagnostic ability and could be useful for many applications including dosimetry assessment. Using ¹⁷⁷Lu, we developed a QSPECT method using a commercially available SPECT/CT system. Methods: Serial SPECT of ¹⁷⁷Lu sources (89–12,400 MBq) were acquired with multiple contiguous energy windows along with a co-registered CT, and were reconstructed using an iterative algorithm with attenuation and scatter correction. Camera sensitivity (based on reconstructed SPECT count rate) and dead-time (based on wide-energy spectrum count rate) were resolved by non-linear curve fit. Utilizing these parameters, a SPECT dataset can be converted to a QSPECT dataset allowing quantitation in Becquerels per cubic centimetre or standardized uptake value (SUV). Validation QSPECT/CT studies were performed on a ¹⁷⁷Lu cylindrical phantom (7 studies) and on 5 patients (6 studies) who were administered a therapeutic dose of [¹⁷⁷Lu]octreotate. Results: The QSPECT sensitivity was 1.08 × 10⁻⁵± 0.02 × 10⁻⁵ s⁻¹ Bq⁻¹. The paralyzing dead-time constant was 0.78 ± 0.03 µs. The measured total activity with QSPECT deviated from the calibrated activity by 5.6 ± 1.9% and 2.6 ± 1.8%, respectively, in phantom and patients. Dead-time count loss up to 11.7% was observed in patient studies. Conclusion: QSPECT has high accuracy both in our phantom model and in clinical practice following [¹⁷⁷Lu]octreotate therapy. This has the potential to yield more accurate dosimetry estimates than planar imaging and facilitate therapeutic response assessment. Validating this method with other radionuclides could open the way for many other research and clinical applications.

Absolute quantification in SPECT

Single-photon emission computed tomography (SPECT) allows the three-dimensional visualization of radioactivity within the human body and is widely used for clinical purposes. In SPECT, image quality is compromised by several factors including photon attenuation, photon scatter, the partial volume effect, and motion artefacts. These variables also confound the capacity of SPECT to quantify the concentration of radioactivity within given volumes of interest in absolute units, e.g. as kilobecquerels per cubic centimetre. In the last decade, considerable technical progress has been achieved in SPECT image reconstruction, involving, in particular, the development of iterative image reconstruction techniques. Furthermore, hybrid cameras integrating a SPECT camera with an X-ray CT scanner have become commercially available. These systems allow the acquisition of SPECT and CT datasets registered to each other with a high anatomical accuracy. First studies have shown that iterative SPECT image reconstruction techniques incorporating information from SPECT/CT image datasets greatly increase the accuracy of SPECT in quantifying radioactivity concentrations in phantoms and also in humans. This new potential of SPECT may improve not only diagnostic accuracy, but also dosimetry for internal radiotherapy.

Quantitative accuracy of clinical 99mTc SPECT/CT using ordered-subset expectation maximization with 3-dimensional resolution recovery, attenuation, and scatter correction

We present a calibration method of a clinical SPECT/CT device for quantitative (99m)Tc SPECT. We use a commercially available reconstruction package including ordered-subset expectation maximization (OSEM) with depth-dependent 3-dimensional resolution recovery (OSEM-3D), CT-based attenuation correction, and scatter correction. We validated the method in phantom studies and applied it to images from patients injected with (99m)Tc-diphosponate.

METHODS

The following 3 steps were performed to derive absolute quantitative values from SPECT reconstructed images. In step 1, we used simulations to characterize the SPECT/CT system and derive emission recovery values for various imaging parameter settings. We simulated spheres of varying diameters and focused on the dependencies of activity estimation errors on structure size and position, pixel size, count density, and reconstruction parameters. In step 2, we cross-calibrated our clinical SPECT/CT system with the well counter using a large cylinder phantom. This step provided the mapping from image counts to kBq/mL. And in step 3, correction factors from steps 1 and 2 were applied to reconstructed images. We used a cylinder phantom with variable-sized spheres for verification of the method. For in vivo validation, SPECT/CT datasets from 16 patients undergoing (99m)Tc-diphosponate SPECT/CT examinations of the pelvis including the bladder were acquired. The radioactivity concentration in the patients' urine served as the gold standard. Mean quantitative accuracy and SEs were calculated.

RESULTS

In the phantom experiments, the mean accuracy in quantifying radioactivity concentration in absolute terms was within 3.6% (SE, 8.0%), with a 95% confidence interval between -19.4% and +12.2%. In the patient studies, the mean accuracy was within 1.1% (SE, 8.4%), with a 95% confidence interval between -15.4% and +17.5%.

CONCLUSION

Current commercially available SPECT/CT technology using OSEM-3D reconstruction, scatter correction, and CT-based attenuation correction allows quantification of (99m)Tc radioactivity concentration in absolute terms within 3.6% in phantoms and 1.1% in patients with a focus on the bladder. This opens up the opportunity of SPECT quantitation entering the routine clinical arena. Still, the imprecision caused by unavoidable measurement errors is a dominant factor for absolute quantitation in a clinical setup.

An investigation of potential sources of artifacts in SPECT-CT myocardial perfusion studies

BACKGROUND

The increased use of hybrid SPECT-CT systems in myocardial perfusion imaging calls for a careful review of protocols that are employed in data acquisition and processing. Our study investigates the cases of potential false perfusion defects that may appear in cardiac images reconstructed with CT-based attenuation correction and high-resolution acquisition matrix.

METHODS

Phantom experiments performed on Infinia-Hawkeye (GE Healthcare) modeling patients with normal hearts and different body sizes and shapes were reconstructed using standard ordered subsets expectation maximization (OSEM) method with 10 subsets and 2-50 iterations. The CT-based attenuation correction (AC) with and without distance-dependent resolution recovery (RR) were employed.

RESULTS

Images reconstructed from scans with a standard thorax phantom did not show any artifacts. However, some images reconstructed from the data with extra water bags (modeling high and non-uniformly distributed attenuation) clearly displayed false perfusion defects. A potential for creating such artifacts was especially high in complex reconstruction cases. In most cases, the severity of defects decreased when reconstructions with more iterations were used.

CONCLUSIONS

In situations with strong and non-uniform attenuation, the reconstruction methods with sophisticated data processing (large matrix size, AC + RR corrections) may require substantially more processing than is currently used and/or recommended.

In vivo tumor grading of prostate cancer using quantitative 111In-capromab pendetide SPECT/CT

We have developed an in vivo method to quantify antibody uptake using (111)In-capromab pendetide SPECT combined with CT (SPECT/CT). Our goal was to evaluate this method for potential grading of prostate tumors.

METHODS

Our phantom experiments focused on the robustness of an advanced iterative reconstruction algorithm that involves corrections for photon attenuation, scatter, and geometric blurring caused by radionuclide collimators. The conversion factors between image values and tracer concentrations (in Bq/mL) were calculated from a uniform phantom filled with an aqueous solution of (111)InCl(3) using the same acquisition protocol and reconstruction parameters as for patient studies. In addition, the spatial resolution of the reconstructed images was measured from a point source phantom. The measured spatial resolution was modeled into a point-spread function, and the point-spread function was implemented in a deconvolution-based partial-volume-error correction algorithm. The recovery capability to correctly estimate true tracer concentrations was tested using prostatelike and bladderlike lesion phantoms fitted in the modified National Electrical Manufacturers Association/International Electrotechnical Commission body phantom. Patients with biopsy-proven prostate cancer (n = 10) who underwent prostatectomy were prospectively enrolled in the preoperative SPECT/CT studies at the San Francisco Veterans Affairs Medical Center. The CT portion of SPECT/CT was used to generate CT-based attenuation maps and as an anatomic localization tool for clinical interpretation. Pathologic Gleason grades were compared with in vivo antibody uptake value (AUV) normalized by injected dose, effective half-life, and injection-scan time difference. AUVs were calculated in each lobe of the prostate gland with cylindric volumes of interest having dimensions of 1.5 cm in both diameter and height.

RESULTS

Reconstructed SPECT images further corrected by the deconvolution-based partial-volume-error correction could recover tracer concentrations up to 90% of true values in measurements of phantom volumes as small as 7.77 mL. From patient studies, there was a statistically significant correlation (rho = 0.71, P = 0.033) between higher AUVs (from either left or right lobe) and higher components of pathologic Gleason scores.

CONCLUSION

Our results strongly indicate potential for noninvasive prostate tumor grading using quantitative (111)In-capromab pendetide SPECT/CT.

Attenuation correction for lung SPECT: evidence of need and validation of an attenuation map derived from the emission data

PURPOSE

The aim of our study was to investigate the importance of attenuation correction (AC) in reconstructed and reprojected images on lung SPECT studies.

METHODS

Simulation studies were undertaken to evaluate the influence of AC on defect-to-normal ratios (D/N), to demonstrate the influence of errors in the correction map values and to detect lung boundaries used for AC. The use of a synthetic map (SM) for AC of the clinical data was also evaluated and the results compared with those obtained with data derived from CT (CTM). Additionally, the role of AC in reprojected SPECT data was assessed and level of noise on the 'planar-like' images was measured.

RESULTS

Phantom studies showed that AC markedly affects the D/N ratio. However, variations in micro values typical of those found in clinical studies resulted in relatively small changes in results. Eroded and dilated conditions did not cause any significant effect on D/N. The level of noise in the reprojected images is reduced in comparison with real planar data. Clinical SPECT/CT data reconstructed with AC using CTM and SM showed an excellent correlation between the two methods.

CONCLUSION

AC improves D/N in lung SPECT studies, thus potentially enhancing the diagnostic capability of the method. The use of a synthetic map for AC is feasible, avoiding the need for an additional procedure and the increased radiation dose involved. Planar-like images generated from reprojected SPECT data are well matched to normal planar images provided AC is performed and attenuation included in the reprojection.

Accuracy of quantitative reconstructions in SPECT/CT imaging

The goal of this study was to determine the quantitative accuracy of our OSEM-APDI reconstruction method based on SPECT/CT imaging for Tc-99m, In-111, I-123, and I-131 isotopes. Phantom studies were performed on a SPECT/low-dose multislice CT system (Infinia-Hawkeye-4 slice, GE Healthcare) using clinical acquisition protocols. Two radioactive sources were centrally and peripherally placed inside an anthropometric Thorax phantom filled with non-radioactive water. Corrections for attenuation, scatter, collimator blurring and collimator septal penetration were applied and their contribution to the overall accuracy of the reconstruction was evaluated. Reconstruction with the most comprehensive set of corrections resulted in activity estimation with error levels of 3-5% for all the isotopes.