1
|
Gu F, Wu Q. Quantitation of dynamic total-body PET imaging: recent developments and future perspectives. Eur J Nucl Med Mol Imaging 2023; 50:3538-3557. [PMID: 37460750 PMCID: PMC10547641 DOI: 10.1007/s00259-023-06299-w] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Accepted: 06/05/2023] [Indexed: 10/04/2023]
Abstract
BACKGROUND Positron emission tomography (PET) scanning is an important diagnostic imaging technique used in disease diagnosis, therapy planning, treatment monitoring, and medical research. The standardized uptake value (SUV) obtained at a single time frame has been widely employed in clinical practice. Well beyond this simple static measure, more detailed metabolic information can be recovered from dynamic PET scans, followed by the recovery of arterial input function and application of appropriate tracer kinetic models. Many efforts have been devoted to the development of quantitative techniques over the last couple of decades. CHALLENGES The advent of new-generation total-body PET scanners characterized by ultra-high sensitivity and long axial field of view, i.e., uEXPLORER (United Imaging Healthcare), PennPET Explorer (University of Pennsylvania), and Biograph Vision Quadra (Siemens Healthineers), further stimulates valuable inspiration to derive kinetics for multiple organs simultaneously. But some emerging issues also need to be addressed, e.g., the large-scale data size and organ-specific physiology. The direct implementation of classical methods for total-body PET imaging without proper validation may lead to less accurate results. CONCLUSIONS In this contribution, the published dynamic total-body PET datasets are outlined, and several challenges/opportunities for quantitation of such types of studies are presented. An overview of the basic equation, calculation of input function (based on blood sampling, image, population or mathematical model), and kinetic analysis encompassing parametric (compartmental model, graphical plot and spectral analysis) and non-parametric (B-spline and piece-wise basis elements) approaches is provided. The discussion mainly focuses on the feasibilities, recent developments, and future perspectives of these methodologies for a diverse-tissue environment.
Collapse
Affiliation(s)
- Fengyun Gu
- School of Mathematics and Physics, North China Electric Power University, 102206, Beijing, China.
- School of Mathematical Sciences, University College Cork, T12XF62, Cork, Ireland.
| | - Qi Wu
- School of Mathematical Sciences, University College Cork, T12XF62, Cork, Ireland
| |
Collapse
|
2
|
Li Y, Hu J, Sari H, Xue S, Ma R, Kandarpa S, Visvikis D, Rominger A, Liu H, Shi K. A deep neural network for parametric image reconstruction on a large axial field-of-view PET. Eur J Nucl Med Mol Imaging 2023; 50:701-714. [PMID: 36326869 DOI: 10.1007/s00259-022-06003-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Accepted: 10/09/2022] [Indexed: 11/06/2022]
Abstract
PURPOSE The PET scanners with long axial field of view (AFOV) having ~ 20 times higher sensitivity than conventional scanners provide new opportunities for enhanced parametric imaging but suffer from the dramatically increased volume and complexity of dynamic data. This study reconstructed a high-quality direct Patlak Ki image from five-frame sinograms without input function by a deep learning framework based on DeepPET to explore the potential of artificial intelligence reducing the acquisition time and the dependence of input function in parametric imaging. METHODS This study was implemented on a large AFOV PET/CT scanner (Biograph Vision Quadra) and twenty patients were recruited with 18F-fluorodeoxyglucose (18F-FDG) dynamic scans. During training and testing of the proposed deep learning framework, the last five-frame (25 min, 40-65 min post-injection) sinograms were set as input and the reconstructed Patlak Ki images by a nested EM algorithm on the vendor were set as ground truth. To evaluate the image quality of predicted Ki images, mean square error (MSE), peak signal-to-noise ratio (PSNR), and structural similarity index measure (SSIM) were calculated. Meanwhile, a linear regression process was applied between predicted and true Ki means on avid malignant lesions and tumor volume of interests (VOIs). RESULTS In the testing phase, the proposed method achieved excellent MSE of less than 0.03%, high SSIM, and PSNR of ~ 0.98 and ~ 38 dB, respectively. Moreover, there was a high correlation (DeepPET: [Formula: see text]= 0.73, self-attention DeepPET: [Formula: see text]=0.82) between predicted Ki and traditionally reconstructed Patlak Ki means over eleven lesions. CONCLUSIONS The results show that the deep learning-based method produced high-quality parametric images from small frames of projection data without input function. It has much potential to address the dilemma of the long scan time and dependency on input function that still hamper the clinical translation of dynamic PET.
Collapse
Affiliation(s)
- Y Li
- College of Biomedical Engineering & Instrument Science, Zhejiang University, Hangzhou, People's Republic of China.,College of Optical Science and Engineering, Zhejiang University, Hangzhou, People's Republic of China
| | - J Hu
- Department of Nuclear Medicine, Inselpital, Bern University Hospital, University of Bern, Bern, Switzerland
| | - H Sari
- Advanced Clinical Imaging Technology, Siemens Healthcare AG, Lausanne, Switzerland
| | - S Xue
- Department of Nuclear Medicine, Inselpital, Bern University Hospital, University of Bern, Bern, Switzerland
| | - R Ma
- Department of Nuclear Medicine, Inselpital, Bern University Hospital, University of Bern, Bern, Switzerland.,Department of Engineering Physics, Tsinghua University, Beijing, China
| | - S Kandarpa
- LaTIM, INSERM, UMR 1101, University of Brest, Brest, France
| | - D Visvikis
- LaTIM, INSERM, UMR 1101, University of Brest, Brest, France
| | - A Rominger
- Department of Nuclear Medicine, Inselpital, Bern University Hospital, University of Bern, Bern, Switzerland
| | - H Liu
- College of Optical Science and Engineering, Zhejiang University, Hangzhou, People's Republic of China.
| | - K Shi
- Department of Nuclear Medicine, Inselpital, Bern University Hospital, University of Bern, Bern, Switzerland.,Computer Aided Medical Procedures and Augmented Reality, Institute of Informatics I16, Technical University of Munich, Munich, Germany
| |
Collapse
|
3
|
Szirmay-Kalos L, Magdics M, Varnyú D. Direct dynamic tomographic reconstruction without explicit blood input function. Biomed Signal Process Control 2023. [DOI: 10.1016/j.bspc.2022.104313] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
4
|
Daube-Witherspoon ME, Pantel AR, Pryma DA, Karp JS. Total-body PET: a new paradigm for molecular imaging. Br J Radiol 2022; 95:20220357. [PMID: 35993615 PMCID: PMC9733603 DOI: 10.1259/bjr.20220357] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Revised: 06/25/2022] [Accepted: 08/12/2022] [Indexed: 11/05/2022] Open
Abstract
Total body (TB) positron emission tomography (PET) instruments have dramatically changed the paradigm of PET clinical and research studies due to their very high sensitivity and capability to image dynamic radiopharmaceutical distributions in the major organs of the body simultaneously. In this manuscript, we review the design of these systems and discuss general challenges and trade-offs to maximize the performance gains of current TB-PET systems. We then describe new concepts and technology that may impact future TB-PET systems. The manuscript summarizes what has been learned from the initial sites with TB-PET and explores potential research and clinical applications of TB-PET. The current generation of TB-PET systems range in axial field-of-view (AFOV) from 1 to 2 m and serve to illustrate the benefits and opportunities of a longer AFOV for various applications in PET. In only a few years of use these new TB-PET systems have shown that they will play an important role in expanding the field of molecular imaging and benefiting clinical practice.
Collapse
Affiliation(s)
| | - Austin R Pantel
- Department of Radiology, University of Pennsylvania, Philadelphia, United States
| | - Daniel A Pryma
- Department of Radiology, University of Pennsylvania, Philadelphia, United States
| | - Joel S Karp
- Department of Radiology, University of Pennsylvania, Philadelphia, United States
| |
Collapse
|
5
|
Automated GMP Production and Preclinical Evaluation of [ 68Ga]Ga-TEoS-DAZA and [ 68Ga]Ga-TMoS-DAZA. Pharmaceutics 2022; 14:pharmaceutics14122695. [PMID: 36559188 PMCID: PMC9783202 DOI: 10.3390/pharmaceutics14122695] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2022] [Revised: 11/28/2022] [Accepted: 11/29/2022] [Indexed: 12/04/2022] Open
Abstract
[68Ga]Ga-TEoS-DAZA and [68Ga]Ga-TMoS-DAZA are two novel radiotracers suitable for functional PET liver imaging. Due to their specific liver uptake and biliary excretion, the tracers may be applied for segmental liver function quantification, gall tree imaging and the differential diagnosis of liver nodules. The purpose of this study was to investigate problems that occurred initially during the development of the GMP compliant synthesis procedure and to evaluate the tracers in a preclinical model. After low radiolabeling yields were attributed to precursor instability at high temperatures, an optimized radiolabeling procedure was established. Quality controls were in accordance with Ph. Eur. requirements and gave compliant results, although the method for the determination of the 68Ga colloid is partially inhibited due to the presence of a radioactive by-product. The determination of logP revealed [68Ga]Ga-TEoS-DAZA (ethoxy bearing) to be more lipophilic than [68Ga]Ga-TMoS-DAZA (methoxy bearing). Accordingly, biodistribution studies in an in ovo model showed a higher liver uptake for [68Ga]Ga-TEoS-DAZA. In dynamic in ovo PET imaging, rapid tracer accumulation in the liver was observed. Similarly, the activity in the intestines rose steadily within the first hour p.i., indicating biliary excretion. As [68Ga]Ga-TEoS-DAZA and [68Ga]Ga-TMoS-DAZA can be prepared according to GMP guidelines, transition into the early clinical phase is now possible.
Collapse
|
6
|
Huh Y, Shrestha UM, Gullberg GT, Seo Y. Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition. Front Cardiovasc Med 2022; 9:871967. [PMID: 35911544 PMCID: PMC9326051 DOI: 10.3389/fcvm.2022.871967] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2022] [Accepted: 06/16/2022] [Indexed: 11/24/2022] Open
Abstract
Purpose Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.
Collapse
Affiliation(s)
- Yoonsuk Huh
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, San Francisco, CA, United States
| | - Uttam M. Shrestha
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, San Francisco, CA, United States
| | - Grant T. Gullberg
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, San Francisco, CA, United States
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA, United States
| | - Youngho Seo
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, San Francisco, CA, United States
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA, United States
- Department of Nuclear Engineering, University of California, Berkeley, Berkeley, CA, United States
- *Correspondence: Youngho Seo,
| |
Collapse
|
7
|
Bevington CWJ, Cheng JC, Sossi V. A 4-D Iterative HYPR Denoising Operator Improves PET Image Quality. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES 2022. [DOI: 10.1109/trpms.2021.3123537] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Connor W. J. Bevington
- Department of Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada
| | - Ju-Chieh Cheng
- Department of Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada
| | - Vesna Sossi
- Department of Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada
| |
Collapse
|
8
|
Lamare F, Bousse A, Thielemans K, Liu C, Merlin T, Fayad H, Visvikis D. PET respiratory motion correction: quo vadis? Phys Med Biol 2021; 67. [PMID: 34915465 DOI: 10.1088/1361-6560/ac43fc] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2021] [Accepted: 12/16/2021] [Indexed: 11/12/2022]
Abstract
Positron emission tomography (PET) respiratory motion correction has been a subject of great interest for the last twenty years, prompted mainly by the development of multimodality imaging devices such as PET/computed tomography (CT) and PET/magnetic resonance imaging (MRI). PET respiratory motion correction involves a number of steps including acquisition synchronization, motion estimation and finally motion correction. The synchronization steps include the use of different external device systems or data driven approaches which have been gaining ground over the last few years. Patient specific or generic motion models using the respiratory synchronized datasets can be subsequently derived and used for correction either in the image space or within the image reconstruction process. Similar overall approaches can be considered and have been proposed for both PET/CT and PET/MRI devices. Certain variations in the case of PET/MRI include the use of MRI specific sequences for the registration of respiratory motion information. The proposed review includes a comprehensive coverage of all these areas of development in field of PET respiratory motion for different multimodality imaging devices and approaches in terms of synchronization, estimation and subsequent motion correction. Finally, a section on perspectives including the potential clinical usage of these approaches is included.
Collapse
Affiliation(s)
- Frederic Lamare
- Nuclear Medicine Department, University Hospital Centre Bordeaux Hospital Group South, ., Bordeaux, Nouvelle-Aquitaine, 33604, FRANCE
| | - Alexandre Bousse
- LaTIM, INSERM UMR1101, Université de Bretagne Occidentale, ., Brest, Bretagne, 29285, FRANCE
| | - Kris Thielemans
- University College London Institute of Nuclear Medicine, UCL Hospital, Tower 5, 235 Euston Road, London, NW1 2BU, UNITED KINGDOM OF GREAT BRITAIN AND NORTHERN IRELAND
| | - Chi Liu
- Department of Diagnostic Radiology, Yale University School of Medicine Department of Radiology and Biomedical Imaging, PO Box 208048, 801 Howard Avenue, New Haven, Connecticut, 06520-8042, UNITED STATES
| | - Thibaut Merlin
- LaTIM, INSERM UMR1101, Universite de Bretagne Occidentale, ., Brest, Bretagne, 29285, FRANCE
| | - Hadi Fayad
- Weill Cornell Medicine - Qatar, ., Doha, ., QATAR
| | - Dimitris Visvikis
- LaTIM, UMR1101, Universite de Bretagne Occidentale, INSERM, Brest, Bretagne, 29285, FRANCE
| |
Collapse
|
9
|
Sanaat A, Mirsadeghi E, Razeghi B, Ginovart N, Zaidi H. Fast dynamic brain PET imaging using stochastic variational prediction for recurrent frame generation. Med Phys 2021; 48:5059-5071. [PMID: 34174787 PMCID: PMC8518550 DOI: 10.1002/mp.15063] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2021] [Revised: 05/30/2021] [Accepted: 06/08/2021] [Indexed: 12/03/2022] Open
Abstract
Purpose We assess the performance of a recurrent frame generation algorithm for prediction of late frames from initial frames in dynamic brain PET imaging. Methods Clinical dynamic 18F‐DOPA brain PET/CT studies of 46 subjects with ten folds cross‐validation were retrospectively employed. A novel stochastic adversarial video prediction model was implemented to predict the last 13 frames (25–90 minutes) from the initial 13 frames (0–25 minutes). The quantitative analysis of the predicted dynamic PET frames was performed for the test and validation dataset using established metrics. Results The predicted dynamic images demonstrated that the model is capable of predicting the trend of change in time‐varying tracer biodistribution. The Bland‐Altman plots reported the lowest tracer uptake bias (−0.04) for the putamen region and the smallest variance (95% CI: −0.38, +0.14) for the cerebellum. The region‐wise Patlak graphical analysis in the caudate and putamen regions for eight subjects from the test and validation dataset showed that the average bias for Ki and distribution volume was 4.3%, 5.1% and 4.4%, 4.2%, (P‐value <0.05), respectively. Conclusion We have developed a novel deep learning approach for fast dynamic brain PET imaging capable of generating the last 65 minutes time frames from the initial 25 minutes frames, thus enabling significant reduction in scanning time.
Collapse
Affiliation(s)
- Amirhossein Sanaat
- Division of Nuclear Medicine and Molecular Imaging, Geneva University Hospital, Geneva, Switzerland
| | - Ehsan Mirsadeghi
- Electrical Engineering Department, Amirkabir University of Technology, Tehran, Iran
| | - Behrooz Razeghi
- Department of Computer Sciences, University of Geneva, Geneva, Switzerland.,School of Engineering and Applied Sciences, Harvard University, Boston, USA
| | - Nathalie Ginovart
- Department of Psychiatry, Geneva University, Geneva, Switzerland.,Department of Basic Neurosciences, Geneva University, Geneva, Switzerland
| | - Habib Zaidi
- Division of Nuclear Medicine and Molecular Imaging, Geneva University Hospital, Geneva, Switzerland.,Geneva University Neurocenter, Geneva University, Geneva, Switzerland.,Department of Nuclear Medicine and Molecular Imaging, University of Groningen, Groningen, Netherlands.,University Medical Center, Groningen, Netherlands.,Department of Nuclear Medicine, University of Southern Denmark, Odense, Denmark
| |
Collapse
|
10
|
Espinós-Morató H, Cascales-Picó D, Vergara M, Hernández-Martínez Á, Benlloch Baviera JM, Rodríguez-Álvarez MJ. Simulation Study of a Frame-Based Motion Correction Algorithm for Positron Emission Imaging. SENSORS 2021; 21:s21082608. [PMID: 33917742 PMCID: PMC8068167 DOI: 10.3390/s21082608] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/26/2021] [Revised: 04/02/2021] [Accepted: 04/04/2021] [Indexed: 11/16/2022]
Abstract
Positron emission tomography (PET) is a functional non-invasive imaging modality that uses radioactive substances (radiotracers) to measure changes in metabolic processes. Advances in scanner technology and data acquisition in the last decade have led to the development of more sophisticated PET devices with good spatial resolution (1–3 mm of full width at half maximum (FWHM)). However, there are involuntary motions produced by the patient inside the scanner that lead to image degradation and potentially to a misdiagnosis. The adverse effect of the motion in the reconstructed image increases as the spatial resolution of the current scanners continues improving. In order to correct this effect, motion correction techniques are becoming increasingly popular and further studied. This work presents a simulation study of an image motion correction using a frame-based algorithm. The method is able to cut the acquired data from the scanner in frames, taking into account the size of the object of study. This approach allows working with low statistical information without losing image quality. The frames are later registered using spatio-temporal registration developed in a multi-level way. To validate these results, several performance tests are applied to a set of simulated moving phantoms. The results obtained show that the method minimizes the intra-frame motion, improves the signal intensity over the background in comparison with other literature methods, produces excellent values of similarity with the ground-truth (static) image and is able to find a limit in the patient-injected dose when some prior knowledge of the lesion is present.
Collapse
Affiliation(s)
- Héctor Espinós-Morató
- Instituto de Instrumentación para Imagen Molecular (i3M), Centro Mixto CSIC—Universitat Politècnica de València, 46022 Valencia, Spain; (D.C.-P.); (M.V.); (Á.H.-M.); (J.M.B.B.); (M.J.R.-Á.)
- Correspondence:
| | - David Cascales-Picó
- Instituto de Instrumentación para Imagen Molecular (i3M), Centro Mixto CSIC—Universitat Politècnica de València, 46022 Valencia, Spain; (D.C.-P.); (M.V.); (Á.H.-M.); (J.M.B.B.); (M.J.R.-Á.)
| | - Marina Vergara
- Instituto de Instrumentación para Imagen Molecular (i3M), Centro Mixto CSIC—Universitat Politècnica de València, 46022 Valencia, Spain; (D.C.-P.); (M.V.); (Á.H.-M.); (J.M.B.B.); (M.J.R.-Á.)
- Department of Imaging and Pathology, Division of Nuclear Medicine, KU Leuven, 3000 Leuven, Belgium
| | - Ángel Hernández-Martínez
- Instituto de Instrumentación para Imagen Molecular (i3M), Centro Mixto CSIC—Universitat Politècnica de València, 46022 Valencia, Spain; (D.C.-P.); (M.V.); (Á.H.-M.); (J.M.B.B.); (M.J.R.-Á.)
| | - José María Benlloch Baviera
- Instituto de Instrumentación para Imagen Molecular (i3M), Centro Mixto CSIC—Universitat Politècnica de València, 46022 Valencia, Spain; (D.C.-P.); (M.V.); (Á.H.-M.); (J.M.B.B.); (M.J.R.-Á.)
| | - María José Rodríguez-Álvarez
- Instituto de Instrumentación para Imagen Molecular (i3M), Centro Mixto CSIC—Universitat Politècnica de València, 46022 Valencia, Spain; (D.C.-P.); (M.V.); (Á.H.-M.); (J.M.B.B.); (M.J.R.-Á.)
| |
Collapse
|
11
|
Cheng JCK, Bevington C, Rahmim A, Klyuzhin I, Matthews J, Boellaard R, Sossi V. Dynamic PET image reconstruction utilizing intrinsic data-driven HYPR4D denoising kernel. Med Phys 2021; 48:2230-2244. [PMID: 33533050 DOI: 10.1002/mp.14751] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2020] [Revised: 10/16/2020] [Accepted: 01/28/2021] [Indexed: 11/07/2022] Open
Abstract
PURPOSE Reconstructed PET images are typically noisy, especially in dynamic imaging where the acquired data are divided into several short temporal frames. High noise in the reconstructed images translates to poor precision/reproducibility of image features. One important role of "denoising" is therefore to improve the precision of image features. However, typical denoising methods achieve noise reduction at the expense of accuracy. In this work, we present a novel four-dimensional (4D) denoised image reconstruction framework, which we validate using 4D simulations, experimental phantom, and clinical patient data, to achieve 4D noise reduction while preserving spatiotemporal patterns/minimizing error introduced by denoising. METHODS Our proposed 4D denoising operator/kernel is based on HighlY constrained backPRojection (HYPR), which is applied either after each update of OSEM reconstruction of dynamic 4D PET data or within the recently proposed kernelized reconstruction framework inspired by kernel methods in machine learning. Our HYPR4D kernel makes use of the spatiotemporal high frequency features extracted from a 4D composite, generated within the reconstruction, to preserve the spatiotemporal patterns and constrain the 4D noise increment of the image estimate. RESULTS Results from simulations, experimental phantom, and patient data showed that the HYPR4D kernel with our proposed 4D composite outperformed other denoising methods, such as the standard OSEM with spatial filter, OSEM with 4D filter, and HYPR kernel method with the conventional 3D composite in conjunction with recently proposed High Temporal Resolution kernel (HYPRC3D-HTR), in terms of 4D noise reduction while preserving the spatiotemporal patterns or 4D resolution within the 4D image estimate. Consequently, the error in outcome measures obtained from the HYPR4D method was less dependent on the region size, contrast, and uniformity/functional patterns within the target structures compared to the other methods. For outcome measures that depend on spatiotemporal tracer uptake patterns such as the nondisplaceable Binding Potential (BPND ), the root mean squared error in regional mean of voxel BPND values was reduced from ~8% (OSEM with spatial or 4D filter) to ~3% using HYPRC3D-HTR and was further reduced to ~2% using our proposed HYPR4D method for relatively small target structures (~10 mm in diameter). At the voxel level, HYPR4D produced two to four times lower mean absolute error in BPND relative to HYPRC3D-HTR. CONCLUSION As compared to conventional methods, our proposed HYPR4D method can produce more robust and accurate image features without requiring any prior information.
Collapse
Affiliation(s)
- Ju-Chieh Kevin Cheng
- Pacific Parkinson's Research Centre, The University of British Columbia, 2215 Wesbrook Mall, Vancouver, BC, V6T 1Z3, Canada.,Department of Physics and Astronomy, The University of British Columbia, 6224 Agricultural Road, Vancouver, BC, V6T 1Z1, Canada
| | - Connor Bevington
- Department of Physics and Astronomy, The University of British Columbia, 6224 Agricultural Road, Vancouver, BC, V6T 1Z1, Canada
| | - Arman Rahmim
- Department of Physics and Astronomy, The University of British Columbia, 6224 Agricultural Road, Vancouver, BC, V6T 1Z1, Canada.,Department of Radiology, University of British Columbia, Vancouver, BC, V5Z 1M9, Canada
| | - Ivan Klyuzhin
- Department of Medicine, Division of Neurology, University of British Columbia, Vancouver, BC, V6T 2B5, Canada
| | - Julian Matthews
- Division of Neuroscience and Experimental Psychology, Wolfson Molecular Imaging Centre, The University of Manchester, Manchester, M20 3LJ, UK
| | - Ronald Boellaard
- Department of Radiology and Nuclear Medicine, VU University Medical Center, De Boelelaan 1117, Amsterdam, 1081 HV, Netherlands.,Department of Nuclear Medicine and Molecular Imaging, University Medical Center Groningen, University of Groningen, Hanzeplein 1, 9713 KC, Groningen, Netherlands
| | - Vesna Sossi
- Department of Physics and Astronomy, The University of British Columbia, 6224 Agricultural Road, Vancouver, BC, V6T 1Z1, Canada
| |
Collapse
|
12
|
Abstract
Total-body PET image reconstruction follows a similar procedure to the image reconstruction process for standard whole-body PET scanners. One unique aspect of total-body imaging is simultaneous coverage of the entire human body, which makes it convenient to perform total-body dynamic PET scans. Therefore, four-dimensional dynamic PET reconstruction and parametric imaging are of great interest in total-body imaging. This article covers some basics of PET image reconstruction and then focuses on three- and four-dimensional PET reconstruction for total-body imaging. Methods for image formation from raw measurements in total-body PET are described. Challenges and opportunities in total-body PET image reconstruction are discussed.
Collapse
Affiliation(s)
- Jinyi Qi
- Department of Biomedical Engineering, University of California, One Shields Avenue, Davis, CA 95616, USA.
| | - Samuel Matej
- Department of Radiology, University of Pennsylvania, 3620 Hamilton Walk, John Morgan Building, Room 156A, Philadelphia, PA 19104-6061, USA
| | - Guobao Wang
- Department of Radiology, University of California Davis Medical Center, Lawrence J. Ellison Ambulatory Care Center Building, Suite 3100, 4860 Y Street, Sacramento, CA 95817, USA
| | - Xuezhu Zhang
- Department of Biomedical Engineering, University of California, One Shields Avenue, Davis, CA 95616, USA
| |
Collapse
|
13
|
Hu J, Panin V, Smith AM, Spottiswoode B, Shah V, CA von Gall C, Baker M, Howe W, Kehren F, Casey M, Bendriem B. Design and Implementation of Automated Clinical Whole Body Parametric PET With Continuous Bed Motion. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES 2020. [DOI: 10.1109/trpms.2020.2994316] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
14
|
Kim K, Gong K, Moon SH, El Fakhri G, Normandin MD, Li Q. Penalized Parametric PET Image Estimation Using Local Linear Fitting. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES 2020. [DOI: 10.1109/trpms.2020.3024123] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
15
|
Petibon Y, Alpert NM, Ouyang J, Pizzagalli DA, Cusin C, Fava M, El Fakhri G, Normandin MD. PET imaging of neurotransmission using direct parametric reconstruction. Neuroimage 2020; 221:117154. [PMID: 32679252 PMCID: PMC7800040 DOI: 10.1016/j.neuroimage.2020.117154] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2020] [Revised: 07/06/2020] [Accepted: 07/07/2020] [Indexed: 11/18/2022] Open
Abstract
Receptor ligand-based dynamic Positron Emission Tomography (PET) permits the measurement of neurotransmitter release in the human brain. For single-scan paradigms, the conventional method of estimating changes in neurotransmitter levels relies on fitting a pharmacokinetic model to activity concentration histories extracted after PET image reconstruction. However, due to the statistical fluctuations of activity concentration data at the voxel scale, parametric images computed using this approach often exhibit low signal-to-noise ratio, impeding characterization of neurotransmitter release. Numerous studies have shown that direct parametric reconstruction (DPR) approaches, which combine image reconstruction and kinetic analysis in a unified framework, can improve the signal-to-noise ratio of parametric mapping. However, there is little experience with DPR in imaging of neurotransmission and the performance of the approach in this application has not been evaluated before in humans. In this report, we present and evaluate a DPR methodology that computes 3-D distributions of ligand transport, binding potential (BPND) and neurotransmitter release magnitude (γ) from a dynamic sequence of PET sinograms. The technique employs the linear simplified reference region model (LSRRM) of Alpert et al. (2003), which represents an extension of the simplified reference region model that incorporates time-varying binding parameters due to radioligand displacement by release of neurotransmitter. Estimation of parametric images is performed by gradient-based optimization of a Poisson log-likelihood function incorporating LSRRM kinetics and accounting for the effects of head movement, attenuation, detector sensitivity, random and scattered coincidences. A 11C-raclopride simulation study showed that the proposed approach substantially reduces the bias and variance of voxel-wise γ estimates as compared to standard methods. Moreover, simulations showed that detection of release could be made more reliable and/or conducted using a smaller sample size using the proposed DPR estimator. Likewise, images of BPND computed using DPR had substantially improved bias and variance properties. Application of the method in human subjects was demonstrated using 11C-raclopride dynamic scans and a reward task, confirming the improved quality of the estimated parametric images using the proposed approach.
Collapse
Affiliation(s)
- Yoann Petibon
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA.
| | - Nathaniel M Alpert
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Jinsong Ouyang
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Diego A Pizzagalli
- Center for Depression, Anxiety & Stress Research, McLean Hospital and Harvard Medical School, Belmont, MA, USA
| | - Cristina Cusin
- Department of Psychiatry, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Maurizio Fava
- Department of Psychiatry, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Georges El Fakhri
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Marc D Normandin
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| |
Collapse
|
16
|
Wang G, Rahmim A, Gunn RN. PET Parametric Imaging: Past, Present, and Future. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES 2020; 4:663-675. [PMID: 33763624 PMCID: PMC7983029 DOI: 10.1109/trpms.2020.3025086] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Positron emission tomography (PET) is actively used in a diverse range of applications in oncology, cardiology, and neurology. The use of PET in the clinical setting focuses on static (single time frame) imaging at a specific time-point post radiotracer injection and is typically considered as semi-quantitative; e.g. standardized uptake value (SUV) measures. In contrast, dynamic PET imaging requires increased acquisition times but has the advantage that it measures the full spatiotemporal distribution of a radiotracer and, in combination with tracer kinetic modeling, enables the generation of multiparametric images that more directly quantify underlying biological parameters of interest, such as blood flow, glucose metabolism, and receptor binding. Parametric images have the potential for improved detection and for more accurate and earlier therapeutic response assessment. Parametric imaging with dynamic PET has witnessed extensive research in the past four decades. In this paper, we provide an overview of past and present activities and discuss emerging opportunities in the field of parametric imaging for the future.
Collapse
Affiliation(s)
- Guobao Wang
- Department of Radiology, University of California Davis Health, Sacramento, CA 95817, USA
| | - Arman Rahmim
- University of British Columbia, Vancouver, BC, Canada
| | | |
Collapse
|
17
|
Vergalasova I, Cai J. A modern review of the uncertainties in volumetric imaging of respiratory-induced target motion in lung radiotherapy. Med Phys 2020; 47:e988-e1008. [PMID: 32506452 DOI: 10.1002/mp.14312] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2019] [Revised: 05/15/2020] [Accepted: 05/26/2020] [Indexed: 12/25/2022] Open
Abstract
Radiotherapy has become a critical component for the treatment of all stages and types of lung cancer, often times being the primary gateway to a cure. However, given that radiation can cause harmful side effects depending on how much surrounding healthy tissue is exposed, treatment of the lung can be particularly challenging due to the presence of moving targets. Careful implementation of every step in the radiotherapy process is absolutely integral for attaining optimal clinical outcomes. With the advent and now widespread use of stereotactic body radiation therapy (SBRT), where extremely large doses are delivered, accurate, and precise dose targeting is especially vital to achieve an optimal risk to benefit ratio. This has largely become possible due to the rapid development of image-guided technology. Although imaging is critical to the success of radiotherapy, it can often be plagued with uncertainties due to respiratory-induced target motion. There has and continues to be an immense research effort aimed at acknowledging and addressing these uncertainties to further our abilities to more precisely target radiation treatment. Thus, the goal of this article is to provide a detailed review of the prevailing uncertainties that remain to be investigated across the different imaging modalities, as well as to highlight the more modern solutions to imaging motion and their role in addressing the current challenges.
Collapse
Affiliation(s)
- Irina Vergalasova
- Department of Radiation Oncology, Rutgers Cancer Institute of New Jersey, Rutgers University, New Brunswick, NJ, USA
| | - Jing Cai
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Kowloon, Hong Kong
| |
Collapse
|
18
|
|
19
|
Direct Parametric Maps Estimation from Dynamic PET Data: An Iterated Conditional Modes Approach. JOURNAL OF HEALTHCARE ENGINEERING 2019; 2018:5942873. [PMID: 30073047 PMCID: PMC6057340 DOI: 10.1155/2018/5942873] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/20/2017] [Revised: 03/27/2018] [Accepted: 05/08/2018] [Indexed: 11/24/2022]
Abstract
We propose and test a novel approach for direct parametric image reconstruction of dynamic PET data. We present a theoretical description of the problem of PET direct parametric maps estimation as an inference problem, from a probabilistic point of view, and we derive a simple iterative algorithm, based on the Iterated Conditional Mode (ICM) framework, which exploits the simplicity of a two-step optimization and the efficiency of an analytic method for estimating kinetic parameters from a nonlinear compartmental model. The resulting method is general enough to be flexible to an arbitrary choice of the kinetic model, and unlike many other solutions, it is capable to deal with nonlinear compartmental models without the need for linearization. We tested its performance on a two-tissue compartment model, including an analytical solution to the kinetic parameters evaluation, based on an auxiliary parameter set, with the aim of reducing computation errors and approximations. The new method is tested on simulated and clinical data. Simulation analysis led to the conclusion that the proposed algorithm gives a good estimation of the kinetic parameters in any noise condition. Furthermore, the application of the proposed method to clinical data gave promising results for further studies.
Collapse
|
20
|
Lingala SG, Guo Y, Bliesener Y, Zhu Y, Lebel RM, Law M, Nayak KS. Tracer kinetic models as temporal constraints during brain tumor DCE-MRI reconstruction. Med Phys 2019; 47:37-51. [PMID: 31663134 PMCID: PMC6980286 DOI: 10.1002/mp.13885] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Revised: 10/17/2019] [Accepted: 10/17/2019] [Indexed: 12/11/2022] Open
Abstract
Purpose To apply tracer kinetic models as temporal constraints during reconstruction of under‐sampled brain tumor dynamic contrast enhanced (DCE) magnetic resonance imaging (MRI). Methods A library of concentration vs time profiles is simulated for a range of physiological kinetic parameters. The library is reduced to a dictionary of temporal bases, where each profile is approximated by a sparse linear combination of the bases. Image reconstruction is formulated as estimation of concentration profiles and sparse model coefficients with a fixed sparsity level. Simulations are performed to evaluate modeling error, and error statistics in kinetic parameter estimation in presence of noise. Retrospective under‐sampling experiments are performed on a brain tumor DCE digital reference object (DRO), and 12 brain tumor in‐vivo 3T datasets. The performances of the proposed under‐sampled reconstruction scheme and an existing compressed sensing‐based temporal finite‐difference (tFD) under‐sampled reconstruction were compared against the fully sampled inverse Fourier Transform‐based reconstruction. Results Simulations demonstrate that sparsity levels of 2 and 3 model the library profiles from the Patlak and extended Tofts‐Kety (ETK) models, respectively. Noise sensitivity analysis showed equivalent kinetic parameter estimation error statistics from noisy concentration profiles, and model approximated profiles. DRO‐based experiments showed good fidelity in recovery of kinetic maps from 20‐fold under‐sampled data. In‐vivo experiments demonstrated reduced bias and uncertainty in kinetic mapping with the proposed approach compared to tFD at under‐sampled reduction factors >= 20. Conclusions Tracer kinetic models can be applied as temporal constraints during brain tumor DCE‐MRI reconstruction. The proposed under‐sampled scheme resulted in model parameter estimates less biased with respect to conventional fully sampled DCE MRI reconstructions and parameter estimation. The approach is flexible, can use nonlinear kinetic models, and does not require tuning of regularization parameters.
Collapse
Affiliation(s)
- Sajan Goud Lingala
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, IA, USA
| | - Yi Guo
- Snap Inc., San Francisco, CA, USA
| | - Yannick Bliesener
- Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, CA, USA
| | | | - R Marc Lebel
- GE Healthcare Applied Sciences Laboratory, Calgary, Canada
| | - Meng Law
- Department of Neuroscience, Monash University, Melbourne, Australia
| | - Krishna S Nayak
- Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, CA, USA
| |
Collapse
|
21
|
Naseri M, Rajabi H, Wang J, Abbasi M, Kalantari F. Simultaneous respiratory motion correction and image reconstruction in 4D-multi pinhole small animal SPECT. Med Phys 2019; 46:5047-5054. [PMID: 31495940 DOI: 10.1002/mp.13807] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2019] [Revised: 08/22/2019] [Accepted: 08/22/2019] [Indexed: 11/05/2022] Open
Abstract
PURPOSE Respiratory motion in the chest region during single photon emission computed tomography (SPECT) is a major degrading factor that reduces the accuracy of image quantification. This effect is more notable when the tumor is very small, or the spatial resolution of the imaging system is less than the respiratory motion amplitude. Small animals imaging systems with sub-millimeter spatial resolution need more attention to the respiratory motion for quantitative studies. We developed a motion-embedded four-dimensional (4D)-multi pinhole SPECT (MPS) reconstruction algorithm for respiratory motion correction. This algorithm makes full use of projection statistics for reconstruction of every individual frame. METHODS The ROBY phantom with small tumors in liver was generated in eight different phases during one respiratory cycle. The MPS projections were modeled using a fast ray tracing method simulating an MPS acquisition. Individual frames were reconstructed and used for motion estimation. The Demons non-rigid registration algorithm was used to calculate deformation vector fields (DVFs) for simultaneous motion correction and image reconstruction. A motion-embedded 4D-MPS method was used to reconstruct images using all the projections and corresponding DVFs, simultaneously. The 4D-MPS reconstructed images were compared to the low-count single frame (LCSF) reconstructed image, the three-dimensional (3D)-MPS images reconstructed using individual frames, and post reconstruction registration (PRR) that aligns all individual phases to a reference frame using Demons-derived DVFs. The tumor volume relative error (TVE), tumor contrast relative error (TCE), and dice index (DI) for 2, 3, and 4 mm liver were calculated and compared for different reconstruction methods. RESULTS For the 4D-MPS reconstruction method, TVE was reduced and DI was higher compared to PRR, 3D-MPS, and LCSF. The extent of the improvement was higher for the small tumor size (i.e. 2 mm). For the biggest tumor in contrast 3 (i.e. 4 mm) TVE for 4D-MPS, PRR, 3D-MPS and, LCSF were 1.33%, 8%, 8%, and 14.67%, respectively. CONCLUSIONS The results suggest that motion-embedded 4D-MPS method is an effective and practical way for respiratory motion correction. It reconstructs high quality gated frames while using all projection data to reconstruct each frame.
Collapse
Affiliation(s)
- Maryam Naseri
- Medical Physics Program, Department of Physics and Astronomy, Louisiana State University, Baton Rouge, LA, USA.,Department of Medical Physics, Faculty of Medicine, Tarbiat Modares University, Tehran, Iran
| | - Hossein Rajabi
- Department of Medical Physics, Faculty of Medicine, Tarbiat Modares University, Tehran, Iran
| | - Jing Wang
- Department of Radiation Oncology, UT Southwestern Medical Center Dallas, Dallas, TX, USA
| | - Mehrshad Abbasi
- Department of Nuclear Medicine, Vali-Asr Hospital, Tehran University of Medical Sciences, Tehran, Iran
| | - Faraz Kalantari
- Department of Radiation Oncology, University of Arkansas for Medical Sciences, Little Rock, AR, USA
| |
Collapse
|
22
|
High Temporal-Resolution Dynamic PET Image Reconstruction Using a New Spatiotemporal Kernel Method. IEEE TRANSACTIONS ON MEDICAL IMAGING 2019; 38:664-674. [PMID: 30222553 PMCID: PMC6422751 DOI: 10.1109/tmi.2018.2869868] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Current clinical dynamic PET has an effective temporal resolution of 5-10 seconds, which can be adequate for traditional compartmental modeling but is inadequate for exploiting the benefit of more advanced tracer kinetic modeling for characterization of diseases (e.g., cancer and heart disease). There is a need to improve dynamic PET to allow fine temporal sampling of 1-2 seconds. However, the reconstruction of these short-time frames from tomographic data is extremely challenging as the count level of each frame is very low and high noise presents in both spatial and temporal domains. Previously, the kernel framework has been developed and demonstrated as a statistically efficient approach to utilizing image prior for low-count PET image reconstruction. Nevertheless, the existing kernel methods mainly explore spatial correlations in the data and only have a limited ability in suppressing temporal noise. In this paper, we propose a new kernel method which extends the previous spatial kernel method to the general spatiotemporal domain. The new kernelized model encodes both spatial and temporal correlations obtained from image prior information and are incorporated into the PET forward projection model to improve themaximumlikelihood(ML) image reconstruction. Computer simulations and an application to real patient scan have shown that the proposed approach can achieve effective noise reduction in both spatial and temporal domains and outperform the spatial kernel method and conventional ML reconstruction method for improving the high temporal-resolution dynamic PET imaging.
Collapse
|
23
|
Rahmim A, Lodge MA, Karakatsanis NA, Panin VY, Zhou Y, McMillan A, Cho S, Zaidi H, Casey ME, Wahl RL. Dynamic whole-body PET imaging: principles, potentials and applications. Eur J Nucl Med Mol Imaging 2018; 46:501-518. [PMID: 30269154 DOI: 10.1007/s00259-018-4153-6] [Citation(s) in RCA: 108] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2018] [Accepted: 08/28/2018] [Indexed: 02/07/2023]
Abstract
PURPOSE In this article, we discuss dynamic whole-body (DWB) positron emission tomography (PET) as an imaging tool with significant clinical potential, in relation to conventional standard uptake value (SUV) imaging. BACKGROUND DWB PET involves dynamic data acquisition over an extended axial range, capturing tracer kinetic information that is not available with conventional static acquisition protocols. The method can be performed within reasonable clinical imaging times, and enables generation of multiple types of PET images with complementary information in a single imaging session. Importantly, DWB PET can be used to produce multi-parametric images of (i) Patlak slope (influx rate) and (ii) intercept (referred to sometimes as "distribution volume"), while also providing (iii) a conventional 'SUV-equivalent' image for certain protocols. RESULTS We provide an overview of ongoing efforts (primarily focused on FDG PET) and discuss potential clinically relevant applications. CONCLUSION Overall, the framework of DWB imaging [applicable to both PET/CT(computed tomography) and PET/MRI (magnetic resonance imaging)] generates quantitative measures that may add significant value to conventional SUV image-derived measures, with limited pitfalls as we also discuss in this work.
Collapse
Affiliation(s)
- Arman Rahmim
- Department of Radiology and Radiological Science, Johns Hopkins University, JHOC Building Room 3245, 601 N. Caroline St, Baltimore, MD, 21287, USA. .,Departments of Radiology and Physics & Astronomy, University of British Columbia, Vancouver, BC, V5Z 1M9, Canada.
| | - Martin A Lodge
- Department of Radiology and Radiological Science, Johns Hopkins University, JHOC Building Room 3245, 601 N. Caroline St, Baltimore, MD, 21287, USA
| | | | | | - Yun Zhou
- Department of Radiology and Radiological Science, Johns Hopkins University, JHOC Building Room 3245, 601 N. Caroline St, Baltimore, MD, 21287, USA
| | - Alan McMillan
- Department of Radiology, University of Wisconsin, Madison, WI, 53705, USA
| | - Steve Cho
- Department of Radiology, University of Wisconsin, Madison, WI, 53705, USA
| | - Habib Zaidi
- Division of Nuclear Medicine and Molecular Imaging, Geneva University Hospital, CH-1211, Geneva, Switzerland
| | | | - Richard L Wahl
- Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, MO, 63110, USA
| |
Collapse
|
24
|
Zaidi H, Alavi A, Naqa IE. Novel Quantitative PET Techniques for Clinical Decision Support in Oncology. Semin Nucl Med 2018; 48:548-564. [PMID: 30322481 DOI: 10.1053/j.semnuclmed.2018.07.003] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Quantitative image analysis has deep roots in the usage of positron emission tomography (PET) in clinical and research settings to address a wide variety of diseases. It has been extensively employed to assess molecular and physiological biomarkers in vivo in healthy and disease states, in oncology, cardiology, neurology, and psychiatry. Quantitative PET allows relating the time-varying activity concentration in tissues/organs of interest and the basic functional parameters governing the biological processes being studied. Yet, quantitative PET is challenged by a number of degrading physical factors related to the physics of PET imaging, the limitations of the instrumentation used, and the physiological status of the patient. Moreover, there is no consensus on the most reliable and robust image-derived PET metric(s) that can be used with confidence in clinical oncology owing to the discrepancies between the conclusions reported in the literature. There is also increasing interest in the use of artificial intelligence based techniques, particularly machine learning and deep learning techniques in a variety of applications to extract quantitative features (radiomics) from PET including image segmentation and outcome prediction in clinical oncology. These novel techniques are revolutionizing clinical practice and are now offering unique capabilities to the clinical molecular imaging community and biomedical researchers at large. In this report, we summarize recent developments and future tendencies in quantitative PET imaging and present example applications in clinical decision support to illustrate its potential in the context of clinical oncology.
Collapse
Affiliation(s)
- Habib Zaidi
- Division of Nuclear Medicine and Molecular Imaging, Geneva University Hospital, Geneva, Switzerland; Geneva Neuroscience Centre, University of Geneva, Geneva, Switzerland; Department of Nuclear Medicine and Molecular Imaging, University of Groningen, Groningen, the Netherlands; Department of Nuclear Medicine, University of Southern Denmark, Odense, Denmark.
| | - Abass Alavi
- Department of Radiology, Hospital of the University of Pennsylvania, Philadelphia, PA
| | - Issam El Naqa
- Department of Radiation Oncology, University of Michigan, Ann Arbor, MI
| |
Collapse
|
25
|
Zuo Y, Qi J, Wang G. Relative Patlak plot for dynamic PET parametric imaging without the need for early-time input function. Phys Med Biol 2018; 63:165004. [PMID: 30020080 DOI: 10.1088/1361-6560/aad444] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
The Patlak graphical method is widely used in parametric imaging for modeling irreversible radiotracer kinetics in dynamic PET. The net influx rate of radiotracer can be determined from the slope of the Patlak plot. The implementation of the standard Patlak method requires the knowledge of full-time input function from the injection time until the scan end time, which presents a challenge for use in the clinic. This paper proposes a new relative Patlak plot method that does not require early-time input function and therefore can be more efficient for parametric imaging. Theoretical analysis proves that the effect of early-time input function is a constant scaling factor on the Patlak slope estimation. Thus, the parametric image of the slope of the relative Patlak plot is related to the parametric image of standard Patlak slope by a global scaling factor. This theoretical finding has been further demonstrated by computer simulation and real patient data. The study indicates that parametric imaging of the relative Patlak slope can be used as a substitute of parametric imaging of standard Patlak slope for tasks that do not require absolute quantification, such as lesion detection and tumor volume segmentation.
Collapse
Affiliation(s)
- Yang Zuo
- Department of Radiology, University of California at Davis, Sacramento, CA 95817, United States of America
| | | | | |
Collapse
|
26
|
Ahnen M, Becker R, Buck A, Casella C, Commichau V, Calafiori DD, Dissertori G, Eleftheriou A, Fischer J, Howard AS, Ito M, Khateri P, Kim J, Lustermann W, Ritzer C, Roser U, Rudin M, Solevi P, Tsoumpas C, Warnock G, Weber B, Wyss M, Zagozdzinska-Bochenek A. Performance Measurements of the SAFIR Prototype Detector With the STiC ASIC Readout. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES 2018. [DOI: 10.1109/trpms.2018.2797484] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
27
|
Merlin T, Visvikis D, Fernandez P, Lamare F. Dynamic PET image reconstruction integrating temporal regularization associated with respiratory motion correction for applications in oncology. ACTA ACUST UNITED AC 2018; 63:045012. [DOI: 10.1088/1361-6560/aaa86a] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
28
|
Lu L, Ma X, Mohy-Ud-Din H, Ma J, Feng Q, Rahmim A, Chen W. Enhancement of dynamic myocardial perfusion PET images based on low-rank plus sparse decomposition. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2018; 154:57-69. [PMID: 29249347 DOI: 10.1016/j.cmpb.2017.10.020] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/17/2017] [Revised: 08/30/2017] [Accepted: 10/16/2017] [Indexed: 06/07/2023]
Abstract
BACKGROUND AND OBJECTIVE The absolute quantification of dynamic myocardial perfusion (MP) PET imaging is challenged by the limited spatial resolution of individual frame images due to division of the data into shorter frames. This study aims to develop a method for restoration and enhancement of dynamic PET images. METHODS We propose that the image restoration model should be based on multiple constraints rather than a single constraint, given the fact that the image characteristic is hardly described by a single constraint alone. At the same time, it may be possible, but not optimal, to regularize the image with multiple constraints simultaneously. Fortunately, MP PET images can be decomposed into a superposition of background vs. dynamic components via low-rank plus sparse (L + S) decomposition. Thus, we propose an L + S decomposition based MP PET image restoration model and express it as a convex optimization problem. An iterative soft thresholding algorithm was developed to solve the problem. Using realistic dynamic 82Rb MP PET scan data, we optimized and compared its performance with other restoration methods. RESULTS The proposed method resulted in substantial visual as well as quantitative accuracy improvements in terms of noise versus bias performance, as demonstrated in extensive 82Rb MP PET simulations. In particular, the myocardium defect in the MP PET images had improved visual as well as contrast versus noise tradeoff. The proposed algorithm was also applied on an 8-min clinical cardiac 82Rb MP PET study performed on the GE Discovery PET/CT, and demonstrated improved quantitative accuracy (CNR and SNR) compared to other algorithms. CONCLUSIONS The proposed method is effective for restoration and enhancement of dynamic PET images.
Collapse
Affiliation(s)
- Lijun Lu
- School of Biomedical Engineering and Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, Guangdong 510515, China.
| | - Xiaomian Ma
- School of Software, Guangdong Food and Drug Vocational College, Guangzhou, Guangdong 510520, China
| | - Hassan Mohy-Ud-Din
- Shaukat Khanum Memorial Cancer Hospital and Research Center, Lahore, Pakistan
| | - Jianhua Ma
- School of Biomedical Engineering and Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, Guangdong 510515, China
| | - Qianjin Feng
- School of Biomedical Engineering and Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, Guangdong 510515, China.
| | - Arman Rahmim
- Department of Radiology, Johns Hopkins University, Baltimore, MD 21287, USA; Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, MD 21218, USA
| | - Wufan Chen
- School of Biomedical Engineering and Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, Guangdong 510515, China.
| |
Collapse
|
29
|
Hutton BF, Erlandsson K, Thielemans K. Advances in clinical molecular imaging instrumentation. Clin Transl Imaging 2018. [DOI: 10.1007/s40336-018-0264-0] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
|
30
|
Cabello J, Ziegler SI. Advances in PET/MR instrumentation and image reconstruction. Br J Radiol 2018; 91:20160363. [PMID: 27376170 PMCID: PMC5966194 DOI: 10.1259/bjr.20160363] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2016] [Revised: 06/26/2016] [Accepted: 06/29/2016] [Indexed: 12/15/2022] Open
Abstract
The combination of positron emission tomography (PET) and MRI has attracted the attention of researchers in the past approximately 20 years in small-animal imaging and more recently in clinical research. The combination of PET/MRI allows researchers to explore clinical and research questions in a wide number of fields, some of which are briefly mentioned here. An important number of groups have developed different concepts to tackle the problems that PET instrumentation poses to the exposition of electromagnetic fields. We have described most of these research developments in preclinical and clinical experiments, including the few commercial scanners available. From the software perspective, an important number of algorithms have been developed to address the attenuation correction issue and to exploit the possibility that MRI provides for motion correction and quantitative image reconstruction, especially parametric modelling of radiopharmaceutical kinetics. In this work, we give an overview of some exemplar applications of simultaneous PET/MRI, together with technological hardware and software developments.
Collapse
Affiliation(s)
- Jorge Cabello
- Nuklearmedizinische Klinik und Poliklinik, Klinikum rechts der Isar, Technische Universität München, Munich, Germany
| | - Sibylle I Ziegler
- Nuklearmedizinische Klinik und Poliklinik, Klinikum rechts der Isar, Technische Universität München, Munich, Germany
| |
Collapse
|
31
|
Germino M, Carson RE. Cardiac-gated parametric images from 82 Rb PET from dynamic frames and direct 4D reconstruction. Med Phys 2017; 45:639-654. [PMID: 29205378 DOI: 10.1002/mp.12710] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2017] [Revised: 11/07/2017] [Accepted: 11/07/2017] [Indexed: 12/23/2022] Open
Abstract
PURPOSE Cardiac perfusion PET data can be reconstructed as a dynamic sequence and kinetic modeling performed to quantify myocardial blood flow, or reconstructed as static gated images to quantify function. Parametric images from dynamic PET are conventionally not gated, to allow use of all events with lower noise. An alternative method for dynamic PET is to incorporate the kinetic model into the reconstruction algorithm itself, bypassing the generation of a time series of emission images and directly producing parametric images. So-called "direct reconstruction" can produce parametric images with lower noise than the conventional method because the noise distribution is more easily modeled in projection space than in image space. In this work, we develop direct reconstruction of cardiac-gated parametric images for 82 Rb PET with an extension of the Parametric Motion compensation OSEM List mode Algorithm for Resolution-recovery reconstruction for the one tissue model (PMOLAR-1T). METHODS PMOLAR-1T was extended to accommodate model terms to account for spillover from the left and right ventricles into the myocardium. The algorithm was evaluated on a 4D simulated 82 Rb dataset, including a perfusion defect, as well as a human 82 Rb list mode acquisition. The simulated list mode was subsampled into replicates, each with counts comparable to one gate of a gated acquisition. Parametric images were produced by the indirect (separate reconstructions and modeling) and direct methods for each of eight low-count and eight normal-count replicates of the simulated data, and each of eight cardiac gates for the human data. For the direct method, two initialization schemes were tested: uniform initialization, and initialization with the filtered iteration 1 result of the indirect method. For the human dataset, event-by-event respiratory motion compensation was included. The indirect and direct methods were compared for the simulated dataset in terms of bias and coefficient of variation as a function of iteration. RESULTS Convergence of direct reconstruction was slow with uniform initialization; lower bias was achieved in fewer iterations by initializing with the filtered indirect iteration 1 images. For most parameters and regions evaluated, the direct method achieved the same or lower absolute bias at matched iteration as the indirect method, with 23%-65% lower noise. Additionally, the direct method gave better contrast between the perfusion defect and surrounding normal tissue than the indirect method. Gated parametric images from the human dataset had comparable relative performance of indirect and direct, in terms of mean parameter values per iteration. Changes in myocardial wall thickness and blood pool size across gates were readily visible in the gated parametric images, with higher contrast between myocardium and left ventricle blood pool in parametric images than gated SUV images. CONCLUSIONS Direct reconstruction can produce parametric images with less noise than the indirect method, opening the potential utility of gated parametric imaging for perfusion PET.
Collapse
Affiliation(s)
- Mary Germino
- Department of Biomedical Engineering, Yale University, P. O. Box 208048, New Haven, CT, 06520-8048, USA
| | - Richard E Carson
- Department of Biomedical Engineering, Yale University, P. O. Box 208048, New Haven, CT, 06520-8048, USA.,Department of Radiology and Biomedical Imaging, Yale University School of Medicine, P. O. Box 208048, New Haven, CT, 06520-8048, USA
| |
Collapse
|
32
|
Beijst C, Kunnen B, Lam MGEH, de Jong HWAM. Technical Advances in Image Guidance of Radionuclide Therapy. J Nucl Med Technol 2017; 45:272-279. [PMID: 29042472 DOI: 10.2967/jnmt.117.190991] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2017] [Accepted: 07/05/2017] [Indexed: 11/16/2022] Open
Abstract
Internal radiation therapy with radionuclides (i.e., radionuclide therapy) owes its success to the many advantages over other, more conventional, treatment options. One distinct advantage of radionuclide therapies is the potential to use (part of) the emitted radiation for imaging of the radionuclide distribution. The combination of diagnostic and therapeutic properties in a set of matched radiopharmaceuticals (sometimes combined in a single radiopharmaceutical) is often referred to as theranostics and allows accurate diagnostic imaging before therapy. The use of imaging benefits treatment planning, dosimetry, and assessment of treatment response. This paper focuses on a selection of advances in imaging technology relevant for image guidance of radionuclide therapy. This involves developments in nuclear imaging modalities, as well as other anatomic and functional imaging modalities. The quality and quantitative accuracy of images used for guidance of radionuclide therapy is continuously being improved, which in turn may improve the therapeutic outcome and efficiency of radionuclide therapies.
Collapse
Affiliation(s)
- Casper Beijst
- Department of Radiology and Nuclear Medicine, UMC Utrecht, Utrecht, The Netherlands; and .,Image Sciences Institute, UMC Utrecht, Utrecht, The Netherlands
| | - Britt Kunnen
- Department of Radiology and Nuclear Medicine, UMC Utrecht, Utrecht, The Netherlands; and.,Image Sciences Institute, UMC Utrecht, Utrecht, The Netherlands
| | - Marnix G E H Lam
- Department of Radiology and Nuclear Medicine, UMC Utrecht, Utrecht, The Netherlands; and
| | - Hugo W A M de Jong
- Department of Radiology and Nuclear Medicine, UMC Utrecht, Utrecht, The Netherlands; and
| |
Collapse
|
33
|
Cui J, Liu X, Wang Y, Liu H. Deep reconstruction model for dynamic PET images. PLoS One 2017; 12:e0184667. [PMID: 28934254 PMCID: PMC5608245 DOI: 10.1371/journal.pone.0184667] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2017] [Accepted: 08/28/2017] [Indexed: 11/18/2022] Open
Abstract
Accurate and robust tomographic reconstruction from dynamic positron emission tomography (PET) acquired data is a difficult problem. Conventional methods, such as the maximum likelihood expectation maximization (MLEM) algorithm for reconstructing the activity distribution-based on individual frames, may lead to inaccurate results due to the checkerboard effect and limitation of photon counts. In this paper, we propose a stacked sparse auto-encoder based reconstruction framework for dynamic PET imaging. The dynamic reconstruction problem is formulated in a deep learning representation, where the encoding layers extract the prototype features, such as edges, so that, in the decoding layers, the reconstructed results are obtained through a combination of those features. The qualitative and quantitative results of the procedure, including the data based on a Monte Carlo simulation and real patient data demonstrates the effectiveness of our method.
Collapse
Affiliation(s)
- Jianan Cui
- State Key Laboratory of Modern Optical Instrumentation, Department of Optical Engineering, Zhejiang University, Hangzhou, China
| | - Xin Liu
- Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzheng, China
| | - Yile Wang
- State Key Laboratory of Modern Optical Instrumentation, Department of Optical Engineering, Zhejiang University, Hangzhou, China
| | - Huafeng Liu
- State Key Laboratory of Modern Optical Instrumentation, Department of Optical Engineering, Zhejiang University, Hangzhou, China
- * E-mail:
| |
Collapse
|
34
|
Cheng JC(K, Matthews J, Sossi V, Anton-Rodriguez J, Salomon A, Boellaard R. Incorporating HYPR de-noising within iterative PET reconstruction (HYPR-OSEM). ACTA ACUST UNITED AC 2017. [DOI: 10.1088/1361-6560/aa7b66] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
35
|
Germino M, Gallezot JD, Yan J, Carson RE. Direct reconstruction of parametric images for brain PET with event-by-event motion correction: evaluation in two tracers across count levels. Phys Med Biol 2017; 62:5344-5364. [PMID: 28504644 PMCID: PMC5783541 DOI: 10.1088/1361-6560/aa731f] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Parametric images for dynamic positron emission tomography (PET) are typically generated by an indirect method, i.e. reconstructing a time series of emission images, then fitting a kinetic model to each voxel time activity curve. Alternatively, 'direct reconstruction', incorporates the kinetic model into the reconstruction algorithm itself, directly producing parametric images from projection data. Direct reconstruction has been shown to achieve parametric images with lower standard error than the indirect method. Here, we present direct reconstruction for brain PET using event-by-event motion correction of list-mode data, applied to two tracers. Event-by-event motion correction was implemented for direct reconstruction in the Parametric Motion-compensation OSEM List-mode Algorithm for Resolution-recovery reconstruction. The direct implementation was tested on simulated and human datasets with tracers [11C]AFM (serotonin transporter) and [11C]UCB-J (synaptic density), which follow the 1-tissue compartment model. Rigid head motion was tracked with the Vicra system. Parametric images of K 1 and distribution volume (V T = K 1/k 2) were compared to those generated by the indirect method by regional coefficient of variation (CoV). Performance across count levels was assessed using sub-sampled datasets. For simulated and real datasets at high counts, the two methods estimated K 1 and V T with comparable accuracy. At lower count levels, the direct method was substantially more robust to outliers than the indirect method. Compared to the indirect method, direct reconstruction reduced regional K 1 CoV by 35-48% (simulated dataset), 39-43% ([11C]AFM dataset) and 30-36% ([11C]UCB-J dataset) across count levels (averaged over regions at matched iteration); V T CoV was reduced by 51-58%, 54-60% and 30-46%, respectively. Motion correction played an important role in the dataset with larger motion: correction increased regional V T by 51% on average in the [11C]UCB-J dataset. Direct reconstruction of dynamic brain PET with event-by-event motion correction is achievable and dramatically more robust to noise in V T images than the indirect method.
Collapse
Affiliation(s)
- Mary Germino
- Department of Biomedical Engineering, Yale University, New Haven, CT, United States of America
| | | | | | | |
Collapse
|
36
|
Petibon Y, Rakvongthai Y, Fakhri GE, Ouyang J. Direct parametric reconstruction in dynamic PET myocardial perfusion imaging: in vivo studies. Phys Med Biol 2017; 62:3539-3565. [PMID: 28379843 PMCID: PMC5739089 DOI: 10.1088/1361-6560/aa6394] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Dynamic PET myocardial perfusion imaging (MPI) used in conjunction with tracer kinetic modeling enables the quantification of absolute myocardial blood flow (MBF). However, MBF maps computed using the traditional indirect method (i.e. post-reconstruction voxel-wise fitting of kinetic model to PET time-activity-curves-TACs) suffer from poor signal-to-noise ratio (SNR). Direct reconstruction of kinetic parameters from raw PET projection data has been shown to offer parametric images with higher SNR compared to the indirect method. The aim of this study was to extend and evaluate the performance of a direct parametric reconstruction method using in vivo dynamic PET MPI data for the purpose of quantifying MBF. Dynamic PET MPI studies were performed on two healthy pigs using a Siemens Biograph mMR scanner. List-mode PET data for each animal were acquired following a bolus injection of ~7-8 mCi of 18F-flurpiridaz, a myocardial perfusion agent. Fully-3D dynamic PET sinograms were obtained by sorting the coincidence events into 16 temporal frames covering ~5 min after radiotracer administration. Additionally, eight independent noise realizations of both scans-each containing 1/8th of the total number of events-were generated from the original list-mode data. Dynamic sinograms were then used to compute parametric maps using the conventional indirect method and the proposed direct method. For both methods, a one-tissue compartment model accounting for spillover from the left and right ventricle blood-pools was used to describe the kinetics of 18F-flurpiridaz. An image-derived arterial input function obtained from a TAC taken in the left ventricle cavity was used for tracer kinetic analysis. For the indirect method, frame-by-frame images were estimated using two fully-3D reconstruction techniques: the standard ordered subset expectation maximization (OSEM) reconstruction algorithm on one side, and the one-step late maximum a posteriori (OSL-MAP) algorithm on the other side, which incorporates a quadratic penalty function. The parametric images were then calculated using voxel-wise weighted least-square fitting of the reconstructed myocardial PET TACs. For the direct method, parametric images were estimated directly from the dynamic PET sinograms using a maximum a posteriori (MAP) parametric reconstruction algorithm which optimizes an objective function comprised of the Poisson log-likelihood term, the kinetic model and a quadratic penalty function. Maximization of the objective function with respect to each set of parameters was achieved using a preconditioned conjugate gradient algorithm with a specifically developed pre-conditioner. The performance of the direct method was evaluated by comparing voxel- and segment-wise estimates of [Formula: see text], the tracer transport rate (ml · min-1 · ml-1), to those obtained using the indirect method applied to both OSEM and OSL-MAP dynamic reconstructions. The proposed direct reconstruction method produced [Formula: see text] maps with visibly lower noise than the indirect method based on OSEM and OSL-MAP reconstructions. At normal count levels, the direct method was shown to outperform the indirect method based on OSL-MAP in the sense that at matched level of bias, reduced regional noise levels were obtained. At lower count levels, the direct method produced [Formula: see text] estimates with significantly lower standard deviation across noise realizations than the indirect method based on OSL-MAP at matched bias level. In all cases, the direct method yielded lower noise and standard deviation than the indirect method based on OSEM. Overall, the proposed direct reconstruction offered a better bias-variance tradeoff than the indirect method applied to either OSEM and OSL-MAP. Direct parametric reconstruction as applied to in vivo dynamic PET MPI data is therefore a promising method for producing MBF maps with lower variance.
Collapse
Affiliation(s)
- Yoann Petibon
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, USA
- Department of Radiology, Harvard Medical School, Boston, USA
| | - Yothin Rakvongthai
- Division of Nuclear Medicine, Department of Radiology, Faculty of Medicine, Chulalongkorn University, Bangkok, Thailand
| | - Georges El Fakhri
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, USA
- Department of Radiology, Harvard Medical School, Boston, USA
| | - Jinsong Ouyang
- Gordon Center for Medical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, USA
- Department of Radiology, Harvard Medical School, Boston, USA
| |
Collapse
|
37
|
Oliver JF, Rafecas M. Modelling Random Coincidences in Positron Emission Tomography by Using Singles and Prompts: A Comparison Study. PLoS One 2016; 11:e0162096. [PMID: 27603143 PMCID: PMC5014417 DOI: 10.1371/journal.pone.0162096] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2016] [Accepted: 08/17/2016] [Indexed: 11/18/2022] Open
Abstract
Random coincidences degrade the image in Positron Emission Tomography, PET. To compensate for their degradation effects, the rate of random coincidences should be estimated. Under certain circumstances, current estimation methods fail to provide accurate results. We propose a novel method, "Singles-Prompts" (SP), that includes the information conveyed by prompt coincidences and models the pile-up. The SP method has the same structure than the well-known "Singles Rate" (SR) approach. Hence, SP can straightforwardly replace SR. In this work, the SP method has been extensively assessed and compared to two conventional methods, SR and the delayed window (DW) method, in a preclinical PET scenario using Monte-Carlo simulations. SP offers accurate estimates for the randoms rates, while SR and DW tend to overestimate the rates (∼10%, and 5%, respectively). With pile-up, the SP method is more robust than SR (but less than DW). At the image level, the contrast is overestimated in SR-corrected images, +16%, while SP produces the correct value. Spill-over is slightly reduced using SP instead of SR. The DW images values are similar to those of SP except for low-statistic scenarios, where DW behaves as if randoms were not compensated for. In particular, the contrast is reduced, -16%. In general, the better estimations of SP translate into better image quality.
Collapse
Affiliation(s)
- Josep F Oliver
- Instituto de Física Corpuscular (IFIC - UV/CSIC), Valencia, Spain
| | - M Rafecas
- Instituto de Física Corpuscular (IFIC - UV/CSIC), Valencia, Spain
| |
Collapse
|
38
|
Karakatsanis NA, Casey ME, Lodge MA, Rahmim A, Zaidi H. Whole-body direct 4D parametric PET imaging employing nested generalized Patlak expectation-maximization reconstruction. Phys Med Biol 2016; 61:5456-85. [PMID: 27383991 DOI: 10.1088/0031-9155/61/15/5456] [Citation(s) in RCA: 58] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Whole-body (WB) dynamic PET has recently demonstrated its potential in translating the quantitative benefits of parametric imaging to the clinic. Post-reconstruction standard Patlak (sPatlak) WB graphical analysis utilizes multi-bed multi-pass PET acquisition to produce quantitative WB images of the tracer influx rate K i as a complimentary metric to the semi-quantitative standardized uptake value (SUV). The resulting K i images may suffer from high noise due to the need for short acquisition frames. Meanwhile, a generalized Patlak (gPatlak) WB post-reconstruction method had been suggested to limit K i bias of sPatlak analysis at regions with non-negligible (18)F-FDG uptake reversibility; however, gPatlak analysis is non-linear and thus can further amplify noise. In the present study, we implemented, within the open-source software for tomographic image reconstruction platform, a clinically adoptable 4D WB reconstruction framework enabling efficient estimation of sPatlak and gPatlak images directly from dynamic multi-bed PET raw data with substantial noise reduction. Furthermore, we employed the optimization transfer methodology to accelerate 4D expectation-maximization (EM) convergence by nesting the fast image-based estimation of Patlak parameters within each iteration cycle of the slower projection-based estimation of dynamic PET images. The novel gPatlak 4D method was initialized from an optimized set of sPatlak ML-EM iterations to facilitate EM convergence. Initially, realistic simulations were conducted utilizing published (18)F-FDG kinetic parameters coupled with the XCAT phantom. Quantitative analyses illustrated enhanced K i target-to-background ratio (TBR) and especially contrast-to-noise ratio (CNR) performance for the 4D versus the indirect methods and static SUV. Furthermore, considerable convergence acceleration was observed for the nested algorithms involving 10-20 sub-iterations. Moreover, systematic reduction in K i % bias and improved TBR were observed for gPatlak versus sPatlak. Finally, validation on clinical WB dynamic data demonstrated the clinical feasibility and superior K i CNR performance for the proposed 4D framework compared to indirect Patlak and SUV imaging.
Collapse
Affiliation(s)
- Nicolas A Karakatsanis
- Division of Nuclear Medicine and Molecular Imaging, School of Medicine, University of Geneva, Geneva, CH-1211, Switzerland
| | | | | | | | | |
Collapse
|
39
|
Kim K, Son YD, Bresler Y, Cho ZH, Ra JB, Ye JC. Dynamic PET reconstruction using temporal patch-based low rank penalty for ROI-based brain kinetic analysis. Phys Med Biol 2016; 60:2019-46. [PMID: 25675392 DOI: 10.1088/0031-9155/60/5/2019] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Dynamic positron emission tomography (PET) is widely used to measure changes in the bio-distribution of radiopharmaceuticals within particular organs of interest over time. However, to retain sufficient temporal resolution, the number of photon counts in each time frame must be limited. Therefore, conventional reconstruction algorithms such as the ordered subset expectation maximization (OSEM) produce noisy reconstruction images, thus degrading the quality of the extracted time activity curves (TACs). To address this issue, many advanced reconstruction algorithms have been developed using various spatio-temporal regularizations. In this paper, we extend earlier results and develop a novel temporal regularization, which exploits the self-similarity of patches that are collected in dynamic images. The main contribution of this paper is to demonstrate that the correlation of patches can be exploited using a low-rank constraint that is insensitive to global intensity variations. The resulting optimization framework is, however, non-Lipschitz and nonconvex due to the Poisson log-likelihood and low-rank penalty terms. Direct application of the conventional Poisson image deconvolution by an augmented Lagrangian (PIDAL) algorithm is, however, problematic due to its large memory requirements, which prevents its parallelization. Thus, we propose a novel optimization framework using the concave-convex procedure (CCCP)
Collapse
Affiliation(s)
- Kyungsang Kim
- Bio Imaging Signal Processing Lab., Department of Bio/Brain Engineering, KAIST, 291 Daehak-ro, Yuseong-gu, Daejeon 305-701, Korea
| | | | | | | | | | | |
Collapse
|
40
|
Shrestha UM, Seo Y, Botvinick EH, Gullberg GT. Image reconstruction in higher dimensions: myocardial perfusion imaging of tracer dynamics with cardiac motion due to deformation and respiration. Phys Med Biol 2015; 60:8275-301. [PMID: 26450115 DOI: 10.1088/0031-9155/60/21/8275] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Myocardial perfusion imaging (MPI) using slow rotating large field of view cameras requires spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration. In vivo, MPI contains additional degrees of freedom involving unavoidable motion of the heart due to quasiperiodic beating and the effects of respiration, which can severely degrade the quality of the images. This work develops a technique for a single photon emission computed tomography (SPECT) that reconstructs the distribution of the radiotracer concentration in the myocardium using a tensor product of different sets of basis functions that approximately describe the spatiotemporal variation of the radiotracer concentration and the motion of the heart. In this study the temporal B-spline basis functions are chosen to reflect the dynamics of the radiotracer, while the intrinsic deformation and the extrinsic motion of the heart are described by a product of a discrete set of Gaussian basis functions. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it deforms due to cardiac beating, and is displaced due to respiratory motion. These results are compared with the conventional 4D-spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. The higher dimensional reconstruction method proposed here improves bias, yet the signal-to-noise ratio (SNR) decreases slightly due to redistribution of the counts over the cardiac-respiratory gates. Additionally, there is a trade-off between the number of gates and the number of projections per gate to achieve high contrast images.
Collapse
Affiliation(s)
- Uttam M Shrestha
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, CA, USA. Structural Biology and Imaging Department, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | | | | | | |
Collapse
|
41
|
Sparse/Low Rank Constrained Reconstruction for Dynamic PET Imaging. PLoS One 2015; 10:e0142019. [PMID: 26540274 PMCID: PMC4634927 DOI: 10.1371/journal.pone.0142019] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2015] [Accepted: 10/15/2015] [Indexed: 11/30/2022] Open
Abstract
In dynamic Positron Emission Tomography (PET), an estimate of the radio activity concentration is obtained from a series of frames of sinogram data taken at ranging in duration from 10 seconds to minutes under some criteria. So far, all the well-known reconstruction algorithms require known data statistical properties. It limits the speed of data acquisition, besides, it is unable to afford the separated information about the structure and the variation of shape and rate of metabolism which play a major role in improving the visualization of contrast for some requirement of the diagnosing in application. This paper presents a novel low rank-based activity map reconstruction scheme from emission sinograms of dynamic PET, termed as SLCR representing Sparse/Low Rank Constrained Reconstruction for Dynamic PET Imaging. In this method, the stationary background is formulated as a low rank component while variations between successive frames are abstracted to the sparse. The resulting nuclear norm and l1 norm related minimization problem can also be efficiently solved by many recently developed numerical methods. In this paper, the linearized alternating direction method is applied. The effectiveness of the proposed scheme is illustrated on three data sets.
Collapse
|
42
|
Loeb R, Navab N, Ziegler SI. Direct Parametric Reconstruction Using Anatomical Regularization for Simultaneous PET/MRI Data. IEEE TRANSACTIONS ON MEDICAL IMAGING 2015; 34:2233-2247. [PMID: 25935030 DOI: 10.1109/tmi.2015.2427777] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Pharmacokinetic analysis of dynamic positron emission tomography (PET) imaging data maps the measured time activity curves to a set of model-specific pharmacokinetic parameters. Voxel-based parameter estimation via curve fitting is conventionally performed indirectly on a sequence of independently reconstructed PET images, leading to high variance and bias in the parametric images. We propose a direct parametric reconstruction algorithm with raw projection data as input that leverages high-resolution anatomical information simultaneously obtained from magnetic resonance (MR) imaging in a PET/MRI scanner for regularization. The reconstruction problem is formulated in a flexible Bayesian framework with Gaussian Markov Random field modeling of activity, parameters, or both simultaneously. MR information is incorporated through a Bowsher-like prior function. Optimization transfer using an expectation-maximization surrogate and a new Bowsher-like penalty surrogate is applied to obtain a voxel-separable algorithm that interleaves a reconstruction with a fitting step. An analytical input function model is used. The algorithm is evaluated on simulated [(18)F]FDG and clinical [(18)F]FET brain data acquired with a Biograph mMR. The results indicate that direct and simultaneously regularized parametric reconstruction increases image quality. Anatomical regularization leads to higher contrast than conventional distance-weighted regularization.
Collapse
|
43
|
Karakatsanis NA, Zhou Y, Lodge MA, Casey ME, Wahl RL, Zaidi H, Rahmim A. Generalized whole-body Patlak parametric imaging for enhanced quantification in clinical PET. Phys Med Biol 2015; 60:8643-73. [PMID: 26509251 DOI: 10.1088/0031-9155/60/22/8643] [Citation(s) in RCA: 56] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
We recently developed a dynamic multi-bed PET data acquisition framework to translate the quantitative benefits of Patlak voxel-wise analysis to the domain of routine clinical whole-body (WB) imaging. The standard Patlak (sPatlak) linear graphical analysis assumes irreversible PET tracer uptake, ignoring the effect of FDG dephosphorylation, which has been suggested by a number of PET studies. In this work: (i) a non-linear generalized Patlak (gPatlak) model is utilized, including a net efflux rate constant kloss, and (ii) a hybrid (s/g)Patlak (hPatlak) imaging technique is introduced to enhance contrast to noise ratios (CNRs) of uptake rate Ki images. Representative set of kinetic parameter values and the XCAT phantom were employed to generate realistic 4D simulation PET data, and the proposed methods were additionally evaluated on 11 WB dynamic PET patient studies. Quantitative analysis on the simulated Ki images over 2 groups of regions-of-interest (ROIs), with low (ROI A) or high (ROI B) true kloss relative to Ki, suggested superior accuracy for gPatlak. Bias of sPatlak was found to be 16-18% and 20-40% poorer than gPatlak for ROIs A and B, respectively. By contrast, gPatlak exhibited, on average, 10% higher noise than sPatlak. Meanwhile, the bias and noise levels for hPatlak always ranged between the other two methods. In general, hPatlak was seen to outperform all methods in terms of target-to-background ratio (TBR) and CNR for all ROIs. Validation on patient datasets demonstrated clinical feasibility for all Patlak methods, while TBR and CNR evaluations confirmed our simulation findings, and suggested presence of non-negligible kloss reversibility in clinical data. As such, we recommend gPatlak for highly quantitative imaging tasks, while, for tasks emphasizing lesion detectability (e.g. TBR, CNR) over quantification, or for high levels of noise, hPatlak is instead preferred. Finally, gPatlak and hPatlak CNR was systematically higher compared to routine SUV values.
Collapse
Affiliation(s)
- Nicolas A Karakatsanis
- Division of Nuclear Medicine and Molecular Imaging, School of Medicine, University of Geneva, Geneva, CH-1211, Switzerland
| | | | | | | | | | | | | |
Collapse
|
44
|
Salavati A, Borofsky S, Boon-Keng TK, Houshmand S, Khiewvan B, Saboury B, Codreanu I, Torigian DA, Zaidi H, Alavi A. Application of partial volume effect correction and 4D PET in the quantification of FDG avid lung lesions. Mol Imaging Biol 2015; 17:140-8. [PMID: 25080325 DOI: 10.1007/s11307-014-0776-6] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
PURPOSE The aim of this study is to assess a software-based method with semiautomated correction for partial volume effect (PVE) to quantify the metabolic activity of pulmonary malignancies in patients who underwent non-gated and respiratory-gated 2-deoxy-2-[(18)F]fluoro-D-glucose (FDG)-positron emission tomography (PET)/x-ray computed tomography(CT). PROCEDURES The study included 106 lesions of 55 lung cancer patients who underwent respiratory-gated FDG-PET/CT for radiation therapy treatment planning. Volumetric PET/CT parameters were determined by using 4D PET/CT and non-gated PET/CT images. We used a semiautomated program employing an adaptive contrast-oriented thresholding algorithm for lesion delineation as well as a lesion-based partial volume effect correction algorithm. We compared respiratory-gated parameters with non-gated parameters by using pairwise comparison and interclass correlation coefficient assessment. In a multivariable regression analysis, we also examined factors, which can affect quantification accuracy, including the size of lesion and the location of tumor. RESULTS This study showed that quantification of volumetric parameters of 4D PET/CT images using an adaptive contrast-oriented thresholding algorithm and 3D lesion-based partial volume correction is feasible. We observed slight increase in FDG uptake by using PET/CT volumetric parameters in comparison of highest respiratory-gated values with non-gated values. After correction for partial volume effect, the mean standardized uptake value (SUVmean) and total lesion glycolysis (TLG) increased substantially (p value <0.001). However, we did not observe a clinically significant difference between partial volume corrected parameters of respiratory-gated and non-gated PET/CT scans. Regression analysis showed that tumor volume was the main predictor of quantification inaccuracy caused by partial volume effect. CONCLUSIONS Based on this study, assessment of volumetric PET/CT parameters and partial volume effect correction for accurate quantification of lung malignant lesions by using respiratory non-gated PET images are feasible and it is comparable to gated measurements. Partial volume correction increased both the respiratory-gated and non-gated values significantly and appears to be the dominant source of quantification error of lung lesions.
Collapse
Affiliation(s)
- Ali Salavati
- Department of Radiology, Hospital of the University of Pennsylvania, 3400 Spruce Street, Philadelphia, PA, USA
| | | | | | | | | | | | | | | | | | | |
Collapse
|
45
|
Mohy-Ud-Din H, Lodge MA, Rahmim A. Quantitative myocardial perfusion PET parametric imaging at the voxel-level. Phys Med Biol 2015. [PMID: 26216052 DOI: 10.1088/0031-9155/60/15/6013] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Quantitative myocardial perfusion (MP) PET has the potential to enhance detection of early stages of atherosclerosis or microvascular dysfunction, characterization of flow-limiting effects of coronary artery disease (CAD), and identification of balanced reduction of flow due to multivessel stenosis. We aim to enable quantitative MP-PET at the individual voxel level, which has the potential to allow enhanced visualization and quantification of myocardial blood flow (MBF) and flow reserve (MFR) as computed from uptake parametric images. This framework is especially challenging for the (82)Rb radiotracer. The short half-life enables fast serial imaging and high patient throughput; yet, the acquired dynamic PET images suffer from high noise-levels introducing large variability in uptake parametric images and, therefore, in the estimates of MBF and MFR. Robust estimation requires substantial post-smoothing of noisy data, degrading valuable functional information of physiological and pathological importance. We present a feasible and robust approach to generate parametric images at the voxel-level that substantially reduces noise without significant loss of spatial resolution. The proposed methodology, denoted physiological clustering, makes use of the functional similarity of voxels to penalize deviation of voxel kinetics from physiological partners. The results were validated using extensive simulations (with transmural and non-transmural perfusion defects) and clinical studies. Compared to post-smoothing, physiological clustering depicted enhanced quantitative noise versus bias performance as well as superior recovery of perfusion defects (as quantified by CNR) with minimal increase in bias. Overall, parametric images obtained from the proposed methodology were robust in the presence of high-noise levels as manifested in the voxel time-activity-curves.
Collapse
Affiliation(s)
- Hassan Mohy-Ud-Din
- Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, MD 21218, USA. Department of Radiology and Radiological Science, Johns Hopkins University, Baltimore, MD 21287, USA. Department of Diagnostic Radiology, Yale University, New Haven, CT 06520, USA
| | | | | |
Collapse
|
46
|
A sinogram warping strategy for pre-reconstruction 4D PET optimization. Med Biol Eng Comput 2015; 54:535-46. [DOI: 10.1007/s11517-015-1339-y] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2014] [Accepted: 06/19/2015] [Indexed: 12/27/2022]
|
47
|
Lee KS, Kim TJ, Pratx G. Single-cell tracking with PET using a novel trajectory reconstruction algorithm. IEEE TRANSACTIONS ON MEDICAL IMAGING 2015; 34:994-1003. [PMID: 25423651 PMCID: PMC4392854 DOI: 10.1109/tmi.2014.2373351] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Virtually all biomedical applications of positron emission tomography (PET) use images to represent the distribution of a radiotracer. However, PET is increasingly used in cell tracking applications, for which the "imaging" paradigm may not be optimal. Here, we investigate an alternative approach, which consists in reconstructing the time-varying position of individual radiolabeled cells directly from PET measurements. As a proof of concept, we formulate a new algorithm for reconstructing the trajectory of one single moving cell directly from list-mode PET data. We model the trajectory as a 3-D B-spline function of the temporal variable and use nonlinear optimization to minimize the mean-square distance between the trajectory and the recorded list-mode coincidence events. Using Monte Carlo simulations (GATE), we show that this new algorithm can track a single source moving within a small-animal PET system with 3 mm accuracy provided that the activity of the cell [Bq] is greater than four times its velocity [mm/s]. The algorithm outperforms conventional ML-EM as well as the "minimum distance" method used for positron emission particle tracking (PEPT). The new method was also successfully validated using experimentally acquired PET data. In conclusion, we demonstrated the feasibility of a new method for tracking a single moving cell directly from PET list-mode data, at the whole-body level, for physiologically relevant activities and velocities.
Collapse
Affiliation(s)
- Keum Sil Lee
- Department of Radiology, Stanford University, CA 94305 USA
| | - Tae Jin Kim
- Department of Radiation Oncology, Stanford University, CA 94305 USA
| | - Guillem Pratx
- Department of Radiation Oncology, Stanford University, CA 94305 USA
| |
Collapse
|
48
|
Chen S, Liu H, Hu Z, Zhang H, Shi P, Chen Y. Simultaneous Reconstruction and Segmentation of Dynamic PET via Low-Rank and Sparse Matrix Decomposition. IEEE Trans Biomed Eng 2015; 62:1784-95. [PMID: 25706502 DOI: 10.1109/tbme.2015.2404296] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Although of great clinical value, accurate and robust reconstruction and segmentation of dynamic positron emission tomography (PET) images are great challenges due to low spatial resolution and high noise. In this paper, we propose a unified framework that exploits temporal correlations and variations within image sequences based on low-rank and sparse matrix decomposition. Thus, the two separate inverse problems, PET image reconstruction and segmentation, are accomplished in a simultaneous fashion. Considering low signal to noise ratio and piece-wise constant assumption of PET images, we also propose to regularize low-rank and sparse matrices with vectorial total variation norm. The resulting optimization problem is solved by augmented Lagrangian multiplier method with variable splitting. The effectiveness of proposed approach is validated on realistic Monte Carlo simulation datasets and the real patient data.
Collapse
|
49
|
Abstract
Image reconstruction from low-count positron emission tomography (PET) projection data is challenging because the inverse problem is ill-posed. Prior information can be used to improve image quality. Inspired by the kernel methods in machine learning, this paper proposes a kernel based method that models PET image intensity in each pixel as a function of a set of features obtained from prior information. The kernel-based image model is incorporated into the forward model of PET projection data and the coefficients can be readily estimated by the maximum likelihood (ML) or penalized likelihood image reconstruction. A kernelized expectation-maximization algorithm is presented to obtain the ML estimate. Computer simulations show that the proposed approach can achieve better bias versus variance trade-off and higher contrast recovery for dynamic PET image reconstruction than the conventional maximum likelihood method with and without post-reconstruction denoising. Compared with other regularization-based methods, the kernel method is easier to implement and provides better image quality for low-count data. Application of the proposed kernel method to a 4-D dynamic PET patient dataset showed promising results.
Collapse
|
50
|
Kazantsev D, M. Thompson W, R. B. Lionheart W, Van Eyndhoven G, P. Kaestner A, J. Dobson K, J. Withers P, D. Lee P. 4D-CT reconstruction with unified spatial-temporal patch-based regularization. ACTA ACUST UNITED AC 2015. [DOI: 10.3934/ipi.2015.9.447] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|