Experimental robustness of Fourier ptychography phase retrieval algorithms

Abstract

Fourier ptychography is a new computational microscopy technique that provides gigapixel-scale intensity and phase images with both wide field-of-view and high resolution. By capturing a stack of low-resolution images under different illumination angles, an inverse algorithm can be used to computationally reconstruct the high-resolution complex field. Here, we compare and classify multiple proposed inverse algorithms in terms of experimental robustness. We find that the main sources of error are noise, aberrations and mis-calibration (i.e. model mis-match). Using simulations and experiments, we demonstrate that the choice of cost function plays a critical role, with amplitude-based cost functions performing better than intensity-based ones. The reason for this is that Fourier ptychography datasets consist of images from both brightfield and darkfield illumination, representing a large range of measured intensities. Both noise (e.g. Poisson noise) and model mis-match errors are shown to scale with intensity. Hence, algorithms that use an appropriate cost function will be more tolerant to both noise and model mis-match. Given these insights, we propose a global Newton's method algorithm which is robust and accurate. Finally, we discuss the impact of procedures for algorithmic correction of aberrations and mis-calibration.

Fourier ptychographic reconstruction using Poisson maximum likelihood and truncated Wirtinger gradient

“…The method guarantees a global optimum, but converges slow which makes it impractical in real applications. Recently, Yeh et al 18. tested different objective functions (intensity based, amplitude based and Poisson maximum likelihood) under the gradient-descent optimization scheme for FP reconstruction.…”

Diffraction tomography with Fourier ptychography

“…We used the well-known alternating projections phase retrieval update. Other solvers based upon convex optimization [55] or alternative gradient descent techniques [56,57] may perform better in the presence of noise. Alternative approximations besides first Born approximation (e.g., Rytov [15]) are also available to simplify the Born series.…”

Single‐shot speckle correlation fluorescence microscopy in thick scattering tissue with image reconstruction priors

“…Perhaps, a joint recovery of both Fourier magnitude and the object itself would further improve reconstruction quality. There is also much to learn from other imaging problems using phase retrieval algorithms .…”