Abstract:In this paper, we propose and investigate the phase retrieval problem with the a priori constraint that the phase is sparse (SPR), which encompasses a number of practical applications, for instance, in characterizing phase-only objects such as microlenses, in phase-contrast microscopy, in optical path difference microscopy, and in Fourier ptychography, where the phase object occupies a small portion of the whole field. The considered problem is strictly more general than the sparse signal recovery problem, which assumes the sparsity of the signal because the sparsity of the signal trivially implies the sparsity of the phase, but the converse is not true. As a result, existing solution algorithms in the literature of sparse signal recovery cannot be applied to SPR and there is an appeal for developing new solution methods for it. In this paper, we propose a new regularization scheme which efficiently captures the sparsity constraint of SPR. The idea behind the proposed approach is to perform a metric projection of the current estimated signal onto the set of all the signals whose phase satisfies the sparsity constraint. The main challenge here is that the latter set is not convex and its associated projector in general does not admit a closed form. One novelty of our analysis is to establish an explicit form of that projector when restricted to those points which are relevant to the solutions of SPR. Note that this result is fundamentally different from the widely known calculation form for projections onto intensity constraint sets. Based on this new result, we propose an efficient solution method, named the sparsity regularization on phase (SROP) algorithm, for the SPR problem in the challenging setting where only one point-spread-function image is given, and we analyze its convergence. The algorithm is the combination of the Gerchberg--Saxton (GS) algorithm with the projection step described above. In view of the GS algorithm being equivalent to the alternating projection for an associated two-set feasibility, the SROP algorithm is shown to be the cyclic projection for an associated three-set feasibility, one of the sets being analyzed in this paper for the first time. Analyzing regularity properties of the involved sets, we obtain convergence results for the SROP algorithm based on our recent convergence theory for the cyclic projection method. Numerical results show clear effectiveness and efficiency of the proposed solution approach for the SPR problem.

A Stochastic Alternating Minimizing Method for Sparse Phase Retrieval

Sample-Efficient Sparse Phase Retrieval Via Stochastic Alternating Minimization

Efficient Recovery of Structured Sparse Signals Via Approximate Message Passing with Structured Spike and Slab Prior

PhaseMax: Stable Guarantees from Noisy Sub-Gaussian Measurements

A Fast and Provable Algorithm for Sparse Phase Retrieval

Phase Retrieval from Fourier Measurements with Masks

Iterative Phase Retrieval with a Sensor Mask.

Phase Retrieval from Quantized Measurements Via Approximate Message Passing

Riemannian Optimization for Phase Retrieval from Masked Fourier Measurements

Phase Retrieval with PhaseLift Algorithm

Large-scale Phase Retrieval Method for Wavefront Reconstruction with Multi-Stage Random Phase Modulation

Sparse Signals Recovery from Noisy Measurements by Orthogonal Matching Pursuit

Subspace Phase Retrieval

GraHTP: A Provable Newton-like Algorithm for Sparse Phase Retrieval

Phase Retrieval with Sparse Phase Constraint

Robust Phase Retrieval by Alternating Minimization

Sparse Affine Sampling: Ambiguity-Free and Efficient Sparse Phase Retrieval

Stochastic IHT With Stochastic Polyak Step-Size for Sparse Signal Recovery

Super Resolution Phase Retrieval for Sparse Signals

Sparse representation utilizing tight frame for phase retrieval

Stochastic Recovery Of Sparse Signals From Random Measurements