
In this paper, we study sparse spike deconvolution over the space of
complexvalued measures when the input measure is a finite sum of Dirac masses.
We introduce a modified version of the Beurling Lasso (BLasso), a semidefinite
program that we refer to as the Concomitant Beurling Lasso (CBLasso). This new
procedure estimates the target measure and the unknown noise level
simultaneously. Contrary to previous estimators in the literature, theory holds
for a tuning parameter that depends only on the sample size, so that it can be
used for unknown noise level problems. Consistent noise level estimation is
standardly proved. As for Radon measure estimation, theoretical guarantees
match the previous stateoftheart results in SuperResolution regarding
minimax prediction and localization. The proofs are based on a bound on the
noise level given by a new tail estimate of the supremum of a stationary
nonGaussian process through the Rice method.

Compressed Sensing (CS) is an appealing framework for applications such as
Magnetic Resonance Imaging (MRI). However, uptodate, the sensing schemes
suggested by CS theories are made of random isolated measurements, which are
usually incompatible with the physics of acquisition. To reflect the physical
constraints of the imaging device, we introduce the notion of blocks of
measurements: the sensing scheme is not a set of isolated measurements anymore,
but a set of groups of measurements which may represent any arbitrary shape
(parallel or radial lines for instance). Structured acquisition with blocks of
measurements are easy to implement, and provide good reconstruction results in
practice. However, very few results exist on the theoretical guarantees of CS
reconstructions in this setting. In this paper, we derive new CS results for
structured acquisitions and signals satisfying a prior structured sparsity. The
obtained results provide a recovery probability of sparse vectors that
explicitly depends on their support. Our results are thus supportdependent and
offer the possibility for flexible assumptions on the sparsity structure.
Moreover, the results are drawingdependent, since we highlight an explicit
dependency between the probability of reconstructing a sparse vector and the
way of choosing the blocks of measurements. Numerical simulations show that the
proposed theory is faithful to experimental observations.

Reducing acquisition time is of fundamental importance in various imaging
modalities. The concept of variable density sampling provides a nice framework
to achieve this. It was justified recently from a theoretical point of view in
the compressed sensing (CS) literature. Unfortunately, the sampling schemes
suggested by current CS theories may not be relevant since they do not take the
acquisition constraints into account (for example, continuity of the
acquisition trajectory in Magnetic Resonance Imaging  MRI). In this paper, we
propose a numerical method to perform variable density sampling with block
constraints. Our main contribution is to propose a new way to draw the blocks
in order to mimic CS strategies based on isolated measurements. The basic idea
is to minimize a tailored dissimilarity measure between a probability
distribution defined on the set of isolated measurements and a probability
distribution defined on a set of blocks of measurements. This problem turns out
to be convex and solvable in high dimension. Our second contribution is to
define an efficient minimization algorithm based on Nesterov's accelerated
gradient descent in metric spaces. We study carefully the choice of the metrics
and of the prox function. We show that the optimal choice may depend on the
type of blocks under consideration. Finally, we show that we can obtain better
MRI reconstruction results using our sampling schemes than standard strategies
such as equiangularly distributed radial lines.