LOFAR Sparse Image Reconstruction

Abstract

Context. The LOFAR Radio Telescope is a giant digital phased array interferometer with multiple antennas gathered in stations placed throughout Europe. As other interferometers, it provides a discrete set of measured Fourier components of the sky brightness. With these samples, recovering the original brightness distribution with aperture synthesis forms an inverse problem that can be solved by different deconvolution and minimization methods. Aims. Recent papers have established a clear link between the discrete nature of radio interferometry measurement and "compressed sensing" theory, which supports sparse recovery methods to reconstruct an image from the measured visibilities. We aimed at the implementation and at the scientific validation of one of these methods. Methods. We evaluated the photometric and resolution performance of the sparse recovery method in the framework of the LOFAR instrument on simulated and real data. Results. We have implemented a sparse recovery method in the standard LOFAR imaging tools, allowing us to compare the reconstructed images from both simulated and real data with images obtained from classical methods such as CLEAN or MS-CLEAN. Conclusions.We show that i) sparse recovery performs as well as CLEAN in recovering the flux of point sources, ii) performs much better on extended objects (the root mean square error is reduced by a factor up to 10), and iii) provides a solution with an effective angular resolution 2-3 times better than the CLEAN map. Applied to a real LOFAR dataset, the sparse recovery has been validated with the correct photometry and realistic recovered structures of Cygnus A, as compared to other methods. Sparse recovery has been implemented as an image recovery method for the LOFAR Radio Telescope and it can be used for other radio interferometers.