Total Variation-Based Phase Retrieval for Poisson Noise Removal

Phase retrieval plays an important role in vast industrial and scientific applications. We consider a noisy phase retrieval problem in which the magnitudes of the Fourier transform (or a general linear transform) of an underling object are corrupted by Poisson noise, since any optical sensors detect photons, and the number of detected photons follows the Poisson distribution. We propose a variational model for phase retrieval based on a total variation regularization as an image prior and maximum a posteriori estimation of a Poisson noise model, which is referred to as “TV-PoiPR”. We also propose an efficient numerical algorithm based on an alternating direction method of multipliers and establish its convergence. Extensive experiments for coded diffraction, holographic, and ptychographic patterns are conducted using both real- and complex-valued images to demonstrate the effectiveness of our proposed methods.

[1]  Babak Hassibi,et al.  Sparse Phase Retrieval: Uniqueness Guarantees and Recovery Algorithms , 2013, IEEE Transactions on Signal Processing.

[2]  A. Fannjiang,et al.  Phase Retrieval With One or Two Diffraction Patterns by Alternating Projections of the Null Vector , 2015 .

[3]  Yonina C. Eldar,et al.  GESPAR: Efficient Phase Retrieval of Sparse Signals , 2013, IEEE Transactions on Signal Processing.

[4]  Raymond H. Chan,et al.  Multilevel algorithm for a Poisson noise removal model with total-variation regularization , 2007 .

[5]  Yonina C. Eldar,et al.  Phase Retrieval via Matrix Completion , 2011, SIAM Rev..

[6]  R. Balan,et al.  On signal reconstruction without phase , 2006 .

[7]  Albert Fannjiang,et al.  Absolute uniqueness of phase retrieval with random illumination , 2011, ArXiv.

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

[9]  Robert W. Harrison,et al.  Phase problem in crystallography , 1993 .

[10]  Tony F. Chan,et al.  A Novel Sparsity Reconstruction Method from Poisson Data for 3D Bioluminescence Tomography , 2012, J. Sci. Comput..

[11]  Alessandro Foi,et al.  Optimal Inversion of the Anscombe Transformation in Low-Count Poisson Image Denoising , 2011, IEEE Transactions on Image Processing.

[12]  Yonina C. Eldar,et al.  Phase Retrieval: An Overview of Recent Developments , 2015, ArXiv.

[13]  J. Xin,et al.  PhaseLiftOff: an Accurate and Stable Phase Retrieval Method Based on Difference of Trace and Frobenius Norms , 2014, 1406.6761.

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

[15]  Thomas J. Asaki,et al.  A Variational Approach to Reconstructing Images Corrupted by Poisson Noise , 2007, Journal of Mathematical Imaging and Vision.

[16]  Heinz H. Bauschke,et al.  Hybrid projection-reflection method for phase retrieval. , 2003, Journal of the Optical Society of America. A, Optics, image science, and vision.

[17]  Raymond H. Chan,et al.  A Two-Stage Image Segmentation Method for Blurry Images with Poisson or Multiplicative Gamma Noise , 2014, SIAM J. Imaging Sci..

[18]  Peter G. Casazza,et al.  Phase retrieval , 2015, SPIE Optical Engineering + Applications.

[19]  L. G. Sodin,et al.  On the ambiguity of the image reconstruction problem , 1979 .

[20]  Chao Yang,et al.  Efficient Algorithms for Ptychographic Phase Retrieval , 2014 .

[21]  S. Marchesini Publisher's Note: “Invited Article: A unified evaluation of iterative projection algorithms for phase retrieval” [Rev. Sci. Instrum. 78, 011301 (2007)] , 2007 .

[22]  Tom Goldstein,et al.  PhaseMax: Convex Phase Retrieval via Basis Pursuit , 2016, IEEE Transactions on Information Theory.

[23]  J. Miao,et al.  Extending X-ray crystallography to allow the imaging of noncrystalline materials, cells, and single protein complexes. , 2008, Annual review of physical chemistry.

[24]  José M. Bioucas-Dias,et al.  Restoration of Poissonian Images Using Alternating Direction Optimization , 2010, IEEE Transactions on Image Processing.

[25]  F. J. Anscombe,et al.  THE TRANSFORMATION OF POISSON, BINOMIAL AND NEGATIVE-BINOMIAL DATA , 1948 .

[26]  Stephen P. Boyd,et al.  Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers , 2011, Found. Trends Mach. Learn..

[27]  Zhi-Quan Luo,et al.  Convergence analysis of alternating direction method of multipliers for a family of nonconvex problems , 2015, ICASSP.

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

[29]  Zhang Fe Phase retrieval from coded diffraction patterns , 2015 .

[30]  M. Lustig,et al.  Compressed Sensing MRI , 2008, IEEE Signal Processing Magazine.

[31]  Stefano Marchesini,et al.  Phase retrieval and saddle-point optimization. , 2006, Journal of the Optical Society of America. A, Optics, image science, and vision.

[32]  Michael P. Friedlander,et al.  Low-Rank Spectral Optimization via Gauge Duality , 2015, SIAM J. Sci. Comput..

[33]  øöö Blockinø Phase retrieval, error reduction algorithm, and Fienup variants: A view from convex optimization , 2002 .

[34]  Martin Burger,et al.  Primal and Dual Bregman Methods with Application to Optical Nanoscopy , 2011, International Journal of Computer Vision.

[35]  A. Walther The Question of Phase Retrieval in Optics , 1963 .

[36]  Yue M. Lu,et al.  Phase Transitions of Spectral Initialization for High-Dimensional Nonconvex Estimation , 2017, Information and Inference: A Journal of the IMA.

[37]  Heinz H. Bauschke,et al.  Phase retrieval, error reduction algorithm, and Fienup variants: a view from convex optimization. , 2002, Journal of the Optical Society of America. A, Optics, image science, and vision.

[38]  H. H. Rachford,et al.  On the numerical solution of heat conduction problems in two and three space variables , 1956 .

[39]  Pengwen Chen,et al.  Fourier Phase Retrieval with a Single Mask by Douglas-Rachford Algorithm , 2015, Applied and computational harmonic analysis.

[40]  Michael K. Ng,et al.  Phase Retrieval from Incomplete Magnitude Information via Total Variation Regularization , 2016, SIAM J. Sci. Comput..

[41]  S. Marchesini,et al.  Alternating projection, ptychographic imaging and phase synchronization , 2014, 1402.0550.

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

[43]  Prateek Jain,et al.  Phase Retrieval Using Alternating Minimization , 2013, IEEE Transactions on Signal Processing.

[44]  Justin Romberg,et al.  Phase Retrieval Meets Statistical Learning Theory: A Flexible Convex Relaxation , 2016, AISTATS.

[45]  Zhiqiang Xu,et al.  Phase Retrieval for Sparse Signals , 2013, ArXiv.

[46]  J. Salmon,et al.  Poisson noise reduction with non-local PCA , 2012, ICASSP.

[47]  S Marchesini,et al.  Invited article: a [corrected] unified evaluation of iterative projection algorithms for phase retrieval. , 2006, The Review of scientific instruments.

[48]  Emmanuel J. Candès,et al.  PhaseLift: Exact and Stable Signal Recovery from Magnitude Measurements via Convex Programming , 2011, ArXiv.

[49]  Rebecca Willett,et al.  A Data-Dependent Weighted LASSO Under Poisson Noise , 2015, IEEE Transactions on Information Theory.

[50]  Chao Yang,et al.  Alternating direction methods for classical and ptychographic phase retrieval , 2012 .

[51]  Xue-Cheng Tai,et al.  Augmented Lagrangian Method, Dual Methods, and Split Bregman Iteration for ROF, Vectorial TV, and High Order Models , 2010, SIAM J. Imaging Sci..

[52]  D. R. Luke Relaxed averaged alternating reflections for diffraction imaging , 2004, math/0405208.

[53]  Yuxin Chen,et al.  Solving Random Quadratic Systems of Equations Is Nearly as Easy as Solving Linear Systems , 2015, NIPS.

[54]  Yonina C. Eldar,et al.  Phase Retrieval with Application to Optical Imaging: A contemporary overview , 2015, IEEE Signal Processing Magazine.

[55]  James V. Burke,et al.  Optical Wavefront Reconstruction: Theory and Numerical Methods , 2002, SIAM Rev..

[56]  Wotao Yin,et al.  Global Convergence of ADMM in Nonconvex Nonsmooth Optimization , 2015, Journal of Scientific Computing.

[57]  Michael Elad,et al.  Sparsity-Based Poisson Denoising With Dictionary Learning , 2013, IEEE Transactions on Image Processing.

[58]  Veit Elser Phase retrieval by iterated projections. , 2003, Journal of the Optical Society of America. A, Optics, image science, and vision.

[59]  Dan Edidin,et al.  An algebraic characterization of injectivity in phase retrieval , 2013, ArXiv.

[60]  Xiaodong Li,et al.  Phase Retrieval via Wirtinger Flow: Theory and Algorithms , 2014, IEEE Transactions on Information Theory.

[61]  Alexandre d'Aspremont,et al.  Phase recovery, MaxCut and complex semidefinite programming , 2012, Math. Program..

[62]  Raymond H. Chan,et al.  Primal-dual algorithms for total variation based image restoration under Poisson noise , 2016 .

[63]  Martin Burger,et al.  Bregman-EM-TV Methods with Application to Optical Nanoscopy , 2009, SSVM.

[64]  J. Sanz Mathematical Considerations for the Problem of Fourier Transform Phase Retrieval from Magnitude , 1985 .

[65]  Ming Yan,et al.  Fast L1–L2 Minimization via a Proximal Operator , 2016, Journal of Scientific Computing.

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

[67]  Shoham Sabach,et al.  Proximal Heterogeneous Block Implicit-Explicit Method and Application to Blind Ptychographic Diffraction Imaging , 2015, SIAM J. Imaging Sci..