Phase Retrieval for Sparse Signals: Uniqueness Conditions

In a variety of fields, in particular those involving imaging and optics, we often measure signals whose phase is missing or has been irremediably distorted. Phase retrieval attempts the recovery of the phase information of a signal from the magnitude of its Fourier transform to enable the reconstruction of the original signal. A fundamental question then is: "Under which conditions can we uniquely recover the signal of interest from its measured magnitudes?" In this paper, we assume the measured signal to be sparse. This is a natural assumption in many applications, such as X-ray crystallography, speckle imaging and blind channel estimation. In this work, we derive a sufficient condition for the uniqueness of the solution of the phase retrieval (PR) problem for both discrete and continuous domains, and for one and multi-dimensional domains. More precisely, we show that there is a strong connection between PR and the turnpike problem, a classic combinatorial problem. We also prove that the existence of collisions in the autocorrelation function of the signal may preclude the uniqueness of the solution of PR. Then, assuming the absence of collisions, we prove that the solution is almost surely unique on 1-dimensional domains. Finally, we extend this result to multi-dimensional signals by solving a set of 1-dimensional problems. We show that the solution of the multi-dimensional problem is unique when the autocorrelation function has no collisions, significantly improving upon a previously known result.

[1]  Solomon W. Golomb,et al.  There Are No Further Counterexamples to S. Piccard's Theorem , 2007, IEEE Transactions on Information Theory.

[2]  Warren D. Smith,et al.  Reconstructing Sets From Interpoint Distances , 2003 .

[3]  Marjorie Senechal A point set puzzle revisited , 2008, Eur. J. Comb..

[4]  Michael C. Hout,et al.  Multidimensional Scaling , 2003, Encyclopedic Dictionary of Archaeology.

[5]  David P. Williamson,et al.  Improved approximation algorithms for maximum cut and satisfiability problems using semidefinite programming , 1995, JACM.

[6]  Bernhard Schölkopf,et al.  Online Multi-frame Blind Deconvolution with Super-resolution and Saturation Correction , 2011 .

[7]  David Sayre X-Ray Crystallography: The Past and Present of the Phase Problem , 2002 .

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

[9]  Zheng Zhang An Exponential Example for a Partial Digest Mapping Algorithm , 1994, J. Comput. Biol..

[10]  Peter G. Casazza,et al.  Equivalence of Reconstruction From the Absolute Value of the Frame Coefficients to a Sparse Representation Problem , 2007, IEEE Signal Processing Letters.

[11]  Gábor Oszlányi,et al.  The charge flipping algorithm. , 2008, Acta crystallographica. Section A, Foundations of crystallography.

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

[13]  Martin Vetterli,et al.  Sparse spectral factorization: Unicity and reconstruction algorithms , 2011, 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[14]  Bernhard G. Bodmann,et al.  Fast algorithms for signal reconstruction without phase , 2007, SPIE Optical Engineering + Applications.

[15]  Ali H. Sayed,et al.  A survey of spectral factorization methods , 2001, Numer. Linear Algebra Appl..

[16]  G. Oszlányi,et al.  Ab initio structure solution by charge flipping. , 2003, Acta crystallographica. Section A, Foundations of crystallography.

[17]  L. Ciric,et al.  A Counterexample to a Theorem of Xu , 1997 .

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

[19]  Ramjee Prasad,et al.  OFDM for Wireless Multimedia Communications , 1999 .

[20]  Babak Hassibi,et al.  Recovery of sparse 1-D signals from the magnitudes of their Fourier transform , 2012, 2012 IEEE International Symposium on Information Theory Proceedings.

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

[22]  Paolo Penna,et al.  Noisy Data Make the Partial Digest Problem NP-hard , 2003, WABI.

[23]  J. Y. Wang,et al.  Modal compensation of atmospheric turbulence phase distortion , 1978 .

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

[25]  Yonina C. Eldar,et al.  Phase Retrieval via Matrix Completion , 2013, SIAM J. Imaging Sci..

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

[27]  A. Oppenheim,et al.  Signal reconstruction from signed Fourier transform magnitude , 1983 .

[28]  Gary S. Bloom,et al.  A Counterexample to a Theorem of S. Piccard , 1977, J. Comb. Theory, Ser. A.

[29]  David Blow,et al.  Direct Methods In Crystallography , 1962 .

[30]  Dan Edidin,et al.  On signal reconstruction from absolute value of frame coefficients , 2005, SPIE Optics + Photonics.

[31]  Sophie Piccard,et al.  Sur les ensembles de distances des emsembles de points d'un espace euclidien , 1940, The Mathematical Gazette.

[32]  Pier Luigi Dragotti,et al.  Tomographic Approach for Sampling Multidimensional Signals with Finite Rate of Innovation , 2006, 2006 International Conference on Image Processing.

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

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

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

[36]  A. L. Patterson A Direct Method for the Determination of the Components of Interatomic Distances in Crystals , 1935 .

[37]  Martin Vetterli,et al.  Frame domain signal processing: Framework and applications , 2010, 2010 IEEE International Conference on Acoustics, Speech and Signal Processing.

[38]  Michael P. Fitz,et al.  Blind estimation of multipath channel parameters: a modal analysis approach , 1999, IEEE Trans. Commun..

[39]  Steven Skiena,et al.  A partial digest approach to restriction site mapping , 1993, ISMB.

[40]  Martin Vetterli,et al.  Fast and Robust Parametric Estimation of Jointly Sparse Channels , 2012, IEEE Journal on Emerging and Selected Topics in Circuits and Systems.