How little data is enough? Phase-diagram analysis of sparsity-regularized X-ray CT

Abstract

We introduce phase-diagram analysis, a standard tool in compressed sensing, to the X-ray CT community as a systematic method for determining how few projections suffice for accurate sparsity-regularized reconstruction. In compressed sensing a phase diagram is a convenient way to study and express certain theoretical relations between sparsity and sufficient sampling. We adapt phase-diagram analysis for empirical use in X-ray CT for which the same theoretical results do not hold. We demonstrate in three case studies the potential of phase-diagram analysis for providing quantitative answers to questions of undersampling: First we demonstrate that there are cases where X-ray CT empirically performs comparable with an optimal compressed sensing strategy, namely taking measurements with Gaussian sensing matrices. Second, we show that, in contrast to what might have been anticipated, taking randomized CT measurements does not lead to improved performance compared to standard structured sampling patterns. Finally, we show preliminary results of how well phase-diagram analysis can predict the sufficient number of projections for accurately reconstructing a large-scale image of a given sparsity by means of total-variation regularization.

arXiv comment24 pages, 13 figures

References in corpus (6)

• Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT
• Observed Universality of Phase Transitions in High-Dimensional Geometry, with Implications for Modern Data Analysis and Signal Processing
• Convex optimization problem prototyping for image reconstruction in computed tomography with the Chambolle-Pock algorithm
• Universality in polytope phase transitions and message passing algorithms
• Near-optimal compressed sensing guarantees for total variation minimization
• Testable uniqueness conditions for empirical assessment of undersampling levels in total variation-regularized x-ray CT

Cited by in corpus (9)

• An algorithm for constrained one-step inversion of spectral CT data
• Accelerated High-Resolution Photoacoustic Tomography via Compressed Sensing
• Learning The Invisible: A Hybrid Deep Learning-Shearlet Framework for Limited Angle Computed Tomography
• Superiorized algorithm for reconstruction of CT images from sparse-view and limited-angle polyenergetic data
• Sparse dynamic tomography. A shearlet-based approach for iodine perfusion in plant stems
• A signal detection model for quantifying over-regularization in non-linear image reconstruction
• Shearlet-based regularization in statistical inverse learning with an application to X-ray tomography
• A Local Fourier Slice Theorem
• Regularization with optimal space-time priors