1,274 research outputs found
Sparse Spikes Deconvolution on Thin Grids
This article analyzes the recovery performance of two popular finite
dimensional approximations of the sparse spikes deconvolution problem over
Radon measures. We examine in a unified framework both the L1 regularization
(often referred to as Lasso or Basis-Pursuit) and the Continuous Basis-Pursuit
(C-BP) methods. The Lasso is the de-facto standard for the sparse
regularization of inverse problems in imaging. It performs a nearest neighbor
interpolation of the spikes locations on the sampling grid. The C-BP method,
introduced by Ekanadham, Tranchina and Simoncelli, uses a linear interpolation
of the locations to perform a better approximation of the infinite-dimensional
optimization problem, for positive measures. We show that, in the small noise
regime, both methods estimate twice the number of spikes as the number of
original spikes. Indeed, we show that they both detect two neighboring spikes
around the locations of an original spikes. These results for deconvolution
problems are based on an abstract analysis of the so-called extended support of
the solutions of L1-type problems (including as special cases the Lasso and
C-BP for deconvolution), which are of an independent interest. They precisely
characterize the support of the solutions when the noise is small and the
regularization parameter is selected accordingly. We illustrate these findings
to analyze for the first time the support instability of compressed sensing
recovery when the number of measurements is below the critical limit (well
documented in the literature) where the support is provably stable
Sampling and Recovery of Pulse Streams
Compressive Sensing (CS) is a new technique for the efficient acquisition of
signals, images, and other data that have a sparse representation in some
basis, frame, or dictionary. By sparse we mean that the N-dimensional basis
representation has just K<<N significant coefficients; in this case, the CS
theory maintains that just M = K log N random linear signal measurements will
both preserve all of the signal information and enable robust signal
reconstruction in polynomial time. In this paper, we extend the CS theory to
pulse stream data, which correspond to S-sparse signals/images that are
convolved with an unknown F-sparse pulse shape. Ignoring their convolutional
structure, a pulse stream signal is K=SF sparse. Such signals figure
prominently in a number of applications, from neuroscience to astronomy. Our
specific contributions are threefold. First, we propose a pulse stream signal
model and show that it is equivalent to an infinite union of subspaces. Second,
we derive a lower bound on the number of measurements M required to preserve
the essential information present in pulse streams. The bound is linear in the
total number of degrees of freedom S + F, which is significantly smaller than
the naive bound based on the total signal sparsity K=SF. Third, we develop an
efficient signal recovery algorithm that infers both the shape of the impulse
response as well as the locations and amplitudes of the pulses. The algorithm
alternatively estimates the pulse locations and the pulse shape in a manner
reminiscent of classical deconvolution algorithms. Numerical experiments on
synthetic and real data demonstrate the advantages of our approach over
standard CS
Super-Resolution in Phase Space
This work considers the problem of super-resolution. The goal is to resolve a
Dirac distribution from knowledge of its discrete, low-pass, Fourier
measurements. Classically, such problems have been dealt with parameter
estimation methods. Recently, it has been shown that convex-optimization based
formulations facilitate a continuous time solution to the super-resolution
problem. Here we treat super-resolution from low-pass measurements in Phase
Space. The Phase Space transformation parametrically generalizes a number of
well known unitary mappings such as the Fractional Fourier, Fresnel, Laplace
and Fourier transforms. Consequently, our work provides a general super-
resolution strategy which is backward compatible with the usual Fourier domain
result. We consider low-pass measurements of Dirac distributions in Phase Space
and show that the super-resolution problem can be cast as Total Variation
minimization. Remarkably, even though are setting is quite general, the bounds
on the minimum separation distance of Dirac distributions is comparable to
existing methods.Comment: 10 Pages, short paper in part accepted to ICASSP 201
Projected gradient descent for non-convex sparse spike estimation
We propose a new algorithm for sparse spike estimation from Fourier
measurements. Based on theoretical results on non-convex optimization
techniques for off-the-grid sparse spike estimation, we present a projected
gradient descent algorithm coupled with a spectral initialization procedure.
Our algorithm permits to estimate the positions of large numbers of Diracs in
2d from random Fourier measurements. We present, along with the algorithm,
theoretical qualitative insights explaining the success of our algorithm. This
opens a new direction for practical off-the-grid spike estimation with
theoretical guarantees in imaging applications
- …