541 research outputs found
Sampling and Reconstruction of Sparse Signals on Circulant Graphs - An Introduction to Graph-FRI
With the objective of employing graphs toward a more generalized theory of
signal processing, we present a novel sampling framework for (wavelet-)sparse
signals defined on circulant graphs which extends basic properties of Finite
Rate of Innovation (FRI) theory to the graph domain, and can be applied to
arbitrary graphs via suitable approximation schemes. At its core, the
introduced Graph-FRI-framework states that any K-sparse signal on the vertices
of a circulant graph can be perfectly reconstructed from its
dimensionality-reduced representation in the graph spectral domain, the Graph
Fourier Transform (GFT), of minimum size 2K. By leveraging the recently
developed theory of e-splines and e-spline wavelets on graphs, one can
decompose this graph spectral transformation into the multiresolution low-pass
filtering operation with a graph e-spline filter, and subsequent transformation
to the spectral graph domain; this allows to infer a distinct sampling pattern,
and, ultimately, the structure of an associated coarsened graph, which
preserves essential properties of the original, including circularity and,
where applicable, the graph generating set.Comment: To appear in Appl. Comput. Harmon. Anal. (2017
Restoration of Poissonian Images Using Alternating Direction Optimization
Much research has been devoted to the problem of restoring Poissonian images,
namely for medical and astronomical applications. However, the restoration of
these images using state-of-the-art regularizers (such as those based on
multiscale representations or total variation) is still an active research
area, since the associated optimization problems are quite challenging. In this
paper, we propose an approach to deconvolving Poissonian images, which is based
on an alternating direction optimization method. The standard regularization
(or maximum a posteriori) restoration criterion, which combines the Poisson
log-likelihood with a (non-smooth) convex regularizer (log-prior), leads to
hard optimization problems: the log-likelihood is non-quadratic and
non-separable, the regularizer is non-smooth, and there is a non-negativity
constraint. Using standard convex analysis tools, we present sufficient
conditions for existence and uniqueness of solutions of these optimization
problems, for several types of regularizers: total-variation, frame-based
analysis, and frame-based synthesis. We attack these problems with an instance
of the alternating direction method of multipliers (ADMM), which belongs to the
family of augmented Lagrangian algorithms. We study sufficient conditions for
convergence and show that these are satisfied, either under total-variation or
frame-based (analysis and synthesis) regularization. The resulting algorithms
are shown to outperform alternative state-of-the-art methods, both in terms of
speed and restoration accuracy.Comment: 12 pages, 12 figures, 2 tables. Submitted to the IEEE Transactions on
Image Processin
A semidiscrete version of the Citti-Petitot-Sarti model as a plausible model for anthropomorphic image reconstruction and pattern recognition
In his beautiful book [66], Jean Petitot proposes a sub-Riemannian model for
the primary visual cortex of mammals. This model is neurophysiologically
justified. Further developments of this theory lead to efficient algorithms for
image reconstruction, based upon the consideration of an associated
hypoelliptic diffusion. The sub-Riemannian model of Petitot and Citti-Sarti (or
certain of its improvements) is a left-invariant structure over the group
of rototranslations of the plane. Here, we propose a semi-discrete
version of this theory, leading to a left-invariant structure over the group
, restricting to a finite number of rotations. This apparently very
simple group is in fact quite atypical: it is maximally almost periodic, which
leads to much simpler harmonic analysis compared to Based upon this
semi-discrete model, we improve on previous image-reconstruction algorithms and
we develop a pattern-recognition theory that leads also to very efficient
algorithms in practice.Comment: 123 pages, revised versio
Structured random measurements in signal processing
Compressed sensing and its extensions have recently triggered interest in
randomized signal acquisition. A key finding is that random measurements
provide sparse signal reconstruction guarantees for efficient and stable
algorithms with a minimal number of samples. While this was first shown for
(unstructured) Gaussian random measurement matrices, applications require
certain structure of the measurements leading to structured random measurement
matrices. Near optimal recovery guarantees for such structured measurements
have been developed over the past years in a variety of contexts. This article
surveys the theory in three scenarios: compressed sensing (sparse recovery),
low rank matrix recovery, and phaseless estimation. The random measurement
matrices to be considered include random partial Fourier matrices, partial
random circulant matrices (subsampled convolutions), matrix completion, and
phase estimation from magnitudes of Fourier type measurements. The article
concludes with a brief discussion of the mathematical techniques for the
analysis of such structured random measurements.Comment: 22 pages, 2 figure
Sampling Theorem and Discrete Fourier Transform on the Riemann Sphere
Using coherent-state techniques, we prove a sampling theorem for Majorana's
(holomorphic) functions on the Riemann sphere and we provide an exact
reconstruction formula as a convolution product of samples and a given
reconstruction kernel (a sinc-type function). We also discuss the effect of
over- and under-sampling. Sample points are roots of unity, a fact which allows
explicit inversion formulas for resolution and overlapping kernel operators
through the theory of Circulant Matrices and Rectangular Fourier Matrices. The
case of band-limited functions on the Riemann sphere, with spins up to , is
also considered. The connection with the standard Euler angle picture, in terms
of spherical harmonics, is established through a discrete Bargmann transform.Comment: 26 latex pages. Final version published in J. Fourier Anal. App
- …