Revisiting the spread spectrum effect in radio interferometric imaging: a sparse variant of the w-projection algorithm
arXiv:1307.3424 · doi:10.1093/mnras/stt1707
Abstract
Next-generation radio interferometric telescopes will exhibit non-coplanar baseline configurations and wide field-of-views, inducing a w-modulation of the sky image, which in turn induces the spread spectrum effect. We revisit the impact of this effect on imaging quality and study a new algorithmic strategy to deal with the associated operator in the image reconstruction process. In previous studies it has been shown that image recovery in the framework of compressed sensing is improved due to the spread spectrum effect, where the w-modulation can act to increase the incoherence between measurement and sparsifying signal representations. For the purpose of computational efficiency, idealised experiments were performed, where only a constant baseline component w in the pointing direction of the telescope was considered. We extend this analysis to the more realistic setting where the w-component varies for each visibility measurement. Firstly, incorporating varying w-components into imaging algorithms is a computational demanding task. We propose a variant of the w-projection algorithm for this purpose, which is based on an adaptive sparsification procedure, and incorporate it in compressed sensing imaging methods. This sparse matrix variant of the w-projection algorithm is generic and adapts to the support of each kernel. Consequently, it is applicable for all types of direction-dependent effects. Secondly, we show that for w-modulation with varying w-components, reconstruction quality is significantly improved compared to the setting where there is no w-modulation (i.e. w=0), reaching levels comparable to the quality of a constant, maximal w-component. This finding confirms that one may seek to optimise future telescope configurations to promote large w-components, thus enhancing the spread spectrum effect and consequently the fidelity of image reconstruction.
12 pages, 9 figures
References in corpus (12)
- Sparsity and Incoherence in Compressive Sampling
- Multi-Scale CLEAN deconvolution of radio synthesis images
- The non-coplanar baselines effect in radio interferometry: The W-Projection algorithm
- Compressed sensing imaging techniques for radio interferometry
- Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
- Correcting direction-dependent gains in the deconvolution of radio interferometric images
- Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
- PURIFY: a new approach to radio-interferometric imaging
- Applying full polarization A-Projection to very wide field of view instruments: An imager for LOFAR
- Compressed sensing for wide-field radio interferometric imaging
- A generalised Measurement Equation and van Cittert-Zernike theorem for wide-field radio astronomical interferometry
- Simulating full-sky interferometric observations
Cited by in corpus (13)
- PURIFY: a new approach to radio-interferometric imaging
- LOFAR Sparse Image Reconstruction
- Scalable splitting algorithms for big-data interferometric imaging in the SKA era
- Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
- Robust sparse image reconstruction of radio interferometric observations with purify
- Uncertainty quantification for radio interferometric imaging: II. MAP estimation
- Probabilistic image reconstruction for radio interferometers
- A Fourier dimensionality reduction model for big data interferometric imaging
- The w-effect in interferometric imaging: from a fast sparse measurement operator to super-resolution
- Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): II. Code and real data proof of concept
- -stacking -projection hybrid algorithm for wide-field interferometric imaging: implementation details and improvements
- A randomised primal-dual algorithm for distributed radio-interferometric imaging
- PURIFY: a new algorithmic framework for next-generation radio-interferometric imaging