ABSTRACT
We present a new formulation of a family of proximity operators that generalize the projector step for phase retrieval. These proximity operators for noisy intensity measurements can replace the classical "noise-free" projection in any projection-based algorithm. They are derived from a maximum-likelihood formulation and admit closed form solutions for both the Gaussian and the Poisson cases. In addition, we extend these proximity operators to under-sampled intensity measurements. To assess their performance, these operators are exploited in a classical Gerchberg-Saxton algorithm. We present numerical experiments showing that the reconstructed complex amplitudes with these proximity operators always perform better than using the classical intensity projector, while their computational overhead is moderate.
ABSTRACT
Astronomical optical interferometers sample the Fourier transform of the intensity distribution of a source at the observation wavelength. Because of rapid perturbations caused by atmospheric turbulence, the phases of the complex Fourier samples (visibilities) cannot be directly exploited. Consequently, specific image reconstruction methods have been devised in the last few decades. Modern polychromatic optical interferometric instruments are now paving the way to multiwavelength imaging. This paper is devoted to the derivation of a spatiospectral (3D) image reconstruction algorithm, coined Polychromatic opticAl INTErferometric Reconstruction software (PAINTER). The algorithm relies on an iterative process, which alternates estimation of polychromatic images and complex visibilities. The complex visibilities are not only estimated from squared moduli and closure phases, but also differential phases, which helps to better constrain the polychromatic reconstruction. Simulations on synthetic data illustrate the efficiency of the algorithm and, in particular, the relevance of injecting a differential phases model in the reconstruction.