How little data is enough? Phase-diagram analysis of sparsity-regularized X-ray computed tomography.

J S Jørgensen, E Y Sidky
Author Information
  1. J S Jørgensen: Department of Applied Mathematics and Computer Science, Technical University of Denmark, Richard Petersens Plads, Kongens Lyngby 2800, Denmark jakj@dtu.dk.
  2. E Y Sidky: Department of Radiology MC-2026, University of Chicago, 5841 South Maryland Avenue, Chicago, IL 60637, USA.

Abstract

We introduce phase-diagram analysis, a standard tool in compressed sensing (CS), to the X-ray computed tomography (CT) community as a systematic method for determining how few projections suffice for accurate sparsity-regularized reconstruction. In CS, 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 comparably with a near-optimal CS 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 with 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.

Keywords

References

  1. IEEE Trans Med Imaging. 2000 May;19(5):548-55 [PMID: 11021698]
  2. Phys Med Biol. 2002 Aug 7;47(15):2599-609 [PMID: 12200927]
  3. Phys Med Biol. 2004 Jun 7;49(11):2463-75 [PMID: 15248590]
  4. Proc Natl Acad Sci U S A. 2005 Jul 5;102(27):9446-51 [PMID: 15976026]
  5. Med Phys. 2007 Nov;34(11):4476-83 [PMID: 18072512]
  6. IEEE Trans Image Process. 1998;7(2):204-21 [PMID: 18267394]
  7. Med Phys. 2008 Feb;35(2):660-3 [PMID: 18383687]
  8. Phys Med Biol. 2008 Sep 7;53(17):4777-807 [PMID: 18701771]
  9. Philos Trans A Math Phys Eng Sci. 2009 Nov 13;367(1906):4273-93 [PMID: 19805445]
  10. Phys Med Biol. 2010 Nov 21;55(22):6575-99 [PMID: 20962368]
  11. Phys Med Biol. 2011 Mar 21;56(6):1545-61 [PMID: 21325707]
  12. Phys Med Biol. 2012 May 21;57(10):3065-91 [PMID: 22538474]
  13. IEEE Trans Med Imaging. 2013 Feb;32(2):460-73 [PMID: 23204282]
  14. Med Phys. 2013 Mar;40(3):031115 [PMID: 23464295]
  15. IEEE Trans Image Process. 2013 Oct;22(10):3941-9 [PMID: 23708808]
  16. J Opt Soc Am A Opt Image Sci Vis. 2014 Jul 1;31(7):1369-94 [PMID: 25121423]
  17. IEEE J Transl Eng Health Med. 2014 Jun 30;2: [PMID: 25401059]
  18. Inverse Probl Imaging (Springfield). 2015 May;9(2):431-446 [PMID: 27019675]

Grants

  1. R01-CA158446/NCI NIH HHS
  2. R01-CA182264/NCI NIH HHS
  3. R01-EB000225/NIBIB NIH HHS
  4. R01-EB018102/NIBIB NIH HHS

MeSH Term

Algorithms
Computer Simulation
Data Compression
Data Interpretation, Statistical
Humans
Models, Statistical
Radiographic Image Enhancement
Radiographic Image Interpretation, Computer-Assisted
Reproducibility of Results
Sample Size
Sensitivity and Specificity
Tomography, X-Ray Computed

Word Cloud

Created with Highcharts 10.0.0analysisphase-diagramX-rayCTsensingCScomputedtomographysparsitysamplingstandardcompressedprojectionssparsity-regularizedreconstructiontheoreticalsufficientresultsdemonstratetakingmeasurementsshowimageregularizationintroducetoolcommunitysystematicmethoddeterminingsufficeaccuratephasediagramconvenientwaystudyexpresscertainrelationsadaptempiricaluseholdthreecasestudiespotentialprovidingquantitativeanswersquestionsundersamplingFirstcasesempiricallyperformscomparablynear-optimalstrategynamelyGaussianmatricesSecondcontrastmightanticipatedrandomizedleadimprovedperformancecomparedstructuredpatternsFinallypreliminarywellcanpredictnumberaccuratelyreconstructinglarge-scalegivenmeanstotal-variationlittledataenough?Phase-diagram

Similar Articles

Cited By