Compressive Phase Retrieval via Generalized Approximate Message Passing

Abstract

In phase retrieval, the goal is to recover a signal $\mathbf{x}\in {\mathbb{C}}^N$ from the magnitudes of linear measurements $\mathbf{Ax}\in {\mathbb{C}}^M$. While recent theory has established that $M\approx 4N$ intensity measurements are necessary and sufficient to recover generic $\mathbf{x}$, there is great interest in reducing the number of measurements through the exploitation of sparse $\mathbf{x}$, which is known as compressive phase retrieval. In this work, we detail a novel, probabilistic approach to compressive phase retrieval based on the generalized approximate message passing (GAMP) algorithm. We then present a numerical study of the proposed PR-GAMP algorithm, demonstrating its excellent phase-transition behavior, robustness to noise, and runtime. Our experiments suggest that approximately $M\ge 2K{\log }_2(N/K)$ intensity measurements suffice to recover $K$-sparse Bernoulli-Gaussian signals for $\mathbf{A}$ with i.i.d Gaussian entries and $K\ll N$. Meanwhile, when recovering a 6k-sparse 65k-pixel grayscale image from 32k randomly masked and blurred Fourier intensity measurements at 30~dB measurement SNR, PR-GAMP achieved an output SNR of no less than 28~dB in all of 100 random trials, with a median runtime of only 7.3 seconds. Compared to the recently proposed CPRL, sparse-Fienup, and GESPAR algorithms, our experiments suggest that PR-GAMP has a superior phase transition and orders-of-magnitude faster runtimes as the sparsity and problem dimensions increase.

References in corpus (3)

• Wide-field, high-resolution Fourier ptychographic microscopy
• Compressive Phase Retrieval via Generalized Approximate Message Passing
• Probabilistic Reconstruction in Compressed Sensing: Algorithms, Phase Diagrams, and Threshold Achieving Matrices

Cited by in corpus (6)

• Compressive Phase Retrieval via Generalized Approximate Message Passing
• Learning-based real-time method to looking through scattering medium beyond the memory effect
• The estimation error of general first order methods
• A general framework for denoising phaseless diffraction measurements
• Fast Compressive Phase Retrieval from Fourier Measurements
• $L^0$-regularized Variational Methods for Sparse Phase Retrieval