Controlled Wavelet Domain Sparsity in X-ray Tomography
This work addresses the challenging problem of regularization parameter selection in tomographic reconstruction for practitioners, offering an automated approach based on known sparsity priors.
The authors propose an automatic method for selecting the soft threshold parameter in wavelet-based variational regularization for X-ray tomography, using a control algorithm to achieve a target sparsity level. The method eliminates manual tuning of the regularization parameter.
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 absolute values of wavelet coefficients. Daubechies, Defrise and De Mol (Comm. Pure Appl. Math. 57) showed that the minimizer of the variational regularization functional can be computed iteratively using a soft thresholding operation. Choosing the soft threshold parameter $μ>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 {\it a priori} known ratio of nonzero versus zero wavelet coefficients in the unknown function.