Total variation regularization for large-scale X-ray tomography

Keijo Hämäläinen, Lauri Harhanen, Andreas Hauptmann, Aki Kallonen, Esa Niemi, Samuli Siltanen

Research output: Contribution to journalArticleScientificpeer-review

16 Citations (Scopus)


A new large-scale computational total variation regularization algorithm is introduced and tested with examples arising from X-ray tomography with sparsely sampled data. The total variation penalty term is discretized using a basis of discontinuous functions. The approach is motivated by discontinuous Galerkin methods and leads to an additional term of the jump part of total variation. The proposed algorithm combines the usage of the jump part with a subgradient descent scheme. A comparison is provided with the gradient-based projected Barzilai-Borwein method which uses a smoothly approximated total variation penalty. The above two methods are examples of total variation regularization algorithms that can be applied to large-scale tomographic problems in reasonable computation time. A comparison between the methods shows that they use roughly equal computational resources and that the new method produces somewhat blockier reconstructions. Although the test problems are two-dimensional, both methods can be applied to three-dimensional situations as well.

Original languageEnglish
Pages (from-to)1-25
Number of pages25
JournalInternational Journal of Tomography and Simulation
Issue number1
Publication statusPublished - 2014
MoE publication typeA1 Journal article-refereed


  • Barzilai-Borwein
  • First-order methods
  • Subgradient descent
  • Total variation
  • X-ray tomography


Dive into the research topics of 'Total variation regularization for large-scale X-ray tomography'. Together they form a unique fingerprint.

Cite this