Low dose dynamic CT myocardial perfusion imaging using a statistical iterative reconstruction method.

PURPOSE Dynamic CT myocardial perfusion imaging has the potential to provide both functional and anatomical information regarding coronary artery stenosis. However, radiation dose can be potentially high due to repeated scanning of the same region. The purpose of this study is to investigate the use of statistical iterative reconstruction to improve parametric maps of myocardial perfusion derived from a low tube current dynamic CT acquisition. METHODS Four pigs underwent high (500 mA) and low (25 mA) dose dynamic CT myocardial perfusion scans with and without coronary occlusion. To delineate the affected myocardial territory, an N-13 ammonia PET perfusion scan was performed for each animal in each occlusion state. Filtered backprojection (FBP) reconstruction was first applied to all CT data sets. Then, a statistical iterative reconstruction (SIR) method was applied to data sets acquired at low dose. Image voxel noise was matched between the low dose SIR and high dose FBP reconstructions. CT perfusion maps were compared among the low dose FBP, low dose SIR and high dose FBP reconstructions. Numerical simulations of a dynamic CT scan at high and low dose (20:1 ratio) were performed to quantitatively evaluate SIR and FBP performance in terms of flow map accuracy, precision, dose efficiency, and spatial resolution. RESULTS Forin vivo studies, the 500 mA FBP maps gave -88.4%, -96.0%, -76.7%, and -65.8% flow change in the occluded anterior region compared to the open-coronary scans (four animals). The percent changes in the 25 mA SIR maps were in good agreement, measuring -94.7%, -81.6%, -84.0%, and -72.2%. The 25 mA FBP maps gave unreliable flow measurements due to streaks caused by photon starvation (percent changes of +137.4%, +71.0%, -11.8%, and -3.5%). Agreement between 25 mA SIR and 500 mA FBP global flow was -9.7%, 8.8%, -3.1%, and 26.4%. The average variability of flow measurements in a nonoccluded region was 16.3%, 24.1%, and 937.9% for the 500 mA FBP, 25 mA SIR, and 25 mA FBP, respectively. In numerical simulations, SIR mitigated streak artifacts in the low dose data and yielded flow maps with mean error <7% and standard deviation <9% of mean, for 30 × 30 pixel ROIs (12.9 × 12.9 mm(2)). In comparison, low dose FBP flow errors were -38% to +258%, and standard deviation was 6%-93%. Additionally, low dose SIR achieved 4.6 times improvement in flow map CNR(2) per unit input dose compared to low dose FBP. CONCLUSIONS SIR reconstruction can reduce image noise and mitigate streaking artifacts caused by photon starvation in dynamic CT myocardial perfusion data sets acquired at low dose (low tube current), and improve perfusion map quality in comparison to FBP reconstruction at the same dose.

[1]  Max A. Viergever,et al.  Noise Reduction in Computed Tomography Scans Using 3-D Anisotropic Hybrid Diffusion With Continuous Switch , 2009, IEEE Transactions on Medical Imaging.

[2]  Jie Tang,et al.  Prior image constrained compressed sensing: implementation and performance evaluation. , 2011, Medical physics.

[3]  U. Schoepf,et al.  Adenosine-stress dynamic myocardial perfusion imaging with second-generation dual-source CT: comparison with conventional catheter coronary angiography and SPECT nuclear myocardial perfusion imaging. , 2012, AJR. American journal of roentgenology.

[4]  Jeffrey A. Fessler.  Penalized weighted least-squares image reconstruction for positron emission tomography , 1994, IEEE Trans. Medical Imaging.

[5]  C. Rochitte,et al.  Additional value of dipyridamole stress myocardial perfusion by 64-row computed tomography in patients with coronary stents. , 2011, Journal of cardiovascular computed tomography.

[6]  Guang-Hong Chen,et al.  Noise spatial nonuniformity and the impact of statistical image reconstruction in CT myocardial perfusion imaging. , 2012, Medical physics.

[7]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .

[8]  S. Schoenberg,et al.  Radiation dose at coronary CT angiography: second-generation dual-source CT versus single-source 64-MDCT and first-generation dual-source CT. , 2011, AJR. American journal of roentgenology.

[9]  Cynthia H McCollough,et al.  A strategy to decrease partial scan reconstruction artifacts in myocardial perfusion CT: phantom and in vivo evaluation. , 2011, Medical physics.

[10]  Yinghua Tao,et al.  Investigation of the potential causes of partial scan artifacts in dynamic CT myocardial perfusion imaging , 2014, Medical Imaging.

[11]  T-Y Lee,et al.  Theoretic Basis and Technical Implementations of CT Perfusion in Acute Ischemic Stroke, Part 1: Theoretic Basis , 2009, American Journal of Neuroradiology.

[12]  D. Lipkin,et al.  Radiation dose of CT coronary angiography in clinical practice: objective evaluation of strategies for dose optimization. , 2012, European journal of radiology.

[13]  D. Berman,et al.  ACCF/ASNC/ACR/AHA/ASE/SCCT/SCMR/SNM 2009 Appropriate Use Criteria for Cardiac Radionuclide Imaging: A Report of the American College of Cardiology Foundation Appropriate Use Criteria Task Force, the American Society of Nuclear Cardiology, the American College of Radiology, the American Heart Associa , 2009, Journal of the American College of Cardiology.

[14]  R. Cury,et al.  Incremental value of adenosine-induced stress myocardial perfusion imaging with dual-source CT at cardiac CT angiography. , 2010, Radiology.

[15]  D. Paech,et al.  A systematic review of the clinical effectiveness of 64-slice or higher computed tomography angiography as an alternative to invasive coronary angiography in the investigation of suspected coronary artery disease , 2011, BMC cardiovascular disorders.

[16]  E. Antman,et al.  ACC/AHA/ASNC guidelines for the clinical use of cardiac radionuclide imaging--executive summary: a report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines (ACC/AHA/ASNC Committee to Revise the 1995 Guidelines for the Clinical Use of Cardiac Radionucl , 2003, Journal of the American College of Cardiology.

[17]  B. van Ginneken,et al.  TIPS bilateral noise reduction in 4D CT perfusion scans produces high-quality cerebral blood flow maps , 2011, Physics in medicine and biology.

[18]  E. Topol,et al.  Our preoccupation with coronary luminology. The dissociation between clinical and angiographic findings in ischemic heart disease. , 1995, Circulation.

[19]  R. Cury,et al.  Adenosine-induced stress myocardial perfusion imaging using dual-source cardiac computed tomography. , 2009, Journal of the American College of Cardiology.

[20]  Jonathan G Goldin,et al.  Assessment of Coronary Artery Disease by Cardiac Computed Tomography: A Scientific Statement From the American Heart Association Committee on Cardiovascular Imaging and Intervention, Council on Cardiovascular Radiology and Intervention, and Committee on Cardiac Imaging, Council on Clinical Cardiolog , 2006, Circulation.

[21]  W. Kalender,et al.  Quantitative Whole Heart Stress Perfusion CT Imaging as Noninvasive Assessment of Hemodynamics in Coronary Artery Stenosis: Preliminary Animal Experience , 2010, Investigative radiology.

[22]  Rainer Raupach,et al.  Wavelet Based Noise Reduction in CT-Images Using Correlation Analysis , 2008, IEEE Transactions on Medical Imaging.

[23]  J. Bartunek,et al.  Quantitative coronary angiography in predicting functional significance of stenoses in an unselected patient cohort. , 1995, Journal of the American College of Cardiology.

[24]  C. Fink,et al.  Adenosine-stress dynamic real-time myocardial perfusion CT and adenosine-stress first-pass dual-energy myocardial perfusion CT for the assessment of acute chest pain: initial results. , 2012, European journal of radiology.

[25]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[26]  Jean-Baptiste Thibault,et al.  A three-dimensional statistical approach to improved image quality for multislice helical CT. , 2007, Medical physics.

[27]  Yinghua Tao,et al.  Radiation dose reduction in computed tomography (CT) using a new implementation of wavelet denoising in low tube current acquisitions , 2011, Medical Imaging.

[28]  Mathias Prokop,et al.  Diagnostic accuracy of 64-slice computed tomography coronary angiography: a prospective, multicenter, multivendor study. , 2008, Journal of the American College of Cardiology.

[29]  Ge Wang,et al.  Reduction of half-scan shading artifact based on full-scan correction. , 2006, Academic radiology.

[30]  U. Schoepf,et al.  Adenosine-Stress Dynamic Myocardial CT Perfusion Imaging: Initial Clinical Experience , 2010, Investigative radiology.

[31]  L. Shepp,et al.  A Statistical Model for Positron Emission Tomography , 1985 .

[32]  Roberto Manduchi,et al.  Bilateral filtering for gray and color images , 1998, Sixth International Conference on Computer Vision (IEEE Cat. No.98CH36271).

[33]  R. Cury,et al.  Dipyridamole stress and rest transmural myocardial perfusion ratio evaluation by 64 detector-row computed tomography. , 2011, Journal of cardiovascular computed tomography.

[34]  S. Reeder,et al.  Reduction of image noise in low tube current dynamic CT myocardial perfusion imaging using HYPR processing: a time-attenuation curve analysis. , 2012, Medical physics.

[35]  Ken D. Sauer,et al.  A unified approach to statistical tomography using coordinate descent optimization , 1996, IEEE Trans. Image Process..

[36]  Kia-Chong Chua,et al.  Stress and rest dynamic myocardial perfusion imaging by evaluation of complete time-attenuation curves with dual-source CT. , 2010, JACC. Cardiovascular imaging.

[37]  Marc Kachelriess,et al.  Partial scan artifact reduction (PSAR) for the assessment of cardiac perfusion in dynamic phase-correlated CT , 2009, 2008 IEEE Nuclear Science Symposium Conference Record.

[38]  장혁재.  Adenosine stress 64- and 256-row detector computed tomography angiography and perfusion imaging: a pilot study evaluating the transmural extent of perfusion abnormalities to predict atherosclerosis causing myocardial ischemia. , 2009 .

[39]  M. Reiser,et al.  Detection of hemodynamically significant coronary artery stenosis: incremental diagnostic value of dynamic CT-based myocardial perfusion imaging. , 2011, Radiology.

[40]  Ken D. Sauer,et al.  A local update strategy for iterative reconstruction from projections , 1993, IEEE Trans. Signal Process..

[41]  J. Lima,et al.  Multidetector computed tomography myocardial perfusion imaging during adenosine stress. , 2006, Journal of the American College of Cardiology.

[42]  Dennis M. Healy,et al.  Wavelet transform domain filters: a spatially selective noise filtration technique , 1994, IEEE Trans. Image Process..

[43]  Kim Mouridsen,et al.  Comparison of 10 Perfusion MRI Parameters in 97 Sub-6-Hour Stroke Patients Using Voxel-Based Receiver Operating Characteristics Analysis , 2009, Stroke.

[44]  A. Arai,et al.  Submillisievert median radiation dose for coronary angiography with a second-generation 320-detector row CT scanner in 107 consecutive patients. , 2013, Radiology.