Testable uniqueness conditions for empirical assessment of undersampling levels in total variation-regularized x-ray CT

Abstract

We study recoverability in fan-beam computed tomography (CT) with sparsity and total variation priors: how many underdetermined linear measurements suffice for recovering images of given sparsity? Results from compressed sensing (CS) establish such conditions for, e.g., random measurements, but not for CT. Recoverability is typically tested by checking whether a computed solution recovers the original. This approach cannot guarantee solution uniqueness and the recoverability decision therefore depends on the optimization algorithm. We propose new computational methods to test recoverability by verifying solution uniqueness conditions. Using both reconstruction and uniqueness testing we empirically study the number of CT measurements sufficient for recovery on new classes of sparse test images. We demonstrate an average-case relation between sparsity and sufficient sampling and observe a sharp phase transition as known from CS, but never established for CT. In addition to assessing recoverability more reliably, we show that uniqueness tests are often the faster option.

arXiv comment18 pages, 7 figures, submitted

References in corpus (7)

• 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
• The Cosparse Analysis Model and Algorithms
• Near-optimal compressed sensing guarantees for total variation minimization
• One condition for solution uniqueness and robustness of both l1-synthesis and l1-analysis minimizations
• Polar Polytopes and Recovery of Sparse Representations
• Phase Transitions and Cosparse Tomographic Recovery of Compound Solid Bodies from Few Projections

Cited by in corpus (4)

• How little data is enough? Phase-diagram analysis of sparsity-regularized X-ray CT
• Sampling Limits for Electron Tomography with Sparsity-exploiting Reconstructions
• Shearlet-based regularization in statistical inverse learning with an application to X-ray tomography
• Quantitative Study on Exact Reconstruction Sampling Condition in Limited-view CT