Numerical methods for phase retrieval

In this work we consider the problem of reconstruction of a signal from the magnitude of its Fourier transform, also known as phase retrieval. The problem arises in many areas of astronomy, crystallography, optics, and coherent diffraction imaging (CDI). Our main goal is to develop an efficient reconstruction method based on continuous optimization techniques. Unlike current reconstruction methods, which are based on alternating projections, our approach leads to a much faster and more robust method. However, all previous attempts to employ continuous optimization methods, such as Newton-type algorithms, to the phase retrieval problem failed. In this work we provide an explanation for this failure, and based on this explanation we devise a sufficient condition that allows development of new reconstruction methods---approximately known Fourier phase. We demonstrate that a rough (up to $\pi/2$ radians) Fourier phase estimate practically guarantees successful reconstruction by any reasonable method. We also present a new reconstruction method whose reconstruction time is orders of magnitude faster than that of the current method-of-choice in phase retrieval---Hybrid Input-Output (HIO). Moreover, our method is capable of successful reconstruction even in the situations where HIO is known to fail. We also extended our method to other applications: Fourier domain holography, and interferometry. Additionally we developed a new sparsity-based method for sub-wavelength CDI. Using this method we demonstrated experimental resolution exceeding several times the physical limit imposed by the diffraction light properties (so called diffraction limit).

[1]  Thomas F. Quatieri,et al.  Recursive phase retrieval using boundary conditions , 1983 .

[2]  Yonina C. Eldar,et al.  Sparsity-Based Reconstruction of SubwavelengthImages from their Optical Far-Field , 2010 .

[3]  R. J.,et al.  Reconstruction of a complex-valued object from the modulus of its Fourier transform using a support constraint , 2002 .

[4]  J. Goodman Introduction to Fourier optics , 1969 .

[5]  B. A. D. H. Brandwood A complex gradient operator and its applica-tion in adaptive array theory , 1983 .

[6]  Yonina C. Eldar,et al.  Beyond bandlimited sampling , 2009, IEEE Signal Processing Magazine.

[7]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[8]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[9]  Alexander E. Kaplan,et al.  Optical physics (A) , 1986 .

[10]  Emmanuel J. Candès,et al.  Near-Optimal Signal Recovery From Random Projections: Universal Encoding Strategies? , 2004, IEEE Transactions on Information Theory.

[11]  Richard A. London,et al.  Femtosecond time-delay X-ray holography , 2007, Nature.

[12]  Irad Yavneh,et al.  Fast Reconstruction Method for Diffraction Imaging , 2009, ISVC.

[13]  A. Bruckstein,et al.  K-SVD : An Algorithm for Designing of Overcomplete Dictionaries for Sparse Representation , 2005 .

[14]  J. L. Harris,et al.  Diffraction and Resolving Power , 1964 .

[15]  G. R. Ringo,et al.  Ion source of high brightness using liquid metal , 1975 .

[16]  Irad Yavneh,et al.  IMAGE RECONSTRUCTION FROM NOISY FOURIER MAGNITUDE WITH PARTIAL PHASE INFORMATION , 2009 .

[17]  Yonina C. Eldar,et al.  Far-Field Microscopy of Sparse Subwavelength Objects , 2010, 1010.0631.

[18]  Paul R. Selvin,et al.  Myosin V Walks Hand-Over-Hand: Single Fluorophore Imaging with 1.5-nm Localization , 2003, Science.

[19]  I. Yavneh,et al.  IMG5 Annual Report|Numerical Solution of Inverse Problems in Optics: Phase Retrieval, Holography, Deblurring, Image Reconstruction From Its Defocused Versions, And Combinations Thereof , 2010 .

[20]  A. Bos Complex gradient and Hessian , 1994 .

[21]  Thomas F. Quatieri,et al.  The importance of boundary conditions in the phase retrieval problem , 1982, ICASSP.

[22]  T. D. Harris,et al.  Breaking the Diffraction Barrier: Optical Microscopy on a Nanometric Scale , 1991, Science.

[23]  J. Magnus,et al.  Matrix Differential Calculus with Applications in Statistics and Econometrics , 1991 .

[24]  Yonina C. Eldar,et al.  Blind Compressed Sensing , 2010, IEEE Transactions on Information Theory.

[25]  Nikolay I Zheludev,et al.  Super-resolution without evanescent waves. , 2008, Nano letters.

[26]  Georg Weidenspointner,et al.  Femtosecond X-ray protein nanocrystallography , 2011, Nature.

[27]  A. Papoulis A new algorithm in spectral analysis and band-limited extrapolation. , 1975 .

[28]  Yonina C. Eldar,et al.  Structured Compressed Sensing: From Theory to Applications , 2011, IEEE Transactions on Signal Processing.

[29]  S. Marchesini,et al.  High-resolution ab initio three-dimensional x-ray diffraction microscopy. , 2005, Journal of the Optical Society of America. A, Optics, image science, and vision.

[30]  R. Millane Multidimensional phase problems , 1996 .

[31]  Thierry Blu,et al.  Sampling signals with finite rate of innovation , 2002, IEEE Trans. Signal Process..

[32]  S. Hädrich,et al.  Lensless diffractive imaging using tabletop coherent high-harmonic soft-X-ray beams. , 2007, Physical review letters.

[33]  Michael A. Saunders,et al.  Atomic Decomposition by Basis Pursuit , 1998, SIAM J. Sci. Comput..

[34]  J. Hajdu,et al.  Potential for biomolecular imaging with femtosecond X-ray pulses , 2000, Nature.

[35]  Yonina C. Eldar,et al.  Super-resolution and reconstruction of sparse sub-wavelength images: erratum , 2010 .

[36]  Microsc. Microanal,et al.  Coherent Diffraction Imaging , 2014 .

[37]  M A Fiddy,et al.  Enforcing irreducibility for phase retrieval in two dimensions. , 1983, Optics letters.

[38]  I. Yavneh,et al.  Signal Reconstruction From The Modulus of its Fourier Transform , 2009 .

[39]  Garth J. Williams,et al.  Three-dimensional mapping of a deformation field inside a nanocrystal , 2006, Nature.

[40]  R A Linke,et al.  Beaming Light from a Subwavelength Aperture , 2002, Science.

[41]  W. Wirtinger Zur formalen Theorie der Funktionen von mehr komplexen Veränderlichen , 1927 .

[42]  Irad Yavneh,et al.  Simultaneous deconvolution and phase retrieval from noisy data , 2010 .

[43]  S. Hell,et al.  Diffraction-unlimited three-dimensional optical nanoscopy with opposing lenses , 2009 .

[44]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[45]  Irad Yavneh,et al.  Approximate Fourier phase information in the phase retrieval problem: what it gives and how to use it. , 2011, Journal of the Optical Society of America. A, Optics, image science, and vision.

[46]  Germund Dahlquist,et al.  Numerical Methods in Scientific Computing: Volume 1 , 2008 .

[47]  R. J.,et al.  Reconstruction of objects having latent reference points , 2004 .

[48]  A. Oppenheim,et al.  Signal reconstruction from phase or magnitude , 1980 .

[49]  P. Prewett,et al.  CORRIGENDUM: Characteristics of a gallium liquid metal field emission ion source , 1980 .

[50]  H. Quiney Coherent diffractive imaging using short wavelength light sources , 2010 .

[51]  Yonina C. Eldar,et al.  Super-resolution and reconstruction of sparse images carried by incoherent light. , 2010, Optics letters.

[52]  J. Zuo,et al.  Atomic Resolution Imaging of a Carbon Nanotube from Diffraction Intensities , 2003, Science.

[53]  R. Gerchberg A practical algorithm for the determination of phase from image and diffraction plane pictures , 1972 .

[54]  J. Miao,et al.  An approach to three-dimensional structures of biomolecules by using single-molecule diffraction images , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[55]  M. Nieto-Vesperinas A Study of the Performance of Nonlinear Least-square Optimization Methods in the Problem of Phase Retrieval , 1986 .

[56]  Irad Yavneh,et al.  Algorithms for phase retrieval with a (rough) phase estimate available: a comparison , 2010 .

[57]  G. Toraldo di Francia,et al.  Super-gain antennas and optical resolving power , 1952 .

[58]  A.V. Oppenheim,et al.  The importance of phase in signals , 1980, Proceedings of the IEEE.

[59]  EDWARD M. HOFSTETTER,et al.  Construction of time-limited functions with specified autocorrelation functions , 1964, IEEE Trans. Inf. Theory.

[60]  J. Miao,et al.  Extending the methodology of X-ray crystallography to allow imaging of micrometre-sized non-crystalline specimens , 1999, Nature.

[61]  Yonina C. Eldar,et al.  Sparsity Based Sub-wavelength Imaging with Partially Incoherent Light via Quadratic Compressed Sensing References and Links , 2022 .

[62]  M. Isaacson,et al.  Development of a 500 Å spatial resolution light microscope: I. light is efficiently transmitted through λ/16 diameter apertures , 1984 .

[63]  Rick P. Millane,et al.  Phase retrieval in crystallography and optics , 1990 .

[64]  M. Zibulevsky,et al.  Sequential Subspace Optimization Method for Large-Scale Unconstrained Problems , 2005 .

[65]  James R. Fienup,et al.  Phase-retrieval stagnation problems and solutions , 1986 .

[66]  M. Hayes The reconstruction of a multidimensional sequence from the phase or magnitude of its Fourier transform , 1982 .

[67]  E.J. Candes,et al.  An Introduction To Compressive Sampling , 2008, IEEE Signal Processing Magazine.

[68]  S. D. Babacan,et al.  Cell imaging beyond the diffraction limit using sparse deconvolution spatial light interference microscopy , 2011, Biomedical optics express.

[69]  D. Sayre Some implications of a theorem due to Shannon , 1952 .

[70]  R. Gerchberg Super-resolution through Error Energy Reduction , 1974 .

[71]  M. Elad,et al.  $rm K$-SVD: An Algorithm for Designing Overcomplete Dictionaries for Sparse Representation , 2006, IEEE Transactions on Signal Processing.

[72]  Keith A. Nugent,et al.  Coherent lensless X-ray imaging , 2010 .

[73]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[74]  J R Fienup,et al.  Phase retrieval algorithms: a comparison. , 1982, Applied optics.

[75]  A. Yagle,et al.  Avoiding phase-retrieval algorithm stagnation using the zeros of the Fourier magnitude , 1989, Sixth Multidimensional Signal Processing Workshop,.

[76]  James R. Fienup,et al.  Iterative Method Applied To Image Reconstruction And To Computer-Generated Holograms , 1979, Optics & Photonics.

[77]  Stéphane Mallat,et al.  Matching pursuits with time-frequency dictionaries , 1993, IEEE Trans. Signal Process..

[78]  R. J.,et al.  Phase retrieval using boundary conditions , 2002 .

[79]  Yonina C Eldar,et al.  Super-resolution and reconstruction of sparse sub-wavelength images. , 2009, Optics express.