2
|
Kim SM, Alessio AM, De Man B, Kinahan PE. Direct Reconstruction of CT-based Attenuation Correction Images for PET with Cluster-Based Penalties. IEEE TRANSACTIONS ON NUCLEAR SCIENCE 2017; 64:959-968. [PMID: 30337765 PMCID: PMC6191195 DOI: 10.1109/tns.2017.2654680] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Extremely low-dose CT acquisitions used for PET attenuation correction have high levels of noise and potential bias artifacts due to photon starvation. This work explores the use of a priori knowledge for iterative image reconstruction of the CT-based attenuation map. We investigate a maximum a posteriori framework with cluster-based multinomial penalty for direct iterative coordinate decent (dICD) reconstruction of the PET attenuation map. The objective function for direct iterative attenuation map reconstruction used a Poisson log-likelihood data fit term and evaluated two image penalty terms of spatial and mixture distributions. The spatial regularization is based on a quadratic penalty. For the mixture penalty, we assumed that the attenuation map may consist of four material clusters: air+background, lung, soft tissue, and bone. Using simulated noisy sinogram data, dICD reconstruction was performed with different strengths of the spatial and mixture penalties. The combined spatial and mixture penalties reduced the RMSE by roughly 2 times compared to a weighted least square and filtered backprojection reconstruction of CT images. The combined spatial and mixture penalties resulted in only slightly lower RMSE compared to a spatial quadratic penalty alone. For direct PET attenuation map reconstruction from ultra-low dose CT acquisitions, the combination of spatial and mixture penalties offers regularization of both variance and bias and is a potential method to reconstruct attenuation maps with negligible patient dose. The presented results, using a best-case histogram suggest that the mixture penalty does not offer a substantive benefit over conventional quadratic regularization and diminishes enthusiasm for exploring future application of the mixture penalty.
Collapse
Affiliation(s)
- Soo Mee Kim
- Department of Radiology, University of Washington, Seattle, WA 98185, USA, telephone: +1-206-543-0236
| | - Adam M Alessio
- Department of Radiology, University of Washington, Seattle, WA 98185, USA, telephone: +1-206-543-0236
| | - Bruno De Man
- Image Reconstruction Laboratory, General Electric Global Research Center, Niskayuna, NY 12309, USA
| | - Paul E Kinahan
- Department of Radiology, University of Washington, Seattle, WA 98185, USA, telephone: +1-206-543-0236
| |
Collapse
|
3
|
Petibon Y, Guehl NJ, Reese TG, Ebrahimi B, Normandin MD, Shoup TM, Alpert NM, El Fakhri G, Ouyang J. Impact of motion and partial volume effects correction on PET myocardial perfusion imaging using simultaneous PET-MR. Phys Med Biol 2017; 62:326-343. [PMID: 27997375 PMCID: PMC5241952 DOI: 10.1088/1361-6560/aa5087] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
PET is an established modality for myocardial perfusion imaging (MPI) which enables quantification of absolute myocardial blood flow (MBF) using dynamic imaging and kinetic modeling. However, heart motion and partial volume effects (PVE) significantly limit the spatial resolution and quantitative accuracy of PET MPI. Simultaneous PET-MR offers a solution to the motion problem in PET by enabling MR-based motion correction of PET data. The aim of this study was to develop a motion and PVE correction methodology for PET MPI using simultaneous PET-MR, and to assess its impact on both static and dynamic PET MPI using 18F-Flurpiridaz, a novel 18F-labeled perfusion tracer. Two dynamic 18F-Flurpiridaz MPI scans were performed on healthy pigs using a PET-MR scanner. Cardiac motion was tracked using a dedicated tagged-MRI (tMR) sequence. Motion fields were estimated using non-rigid registration of tMR images and used to calculate motion-dependent attenuation maps. Motion correction of PET data was achieved by incorporating tMR-based motion fields and motion-dependent attenuation coefficients into image reconstruction. Dynamic and static PET datasets were created for each scan. Each dataset was reconstructed as (i) Ungated, (ii) Gated (end-diastolic phase), and (iii) Motion-Corrected (MoCo), each without and with point spread function (PSF) modeling for PVE correction. Myocardium-to-blood concentration ratios (MBR) and apparent wall thickness were calculated to assess image quality for static MPI. For dynamic MPI, segment- and voxel-wise MBF values were estimated by non-linear fitting of a 2-tissue compartment model to tissue time-activity-curves. MoCo and Gating respectively decreased mean apparent wall thickness by 15.1% and 14.4% and increased MBR by 20.3% and 13.6% compared to Ungated images (P < 0.01). Combined motion and PSF correction (MoCo-PSF) yielded 30.9% (15.7%) lower wall thickness and 82.2% (20.5%) higher MBR compared to Ungated data reconstructed without (with) PSF modeling (P < 0.01). For dynamic PET, mean MBF across all segments were comparable for MoCo (0.72 ± 0.21 ml/min/ml) and Gating (0.69 ± 0.18 ml/min/ml). Ungated data yielded significantly lower mean MBF (0.59 ± 0.16 ml/min/ml). Mean MBF for MoCo-PSF was 0.80 ± 0.22 ml/min/ml, which was 37.9% (25.0%) higher than that obtained from Ungated data without (with) PSF correction (P < 0.01). The developed methodology holds promise to improve the image quality and sensitivity of PET MPI studies performed using PET-MR.
Collapse
Affiliation(s)
- Yoann Petibon
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, MA 02114
- Department of Radiology, Harvard Medical School, Boston, MA 02115
| | - Nicolas J. Guehl
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, MA 02114
| | - Timothy G. Reese
- Department of Radiology, Harvard Medical School, Boston, MA 02115
- Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129
| | - Behzad Ebrahimi
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, MA 02114
- Department of Radiology, Harvard Medical School, Boston, MA 02115
| | - Marc D. Normandin
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, MA 02114
- Department of Radiology, Harvard Medical School, Boston, MA 02115
| | - Timothy M. Shoup
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, MA 02114
- Department of Radiology, Harvard Medical School, Boston, MA 02115
| | - Nathaniel M. Alpert
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, MA 02114
- Department of Radiology, Harvard Medical School, Boston, MA 02115
| | - Georges El Fakhri
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, MA 02114
- Department of Radiology, Harvard Medical School, Boston, MA 02115
| | - Jinsong Ouyang
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, MA 02114
- Department of Radiology, Harvard Medical School, Boston, MA 02115
| |
Collapse
|
4
|
Mehranian A, Arabi H, Zaidi H. Vision 20/20: Magnetic resonance imaging-guided attenuation correction in PET/MRI: Challenges, solutions, and opportunities. Med Phys 2016; 43:1130-55. [PMID: 26936700 DOI: 10.1118/1.4941014] [Citation(s) in RCA: 91] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
Attenuation correction is an essential component of the long chain of data correction techniques required to achieve the full potential of quantitative positron emission tomography (PET) imaging. The development of combined PET/magnetic resonance imaging (MRI) systems mandated the widespread interest in developing novel strategies for deriving accurate attenuation maps with the aim to improve the quantitative accuracy of these emerging hybrid imaging systems. The attenuation map in PET/MRI should ideally be derived from anatomical MR images; however, MRI intensities reflect proton density and relaxation time properties of biological tissues rather than their electron density and photon attenuation properties. Therefore, in contrast to PET/computed tomography, there is a lack of standardized global mapping between the intensities of MRI signal and linear attenuation coefficients at 511 keV. Moreover, in standard MRI sequences, bones and lung tissues do not produce measurable signals owing to their low proton density and short transverse relaxation times. MR images are also inevitably subject to artifacts that degrade their quality, thus compromising their applicability for the task of attenuation correction in PET/MRI. MRI-guided attenuation correction strategies can be classified in three broad categories: (i) segmentation-based approaches, (ii) atlas-registration and machine learning methods, and (iii) emission/transmission-based approaches. This paper summarizes past and current state-of-the-art developments and latest advances in PET/MRI attenuation correction. The advantages and drawbacks of each approach for addressing the challenges of MR-based attenuation correction are comprehensively described. The opportunities brought by both MRI and PET imaging modalities for deriving accurate attenuation maps and improving PET quantification will be elaborated. Future prospects and potential clinical applications of these techniques and their integration in commercial systems will also be discussed.
Collapse
Affiliation(s)
- Abolfazl Mehranian
- Division of Nuclear Medicine and Molecular Imaging, Geneva University Hospital, Geneva CH-1211, Switzerland
| | - Hossein Arabi
- Division of Nuclear Medicine and Molecular Imaging, Geneva University Hospital, Geneva CH-1211, Switzerland
| | - Habib Zaidi
- Division of Nuclear Medicine and Molecular Imaging, Geneva University Hospital, Geneva CH-1211, Switzerland; Geneva Neuroscience Centre, University of Geneva, Geneva CH-1205, Switzerland; and Department of Nuclear Medicine and Molecular Imaging, University of Groningen, Groningen 9700 RB, Netherlands
| |
Collapse
|
5
|
Mehranian A, Zaidi H. Joint Estimation of Activity and Attenuation in Whole-Body TOF PET/MRI Using Constrained Gaussian Mixture Models. IEEE TRANSACTIONS ON MEDICAL IMAGING 2015; 34:1808-1821. [PMID: 25769148 DOI: 10.1109/tmi.2015.2409157] [Citation(s) in RCA: 67] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
It has recently been shown that the attenuation map can be estimated from time-of-flight (TOF) PET emission data using joint maximum likelihood reconstruction of attenuation and activity (MLAA). In this work, we propose a novel MRI-guided MLAA algorithm for emission-based attenuation correction in whole-body PET/MR imaging. The algorithm imposes MR spatial and CT statistical constraints on the MLAA estimation of attenuation maps using a constrained Gaussian mixture model (GMM) and a Markov random field smoothness prior. Dixon water and fat MR images were segmented into outside air, lung, fat and soft-tissue classes and an MR low-intensity (unknown) class corresponding to air cavities, cortical bone and susceptibility artifacts. The attenuation coefficients over the unknown class were estimated using a mixture of four Gaussians, and those over the known tissue classes using unimodal Gaussians, parameterized over a patient population. To eliminate misclassification of spongy bones with surrounding tissues, and thus include them in the unknown class, we heuristically suppressed fat in water images and also used a co-registered bone probability map. The proposed MLAA-GMM algorithm was compared with the MLAA algorithms proposed by Rezaei and Salomon using simulation and clinical studies with two different tracer distributions. The results showed that our proposed algorithm outperforms its counterparts in suppressing the cross-talk and scaling problems of activity and attenuation and thus produces PET images of improved quantitative accuracy. It can be concluded that the proposed algorithm effectively exploits the MR information and can pave the way toward accurate emission-based attenuation correction in TOF PET/MRI.
Collapse
|
6
|
Sun T, Wu TH, Wang SJ, Yang BH, Wu NY, Mok GSP. Low dose interpolated average CT for thoracic PET/CT attenuation correction using an active breathing controller. Med Phys 2013; 40:102507. [DOI: 10.1118/1.4820976] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|