Chi Liu
Yale University
Network
Latest external collaboration on country level. Dive into details by clicking on the dots.
Publication
Featured researches published by Chi Liu.
Physics in Medicine and Biology | 2009
Chi Liu; Larry Pierce; Adam M. Alessio; Paul E. Kinahan
Our aim is to investigate the impact of respiratory motion on tumor quantification and delineation in static PET/CT imaging using a population of patient respiratory traces. A total of 1295 respiratory traces acquired during whole body PET/CT imaging were classified into three types according to the qualitative shape of their signal histograms. Each trace was scaled to three diaphragm motion amplitudes (6 mm, 11 mm and 16 mm) to drive a whole body PET/CT computer simulation that was validated with a physical phantom experiment. Three lung lesions and one liver lesion were simulated with diameters of 1 cm and 2 cm. PET data were reconstructed using the OS-EM algorithm with attenuation correction using CT images at the end-expiration phase and respiratory-averaged CT. The errors of the lesion maximum standardized uptake values (SUV(max)) and lesion volumes between motion-free and motion-blurred PET/CT images were measured and analyzed. For respiration with 11 mm diaphragm motion and larger quiescent period fraction, respiratory motion can cause a mean lesion SUV(max) underestimation of 28% and a mean lesion volume overestimation of 130% in PET/CT images with 1 cm lesions. The errors of lesion SUV(max) and volume are larger for patient traces with larger motion amplitudes. Smaller lesions are more sensitive to respiratory motion than larger lesions for the same motion amplitude. Patient respiratory traces with relatively larger quiescent period fraction yield results less subject to respiratory motion than traces with long-term amplitude variability. Mismatched attenuation correction due to respiratory motion can cause SUV(max) overestimation for lesions in the lower lung region close to the liver dome. Using respiratory-averaged CT for attenuation correction yields smaller mismatch errors than those using end-expiration CT. Respiratory motion can have a significant impact on static oncological PET/CT imaging where SUV and/or volume measurements are important. The impact is highly dependent upon motion amplitude, lesion location and size, attenuation map and respiratory pattern. To overcome the motion effect, motion compensation techniques may be necessary in clinical practice to improve the tumor quantification for determining the response to therapy or for radiation treatment planning.
Medical Physics | 2010
Chi Liu; Adam M. Alessio; Larry Pierce; Kris Thielemans; Scott D. Wollenweber; Alexander Ganin; Paul E. Kinahan
PURPOSE To minimize respiratory motion artifacts, this work proposes quiescent period gating (QPG) methods that extract PET data from the end-expiration quiescent period and form a single PET frame with reduced motion and improved signal-to-noise properties. METHODS Two QPG methods are proposed andevaluated. Histogram-based quiescent period gating (H-QPG) extracts a fraction of PET data determined by a window of the respiratory displacement signal histogram. Cycle-based quiescent period gating (C-QPG) extracts data with a respiratory displacement signal below a specified threshold of the maximum amplitude of each individual respiratory cycle. Performances of both QPG methods were compared to ungated and five-bin phase-gated images across 21 FDG-PET/CT patient data sets containing 31 thorax and abdomen lesions as well as with computer simulations driven by 1295 different patient respiratory traces. Image quality was evaluated in terms of the lesion SUV(max) and the fraction of counts included in each gate as a surrogate for image noise. RESULTS For all the gating methods, image noise artifactually increases SUV(max) when the fraction of counts included in each gate is less than 50%. While simulation data show that H-QPG is superior to C-QPG, the H-QPG and C-QPG methods lead to similar quantification-noise tradeoffs in patient data. Compared to ungated images, both QPG methods yield significantly higher lesion SUV(max). Compared to five-bin phase gating, the QPG methods yield significantly larger fraction of counts with similar SUV(max) improvement. Both QPG methods result in increased lesion SUV(max) for patients whose lesions have longer quiescent periods. CONCLUSIONS Compared to ungated and phase-gated images, the QPG methods lead to images with less motion blurring and an improved compromise between SUV(max) and fraction of counts. The QPG methods for respiratory motion compensation could effectively improve tumor quantification with minimal noise increase.
Medical Physics | 2011
Lawrence R. MacDonald; Steve Kohlmyer; Chi Liu; Thomas K. Lewellen; Paul E. Kinahan
PURPOSE The goal of this work was to investigate the effects of MRI surface coils on attenuation-corrected PET emission data. The authors studied the cases where either an MRI or a CT scan would be used to provide PET attenuation correction (AC). Combined MR/PET scanners that use the MRI for PET AC (MR-AC) face the challenge of absent surface coils in MR images and thus cannot directly account for attenuation in the coils. Combining MR and PET images could be achieved by transporting the subject on a stereotactically registered table between independent MRI and PET scanners. In this case, conventional PET CT-AC methods could be used. A challenge here is that high atomic number materials within MR coils cause artifacts in CT images and CT based AC is typically not validated for coil materials. METHODS The authors evaluated PET artifacts when MR coils were absent from AC data (MR-AC), or when coil attenuation was measured by CT scanning (CT-AC). They scanned PET phantoms with MR surface coils on a clinical PET/CT system and used CT-AC to reconstruct PET data. The authors then omitted the coil from the CT-AC image to mimic the MR-AC scenario. Images were acquired using cylinder and anthropomorphic phantoms. They evaluated and compared the following five scenarios: (1) A uniform cylinder phantom and head coil scanned and reconstructed using CT-AC; (2) similar emission data (with head coil present) were reconstructed without the head coil in the AC data; (3) the same cylinder scanned without the head coil present (reference scan); (4) a PET torso phantom with a full MR torso coil present in both PET and CT; (5) only half of the separable torso coil present in the PET/CT acquisition. The authors also performed analytic simulations of the first three scenarios. RESULTS Streak artifacts were present in CT images containing MR surface coils due to metal components. These artifacts persisted after the CT images were converted for PET AC. The artifacts were significantly reduced when half of the separable coil was removed during the scan. CT scans tended to over-estimate the linear attenuation coefficient (micro) of the metal components when using conventional methods for converting from CT number to micro(511 keV). Artifacts were visible outside the phantom in some of the PET emission images, corresponding to the MRI coil geometry. However, only subtle artifacts were apparent in the emission images inside the phantoms. On the other hand, the PET emission image quantitative accuracy was significantly affected: the activity was underestimated by 19% when AC did not include the head coil, and overestimated by 28% when the CT-AC included the head coil. CONCLUSIONS The presence of MR coils during PET or PET/CT scanning can cause subtle artifacts and potentially important quantification errors. Alternative CT techniques that mitigate artifacts should be used to improve AC accuracy. When possible, removing segments of an MR coil prior to the PET/CT exam is recommended. Further, MR coils could be redesigned to reduce artifacts by rearranging placement of the most attenuating materials.
Medical Physics | 2011
Chi Liu; Adam M. Alessio; Paul E. Kinahan
PURPOSE We present a method to correct respiratory motion blurring in PET/CT imaging using internal-external (INTEX) motion correlation. The internal motion of a known tumor is derived from respiratory-gated PET images; this internal motion is then correlated with external respiratory signals to determine the complete information of tumor motion during the scan. METHODS For each PET/CT data, PET listmode data were phase-gated into five bins and reconstructed. The centroid of a targeted tumor in each bin was determined and correlated with the corresponding mean displacement of externally monitored respiratory motion signal. Based on this correlation, the external motion signal was converted into internal tumor motion information in the superior-inferior direction. Then, the PET listmode data were binned sequentially to multiple 1-s sinograms. According to the converted internal tumor motion signal, each 1-s sinogram was registered to a reference frame, which best matched the helical CT attenuation map based on consistency conditions. The registered sinograms were summed and reconstructed to form an image, corrected for the motion of the specific tumor. In this study, the proposed INTEX method was evaluated with phantom and patient studies in terms of tracer concentration and volume. RESULTS The INTEX method effectively recovered the tracer concentration to the level of the stationary scan data in the phantom experiment. In the patient study, the INTEX method yielded a (17 +/- 22)% tumor volume decrease and a (10 +/- 10)% tumor SUVmax increase compared to non-gated images. CONCLUSIONS The proposed INTEX method reduces respiratory motion degradation of PET tumor quantification and delineation in an effective manner. This can be used to improve the assessment of response to therapy for a known tumor by minimizing residual motion and matching the attenuation correction, without increasing image noise.
The Journal of Nuclear Medicine | 2016
Yunhan Yu; Chung Chan; Tianyu Ma; Yaqiang Liu; Jean Dominique Gallezot; Mika Naganawa; Olivia J. Kelada; Mary Germino; Albert J. Sinusas; Richard E. Carson; Chi Liu
Existing respiratory motion-correction methods are applied only to static PET imaging. We have previously developed an event-by-event respiratory motion-correction method with correlations between internal organ motion and external respiratory signals (INTEX). This method is uniquely appropriate for dynamic imaging because it corrects motion for each time point. In this study, we applied INTEX to human dynamic PET studies with various tracers and investigated the impact on kinetic parameter estimation. Methods: The use of 3 tracers—a myocardial perfusion tracer, 82Rb (n = 7); a pancreatic β-cell tracer, 18F-FP(+)DTBZ (n = 4); and a tumor hypoxia tracer, 18F-fluoromisonidazole (18F-FMISO) (n = 1)—was investigated in a study of 12 human subjects. Both rest and stress studies were performed for 82Rb. The Anzai belt system was used to record respiratory motion. Three-dimensional internal organ motion in high temporal resolution was calculated by INTEX to guide event-by-event respiratory motion correction of target organs in each dynamic frame. Time–activity curves of regions of interest drawn based on end-expiration PET images were obtained. For 82Rb studies, K1 was obtained with a 1-tissue model using a left-ventricle input function. Rest–stress myocardial blood flow (MBF) and coronary flow reserve (CFR) were determined. For 18F-FP(+)DTBZ studies, the total volume of distribution was estimated with arterial input functions using the multilinear analysis 1 method. For the 18F-FMISO study, the net uptake rate Ki was obtained with a 2-tissue irreversible model using a left-ventricle input function. All parameters were compared with the values derived without motion correction. Results: With INTEX, K1 and MBF increased by 10% ± 12% and 15% ± 19%, respectively, for 82Rb stress studies. CFR increased by 19% ± 21%. For studies with motion amplitudes greater than 8 mm (n = 3), K1, MBF, and CFR increased by 20% ± 12%, 30% ± 20%, and 34% ± 23%, respectively. For 82Rb rest studies, INTEX had minimal effect on parameter estimation. The total volume of distribution of 18F-FP(+)DTBZ and Ki of 18F-FMISO increased by 17% ± 6% and 20%, respectively. Conclusion: Respiratory motion can have a substantial impact on dynamic PET in the thorax and abdomen. The INTEX method using continuous external motion data substantially changed parameters in kinetic modeling. More accurate estimation is expected with INTEX.
Medical Physics | 2015
Peng Fan; Brian F. Hutton; Maria Holstensson; Michael Ljungberg; P. Hendrik Pretorius; Rameshwar Prasad; Tianyu Ma; Yaqiang Liu; Shi Wang; Stephanie Thorn; Mitchel R. Stacy; Albert J. Sinusas; Chi Liu
PURPOSE The energy spectrum for a cadmium zinc telluride (CZT) detector has a low energy tail due to incomplete charge collection and intercrystal scattering. Due to these solid-state detector effects, scatter would be overestimated if the conventional triple-energy window (TEW) method is used for scatter and crosstalk corrections in CZT-based imaging systems. The objective of this work is to develop a scatter and crosstalk correction method for (99m)Tc/(123)I dual-radionuclide imaging for a CZT-based dedicated cardiac SPECT system with pinhole collimators (GE Discovery NM 530c/570c). METHODS A tailing model was developed to account for the low energy tail effects of the CZT detector. The parameters of the model were obtained using (99m)Tc and (123)I point source measurements. A scatter model was defined to characterize the relationship between down-scatter and self-scatter projections. The parameters for this model were obtained from Monte Carlo simulation using SIMIND. The tailing and scatter models were further incorporated into a projection count model, and the primary and self-scatter projections of each radionuclide were determined with a maximum likelihood expectation maximization (MLEM) iterative estimation approach. The extracted scatter and crosstalk projections were then incorporated into MLEM image reconstruction as an additive term in forward projection to obtain scatter- and crosstalk-corrected images. The proposed method was validated using Monte Carlo simulation, line source experiment, anthropomorphic torso phantom studies, and patient studies. The performance of the proposed method was also compared to that obtained with the conventional TEW method. RESULTS Monte Carlo simulations and line source experiment demonstrated that the TEW method overestimated scatter while their proposed method provided more accurate scatter estimation by considering the low energy tail effect. In the phantom study, improved defect contrasts were observed with both correction methods compared to no correction, especially for the images of (99m)Tc in dual-radionuclide imaging where there is heavy contamination from (123)I. In this case, the nontransmural defect contrast was improved from 0.39 to 0.47 with the TEW method and to 0.51 with their proposed method and the transmural defect contrast was improved from 0.62 to 0.74 with the TEW method and to 0.73 with their proposed method. In the patient study, the proposed method provided higher myocardium-to-blood pool contrast than that of the TEW method. Similar to the phantom experiment, the improvement was the most substantial for the images of (99m)Tc in dual-radionuclide imaging. In this case, the myocardium-to-blood pool ratio was improved from 7.0 to 38.3 with the TEW method and to 63.6 with their proposed method. Compared to the TEW method, the proposed method also provided higher count levels in the reconstructed images in both phantom and patient studies, indicating reduced overestimation of scatter. Using the proposed method, consistent reconstruction results were obtained for both single-radionuclide data with scatter correction and dual-radionuclide data with scatter and crosstalk corrections, in both phantom and human studies. CONCLUSIONS The authors demonstrate that the TEW method leads to overestimation in scatter and crosstalk for the CZT-based imaging system while the proposed scatter and crosstalk correction method can provide more accurate self-scatter and down-scatter estimations for quantitative single-radionuclide and dual-radionuclide imaging.
EJNMMI research | 2016
Mary Germino; Jim Ropchan; Tim Mulnix; Kathryn Fontaine; Nabeel Nabulsi; Eric Ackah; Herman H Feringa; Albert J. Sinusas; Chi Liu; Richard E. Carson
BackgroundWe quantified myocardial blood flow with 82Rb PET using parameters of the generalized Renkin-Crone model estimated from 82Rb and 15O-water images reconstructed with time-of-flight and point spread function modeling. Previous estimates of rubidium extraction have used older-generation scanners without time-of-flight or point spread function modeling. We validated image-derived input functions with continuously collected arterial samples.MethodsNine healthy subjects were scanned at rest and under pharmacological stress on the Siemens Biograph mCT with 82Rb and 15O-water PET, undergoing arterial blood sampling with each scan. Image-derived input functions were estimated from the left ventricle cavity and corrected with tracer-specific population-based scale factors determined from arterial data. Kinetic parametric images were generated from the dynamic PET images by fitting the one-tissue compartment model to each voxel’s time activity curve. Mean myocardial blood flow was determined from each subject’s 15O-water k2 images. The parameters of the generalized Renkin-Crone model were estimated from these water-based flows and mean myocardial 82Rb K1 estimates.ResultsImage-derived input functions showed improved agreement with arterial measurements after a scale correction. The Renkin-Crone model fit (a = 0.77, b = 0.39) was similar to those previously published, though b was lower.ConclusionsWe have presented parameter estimates for the generalized Renkin-Crone model of extraction for 82Rb PET using human 82Rb and 15O-water PET from high-resolution images using a state-of-the-art time-of-flight-capable scanner. These results provide a state-of-the-art methodology for myocardial blood flow measurement with 82Rb PET.
Medical Physics | 2013
Chung Chan; Xiao Jin; Edward K. Fung; Mika Naganawa; Tim Mulnix; Richard E. Carson; Chi Liu
PURPOSE Respiratory motion during PET∕CT imaging can cause substantial image blurring and underestimation of tracer concentration for both static and dynamic studies. In this study, the authors developed an event-by-event respiratory motion correction method that used three-dimensional internal-one-dimensional external motion correlation (INTEX3D) in listmode reconstruction. The authors aim to fully correct for organ/tumor-specific rigid motion caused by respiration using all detected events to eliminate both intraframe and interframe motion, and investigate the quantitative improvement in static and dynamic imaging. METHODS The positional translation of an internal organ or tumor during respiration was first determined from the reconstructions of multiple phase-gated images. A level set (active contour) method was used to segment the targeted internal organs/tumors whose centroids were determined. The mean displacement of the external respiratory signal acquired by the Anzai system that corresponded to each phase-gated frame was determined. Three linear correlations between the 1D Anzai mean displacements and the 3D centroids of the internal organ/tumor were established. The 3D internal motion signal with high temporal resolution was then generated by applying each of the three correlation functions to the entire Anzai trace (40 Hz) to guide event-by-event motion correction in listmode reconstruction. The reference location was determined as the location where CT images were acquired to facilitate phase-matched attenuation correction and anatomical-based postfiltering. The proposed method was evaluated with a NEMA phantom driven by a QUASAR respiratory motion platform, and human studies with two tracers: pancreatic beta cell tracer [(18)F]FP(+)DTBZ and tumor hypoxia tracer [(18)F]fluoromisonidazole (FMISO). An anatomical-based postreconstruction filter was applied to the motion-corrected images to reduce noise while preserving quantitative accuracy and organ boundaries in the patient studies. RESULTS The INTEX3D method yielded an increase of 5%-9% and 32%-40% in contrast recovery coefficient on the hot spheres in the NEMA phantom, compared to the reconstructions with only 1D motion correction (INTEX1D) and no motion correction, respectively. The proposed method also increased the mean activities of the pancreas and kidney by 9.3% and 11.2%, respectively, across three subjects in the FPDTBZ studies, and the average lesion-to-blood ratio by 20% across three lesions in the FMISO study, compared to the reconstructions without motion correction. In addition, the proposed method reduced intragate motion as compared to phase-gated images. The application of the anatomical-based postreconstruction filter further reduced noise in the background by >50% compared to reconstructions without postfiltering, while preserving quantitative accuracy and organ boundaries. Finally, the measurements of the time-activity curves from a subject with FPDTBZ showed that INTEX3D yielded 18% and 11% maximum increases in tracer concentration in the pancreas and kidney cortex, respectively. CONCLUSIONS These results suggest that the proposed method can effectively compensate for both intragate and intergate respiratory motion while preserving all the counts, and is applicable to dynamic studies.
The Journal of Nuclear Medicine | 2014
Chi Liu; Albert J. Sinusas
Quantification of absolute myocardial flow and flow reserve using dynamic PET imaging has now been more widely applied clinically to detect microvascular disease, improve diagnostic accuracy for detection of coronary artery disease, and improve prognostic implications of stress perfusion imaging in patients with heart disease (1). PET imaging offers some distinct advantages over SPECT imaging for the assessment of absolute myocardial flow, including higher count sensitivity and higher spatial resolution than SPECT; well-established approaches for attenuation correction for all PET imaging systems; and approved radiotracers with favorable kinetics for the assessment of absolute myocardial perfusion. The clinical quantification of myocardial flow with PET has been facilitated by the availability of several Food and Drug Administration–approved quantitative software programs that are commercially available. These programs have been shown to be highly reproducible and comparable for estimation of myocardial flow reserve in patients, although some differences exist among these programs for estimation of absolute myocardial flow (2). Although the availability of hybrid PET/CT scanners is growing,
Physics in Medicine and Biology | 2017
Silin Ren; Xiao Jin; Chung Chan; Yiqiang Jian; Tim Mulnix; Chi Liu; Richard E. Carson
Data-driven respiratory gating techniques were developed to correct for respiratory motion in PET studies, without the help of external motion tracking systems. Due to the greatly increased image noise in gated reconstructions, it is desirable to develop a data-driven event-by-event respiratory motion correction method. In this study, using the Centroid-of-distribution (COD) algorithm, we established a data-driven event-by-event respiratory motion correction technique using TOF PET list-mode data, and investigated its performance by comparing with an external system-based correction method. Ten human scans with the pancreatic β-cell tracer 18F-FP-(+)-DTBZ were employed. Data-driven respiratory motions in superior-inferior (SI) and anterior-posterior (AP) directions were first determined by computing the centroid of all radioactive events during each short time frame with further processing. The Anzai belt system was employed to record respiratory motion in all studies. COD traces in both SI and AP directions were first compared with Anzai traces by computing the Pearson correlation coefficients. Then, respiratory gated reconstructions based on either COD or Anzai traces were performed to evaluate their relative performance in capturing respiratory motion. Finally, based on correlations of displacements of organ locations in all directions and COD information, continuous 3D internal organ motion in SI and AP directions was calculated based on COD traces to guide event-by-event respiratory motion correction in the MOLAR reconstruction framework. Continuous respiratory correction results based on COD were compared with that based on Anzai, and without motion correction. Data-driven COD traces showed a good correlation with Anzai in both SI and AP directions for the majority of studies, with correlation coefficients ranging from 63% to 89%. Based on the determined respiratory displacements of pancreas between end-expiration and end-inspiration from gated reconstructions, there was no significant difference between COD-based and Anzai-based methods. Finally, data-driven COD-based event-by-event respiratory motion correction yielded comparable results to that based on Anzai respiratory traces, in terms of contrast recovery and reduced motion-induced blur. Data-driven event-by-event respiratory motion correction using COD showed significant image quality improvement compared with reconstructions with no motion correction, and gave comparable results to the Anzai-based method.