Abstract:The phase retrieval problem, where one aims to recover a complex-valued image from far-field intensity measurements, is a classic problem encountered in a range of imaging applications. Modern phase retrieval approaches usually rely on gradient descent methods in a nonlinear minimization framework. Calculating closed-form gradients for use in these methods is tedious work, and formulating second order derivatives is even more laborious. Additionally, second order techniques often require the storage and inversion of large matrices of partial derivatives, with memory requirements that can be prohibitive for data-rich imaging modalities. We use a reverse-mode automatic differentiation (AD) framework to implement an efficient matrix-free version of the Levenberg-Marquardt (LM) algorithm, a longstanding method that finds popular use in nonlinear least-square minimization problems but which has seen little use in phase retrieval. Furthermore, we extend the basic LM algorithm so that it can be applied for general constrained optimization problems beyond just the least-square applications. Since we use AD, we only need to specify the physics-based forward model for a specific imaging application; the derivative terms are calculated automatically through matrix-vector products, without explicitly forming any large Jacobian or Gauss-Newton matrices. We demonstrate that this algorithm can be used to solve both the unconstrained ptychographic object retrieval problem and the constrained "blind" ptychographic object and probe retrieval problems, under both the Gaussian and Poisson noise models, and that this method outperforms best-in-class first-order ptychographic reconstruction methods: it provides excellent convergence guarantees with (in many cases) a superlinear rate of convergence, all with a computational cost comparable to, or lower than, the tested first-order algorithms.

Alternating Direction Methods For Classical And Ptychographic Phase Retrieval

Efficient Image Completion Method Based on Alternating Direction Theory.

Iterative Algorithms for Ptychographic Phase Retrieval

Fast and Robust Phase Retrieval for Masked Coherent Diffractive Imaging

Simultaneous Reconstruction of Phase and Amplitude for Wavefront Measurements Based on Nonlinear Optimization Algorithms

Fast Autofocusing Strategy for Phase Retrieval Based on Statistical Gradient Optimization

Adaptive Selection of Illumination Based on 3d Fourier Ptychography

Application of Vector Phase Retrieval Model in High Numerical Aperture Optical Systems

Iterative Phase Retrieval with a Sensor Mask.

A matrix-free Levenberg-Marquardt algorithm for efficient ptychographic phase retrieval

An adaptive noise-blind-separation algorithm for ptychography

Ptychographic phase-retrieval by proximal algorithms

Cross-iteration Multi-Step Optimization Strategy for Three-Dimensional Intensity Position Correction in Phase Diverse Phase Retrieval.

A Stochastic ADMM Algorithm for Large-Scale Ptychography with Weighted Difference of Anisotropic and Isotropic Total Variation

Extrapolative Phase Retrieval Based on a Hybrid of PhaseCut and Alternating Projection Techniques

Differential Phase Retrieval Based on Phase Modulating for Wavefront Detection

Admissible Measurements and Robust Algorithms for Ptychography

Practical Phase Retrieval Using Double Deep Image Priors

Set projection algorithms for blind ptychographic phase retrieval

Phase Retrieval with One or Two Diffraction Patterns by Alternating Projection with Null Initialization

Two Noise-Robust Axial Scanning Multi-Image Phase Retrieval Algorithms Based on Pauta Criterion and Smoothness Constraint.