42 research outputs found
A Cheeger Inequality for the Graph Connection Laplacian
The O(d) Synchronization problem consists of estimating a set of unknown
orthogonal transformations O_i from noisy measurements of a subset of the
pairwise ratios O_iO_j^{-1}. We formulate and prove a Cheeger-type inequality
that relates a measure of how well it is possible to solve the O(d)
synchronization problem with the spectra of an operator, the graph Connection
Laplacian. We also show how this inequality provides a worst case performance
guarantee for a spectral method to solve this problem.Comment: To appear in the SIAM Journal on Matrix Analysis and Applications
(SIMAX
Orthogonal Matrix Retrieval in Cryo-Electron Microscopy
In single particle reconstruction (SPR) from cryo-electron microscopy
(cryo-EM), the 3D structure of a molecule needs to be determined from its 2D
projection images taken at unknown viewing directions. Zvi Kam showed already
in 1980 that the autocorrelation function of the 3D molecule over the rotation
group SO(3) can be estimated from 2D projection images whose viewing directions
are uniformly distributed over the sphere. The autocorrelation function
determines the expansion coefficients of the 3D molecule in spherical harmonics
up to an orthogonal matrix of size for each
. In this paper we show how techniques for solving the phase
retrieval problem in X-ray crystallography can be modified for the cryo-EM
setup for retrieving the missing orthogonal matrices. Specifically, we present
two new approaches that we term Orthogonal Extension and Orthogonal
Replacement, in which the main algorithmic components are the singular value
decomposition and semidefinite programming. We demonstrate the utility of these
approaches through numerical experiments on simulated data.Comment: Modified introduction and summary. Accepted to the IEEE International
Symposium on Biomedical Imagin
Phase retrieval from power spectra of masked signals
In diffraction imaging, one is tasked with reconstructing a signal from its
power spectrum. To resolve the ambiguity in this inverse problem, one might
invoke prior knowledge about the signal, but phase retrieval algorithms in this
vein have found limited success. One alternative is to create redundancy in the
measurement process by illuminating the signal multiple times, distorting the
signal each time with a different mask. Despite several recent advances in
phase retrieval, the community has yet to construct an ensemble of masks which
uniquely determines all signals and admits an efficient reconstruction
algorithm. In this paper, we leverage the recently proposed polarization method
to construct such an ensemble. We also present numerical simulations to
illustrate the stability of the polarization method in this setting. In
comparison to a state-of-the-art phase retrieval algorithm known as PhaseLift,
we find that polarization is much faster with comparable stability.Comment: 18 pages, 3 figure
Stable optimizationless recovery from phaseless linear measurements
We address the problem of recovering an n-vector from m linear measurements
lacking sign or phase information. We show that lifting and semidefinite
relaxation suffice by themselves for stable recovery in the setting of m = O(n
log n) random sensing vectors, with high probability. The recovery method is
optimizationless in the sense that trace minimization in the PhaseLift
procedure is unnecessary. That is, PhaseLift reduces to a feasibility problem.
The optimizationless perspective allows for a Douglas-Rachford numerical
algorithm that is unavailable for PhaseLift. This method exhibits linear
convergence with a favorable convergence rate and without any parameter tuning