91,811 research outputs found
Communications-Inspired Projection Design with Application to Compressive Sensing
We consider the recovery of an underlying signal x \in C^m based on
projection measurements of the form y=Mx+w, where y \in C^l and w is
measurement noise; we are interested in the case l < m. It is assumed that the
signal model p(x) is known, and w CN(w;0,S_w), for known S_W. The objective is
to design a projection matrix M \in C^(l x m) to maximize key
information-theoretic quantities with operational significance, including the
mutual information between the signal and the projections I(x;y) or the Renyi
entropy of the projections h_a(y) (Shannon entropy is a special case). By
capitalizing on explicit characterizations of the gradients of the information
measures with respect to the projections matrix, where we also partially extend
the well-known results of Palomar and Verdu from the mutual information to the
Renyi entropy domain, we unveil the key operations carried out by the optimal
projections designs: mode exposure and mode alignment. Experiments are
considered for the case of compressive sensing (CS) applied to imagery. In this
context, we provide a demonstration of the performance improvement possible
through the application of the novel projection designs in relation to
conventional ones, as well as justification for a fast online projections
design method with which state-of-the-art adaptive CS signal recovery is
achieved.Comment: 25 pages, 7 figures, parts of material published in IEEE ICASSP 2012,
submitted to SIIM
Heterogeneous multireference alignment: a single pass approach
Multireference alignment (MRA) is the problem of estimating a signal from
many noisy and cyclically shifted copies of itself. In this paper, we consider
an extension called heterogeneous MRA, where signals must be estimated, and
each observation comes from one of those signals, unknown to us. This is a
simplified model for the heterogeneity problem notably arising in cryo-electron
microscopy. We propose an algorithm which estimates the signals without
estimating either the shifts or the classes of the observations. It requires
only one pass over the data and is based on low-order moments that are
invariant under cyclic shifts. Given sufficiently many measurements, one can
estimate these invariant features averaged over the signals. We then design
a smooth, non-convex optimization problem to compute a set of signals which are
consistent with the estimated averaged features. We find that, in many cases,
the proposed approach estimates the set of signals accurately despite
non-convexity, and conjecture the number of signals that can be resolved as
a function of the signal length is on the order of .Comment: 6 pages, 3 figure
Bispectrum Inversion with Application to Multireference Alignment
We consider the problem of estimating a signal from noisy
circularly-translated versions of itself, called multireference alignment
(MRA). One natural approach to MRA could be to estimate the shifts of the
observations first, and infer the signal by aligning and averaging the data. In
contrast, we consider a method based on estimating the signal directly, using
features of the signal that are invariant under translations. Specifically, we
estimate the power spectrum and the bispectrum of the signal from the
observations. Under mild assumptions, these invariant features contain enough
information to infer the signal. In particular, the bispectrum can be used to
estimate the Fourier phases. To this end, we propose and analyze a few
algorithms. Our main methods consist of non-convex optimization over the smooth
manifold of phases. Empirically, in the absence of noise, these non-convex
algorithms appear to converge to the target signal with random initialization.
The algorithms are also robust to noise. We then suggest three additional
methods. These methods are based on frequency marching, semidefinite relaxation
and integer programming. The first two methods provably recover the phases
exactly in the absence of noise. In the high noise level regime, the invariant
features approach for MRA results in stable estimation if the number of
measurements scales like the cube of the noise variance, which is the
information-theoretic rate. Additionally, it requires only one pass over the
data which is important at low signal-to-noise ratio when the number of
observations must be large
The dependence of intrinsic alignment of galaxies on wavelength using KiDS and GAMA
The outer regions of galaxies are more susceptible to the tidal interactions
that lead to intrinsic alignments of galaxies. The resulting alignment signal
may therefore depend on the passband if the colours of galaxies vary spatially.
To quantify this, we measured the shapes of galaxies with spectroscopic
redshifts from the GAMA survey using deep gri imaging data from the KiloDegree
Survey. The performance of the moment-based shape measurement algorithm DEIMOS
was assessed using dedicated image simulations, which showed that the
ellipticities could be determined with an accuracy better than 1% in all bands.
Additional tests for potential systematic errors did not reveal any issues. We
measure a significant difference of the alignment signal between the g,r and
i-band observations. This difference exceeds the amplitude of the linear
alignment model on scales below 2 Mpc/h. Separating the sample into
central/satellite and red/blue galaxies, we find that that the difference is
dominated by red satellite galaxies.Comment: 16 pages, 13 figures, accepted, to appear in A&
- …