Multiresolution Parameter Choice Method for Total Variation Regularized Tomography

A computational method is introduced for choosing the regularization parameter for total variation (TV) regularization. A partial understanding of the properties of the method is provided by rigorously proving that the TV norms of the reconstructions converge with any choice of regularization parameter. The computational approach is based on computing reconstructions at a few different resolutions and various values of regularization parameter. The chosen parameter is the smallest one resulting in approximately discretization-invariant TV norms of the reconstructions. The method is tested with simulated and experimental X-ray tomography data and compared to the S-curve method. The results are comparable to those of the S-curve method. However, the S-curve method needs quantitative a priori information about the expected sparsity (TV norm) of the unknown, while the proposed method does not require such input parameters.