1
|
Wang H, Huang Z, Zhang Q, Gao D, OuYang Z, Liang D, Liu X, Yang Y, Zheng H, Hu Z. Technical note: A preliminary study of dual-tracer PET image reconstruction guided by FDG and/or MR kernels. Med Phys 2021; 48:5259-5271. [PMID: 34252216 DOI: 10.1002/mp.15089] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2020] [Revised: 06/22/2021] [Accepted: 06/23/2021] [Indexed: 11/09/2022] Open
Abstract
PURPOSE Clinically, single radiotracer positron emission tomography (PET) imaging is a commonly used examination method; however, since each radioactive tracer reflects the information of only one kind of cell, it easily causes false negatives or false positives in disease diagnosis. Therefore, reasonably combining two or more radiotracers is recommended to improve the accuracy of diagnosis and the sensitivity and specificity of the disease when conditions permit. METHODS This paper proposes incorporating 18 F-fluorodeoxyglucose (FDG) as a higher-quality PET image to guide the reconstruction of other lower-count 11 C-methionine (MET) PET datasets to compensate for the lower image quality by a popular kernel algorithm. Specifically, the FDG prior is needed to extract kernel features, and these features were used to build a kernel matrix using a k-nearest-neighbor (kNN) search for MET image reconstruction. We created a 2-D brain phantom to validate the proposed method by simulating sinogram data containing Poisson random noise and quantitatively compared the performance of the proposed FDG-guided kernelized expectation maximization (KEM) method with the performance of Gaussian and non-local means (NLM) smoothed maximum likelihood expectation maximization (MLEM), MR-guided KEM, and multi-guided-S KEM algorithms. Mismatch experiments between FDG/MR and MET data were also carried out to investigate the outcomes of possible clinical situations. RESULTS In the simulation study, the proposed method outperformed the other algorithms by at least 3.11% in the signal-to-noise ratio (SNR) and 0.68% in the contrast recovery coefficient (CRC), and it reduced the mean absolute error (MAE) by 8.07%. Regarding the tumor in the reconstructed image, the proposed method contained more pathological information. Furthermore, the proposed method was still superior to the MR-guided KEM method in the mismatch experiments. CONCLUSIONS The proposed FDG-guided KEM algorithm can effectively utilize and compensate for the tissue metabolism information obtained from dual-tracer PET to maximize the advantages of PET imaging.
Collapse
Affiliation(s)
- Haiyan Wang
- Lauterbur Research Center for Biomedical Imaging, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,Shenzhen College of Advanced Technology, University of Chinese Academy of Sciences, Beijing, China
| | - Zhenxing Huang
- Lauterbur Research Center for Biomedical Imaging, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
| | - Qiyang Zhang
- Lauterbur Research Center for Biomedical Imaging, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
| | - Dongfang Gao
- Lauterbur Research Center for Biomedical Imaging, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
| | - Zhanglei OuYang
- Lauterbur Research Center for Biomedical Imaging, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
| | - Dong Liang
- Lauterbur Research Center for Biomedical Imaging, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,Chinese Academy of Sciences Key Laboratory of Health Informatics, Shenzhen, China
| | - Xin Liu
- Lauterbur Research Center for Biomedical Imaging, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,Chinese Academy of Sciences Key Laboratory of Health Informatics, Shenzhen, China
| | - Yongfeng Yang
- Lauterbur Research Center for Biomedical Imaging, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,Chinese Academy of Sciences Key Laboratory of Health Informatics, Shenzhen, China
| | - Hairong Zheng
- Lauterbur Research Center for Biomedical Imaging, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,Chinese Academy of Sciences Key Laboratory of Health Informatics, Shenzhen, China
| | - Zhanli Hu
- Lauterbur Research Center for Biomedical Imaging, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,Chinese Academy of Sciences Key Laboratory of Health Informatics, Shenzhen, China
| |
Collapse
|
2
|
Arridge SR, Ehrhardt MJ, Thielemans K. (An overview of) Synergistic reconstruction for multimodality/multichannel imaging methods. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2021; 379:20200205. [PMID: 33966461 DOI: 10.1098/rsta.2020.0205] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Imaging is omnipresent in modern society with imaging devices based on a zoo of physical principles, probing a specimen across different wavelengths, energies and time. Recent years have seen a change in the imaging landscape with more and more imaging devices combining that which previously was used separately. Motivated by these hardware developments, an ever increasing set of mathematical ideas is appearing regarding how data from different imaging modalities or channels can be synergistically combined in the image reconstruction process, exploiting structural and/or functional correlations between the multiple images. Here we review these developments, give pointers to important challenges and provide an outlook as to how the field may develop in the forthcoming years. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.
Collapse
Affiliation(s)
- Simon R Arridge
- Department of Computer Science, University College London, London, UK
| | - Matthias J Ehrhardt
- Department of Mathematical Sciences, University of Bath, Bath, UK
- Institute for Mathematical Innovation, University of Bath, Bath, UK
| | - Kris Thielemans
- Institute of Nuclear Medicine, University College London, London, UK
| |
Collapse
|
3
|
Munoz C, Ellis S, Nekolla SG, Kunze KP, Vitadello T, Neji R, Botnar RM, Schnabel JA, Reader AJ, Prieto C. MR-guided motion-corrected PET image reconstruction for cardiac PET-MR. J Nucl Med 2021; 62:jnumed.120.254235. [PMID: 34049978 PMCID: PMC8612202 DOI: 10.2967/jnumed.120.254235] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2020] [Revised: 03/09/2021] [Accepted: 03/09/2021] [Indexed: 11/16/2022] Open
Abstract
Simultaneous PET-MR imaging has shown potential for the comprehensive assessment of myocardial health from a single examination. Furthermore, MR-derived respiratory motion information has been shown to improve PET image quality by incorporating this information into the PET image reconstruction. Separately, MR-based anatomically guided PET image reconstruction has been shown to perform effective denoising, but this has been so far demonstrated mainly in brain imaging. To date the combined benefits of motion compensation and anatomical guidance have not been demonstrated for myocardial PET-MR imaging. This work addresses this by proposing a single cardiac PET-MR image reconstruction framework which fully utilises MR-derived information to allow both motion compensation and anatomical guidance within the reconstruction. Methods: Fifteen patients underwent a 18F-FDG cardiac PET-MR scan with a previously introduced acquisition framework. The MR data processing and image reconstruction pipeline produces respiratory motion fields and a high-resolution respiratory motion-corrected MR image with good tissue contrast. This MR-derived information was then included in a respiratory motion-corrected, cardiac-gated, anatomically guided image reconstruction of the simultaneously acquired PET data. Reconstructions were evaluated by measuring myocardial contrast and noise and compared to images from several comparative intermediate methods using the components of the proposed framework separately. Results: Including respiratory motion correction, cardiac gating, and anatomical guidance significantly increased contrast. In particular, myocardium-to-blood pool contrast increased by 143% on average (p<0.0001) compared to conventional uncorrected, non-guided PET images. Furthermore, anatomical guidance significantly reduced image noise compared to non-guided image reconstruction by 16.1% (p<0.0001). Conclusion: The proposed framework for MR-derived motion compensation and anatomical guidance of cardiac PET data was shown to significantly improve image quality compared to alternative reconstruction methods. Each component of the reconstruction pipeline was shown to have a positive impact on the final image quality. These improvements have the potential to improve clinical interpretability and diagnosis based on cardiac PET-MR images.
Collapse
Affiliation(s)
- Camila Munoz
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Sam Ellis
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Stephan G. Nekolla
- Nuklearmedizinische Klinik und Poliklinik, Technische Technical University of Munich, Munich, Germany
| | - Karl P. Kunze
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
- MR Research Collaborations, Siemens Healthcare, Frimley, United Kingdom
| | - Teresa Vitadello
- Department of Internal Medicine I, University Hospital Rechts der Isar, School of Medicine, Technical University of Munich, Munich, Germany; and
| | - Radhouene Neji
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
- MR Research Collaborations, Siemens Healthcare, Frimley, United Kingdom
| | - Rene M. Botnar
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
- Escuela de Ingeniería, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Julia A. Schnabel
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Andrew J. Reader
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Claudia Prieto
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
- Escuela de Ingeniería, Pontificia Universidad Católica de Chile, Santiago, Chile
| |
Collapse
|
4
|
Marquis H, Deidda D, Gillman A, Willowson KP, Gholami Y, Hioki T, Eslick E, Thielemans K, Bailey DL. Theranostic SPECT reconstruction for improved resolution: application to radionuclide therapy dosimetry. EJNMMI Phys 2021; 8:16. [PMID: 33598750 PMCID: PMC7889770 DOI: 10.1186/s40658-021-00362-x] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2020] [Accepted: 02/02/2021] [Indexed: 12/20/2022] Open
Abstract
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 68Ga/177Lu IEC phantom study and a clinical example using 64Cu/67Cu. 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.
Collapse
Affiliation(s)
- H Marquis
- Sydney Vital Translational Cancer Research Centre, Sydney, Australia.,Institute of Medical Physics, University of Sydney, Sydney, Australia
| | - D Deidda
- National Physical Laboratory, Teddington, UK
| | - A Gillman
- Australian e-Health Research Centre, CSIRO, Brisbane, Australia
| | - K P Willowson
- Department of Nuclear Medicine, Royal North Shore Hospital, Sydney, Australia
| | - Y Gholami
- Sydney Vital Translational Cancer Research Centre, Sydney, Australia.,Institute of Medical Physics, University of Sydney, Sydney, Australia.,Faculty of Medicine and Health, University of Sydney, Sydney, Australia
| | - T Hioki
- Institute of Medical Physics, University of Sydney, Sydney, Australia
| | - E Eslick
- Department of Nuclear Medicine, Royal North Shore Hospital, Sydney, Australia
| | - K Thielemans
- Institute of Nuclear Medicine, University College London, London, UK
| | - D L Bailey
- Sydney Vital Translational Cancer Research Centre, Sydney, Australia. .,Department of Nuclear Medicine, Royal North Shore Hospital, Sydney, Australia. .,Faculty of Medicine and Health, University of Sydney, Sydney, Australia.
| |
Collapse
|
5
|
Reader AJ, Ellis S. Bootstrap-Optimised Regularised Image Reconstruction for Emission Tomography. IEEE TRANSACTIONS ON MEDICAL IMAGING 2020; 39:2163-2175. [PMID: 31944935 PMCID: PMC7273977 DOI: 10.1109/tmi.2019.2956878] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/13/2019] [Revised: 11/22/2019] [Accepted: 11/24/2019] [Indexed: 06/10/2023]
Abstract
In emission tomography, iterative image reconstruction from noisy measured data usually results in noisy images, and so regularisation is often used to compensate for noise. However, in practice, an appropriate, automatic and precise specification of the strength of regularisation for image reconstruction from a given noisy measured dataset remains unresolved. Existing approaches are either empirical approximations with no guarantee of generalisation, or else are computationally intensive cross-validation methods requiring full reconstructions for a limited set of preselected regularisation strengths. In contrast, we propose a novel methodology embedded within iterative image reconstruction, using one or more bootstrapped replicates of the measured data for precise optimisation of the regularisation. The approach uses a conventional unregularised iterative update of a current image estimate from the noisy measured data, and then also uses the bootstrap replicate to obtain a bootstrap update of the current image estimate. The method then seeks the regularisation hyperparameters which, when applied to the bootstrap update of the image, lead to a best fit of the regularised bootstrap update to the conventional measured data update. This corresponds to estimating the degree of regularisation needed in order to map the noisy update to a model of the mean of an ensemble of noisy updates. For a given regularised objective function (e.g. penalised likelihood), no hyperparameter selection or tuning is required. The method is demonstrated for positron emission tomography (PET) data at different noise levels, and delivers near-optimal reconstructions (in terms of reconstruction error) without any knowledge of the ground truth, nor any form of training data.
Collapse
Affiliation(s)
- Andrew J. Reader
- School of Biomedical Engineering and Imaging SciencesKing’s College London, King’s Health Partners, St Thomas’ HospitalLondonSE1 7EHU.K.
| | - Sam Ellis
- School of Biomedical Engineering and Imaging SciencesKing’s College London, King’s Health Partners, St Thomas’ HospitalLondonSE1 7EHU.K.
| |
Collapse
|
6
|
Bland J, Mehranian A, Belzunce MA, Ellis S, da Costa‐Luis C, McGinnity CJ, Hammers A, Reader AJ. Intercomparison of MR-informed PET image reconstruction methods. Med Phys 2019; 46:5055-5074. [PMID: 31494961 PMCID: PMC6899618 DOI: 10.1002/mp.13812] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2019] [Revised: 08/23/2019] [Accepted: 08/23/2019] [Indexed: 12/28/2022] Open
Abstract
PURPOSE Numerous image reconstruction methodologies for positron emission tomography (PET) have been developed that incorporate magnetic resonance (MR) imaging structural information, producing reconstructed images with improved suppression of noise and reduced partial volume effects. However, the influence of MR structural information also increases the possibility of suppression or bias of structures present only in the PET data (PET-unique regions). To address this, further developments for MR-informed methods have been proposed, for example, through inclusion of the current reconstructed PET image, alongside the MR image, in the iterative reconstruction process. In this present work, a number of kernel and maximum a posteriori (MAP) methodologies are compared, with the aim of identifying methods that enable a favorable trade-off between the suppression of noise and the retention of unique features present in the PET data. METHODS The reconstruction methods investigated were: the MR-informed conventional and spatially compact kernel methods, referred to as KEM and KEM largest value sparsification (LVS) respectively; the MR-informed Bowsher and Gaussian MR-guided MAP methods; and the PET-MR-informed hybrid kernel and anato-functional MAP methods. The trade-off between improving the reconstruction of the whole brain region and the PET-unique regions was investigated for all methods in comparison with postsmoothed maximum likelihood expectation maximization (MLEM), evaluated in terms of structural similarity index (SSIM), normalized root mean square error (NRMSE), bias, and standard deviation. Both simulated BrainWeb (10 noise realizations) and real [18 F] fluorodeoxyglucose (FDG) three-dimensional datasets were used. The real [18 F]FDG dataset was augmented with simulated tumors to allow comparison of the reconstruction methodologies for the case of known regions of PET-MR discrepancy and evaluated at full counts (100%) and at a reduced (10%) count level. RESULTS For the high-count simulated and real data studies, the anato-functional MAP method performed better than the other methods under investigation (MR-informed, PET-MR-informed and postsmoothed MLEM), in terms of achieving the best trade-off for the reconstruction of the whole brain and PET-unique regions, assessed in terms of the SSIM, NRMSE, and bias vs standard deviation. The inclusion of PET information in the anato-functional MAP method enables the reconstruction of PET-unique regions to attain similarly low levels of bias as unsmoothed MLEM, while moderately improving the whole brain image quality for low levels of regularization. However, for low count simulated datasets the anato-functional MAP method performs poorly, due to the inclusion of noisy PET information in the regularization term. For the low counts simulated dataset, KEM LVS and to a lesser extent, HKEM performed better than the other methods under investigation in terms of achieving the best trade-off for the reconstruction of the whole brain and PET-unique regions, assessed in terms of the SSIM, NRMSE, and bias vs standard deviation. CONCLUSION For the reconstruction of noisy data, multiple MR-informed methods produce favorable whole brain vs PET-unique region trade-off in terms of the image quality metrics of SSIM and NRMSE, comfortably outperforming the whole image denoising of postsmoothed MLEM.
Collapse
Affiliation(s)
- James Bland
- School of Biomedical Engineering and Imaging SciencesKing's College LondonSt Thomas' HospitalLondonSE1 7EHUK
| | - Abolfazl Mehranian
- School of Biomedical Engineering and Imaging SciencesKing's College LondonSt Thomas' HospitalLondonSE1 7EHUK
| | - Martin A. Belzunce
- School of Biomedical Engineering and Imaging SciencesKing's College LondonSt Thomas' HospitalLondonSE1 7EHUK
| | - Sam Ellis
- School of Biomedical Engineering and Imaging SciencesKing's College LondonSt Thomas' HospitalLondonSE1 7EHUK
| | - Casper da Costa‐Luis
- School of Biomedical Engineering and Imaging SciencesKing's College LondonSt Thomas' HospitalLondonSE1 7EHUK
| | - Colm J. McGinnity
- King's College London & Guy's and St Thomas' PET CentreSt Thomas' HospitalLondonSE1 7EHUK
| | - Alexander Hammers
- King's College London & Guy's and St Thomas' PET CentreSt Thomas' HospitalLondonSE1 7EHUK
| | - Andrew J. Reader
- School of Biomedical Engineering and Imaging SciencesKing's College LondonSt Thomas' HospitalLondonSE1 7EHUK
| |
Collapse
|
7
|
Tahaei MS, Reader AJ, Collins DL. Two novel PET image restoration methods guided by PET-MR kernels: Application to brain imaging. Med Phys 2019; 46:2085-2102. [PMID: 30710342 DOI: 10.1002/mp.13418] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2018] [Revised: 12/27/2018] [Accepted: 01/18/2019] [Indexed: 11/06/2022] Open
Abstract
PURPOSE Post-reconstruction positron emission tomography (PET) image restoration methods that take advantage of available anatomical information can play an important role in accurate quantification of PET images. However, when using anatomical information, the resulting PET image may lose resolution in certain regions where the anatomy does not agree with the change in functional activity. In this work, this problem is addressed by using both magnetic resonance (MR) and filtered PET images to guide the denoising process. METHODS In this work, two novel post-reconstruction methods for restoring PET images using the subject's registered T1-weighted MR image are proposed. The first method is based on a representation of the image using basis functions extracted from T1-weighted MR and filtered PET image. The coefficients for these basis functions are estimated using a sparsity-penalized least squares objective function. The second method is a noniterative fast method that uses guided kernel filtering in combination with twicing to restore the noisy PET image. When applied after conventional PVE correction, these methods can be considered as voxel-based MR-guided partial volume effect (PVE) correction methods. RESULTS Using simulation analyses of [ 18 F]FDG PET images of the brain with lesions, the proposed methods are compared to other denoising methods through different figures of merit. The results show promising improvements in image quality as well as reduction in bias and variance of the lesions. We also show the application of the proposed methods on real [ 18 F]FDG data. CONCLUSION Two methods for restoring PET images were proposed. The methods were evaluated on simulation and real brain images. Most MR-guided PVE correction methods are only based on segmented T1-weighted images and their accuracy is very sensitive to segmentation errors, especially in regions of abnormalities and lesions. However, both proposed methods can use the T1-weighted image without segmentation. The simplicity and the very low computational cost of the second method make it suitable for clinical applications and large data studies. The proposed methods can be naturally extended to PVE correction and denoising of other functional modalities using corresponding anatomical information.
Collapse
Affiliation(s)
- Marzieh S Tahaei
- McConnell Brain Imaging Centre, Montreal Neurological Institute, Montreal, Canada.,Department of Biomedical Engineering, McGill University, Montreal, Canada
| | - Andrew J Reader
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, St. Thomas' Hospital, King's College London, London, UK
| | - D Louis Collins
- McConnell Brain Imaging Centre, Montreal Neurological Institute, Montreal, Canada.,Department of Biomedical Engineering, McGill University, Montreal, Canada
| |
Collapse
|