No Arabic abstract
We consider the recovery of square-integrable signals from discrete, equidistant samples of their Gabor transform magnitude and show that, in general, signals can not be recovered from such samples. In particular, we show that for any lattice, one can construct functions in $L^2(mathbb{R})$ which do not agree up to global phase but whose Gabor transform magnitudes sampled on the lattice agree. These functions can be constructed to be either real-valued or complex-valued and have good concentration in both time and frequency.
We establish novel uniqueness results for the Gabor phase retrieval problem: If $mathcal{G} : L^2(mathbb{R}) to L^2(mathbb{R}^2)$ denotes the Gabor transform then every $f in L^4[-c/2,c/2]$ is determined up to a global phase by the values $|mathcal{G}f(x,omega)|$ where $(x,omega)$ are points on the lattice $mathbb{Z} times (2c)^{-1}mathbb{Z}$. This for the first time shows that compactly supported functions can be uniquely reconstructed from lattice samples of their spectrogram. Moreover, by making use of recent developments related to sampling in shift-invariant spaces by Grochenig, Romero and Stockler, we prove analogous uniqueness results for functions in general shift-invariant spaces with Gaussian generator. Generalizations to nonuniform lattices are also presented. The results are based on a combination of certain Muntz-type theorems as well as sampling theorems in shift-invariant spaces.
In recent work [P. Grohs and M. Rathmair. Stable Gabor Phase Retrieval and Spectral Clustering. Communications on Pure and Applied Mathematics (2018)] the instabilities of the Gabor phase retrieval problem, i.e., the problem of reconstructing a function $f$ from its spectrogram $|mathcal{G}f|$, where $$ mathcal{G}f(x,y)=int_{mathbb{R}^d} f(t) e^{-pi|t-x|^2} e^{-2pi i tcdot y} dt, quad x,yin mathbb{R}^d, $$ have been completely classified in terms of the disconnectedness of the spectrogram. These findings, however, were crucially restricted to the onedimensional case ($d=1$) and therefore not relevant for many practical applications. In the present paper we not only generalize the aforementioned results to the multivariate case but also significantly improve on them. Our new results have comprehensive implications in various applications such as ptychography, a highly popular method in coherent diffraction imaging.
The problem of reconstructing a function from the magnitudes of its frame coefficients has recently been shown to be never uniformly stable in infinite-dimensional spaces [5]. This result also holds for frames that are possibly continuous [2]. On the other hand, the problem is always stable in finite-dimensional settings. A prominent example of such a phase retrieval problem is the recovery of a signal from the modulus of its Gabor transform. In this paper, we study Gabor phase retrieval and ask how the stability degrades on a natural family of finite-dimensional subspaces of the signal domain $L^2(mathbb{R})$. We prove that the stability constant scales at least quadratically exponentially in the dimension of the subspaces. Our construction also shows that typical priors such as sparsity or smoothness promoting penalties do not constitute regularization terms for phase retrieval.
We consider the problem of reconstructing the missing phase information from spectrogram data $|mathcal{G} f|,$ with $$ mathcal{G}f(x,y)=int_mathbb{R} f(t) e^{-pi(t-x)^2}e^{-2pi i t y}dt, $$ the Gabor transform of a signal $fin L^2(mathbb{R})$. More specifically, we are interested in domains $Omegasubseteq mathbb{R}^2$, which allow for stable local reconstruction, that is $$ |mathcal{G}g| approx |mathcal{G}f| quad text{in} ~Omega quadLongrightarrow quad exists tauinmathbb{T}:quad mathcal{G}g approx taumathcal{G}f quad text{in} ~Omega. $$ In recent work [P. Grohs and M. Rathmair. Stable Gabor Phase Retrieval and Spectral Clustering. Comm. Pure Appl. Math. (2019)] and [P. Grohs and M. Rathmair. Stable Gabor phase retrieval for multivariate functions. J. Eur. Math. Soc. (2021)] we established a characterization of the stability of this phase retrieval problem in terms of the connectedness of the observed measurements. The main downside of the aforementioned results is that the similarity of two spectrograms is measured w.r.t. a first order weighted Sobolev norm. In this article we remove this flaw and essentially show that the Sobolev norm may be replaced by the $L^2-$norm. Using this result allows us to show that it suffices to sample the spectrogram on suitable discrete sampling sets -- a property of crucial importance for practical applications.
We study the phase reconstruction of signals $f$ belonging to complex Gaussian shift-invariant spaces $V^infty(varphi)$ from spectrogram measurements $|mathcal{G}f(X)|$ where $mathcal{G}$ is the Gabor transform and $X subseteq mathbb{R}^2$. An explicit reconstruction formula will demonstrate that such signals can be recovered from measurements located on parallel lines in the time-frequency plane by means of a Riesz basis expansion. Moreover, connectedness assumptions on $|f|$ result in stability estimates in the situation where one aims to reconstruct $f$ on compact intervals. Driven by a recent observation that signals in Gaussian shift-invariant spaces are determined by lattice measurements [Grohs, P., Liehr, L., Injectivity of Gabor phase retrieval from lattice measurements, arXiv:2008.07238] we prove a sampling result on the stable approximation from finitely many spectrogram samples. The resulting algorithm provides a non-iterative, provably stable and convergent approximation technique. In addition, it constitutes a method of approximating signals in function spaces beyond $V^infty(varphi)$, such as Paley-Wiener spaces.