Controlled Wavelet Domain Sparsity for X-ray Tomography

Zenith Purisha, Juho Rimpeläinen, Tatiana Bubba, Samuli Siltanen

Research output: Contribution to journalArticlepeer-review

13 Citations (SciVal)


Tomographic reconstruction is an ill-posed inverse problem that calls for regularization. One possibility is to require sparsity of the unknown in an orthonormal wavelet basis. This, in turn, can be achieved by variational regularization, where the penalty term is the sum of the absolute values of the wavelet coefficients. The primal-dual fixed point algorithm showed that the minimizer of the variational regularization functional can be computed iteratively using a soft-thresholding operation. Choosing the soft-thresholding parameter $ \newcommand{\m}{{\boldsymbol m}} \mu>0$ is analogous to the notoriously difficult problem of picking the optimal regularization parameter in Tikhonov regularization. Here, a novel automatic method is introduced for choosing μ, based on a control algorithm driving the sparsity of the reconstruction to an a priori known ratio of nonzero versus zero wavelet coefficients in the unknown.
Original languageEnglish
Article number014002
JournalMeasurement Science and Technology
Issue number1
Publication statusPublished - 14 Dec 2017


Dive into the research topics of 'Controlled Wavelet Domain Sparsity for X-ray Tomography'. Together they form a unique fingerprint.

Cite this