407 research outputs found
Undersampled Phase Retrieval with Outliers
We propose a general framework for reconstructing transform-sparse images
from undersampled (squared)-magnitude data corrupted with outliers. This
framework is implemented using a multi-layered approach, combining multiple
initializations (to address the nonconvexity of the phase retrieval problem),
repeated minimization of a convex majorizer (surrogate for a nonconvex
objective function), and iterative optimization using the alternating
directions method of multipliers. Exploiting the generality of this framework,
we investigate using a Laplace measurement noise model better adapted to
outliers present in the data than the conventional Gaussian noise model. Using
simulations, we explore the sensitivity of the method to both the
regularization and penalty parameters. We include 1D Monte Carlo and 2D image
reconstruction comparisons with alternative phase retrieval algorithms. The
results suggest the proposed method, with the Laplace noise model, both
increases the likelihood of correct support recovery and reduces the mean
squared error from measurements containing outliers. We also describe exciting
extensions made possible by the generality of the proposed framework, including
regularization using analysis-form sparsity priors that are incompatible with
many existing approaches.Comment: 11 pages, 9 figure
Robust Subspace Learning: Robust PCA, Robust Subspace Tracking, and Robust Subspace Recovery
PCA is one of the most widely used dimension reduction techniques. A related
easier problem is "subspace learning" or "subspace estimation". Given
relatively clean data, both are easily solved via singular value decomposition
(SVD). The problem of subspace learning or PCA in the presence of outliers is
called robust subspace learning or robust PCA (RPCA). For long data sequences,
if one tries to use a single lower dimensional subspace to represent the data,
the required subspace dimension may end up being quite large. For such data, a
better model is to assume that it lies in a low-dimensional subspace that can
change over time, albeit gradually. The problem of tracking such data (and the
subspaces) while being robust to outliers is called robust subspace tracking
(RST). This article provides a magazine-style overview of the entire field of
robust subspace learning and tracking. In particular solutions for three
problems are discussed in detail: RPCA via sparse+low-rank matrix decomposition
(S+LR), RST via S+LR, and "robust subspace recovery (RSR)". RSR assumes that an
entire data vector is either an outlier or an inlier. The S+LR formulation
instead assumes that outliers occur on only a few data vector indices and hence
are well modeled as sparse corruptions.Comment: To appear, IEEE Signal Processing Magazine, July 201
Retrieving shallow shear-wave velocity profiles from 2D seismic-reflection data with severely aliased surface waves
The inversion of surface-wave phase-velocity dispersion curves provides a reliable method to derive near-surface shear-wave velocity profiles. In this work, we invert phase-velocity dispersion curves estimated from 2D seismic-reflection data. These data cannot be used to image the first 50 m with seismic-reflection processing techniques due to the presence of indistinct first breaks and significant NMO-stretching of the shallow reflections. A surface-wave analysis was proposed to derive information about the near surface in order to complement the seismic-reflection stacked sections, which are satisfactory for depths between 50 and 700 m. In order to perform the analysis, we had to overcome some problems, such as the short acquisition time and the large receiver spacing, which resulted in severe spatial aliasing. The analysis consists of spatial partitioning of each line in segments, picking of the phase-velocity dispersion curves for each segment in the f-k domain, and inversion of the picked curves using the neighborhood algorithm. The spatial aliasing is successfully circumvented by continuously tracking the surface-wave modal curves in the f-k domain. This enables us to sample the curves up to a frequency of 40 Hz, even though most components beyond 10 Hz are spatially aliased. The inverted 2D VS sections feature smooth horizontal layers, and a sensitivity analysis yields a penetration depth of 20–25 m. The results suggest that long profiles may be more efficiently surveyed by using a large receiver separation and dealing with the spatial aliasing in the described way, rather than ensuring that no spatially aliased surface waves are acquired.Fil: Onnis, Luciano Emanuel. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Física; ArgentinaFil: Osella, Ana Maria. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Física; ArgentinaFil: Carcione, Jose M.. Istituto Nazionale di Oceanografia e di Geofisica Sperimentale; Itali
Combining ptychographical algorithms with the Hybrid Input-Output (HIO) algorithm
AbstractIn this article we combine the well-known Ptychographical Iterative Engine (PIE) with the Hybrid Input-Output (HIO) algorithm. The important insight is that the HIO feedback function should be kept strictly separate from the reconstructed object, which is done by introducing a separate feedback function per probe position. We have also combined HIO with floating PIE (fPIE) and extended PIE (ePIE). Simulations indicate that the combined algorithm performs significantly better in many situations. Although we have limited our research to a combination with HIO, the same insight can be used to combine ptychographical algorithms with any phase retrieval algorithm that uses a feedback function
- …