Attenuation map estimation with SPECT emission data only

Deep-learning-based methods of attenuation correction for SPECT and PET

Attenuation correction (AC) is essential for quantitative analysis and clinical diagnosis of single-photon emission computed tomography (SPECT) and positron emission tomography (PET). In clinical practice, computed tomography (CT) is utilized to generate attenuation maps (μ-maps) for AC of hybrid SPECT/CT and PET/CT scanners. However, CT-based AC methods frequently produce artifacts due to CT artifacts and misregistration of SPECT-CT and PET-CT scans. Segmentation-based AC methods using magnetic resonance imaging (MRI) for PET/MRI scanners are inaccurate and complicated since MRI does not contain direct information of photon attenuation. Computational AC methods for SPECT and PET estimate attenuation coefficients directly from raw emission data, but suffer from low accuracy, cross-talk artifacts, high computational complexity, and high noise level. The recently evolving deep-learning-based methods have shown promising results in AC of SPECT and PET, which can be generally divided into two categories: indirect and direct strategies. Indirect AC strategies apply neural networks to transform emission, transmission, or MR images into synthetic μ-maps or CT images which are then incorporated into AC reconstruction. Direct AC strategies skip the intermediate steps of generating μ-maps or CT images and predict AC SPECT or PET images from non-attenuation-correction (NAC) SPECT or PET images directly. These deep-learning-based AC methods show comparable and even superior performance to non-deep-learning methods. In this article, we first discussed the principles and limitations of non-deep-learning AC methods, and then reviewed the status and prospects of deep-learning-based methods for AC of SPECT and PET.

Joint Estimation of Metal Density and Attenuation Maps with Pencil Beam XFET

X-ray fluorescence emission tomography (XFET) is an emerging imaging modality that images the spatial distribution of metal without requiring biochemical modification or radioactivity. This work investigates the joint estimation of metal and attenuation maps with a pencil-beam XFET system that allows for direct metal measurement in the absence of attenuation. Using singular value decomposition on a simplified imaging model, we show that reconstructing metal and attenuation voxels far from the detector is an ill-conditioned problem. Using simulated data, we develop and compare two image reconstruction methods for joint estimation. The first method alternates between updating the attenuation map with a separable paraboloidal surrogates algorithm and updating the metal map with a closed-form solution. The second method performs simultaneous joint estimation with conjugate gradients based on a linearized imaging model. The alternating approach outperforms the linearized approach for iron and gold numerical phantom reconstructions. Reconstructing an (8 cm)3 object containing gold concentrations of 5 mg/cm3 and an unknown beam attenuation map using the alternating approach yields an accurate gold map (NRMSE = 0.19) and attenuation map (NRMSE = 0.14). This simulation demonstrates an accurate joint reconstruction of metal and attenuation maps, from emission data, without previous knowledge of any attenuation map.

Attenuation Correction for Dedicated Cardiac SPECT Imaging Without Using Transmission Data

Objectives

Attenuation correction (AC) using transmission scanning-like computed tomography (CT) is the standard method to increase the accuracy of cardiac single-photon emission computed tomography (SPECT) images. Recently developed dedicated cardiac SPECT do not support CT, and thus, scans on these systems are vulnerable to attenuation artifacts. This study presented a new method for generating an attenuation map directly from emission data by segmentation of precisely non-rigid registration extended cardiac-torso (XCAT)-digital phantom with cardiac SPECT images.

Methods

In-house developed non-rigid registration algorithm automatically aligns the XCAT- phantom with cardiac SPECT image to precisely segment the contour of organs. Pre-defined attenuation coefficients for given photon energies were assigned to generate attenuation maps. The CT-based attenuation maps were used for validation with which cardiac SPECT/CT data of 38 patients were included. Segmental myocardial counts of a 17-segment model from these databases were compared based on the basis of the paired t-test.

Results

The mean, and standard deviation of the mean square error and structural similarity index measure of the female stress phase between the proposed attenuation maps and the CT attenuation maps were 6.99±1.23% and 92±2.0%, of the male stress were 6.87±3.8% and 96±1.0%. Proposed attenuation correction and computed tomography based attenuation correction average myocardial perfusion count was significantly higher than that in non-AC in the mid-inferior, mid-lateral, basal-inferior, and lateral regions (p<0.001).

Conclusion

The proposed attenuation maps showed good agreement with the CT-based attenuation map. Therefore, it is feasible to enable AC for a dedicated cardiac SPECT or SPECT standalone scanners.

Convolutional neural network based attenuation correction for123I-FP-CIT SPECT with focused striatum imaging

SPECT imaging with¹²³I-FP-CIT is used for diagnosis of neurodegenerative disorders like Parkinson's disease. Attenuation correction (AC) can be useful for quantitative analysis of¹²³I-FP-CIT SPECT. Ideally, AC would be performed based on attenuation maps (μ-maps) derived from perfectly registered CT scans. Suchμ-maps, however, are most times not available and possible errors in image registration can induce quantitative inaccuracies in AC corrected SPECT images. Earlier, we showed that a convolutional neural network (CNN) based approach allows to estimate SPECT-alignedμ-maps for full brain perfusion imaging using only emission data. Here we investigate the feasibility of similar CNN methods for axially focused¹²³I-FP-CIT scans. We tested our approach on a high-resolution multi-pinhole prototype clinical SPECT system in a Monte Carlo simulation study. Three CNNs that estimateμ-maps in a voxel-wise, patch-wise and image-wise manner were investigated. As the added value of AC on clinical¹²³I-FP-CIT scans is still debatable, the impact of AC was also reported to check in which cases CNN based AC could be beneficial. AC using the ground truthμ-maps (GT-AC) and CNN estimatedμ-maps (CNN-AC) were compared with the case when no AC was done (No-AC). Results show that the effect of using GT-AC versus CNN-AC or No-AC on striatal shape and symmetry is minimal. Specific binding ratios (SBRs) from localized regions show a deviation from GT-AC≤2.5% for all three CNN-ACs while No-AC systematically underestimates SBRs by 13.1%. A strong correlation (r≥0.99) was obtained between GT-AC based SBRs and SBRs from CNN-ACs and No-AC. Absolute quantification (in kBq ml^-1) shows a deviation from GT-AC within 2.2% for all three CNN-ACs and of 71.7% for No-AC. To conclude, all three CNNs show comparable performance in accurateμ-map estimation and¹²³I-FP-CIT quantification. CNN-estimatedμ-map can be a promising substitute for CT-basedμ-map.

A physics and learning-based transmission-less attenuation compensation method for SPECT

Attenuation compensation (AC) is a pre-requisite for reliable quantification and beneficial for visual interpretation tasks in single-photon emission computed tomography (SPECT). Typical AC methods require the availability of an attenuation map, which is obtained using a transmission scan, such as a CT scan. This has several disadvantages such as increased radiation dose, higher costs, and possible misalignment between SPECT and CT scans. Also, often a CT scan is unavailable. In this context, we and others are showing that scattered photons in SPECT contain information to estimate the attenuation distribution. To exploit this observation, we propose a physics and learning-based method that uses the SPECT emission data in the photopeak and scatter windows to perform transmission-less AC in SPECT. The proposed method uses data acquired in the scatter window to reconstruct an initial estimate of the attenuation map using a physics-based approach. A convolutional neural network is then trained to segment this initial estimate into different regions. Pre-defined attenuation coefficients are assigned to these regions, yielding the reconstructed attenuation map, which is then used to reconstruct the activity distribution using an ordered subsets expectation maximization (OSEM)-based reconstruction approach. We objectively evaluated the performance of this method using highly realistic simulation studies conducted on the clinically relevant task of detecting perfusion defects in myocardial perfusion SPECT. Our results showed no statistically significant differences between the performance achieved using the proposed method and that with the true attenuation maps. Visually, the images reconstructed using the proposed method looked similar to those with the true attenuation map. Overall, these results provide evidence of the capability of the proposed method to perform transmission-less AC and motivate further evaluation.

Fisher information analysis of list-mode SPECT emission data for joint estimation of activity and attenuation distribution

The potential to perform attenuation and scatter compensation (ASC) in single-photon emission computed tomography (SPECT) imaging without a separate transmission scan is highly significant. In this context, attenuation in SPECT is primarily due to Compton scattering, where the probability of Compton scatter is proportional to the attenuation coefficient of the tissue and the energy of the scattered photon and the scattering angle are related. Based on this premise, we investigated whether the SPECT scattered-photon data acquired in list-mode (LM) format and including the energy information can be used to estimate the attenuation map. For this purpose, we propose a Fisher-information-based method that yields the Cramer-Rao bound (CRB) for the task of jointly estimating the activity and attenuation distribution using only the SPECT emission data. In the process, a path-based formalism to process the LM SPECT emission data, including the scattered-photon data, is proposed. The Fisher information method was implemented on NVIDIA graphics processing units (GPU) for acceleration. The method was applied to analyze the information content of SPECT LM emission data, which contains up to first-order scattered events, in a simulated SPECT system with parameters modeling a clinical system using realistic computational studies with 2-D digital synthetic and anthropomorphic phantoms. The method was also applied to LM data containing up to second-order scatter for a synthetic phantom. Experiments with anthropomorphic phantoms simulated myocardial perfusion and dopamine transporter (DaT)-Scan SPECT studies. The results show that the CRB obtained for the attenuation and activity coefficients was typically much lower than the true value of these coefficients. An increase in the number of detected photons yielded lower CRB for both the attenuation and activity coefficients. Further, we observed that systems with better energy resolution yielded a lower CRB for the attenuation coefficient. Overall, the results provide evidence that LM SPECT emission data, including the scattered photons, contains information to jointly estimate the activity and attenuation coefficients.

Deep learning-based attenuation map generation for myocardial perfusion SPECT

PURPOSE

Attenuation correction using CT transmission scanning increases the accuracy of single-photon emission computed tomography (SPECT) and enables quantitative analysis. Current existing SPECT-only systems normally do not support transmission scanning and therefore scans on these systems are susceptible to attenuation artifacts. Moreover, the use of CT scans also increases radiation dose to patients and significant artifacts can occur due to the misregistration between the SPECT and CT scans as a result of patient motion. The purpose of this study is to develop an approach to estimate attenuation maps directly from SPECT emission data using deep learning methods.

METHODS

Both photopeak window and scatter window SPECT images were used as inputs to better utilize the underlying attenuation information embedded in the emission data. The CT-based attenuation maps were used as labels with which cardiac SPECT/CT images of 65 patients were included for training and testing. We implemented and evaluated deep fully convolutional neural networks using both standard training and training using an adversarial strategy.

RESULTS

The synthetic attenuation maps were qualitatively and quantitatively consistent with the CT-based attenuation map. The globally normalized mean absolute error (NMAE) between the synthetic and CT-based attenuation maps were 3.60% ± 0.85% among the 25 testing subjects. The SPECT reconstructed images corrected using the CT-based attenuation map and synthetic attenuation map are highly consistent. The NMAE between the reconstructed SPECT images that were corrected using the synthetic and CT-based attenuation maps was 0.26% ± 0.15%, whereas the localized absolute percentage error was 1.33% ± 3.80% in the left ventricle (LV) myocardium and 1.07% ± 2.58% in the LV blood pool.

CONCLUSION

We developed a deep convolutional neural network to estimate attenuation maps for SPECT directly from the emission data. The proposed method is capable of generating highly reliable attenuation maps to facilitate attenuation correction for SPECT-only scanners for myocardial perfusion imaging.

Attenuation correction in emission tomography using the emission data--A review

The problem of attenuation correction (AC) for quantitative positron emission tomography (PET) had been considered solved to a large extent after the commercial availability of devices combining PET with computed tomography (CT) in 2001; single photon emission computed tomography (SPECT) has seen a similar development. However, stimulated in particular by technical advances toward clinical systems combining PET and magnetic resonance imaging (MRI), research interest in alternative approaches for PET AC has grown substantially in the last years. In this comprehensive literature review, the authors first present theoretical results with relevance to simultaneous reconstruction of attenuation and activity. The authors then look back at the early history of this research area especially in PET; since this history is closely interwoven with that of similar approaches in SPECT, these will also be covered. We then review algorithmic advances in PET, including analytic and iterative algorithms. The analytic approaches are either based on the Helgason-Ludwig data consistency conditions of the Radon transform, or generalizations of John's partial differential equation; with respect to iterative methods, we discuss maximum likelihood reconstruction of attenuation and activity (MLAA), the maximum likelihood attenuation correction factors (MLACF) algorithm, and their offspring. The description of methods is followed by a structured account of applications for simultaneous reconstruction techniques: this discussion covers organ-specific applications, applications specific to PET/MRI, applications using supplemental transmission information, and motion-aware applications. After briefly summarizing SPECT applications, we consider recent developments using emission data other than unscattered photons. In summary, developments using time-of-flight (TOF) PET emission data for AC have shown promising advances and open a wide range of applications. These techniques may both remedy deficiencies of purely MRI-based AC approaches in PET/MRI and improve standalone PET imaging.

Joint reconstruction of activity and attenuation map using LM SPECT emission data

Attenuation and scatter correction in single photon emission computed tomography (SPECT) imaging often requires a computed tomography (CT) scan to compute the attenuation map of the patient. This results in increased radiation dose for the patient, and also has other disadvantages such as increased costs and hardware complexity. Attenuation in SPECT is a direct consequence of Compton scattering, and therefore, if the scattered photon data can give information about the attenuation map, then the CT scan may not be required. In this paper, we investigate the possibility of joint reconstruction of the activity and attenuation map using list-mode (LM) SPECT emission data, including the scattered-photon data. We propose a path-based formalism to process scattered-photon data. Following this, we derive analytic expressions to compute the Cramér-Rao bound (CRB) of the activity and attenuation map estimates, using which, we can explore the fundamental limit of information-retrieval capacity from LM SPECT emission data. We then suggest a maximum-likelihood (ML) scheme that uses the LM emission data to jointly reconstruct the activity and attenuation map. We also propose an expectation-maximization (EM) algorithm to compute the ML solution.