Continuous-Domain Solutions of Linear Inverse Problems With Tikhonov Versus Generalized TV Regularization

We consider one-dimensional (1-D) linear inverse problems that are formulated in the continuous domain. The object of recovery is a function that is assumed to minimize a convex objective functional. The solutions are constrained by imposing a continuous-domain regularization. We derive the parametric form of the solution (representer theorems) for Tikhonov (quadratic) and generalized total-variation (gTV) regularizations. We show that, in both cases, the solutions are splines that are intimately related to the regularization operator. In the Tikhonov case, the solution is smooth and constrained to live in a fixed subspace that depends on the measurement operator. By contrast, the gTV regularization results in a sparse solution composed of only a few dictionary elements that are upper-bounded by the number of measurements and independent of the measurement operator. Our findings for the gTV regularization resonates with the minimization of the $\ell _1$ norm, which is its discrete counterpart and also produces sparse solutions. Finally, we find the experimental solutions for some measurement models in one dimension. We discuss the special case when the gTV regularization results in multiple solutions and devise an algorithm to find an extreme point of the solution set which is guaranteed to be sparse.

[1]  Robert D. Nowak,et al.  An EM algorithm for wavelet-based image restoration , 2003, IEEE Trans. Image Process..

[2]  Gongguo Tang,et al.  Sparse recovery over continuous dictionaries-just discretize , 2013, 2013 Asilomar Conference on Signals, Systems and Computers.

[3]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[4]  Mário A. T. Figueiredo,et al.  Gradient Projection for Sparse Reconstruction: Application to Compressed Sensing and Other Inverse Problems , 2007, IEEE Journal of Selected Topics in Signal Processing.

[5]  Didier Henrion,et al.  Exact Solutions to Super Resolution on Semi-Algebraic Domains in Higher Dimensions , 2015, IEEE Transactions on Information Theory.

[6]  Thierry Blu,et al.  Generalized sampling: a variational approach .I. Theory , 2002, IEEE Trans. Signal Process..

[7]  Gabriel Peyré,et al.  Support Recovery for Sparse Super-Resolution of Positive Measures , 2017 .

[8]  Mario Bertero,et al.  Introduction to Inverse Problems in Imaging , 1998 .

[9]  Thierry Blu,et al.  Generalized smoothing splines and the optimal discretization of the Wiener filter , 2005, IEEE Transactions on Signal Processing.

[10]  I. Csiszár Why least squares and maximum entropy? An axiomatic approach to inference for linear inverse problems , 1991 .

[11]  G. Peyré,et al.  Sparse regularization on thin grids I: the Lasso , 2017 .

[12]  Holger Rauhut,et al.  A Mathematical Introduction to Compressive Sensing , 2013, Applied and Numerical Harmonic Analysis.

[13]  R. Tibshirani The Lasso Problem and Uniqueness , 2012, 1206.0313.

[14]  Emmanuel J. Candès,et al.  Towards a Mathematical Theory of Super‐resolution , 2012, ArXiv.

[15]  Michael Unser,et al.  A representer theorem for deep neural networks , 2018, J. Mach. Learn. Res..

[16]  A. Chambolle,et al.  On the convergence of the iterates of "FISTA" , 2015 .

[17]  M. Reed Methods of Modern Mathematical Physics. I: Functional Analysis , 1972 .

[18]  M. Urner Scattered Data Approximation , 2016 .

[19]  K. Bredies,et al.  Inverse problems in spaces of measures , 2013 .

[20]  Emmanuel J. Candès,et al.  Super-Resolution from Noisy Data , 2012, Journal of Fourier Analysis and Applications.

[21]  Joseph W. Jerome,et al.  Spline solutions to L1 extremal problems in one and several variables , 1975 .

[22]  P. Weiss,et al.  Exact solutions of infinite dimensional total-variation regularized problems , 2017, Information and Inference: A Journal of the IMA.

[23]  D. Donoho,et al.  Sparse MRI: The application of compressed sensing for rapid MR imaging , 2007, Magnetic resonance in medicine.

[24]  Arthur E. Hoerl,et al.  Ridge Regression: Biased Estimation for Nonorthogonal Problems , 2000, Technometrics.

[25]  G. Wahba Support vector machines, reproducing kernel Hilbert spaces, and randomized GACV , 1999 .

[26]  R. Cooke Real and Complex Analysis , 2011 .

[27]  E. Polak Introduction to linear and nonlinear programming , 1973 .

[28]  Anthony Widjaja,et al.  Learning with Kernels: Support Vector Machines, Regularization, Optimization, and Beyond , 2003, IEEE Transactions on Neural Networks.

[29]  Michael Unser,et al.  An Introduction to Sparse Stochastic Processes , 2014 .

[30]  Michael Unser,et al.  Gaussian and Sparse Processes Are Limits of Generalized Poisson Processes , 2017, Applied and Computational Harmonic Analysis.

[31]  Stephan Didas,et al.  Splines in Higher Order TV Regularization , 2006, International Journal of Computer Vision.

[32]  G. Wahba Spline models for observational data , 1990 .

[33]  Michael Unser,et al.  Stochastic Models for Sparse and Piecewise-Smooth Signals , 2011, IEEE Transactions on Signal Processing.

[34]  Thierry Blu,et al.  Generalized sampling: a variational approach .II. Applications , 2002, IEEE Trans. Signal Process..

[35]  A. Y. Bezhaev,et al.  Variational Theory of Splines , 2001 .

[36]  Andrea Braides Gamma-Convergence for Beginners , 2002 .

[37]  M. Unser Generalized Sampling : A Variational Approach — Part I : Theory , 2001 .

[38]  A. Chambolle,et al.  Geometric properties of solutions to the total variation denoising problem , 2016, 1602.00087.

[39]  M. Unser,et al.  Generalized Sampling : A Variational Approach . Part II : Applications , 2002 .

[40]  S. Frick,et al.  Compressed Sensing , 2014, Computer Vision, A Reference Guide.

[41]  R. Tibshirani,et al.  Discussion: The Dantzig selector: Statistical estimation when p is much larger than n , 2007, 0803.3126.

[42]  R. Tibshirani Regression Shrinkage and Selection via the Lasso , 1996 .

[43]  Gabriel Peyré,et al.  Exact Support Recovery for Sparse Spikes Deconvolution , 2013, Foundations of Computational Mathematics.

[44]  lexander,et al.  THE GENERALIZED SIMPLEX METHOD FOR MINIMIZING A LINEAR FORM UNDER LINEAR INEQUALITY RESTRAINTS , 2012 .

[45]  Parikshit Shah,et al.  Compressed Sensing Off the Grid , 2012, IEEE Transactions on Information Theory.

[46]  Bernhard Schölkopf,et al.  A Generalized Representer Theorem , 2001, COLT/EuroCOLT.

[47]  E. Caianiello Functional analysis and optimization , 1966 .

[48]  S. Geer,et al.  Locally adaptive regression splines , 1997 .

[49]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[50]  A Tikhonov,et al.  Solution of Incorrectly Formulated Problems and the Regularization Method , 1963 .

[51]  Pierre Vandergheynst,et al.  Compressed Sensing and Redundant Dictionaries , 2007, IEEE Transactions on Information Theory.

[52]  E. Candès,et al.  Sparsity and incoherence in compressive sampling , 2006, math/0611957.

[53]  Michael Unser,et al.  Splines Are Universal Solutions of Linear Inverse Problems with Generalized TV Regularization , 2016, SIAM Rev..