Chen, Donghui3; Kilmer, Misha E.4; Hansen, Per Christian1
1 Department of Applied Mathematics and Computer Science, Technical University of Denmark2 Scientific Computing, Department of Applied Mathematics and Computer Science, Technical University of Denmark3 Southwestern University of Finance and Economics4 Tufts University
In many inverse problems it is essential to use regularization methods that preserve edges in the reconstructions, and many reconstruction models have been developed for this task, such as the Total Variation (TV) approach. The associated algorithms are complex and require a good knowledge of large-scale optimization algorithms, and they involve certain tolerances that the user must choose. We present a simpler approach that relies only on standard computational building blocks in matrix computations, such as orthogonal transformations, preconditioned iterative solvers, Kronecker products, and the discrete cosine transform.hence the term "plug-and-play" . We do not attempt to improve on TV reconstructions, but rather provide an easy-to-use approach to computing reconstructions with similar properties.
Electronic Transactions on Numerical Analysis, 2014, Vol 41, p. 465-477