Fourier Phase Retrieval: Uniqueness and Algorithms Arxiv:1705.09590V3 [Cs.IT] 6 Nov 2017
Total Page:16
File Type:pdf, Size:1020Kb
Fourier Phase Retrieval: Uniqueness and Algorithms Tamir Bendory∗1, Robert Beinerty2 and Yonina C. Eldarz3 1The Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA 2Institute of Mathematics and Scientific Computing, University of Graz, Heinrichstraße 36, 8010 Graz, Austria 3The Andrew and Erna Viterbi Faculty of Electrical Engineering, Technion - Israel Institute of Technology, Haifa, Israel Abstract The problem of recovering a signal from its phaseless Fourier transform measurements, called Fourier phase retrieval, arises in many applications in engineering and science. Fourier phase retrieval poses fundamental theo- retical and algorithmic challenges. In general, there is no unique mapping between a one-dimensional signal and its Fourier magnitude and therefore the problem is ill-posed. Additionally, while almost all multidimensional signals are uniquely mapped to their Fourier magnitude, the performance of existing algorithms is generally not well-understood. In this chapter we survey methods to guarantee uniqueness in Fourier phase retrieval. We then present different algorithmic approaches to retrieve the signal in practice. We conclude by outlining some of the main open questions in this field. arXiv:1705.09590v3 [cs.IT] 6 Nov 2017 ∗[email protected] [email protected]. The Institute of Mathematics and Scientific Computing is a mem- ber of NAWI Graz (http://www.nawigraz.at). The author is supported by the Austrian Science Fund (FWF) within the project P 28858. [email protected]. The author is supported by the European Unions Horizon 2020 research and innovation program under grant agreement no. 646804-ERC-COG-BNYQ, and from the Israel Science Foundation under Grant no. 335/14. 1 1 Introduction The task of recovering a signal from its Fourier transform magnitude, called Fourier phase retrieval, arises in many areas in engineering and science. The problem has a rich history, tracing back to 1952 [110]. Important ex- amples for Fourier phase retrieval naturally appear in many optical settings since optical sensors, such as a charge-coupled device (CCD) and the hu- man eye, are insensitive to phase information of the light wave. A typi- cal example is coherent diffraction imaging (CDI) which is used in a va- riety of imaging techniques [24, 33, 34, 92, 105, 109]. In CDI, an object is illuminated with a coherent electro-magnetic wave and the far-field in- tensity diffraction pattern is measured. This pattern is proportional to the object’s Fourier transform and therefore the measured data is proportional to its Fourier magnitude. Phase retrieval also played a key role in the de- velopment of the DNA double helix model [55]. This discovery awarded Watson, Crick and Wilkins the Nobel prize in Physiology or Medicine in 1962 [1]. Additional examples for applications in which Fourier phase re- trieval appear are X-ray crystallography, speech recognition, blind channel estimation, astronomy, computational biology, alignment and blind decon- volution [10, 19, 54, 62, 100, 114, 125, 132, 137]. Fourier phase retrieval has been a long-standing problem since it raises difficult challenges. In general, there is no unique mapping between a one- dimensional signal and its Fourier magnitude and therefore the problem is ill-posed. Additionally, while almost all multidimensional signals are uniquely mapped to their Fourier magnitude, the performance and stabil- ity of existing algorithms is generally not well-understood. In particular, it is not clear when given methods recover the true underlying signal. To simplify the mathematical analysis, in recent years attention has been de- voted to a family of related problems, frequently called generalized phase retrieval. This refers to the setting in which the measurements are the phase- less inner products of the signal with known vectors. Particularly, the major- ity of works studied inner products with random vectors. Based on proba- bilistic considerations, a variety of convex and non-convex algorithms were suggested, equipped with stability guarantees from near-optimal number of measurements; see [3, 5, 30, 37, 44, 58, 124, 131, 133] to name a few works along these lines. Here, we focus on the original Fourier phase retrieval problem and study it in detail. We begin by considering the ambiguities of Fourier phase re- trieval [12, 15, 16]. We show that while in general a one-dimensional sig- nal cannot be determined from its Fourier magnitude, there are several ex- ceptional cases, such as minimum phase [66] and sparse signals [73, 101]. For general signals, one can guarantee uniqueness by taking multiple mea- surements, each one with a different mask. This setup is called masked 2 phase retrieval [30, 69] and has several interesting special cases, such as the short-time Fourier transform (STFT) phase retrieval [22, 45, 71] and vecto- rial phase retrieval [82, 102, 103, 104]. For all aforementioned setups, we present algorithms and discuss their properties. We also study the closely– related Frequency-Resolved Optical Gating (FROG) methods [23] and mul- tidimensional Fourier phase retrieval [78]. The outline of this chapter is as follows. In Section 2 we formulate the Fourier phase retrieval problem. We also introduce several of its variants, such as masked Fourier phase retrieval and STFT phase retrieval. In Section 3 we discuss the fundamental problem of uniqueness. Namely, conditions under which there is a unique mapping between a signal and its phaseless measurements. Section 4 is devoted to different algorithmic approaches to recover a signal from its phaseless measurements. Section 5 concludes the chapter and outlines some open questions. We hope that highlighting the gaps in the theory of phase retrieval will motivate more research on these issues. 2 Problem formulation In this section, we formulate the Fourier phase retrieval problem and intro- duce notation. N Let x 2 C be the underlying signal we wish to recover. In Fourier phase retrieval the measurements are given by 2 N−1 −2p jkn=N˜ (2.1) y[k] = ∑ x[n]e ; k = 0;:::;K − 1: n=0 Unless otherwise mentioned, we consider the over-sampled Fourier trans- form, i.e., N˜ = K = 2N −1, since in this case the acquired data is equivalent to the auto-correlation of x as explained in Section 3.1. We refer to this case as the classical phase retrieval problem. As will be discussed in the next section, in general the classical phase retrieval problem is ill-posed. Nev- ertheless, some special structures may impose uniqueness. Two important examples are sparse signals obeying a non-periodic support [73, 101] and minimum phase signals [66]; see Section 3.2. For general signals, a popular method to guarantee a unique mapping between the signal and its phaseless Fourier measurements is by utilizing several masks to introduce redundancy in the acquired data. In this case, the measurements are given by 2 N−1 −2p jkn=N˜ y[m;k] = ∑ x[n]dm[n]e ; k = 0;:::;K − 1; m = 0;:::;M − 1; n=0 (2.2) 3 Figure 2.1: An illustration of a typical masked phase retrieval setup (courtesy of [30]). where dm are M known masks. In matrix notation, this model can be written as ∗ 2 y[m;k] = j fk Dmxj ; k = 0;:::;K − 1; m = 0;:::;M − 1; (2.3) ∗ K×N N×N where fk is the kth row of the DFT matrix F 2 C and Dm 2 C is a diagonal matrix that contains the entries of the mth mask. Classical phase retrieval is a special case in which M = 1 and D0 = IN is the identity matrix. There are several experimental techniques to generate masked Fourier measurements in optical setups [30]. One method is to insert a mask or a phase plate after the object [84]. Another possibility is to modulate the illu- minating beam by an optical grating [85]. A third alternative is oblique il- lumination by illuminating beams hitting the object at specified angles [49]. An illustration of a masked phase retrieval setup is shown in Fig. 2.1. An interesting special case of masked phase retrieval is signal recon- struction from phaseless STFT measurements. Here, all masks are transla- tions of a reference mask, i.e., dm[n] = d[mL−n], where L is a parameter that determines the overlapping factor between adjacent windows. Explicitly, the STFT phase retrieval problem takes on the form 2 N−1 −2p jkn=N˜ y[m;k] = ∑ x[n]d[mL − n]e ; n=0 (2.4) k = 0;:::;K − 1; m = 0;:::;dN=Le − 1: The reference mask d is referred to as STFT window. We denote the length of the STFT window by W, namely, d[n] = 0 for n = W;:::;N − 1 for some W ≤ N. 4 Figure 2.2: An illustration of a conventional ptychography setup (courtesy of [118]). The problem of recovering a signal from its STFT magnitude arises in several applications in optics and speech processing. Particularly, it serves as the model of a popular variant of an ultra-short laser pulse measurement technique called Frequency-Resolved Optical Gating (FROG) which is in- troduced in Section 3.5 (the variant is referred to as X-FROG) [23, 127]. An- other application is ptychography in which a moving probe (pinhole) is used to sense multiple diffraction measurements [87, 90, 106]. An illustration of a conventional ptychography setup is given in Fig. 2.2. A closely-related problem is Fourier ptychography [137]. The next section is devoted to the question of uniqueness. Namely, un- der what conditions on the signal x and the masks dm there exists a unique mapping between the signal and the phaseless measurements. In Section 4 we survey different algorithmic approaches to recover the underlying signal x from the acquired data. 3 Uniqueness guarantees The aim of this section is to survey several approaches to ensure uniqueness of the discrete phase retrieval problem.