Positron Emission tomography (PET) is a major molecular imaging tool in oncology, with applications ranging from diagnosis and staging to patient management. Despite the broad use of PET in the clinical environment, there is no quantitative PET imaging method available for routine clinical practice. The currently used static scan can provide a semi-quantitative measurement, standardized uptake value (SUV), for a whole body scan. However, it completely ignores the dynamic nature of radiopharmaceutical kinetics. The popular semi-quantitative dual time point method can approximate the kinetic differences at two time points by comparing activities but usually requires an extended waiting time for the second scan. The multiple time point method can calculate the net influx rate but still requires long scan duration and makes a whole body scan infeasible. The challenge of a quantitative whole body dynamic PET scan lies in how to estimate the quantitative functional values, such as net flux rate, using data from a short acquisition period, and how to accelerate the computation to make it practical in a clinical setting. We address this challenge by developing and optimizing a novel data analysis method and implementing it using a high performance computing tool. We take advantage of the linearity of Patlak graphic analysis to model the tracer activity in each voxel as a linear combination of the blood input function and its integral, weighted by the Patlak parameters including net influx rate. In addition, we derive a simplified model of the blood input function, based on the same assumptions used to derive Patlak parameters from the kinetic compartment model. We then estimate the Patlak parameters and the parameters in the blood input function in a penalized maximum likelihood estimation framework using the list mode data and its associated inhomogeneous Poisson statistical model. We also theoretically analyze the performance of our Patlak estimator in terms of noise, resolution and signal-to-noise ratio (SNR), and use the results to guide us in optimizing the scan duration and any movement of imaging bed to achieve the best SNR. The advanced estimation algorithm, along with an accurate imaging system model, can robustly compute the net influx rate using the list mode data in a short acquisition without a measured blood input function, and make whole body dynamic scans practical. Our algorithms will be implemented on an Nvidia Tesla GPU (graphics processing unit) based workstation, a new computing tool that provides computational power previously available only on a mini super computer. We will further accelerate our algorithms using a combination of efficient representation of the list mode data and the system matrix. We will evaluate the performance of the proposed method and compare it with SUV, the dual time point method, and the traditional Patlak method, using simulated and clinical data. We will use a range of performance metrics, including region of interest (ROI) bias, ROI variance, lesion detectability, and computer and human observers. This project will eventually provide a quantitative dynamic whole body PET imaging protocol that can potentially improve the sensitivity and specificity of PET imaging in oncology.
Positron Emission Tomography (PET) has been widely used in cancer diagnosis, staging, treatment planning, management and evaluation. However, its potential is not yet fully realized, in part because we are not able to take full advantage of the dynamic information that can be collected by the PET scanner. In this project we will develop a new approach to the acquisition and analysis of PET data that will allow us for the first time to scan the whole body of the patient and produce quantitative estimates of PET tracer uptake from dynamically acquired data. These measures may be more sensitive indicators of the presence and metabolic activity of tumors, so that their use would lead to improved detection, staging and monitoring of primary and metastatic tumors.
|Wang, Mengdie; Guo, Ning; Hu, Guangshu et al. (2016) A novel approach to assess the treatment response using Gaussian random field in PET. Med Phys 43:833-42|
|Hu, Chenhui; Sepulcre, Jorge; Johnson, Keith A et al. (2016) Matched signal detection on graphs: Theory and application to brain imaging data classification. Neuroimage 125:587-600|
|Lorsakul, Auranuch; Fakhri, Georges El; Worstell, William et al. (2016) Numerical observer for atherosclerotic plaque classification in spectral computed tomography. J Med Imaging (Bellingham) 3:035501|
|Hu, Chenhui; Cheng, Lin; Sepulcre, Jorge et al. (2015) A spectral graph regression model for learning brain connectivity of Alzheimer's disease. PLoS One 10:e0128136|
|Dutta, Joyita; Huang, Chuan; Li, Quanzheng et al. (2015) Pulmonary imaging using respiratory motion compensated simultaneous PET/MR. Med Phys 42:4227-40|
|Kim, Kyungsang; Ye, Jong Chul; Worstell, William et al. (2015) Sparse-view spectral CT reconstruction using spectral patch-based low-rank penalty. IEEE Trans Med Imaging 34:748-60|
|Zhu, Wentao; Li, Quanzheng; Bai, Bing et al. (2014) Patlak image estimation from dual time-point list-mode PET data. IEEE Trans Med Imaging 33:913-24|
|Lin, Yanguang; Haldar, Justin P; Li, Quanzheng et al. (2014) Sparsity Constrained Mixture Modeling for the Estimation of Kinetic Parameters in Dynamic PET. IEEE Trans Med Imaging 33:173-85|
|Lorsakul, Auranuch; Li, Quanzheng; Trott, Cathryn M et al. (2014) 4D numerical observer for lesion detection in respiratory-gated PET. Med Phys 41:102504|
|Bai, Bing; Lin, Yanguang; Zhu, Wentao et al. (2014) MAP reconstruction for Fourier rebinned TOF-PET data. Phys Med Biol 59:925-49|
Showing the most recent 10 out of 20 publications