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.

Sparse Phase Retrieval Via Truncated Amplitude Flow

Sparse Phase Retrieval Via Iteratively Reweighted Amplitude Flow

Approximate Message Passing with Nearest Neighbor Sparsity Pattern Learning

A Fast and Provable Algorithm for Sparse Phase Retrieval

Phase Retrieval Via Reweighted Amplitude Flow

Phase Retrieval via Smooth Amplitude Flow

Sampling Complexity on Phase Retrieval from Masked Fourier Measurements Via Wirtinger Flow

Sparse Signal Phase Retrieval for Phaseless Short-Time Fourier Transform Measurement Based on Local Search

Non-propagation Fast Phase Diverse Phase Retrieval for Wavefront Measurement with Generalized FFT-based Basis Function.

PhaseMax: Stable Guarantees from Noisy Sub-Gaussian Measurements

Phase Retrieval from Quantized Measurements Via Approximate Message Passing

Phase Retrieval with Sparse Phase Constraint

Sparse representation utilizing tight frame for phase retrieval

Subspace Phase Retrieval

Phase Retrieval for Signals with Block Sparsity Using BOMP: Algorithms and Recovery Guarantees

Frequency Subspace Amplitude Flow for Phase Retrieval.

Truncated amplitude flow with coded diffraction patterns

A Stochastic Alternating Minimizing Method for Sparse Phase Retrieval

An Alternating Optimization Approach For Phase Retrieval

Phase Retrieval of Block-Sparsity Via Convex Optimization.

Sample-Efficient Sparse Phase Retrieval Via Stochastic Alternating Minimization