Sparse spikes super-resolution on thin grids II: the continuous basis pursuit

This article analyzes the performance of the Continuous Basis Pursuit (C-BP) method for sparse super-resolution. The C-BP has been recently proposed by Ekanadham, Tranchina and Simoncelli as a refined discretization scheme for the recovery of spikes in inverse problems regularization. One of the most well known discretization scheme, the Basis Pursuit (BP, also known as Lasso) makes use of a finite dimensional l1 norm on a grid. In contrast, the C-BP uses a linear interpolation of the spikes positions to enable the recovery of spikes between grid points. When the sought-after solution is constrained to be positive, a remarkable feature of this approach is that it retains the convexity of the initial l1 problem. The present paper shows how the C-BP is able to recover the spikes locations with sub-grid accuracy in the favorable case. We also prove that this regime generally breaks when the grid is too thin, and we describe precisely the artifacts that appear: each spike is approximated by a pair of Dirac masses. We show numerical illustrations of these phenomena, and evaluate numerically the validity of the technical assumptions of our analysis.

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

[2]  Michael A. Saunders,et al.  Atomic Decomposition by Basis Pursuit , 1998, SIAM J. Sci. Comput..

[3]  Jean-Jacques Fuchs,et al.  On sparse representations in arbitrary redundant bases , 2004, IEEE Transactions on Information Theory.

[4]  Oskar Maria Baksalary,et al.  Particular formulae for the Moore-Penrose inverse of a columnwise partitioned matrix , 2007 .

[5]  R. Tibshirani,et al.  Regression shrinkage and selection via the lasso: a retrospective , 2011 .

[6]  Heinz H. Bauschke,et al.  Convex Analysis and Monotone Operator Theory in Hilbert Spaces , 2011, CMS Books in Mathematics.

[7]  Eero P. Simoncelli,et al.  Recovery of Sparse Translation-Invariant Signals With Continuous Basis Pursuit , 2011, IEEE Transactions on Signal Processing.

[8]  Yohann de Castro,et al.  Exact Reconstruction using Beurling Minimal Extrapolation , 2011, 1103.4951.

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

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

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

[12]  Gongguo Tang Atomic Decomposition of Mixtures of Translation-Invariant Signals , 2013 .

[13]  Jean-Christophe Pesquet,et al.  A constrained optimization approach for complex sparse perturbed models , 2013 .

[14]  Carlos Fernandez-Granda Support detection in super-resolution , 2013, ArXiv.

[15]  F. Gamboa,et al.  Spike detection from inaccurate samplings , 2013, 1301.5873.

[16]  Eero P. Simoncelli,et al.  Journal of Neuroscience Methods , 2022 .

[17]  Bastian Goldlücke,et al.  Variational Analysis , 2014, Computer Vision, A Reference Guide.

[18]  Benjamin Recht,et al.  Superresolution without separation , 2015, 2015 IEEE 6th International Workshop on Computational Advances in Multi-Sensor Adaptive Processing (CAMSAP).

[19]  Gabriel Peyré,et al.  Support Recovery for Sparse Deconvolution of Positive Measures , 2015, ArXiv.

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

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

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