Abstract

Iterative phase retrieval methods based on the Gerchberg–Saxton (GS) or Fienup algorithm typically show stagnation artifacts even after a large number of iterations. We introduce a complexity parameter $\zeta$ that can be computed directly from the Fourier magnitude data and provides a measure of fluctuations in the desired phase retrieval solution. It is observed that when initiated with a constant or a uniformly random phase map, the complexity of the Fienup solution containing stagnation artifacts stabilizes at a numerical value that is higher than $\zeta$. We propose a modified Fienup algorithm that uses a controlled sparsity-enhancing step such that in every iteration the complexity of the resulting guess solution is explicitly made close to $\zeta$. This approach, which we refer to as complexity-guided phase retrieval, is seen to provide an artifact-free phase retrieval solution within a few hundred iterations. Numerical illustrations are provided for both amplitude as well as phase objects with and without Poisson noise introduced in the Fourier intensity data. The complexity-guidance concept may potentially be combined with a variety of phase retrieval algorithms and can enable several practical applications.

© 2019 Optical Society of America

Complexity-guided Fourier phase retrieval from noisy data

Mansi Butola, Sunaina Rajora, and Kedar Khare
J. Opt. Soc. Am. A 38(4) 488-497 (2021)

Sparsity-assisted solution to the twin image problem in phase retrieval

Charu Gaur, Baranidharan Mohan, and Kedar Khare
J. Opt. Soc. Am. A 32(11) 1922-1927 (2015)

Hybrid projection–reflection method for phase retrieval

Heinz H. Bauschke, Patrick L. Combettes, and D. Russell Luke
J. Opt. Soc. Am. A 20(6) 1025-1034 (2003)

Further study on the convergence property of the hybrid input–output algorithm used for phase retrieval

Hiroaki Takajo, Tohru Takahashi, and Takao Shizuma
J. Opt. Soc. Am. A 16(9) 2163-2168 (1999)

Simple constraint for phase retrieval with high efficiency

Haifeng He
J. Opt. Soc. Am. A 23(3) 550-556 (2006)