Optimal arrays for compressed sensing in snapshot-mode radio interferometry

Stanford University, 450 Serra Mall, Stanford, CA 94305, USA
e-mail: clarafj@stanford.edu

Received: 11 January 2013
Accepted: 19 June 2013

Abstract

Context. Radio interferometry has always faced the problem of incomplete sampling of the Fourier plane. A possible remedy can be found in the promising new theory of compressed sensing (CS), which allows for the accurate recovery of sparse signals from sub-Nyquist sampling given certain measurement conditions.

Aims. We provide an introductory assessment of optimal arrays for CS in snapshot-mode radio interferometry, using orthogonal matching pursuit (OMP), a widely used CS recovery algorithm similar in some respects to CLEAN. We focus on comparing centrally condensed (specifically, Gaussian) arrays to uniform arrays, and randomized arrays to deterministic arrays such as the VLA.

Methods. The theory of CS is grounded in a) sparse representation of signals and b) measurement matrices of low coherence. We calculate the mutual coherence of measurement matrices as a theoretical indicator of arrays’ suitability for OMP, based on the recovery error bounds in Donoho et al. (2006, IEEE Trans. Inform. Theory, 52, 1289). OMP reconstructions of both point and extended objects are also run from simulated incomplete data. Optimal arrays are considered for objects represented in 1) the natural pixel basis and 2) the block discrete cosine transform (BDCT).

Results. We find that reconstructions of the pixel representation perform best with the uniform random array, while reconstructions of the BDCT representation perform best with normal random arrays. Slight randomization to the VLA also improves it dramatically for CS recovery with the pixel basis.

Conclusions. In the pixel basis, array design for CS reflects known principles of array design for small numbers of antennas, namely of randomness and uniform distribution. Differing results with the BDCT, however, emphasize the need to study how sparsifying bases affect array design before CS can be optimized for radio interferometry.