Nonconvex prior image constrained compressed sensing (NCPICCS): Theory and simulations on perfusion CT

被引:66
作者
Ramirez-Giraldo, J. C. [1 ,2 ]
Trzasko, J. [3 ]
Leng, S. [1 ]
Yu, L. [1 ]
Manduca, A. [2 ,3 ]
McCollough, C. H. [1 ]
机构
[1] Mayo Clin, CT Clin Innovat Ctr, Dept Radiol, Rochester, MN 55905 USA
[2] Mayo Clin, Dept Physiol & Biomed Engn, Rochester, MN 55905 USA
[3] Mayo Clin, Ctr Adv Imaging Res, Rochester, MN 55905 USA
关键词
image reconstruction; dynamic CT; perfusion CT; compressed sensing; radiation dose reduction; BEAM COMPUTED-TOMOGRAPHY; HIGH TEMPORAL RESOLUTION; DUAL-ENERGY CT; RADIATION-EXPOSURE; RECONSTRUCTION; PICCS; MRI; FEASIBILITY; ALGORITHMS;
D O I
10.1118/1.3560878
中图分类号
R8 [特种医学]; R445 [影像诊断学];
学科分类号
1002 ; 100207 ; 1009 ;
摘要
Purpose: To present and evaluate a new image reconstruction method for dynamic CT based on a nonconvex prior image constrained compressed sensing (NCPICCS) algorithm. The authors systematically compared the undersampling potential, functional information recovery, and solution convergence speed of four compressed sensing (CS) based image reconstruction methods using perfusion CT data: Standard l(1)-based CS, nonconvex CS (NCCS), and l(1)-based and nonconvex CS, including an additional constraint based on a prior image (PICCS and NCPICCS, respectively). Methods: The Shepp-Logan phantom was modified such that its uppermost ellipses changed attenuation through time, simulating both an arterial input function (AIF) and a homogeneous tissue perfusion region. Data were simulated with and without Poisson noise added to the projection data and subsequently reconstructed with all four CS-based methods at four levels of undersampling: 20, 12, 6, and 4 projections. Root mean squared (RMS) error of reconstructed images and recovered time attenuation curves (TACs) were assessed as well as convergence speed. The performance of both PICCS and NCPICCS methods were also evaluated using a kidney perfusion animal experiment data set. Results: All four CS-based methods were able to reconstruct the phantoms with 20 projections, with similar results on the RMS error of the recovered TACs. NCCS allowed accurate reconstructions with as few as 12 projections, PICCS with as few as six projections, and NCPICCS with as few as four projections. These results were consistent for noise-free and noisy data. NCPICCS required the fewest iterations to converge across all simulation conditions, followed by PICCS, NCCS, and then CS. On animal data, at the lowest level of undersampling tested (16 projections), the image quality of NCPICCS was better than PICCS with fewer streaking artifacts, while the TAC accuracy on the selected region of interest was comparable. Conclusions: The authors have presented a novel method for image reconstruction using highly undersampled dynamic CT data. The NCPICCS method takes advantage of the information provided by a prior image, as in PICCS, but employs a more general nonconvex sparsity measure [such as the l(p)-norm (0 <= p <= 1)] rather than the conventional convex l(1)-norm. Despite the lack of guarantees of a globally optimal solution, the proposed nonconvex extension of PICCS consistently allowed for image reconstruction from fewer samples than the analogous l(1)-based PICCS method. Both nonconvex sparsity measures as well as prior image information (when available) significantly reduced the number of iterations required for convergence, potentially providing computational advantages for practical implementation of CS-based image reconstruction techniques. (C) 2011 American Association of Physicists in Medicine. [DOI: 10.1118/1.3560878]
引用
收藏
页码:2157 / 2167
页数:11
相关论文
共 44 条
[21]   Strategies for Reducing Radiation Dose in CT [J].
McCollough, Cynthia H. ;
Primak, Andrew N. ;
Braun, Natalie ;
Kofler, James ;
Yu, Lifeng ;
Christner, Jodie .
RADIOLOGIC CLINICS OF NORTH AMERICA, 2009, 47 (01) :27-+
[22]   Perfusion CT: a worthwhile enhancement? [J].
Miles, KA ;
Griffiths, MR .
BRITISH JOURNAL OF RADIOLOGY, 2003, 76 (904) :220-231
[23]   Perfusion measurements by micro-CT using prior image constrained compressed sensing (PICCS): initial phantom results [J].
Nett, Brian E. ;
Brauweiler, Robert ;
Kalender, Willi ;
Rowley, Howard ;
Chen, Guang-Hong .
PHYSICS IN MEDICINE AND BIOLOGY, 2010, 55 (08) :2333-2350
[24]   Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction? [J].
Pan, Xiaochuan ;
Sidky, Emil Y. ;
Vannier, Michael .
INVERSE PROBLEMS, 2009, 25 (12)
[25]   Three-Dimensional Imaging Properties of Rotation-Free Square and Hexagonal Micro-CT Systems [J].
Quan, Enzhuo Michelle ;
Lalush, David S. .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 2010, 29 (03) :916-923
[26]  
RAMIREZGIRALDO J, 2010, P SPIE MED IM C C, V7622
[27]   3-DIMENSIONAL IMAGING OF HEART, LUNGS, AND CIRCULATION [J].
RITMAN, EL ;
KINSEY, JH ;
ROBB, RA ;
GILBERT, BK ;
HARRIS, LD ;
WOOD, EH .
SCIENCE, 1980, 210 (4467) :273-280
[28]   Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization [J].
Sidky, Emil Y. ;
Pan, Xiaochuan .
PHYSICS IN MEDICINE AND BIOLOGY, 2008, 53 (17) :4777-4807
[29]   Image reconstruction from few views by non-convex optimization [J].
Sidky, Emil Y. ;
Chartrand, Rick ;
Pan, Xiaochuan .
2007 IEEE NUCLEAR SCIENCE SYMPOSIUM CONFERENCE RECORD, VOLS 1-11, 2007, :3526-+
[30]  
Sidky EY, 2006, J X-RAY SCI TECHNOL, V14, P119