Optical interferometers provide multiple wavelength measurements. In order to fully exploit the spectral and spatial resolution of these instruments, new algorithms for image reconstruction have to be developed. Early attempts to deal with multichromatic interferometric data have consisted in recovering a gray image of the object or independent monochromatic images in some spectral bandwidths. The main challenge is now to recover the full three-dimensional (spatiospectral) brightness distribution of the astronomical target given all the available data. We describe an approach to implement multiwavelength image reconstruction in the case where the observed scene is a collection of point-like sources. We show the gain in image quality (both spatially and spectrally) achieved by globally taking into account all the data instead of dealing with independent spectral slices. This is achieved thanks to a regularization that favors spatial sparsity and spectral grouping of the sources. Since the objective function is not differentiable, we had to develop a specialized optimization algorithm that also accounts for non-negativity of the brightness distribution.
© 2013 Optical Society of America
Original Manuscript: September 12, 2012
Revised Manuscript: November 30, 2012
Manuscript Accepted: December 3, 2012
Published: January 8, 2013
Éric Thiébaut, Ferréol Soulez, and Loïc Denis, "Exploiting spatial sparsity for multiwavelength imaging in optical interferometry," J. Opt. Soc. Am. A 30, 160-170 (2013)