407 research outputs found

    Undersampled Phase Retrieval with Outliers

    Full text link
    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

    Full text link
    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

    Get PDF
    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

    Get PDF
    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
    corecore