Phase Retrieval and Hilbert Integral Equations – Beyond Minimum-Phase

Abstract

The Fourier transform (spectrum) of a signal is a complex function and is characterized by the magnitude and phase spectra. Phase retrieval is the reconstruction of the phase spectrum from the measurements of the magnitude spectrum. Such problems are encountered in imaging modalities such as X-ray crystallography, frequency-domain optical coherence tomography (FDOCT), quantitative phase microscopy, digital holography, etc., where only the magnitudes of the wavefront are detected by the sensors. The phase retrieval problem is ill-posed in general, since an in nite number of signals can have the same magnitude spectrum. Typical phase retrieval techniques rely on certain prior knowledge about the signal, such as its support or sparsity, to reconstruct the signal. A classical result in phase retrieval is that minimum-phase signals have log-magnitude and phase spectra that satisfy the Hilbert integral equations, thus facilitating exact phase retrieval. In this thesis, we demonstrate that there exist larger classes of signals beyond minimum-phase signals, for which exact phase retrieval is possible. We generalize Hilbert integral equations to 2-D, and also introduce a variant that we call the composite Hilbert transform in the context of 2-D periodic signals. Our first extension pertains to a particular type of parametric modelling of 2-D signals. While 1-D minimum-phase signals have a parametric representation, in terms of poles and zeros, there exists no such 2-D counterpart. We introduce a new class of parametric 2-D signals that possess the exact phase retrieval property, that is, their magnitude spectrum completely characterizes the signal. Starting from the magnitude spectrum, a sequence of non-linear operations lead us to a sum-of-exponentials signal, from which the parameters are computed employing concepts from high-resolution spectral estimation such as the annihilating filter and algebraically coupled matrix-pencil methods. We demonstrate that, for this new class of signals, our method outperforms existing techniques even in the presence of noise. Our second extension is to continuous-domain signals that lie in a principal shift-invariant space spanned by a known basis. Such signals are characterized by the basis combining coefficients. These signals need not be minimum-phase, but certain conditions on the coefficients lead to exact phase retrieval of the continuous-domain signal. In particular, we introduce the concept of causal, delta dominant (CDD) sequences, and show that such signals are characterized by their magnitude spectra. This condition pertains to the time/spatial-domain description of the signal, in contrast to the minimum-phase condition, which is described in the spectral domain. We show that there exist CDD sequences that are not minimum-phase, and vice versa. However, finite-length CDD sequences are always minimum-phase. Our method reconstructs the signal from the magnitude spectrum up to ma-chine precision. We thus have a class of continuous-domain signals that are neither causal nor minimum phase, and yet allow for exact phase retrieval. The shift-invariant structure is applicable to modelling signals encountered in imaging modalities such as FDOCT. We next present an application of 2-D phase retrieval to continuous-domain CDD signals in the context of quantiative phase microscopy. We develop sufficient conditions on the interfering reference wave for exact phase retrieval from magnitude measurements. In particular, we show that when the reference wave is a plane wave with magnitude greater that the intensity of the object wave, and when the carrier frequency is larger than the band-width of the object wave, we can reconstruct the object wave exactly. We demonstrate high-resolution reconstruction of our method on USAF target images. Our final and perhaps the most unifying contribution is in developing Hilbert integral equations for 2-D first-quadrant signals and in introducing the notion of generalized minimum-phase signals for both 1-D and 2-D signals. For 2-D continuous-domain, first-quadrant signals, we establish partial Hilbert transform relations between the real and imaginary parts of the spectrum. In the context of 2-D discrete-domain signals, we show that the partial Hilbert transform does not suffice and introduce the notion of composite Hilbert transform and establish the integral equations. We then introduce four classes of signals (combinations of 1-D/2-D and continuous/discrete-domain) that we call generalized minimum-phase signals, which satisfy corresponding Hilbert integral equations between log-magnitude and phase spectra, hence facilitating exact phase retrieval. This class of generalized minimum-phase signals subsumes the well known class of minimum-phase signals. We further show that, akin to minimum-phase signals, these signals also have stable inverses, which are also generalized minimum-phase signals.