Adaptive Sparsifying Transforms for Iterative Tomographic Reconstruction

A major challenge in computed tomography imaging is to obtain high-quality images from low-dose measurements. Key to this goal are computationally efficient reconstruction algorithms combined with detailed signal models. We show that the recently introduced adaptive sparsifying transform (AST) signal model provides superior reconstructions from low-dose data at significantly lower cost than competing dictionary learning methods. We further accelerate this technique for tomography by utilizing the Linearized Alternating Direction Method of Multipliers (L-ADMM) to remove the need to solve an expensive least-squares problem that requires computing multiple forward and backward projections. Numerical experiments on data from clinical CT images show that adaptive sparsifying transform regularization outperforms total-variation and dictionary learning methods, and combining our regularizer with L-ADMM provides for faster reconstructions than standard ADMM.

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

[2]  Jeffrey A. Fessler,et al.  Conjugate-gradient preconditioning methods for shift-variant PET image reconstruction , 1999, IEEE Trans. Image Process..

[3]  Guillermo Sapiro,et al.  Sparse representations for limited data tomography , 2008, 2008 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[4]  Stanley Osher,et al.  A Unified Primal-Dual Algorithm Framework Based on Bregman Iteration , 2010, J. Sci. Comput..

[5]  Yoram Bresler,et al.  Learning doubly sparse transforms for image representation , 2012, 2012 19th IEEE International Conference on Image Processing.

[6]  Jeffrey A. Fessler,et al.  A Splitting-Based Iterative Algorithm for Accelerated Statistical X-Ray CT Reconstruction , 2012, IEEE Transactions on Medical Imaging.

[7]  Lei Zhang,et al.  Low-Dose X-ray CT Reconstruction via Dictionary Learning , 2012, IEEE Transactions on Medical Imaging.

[8]  Michael Elad,et al.  Analysis K-SVD: A Dictionary-Learning Algorithm for the Analysis Sparse Model , 2013, IEEE Transactions on Signal Processing.

[9]  Yoram Bresler,et al.  Learning Sparsifying Transforms , 2013, IEEE Transactions on Signal Processing.

[10]  Klaus Diepold,et al.  Analysis Operator Learning and its Application to Image Reconstruction , 2012, IEEE Transactions on Image Processing.

[11]  Yoram Bresler,et al.  Closed-form solutions within sparsifying transform learning , 2013, 2013 IEEE International Conference on Acoustics, Speech and Signal Processing.

[12]  Aleksandra Pizurica,et al.  Iterative CT Reconstruction Using Shearlet-Based Regularization , 2013, IEEE Transactions on Nuclear Science.

[13]  Yoram Bresler,et al.  Sparsifying transform learning for Compressed Sensing MRI , 2013, 2013 IEEE 10th International Symposium on Biomedical Imaging.

[14]  Rémi Gribonval,et al.  Constrained Overcomplete Analysis Operator Learning for Cosparse Signal Modelling , 2012, IEEE Transactions on Signal Processing.

[15]  Yoram Bresler,et al.  Tomographic reconstruction with adaptive sparsifying transforms , 2014, 2014 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).