Sparse Stochastic Processes and Discretization of Linear Inverse Problems
E. Bostan, U.S. Kamilov, M. Nilchian, M. Unser
IEEE Transactions on Image Processing, vol. 22, no. 7, pp. 2699–2710, July 2013.
We present a novel statistically-based discretization paradigm and derive a class of maximum a posteriori (MAP) estimators for solving ill-conditioned linear inverse problems. We are guided by the theory of sparse stochastic processes, which specifies continuous-domain signals as solutions of linear stochastic differential equations. Accordingly, we show that the class of admissible priors for the discretized version of the signal is confined to the family of infinitely divisible distributions. Our estimators not only cover the well-studied methods of Tikhonov and ℓ1-type regularizations as particular cases, but also open the door to a broader class of sparsity-promoting regularization schemes that are typically nonconvex. We provide an algorithm that handles the corresponding nonconvex problems and illustrate the use of our formalism by applying it to deconvolution, magnetic resonance imaging, and X-ray tomographic reconstruction problems. Finally, we compare the performance of estimators associated with models of increasing sparsity.
@ARTICLE(http://bigwww.epfl.ch/publications/bostan1303.html, AUTHOR="Bostan, E. and Kamilov, U.S. and Nilchian, M. and Unser, M.", TITLE="Sparse Stochastic Processes and Discretization of Linear Inverse Problems", JOURNAL="{IEEE} Transactions on Image Processing", YEAR="2013", volume="22", number="7", pages="2699--2710", month="July", note="")