A Fourier dimensionality reduction model for big data interferometric imaging

Abstract

Data dimensionality reduction in radio interferometry can provide savings of computational resources for image reconstruction through reduced memory footprints and lighter computations per iteration, which is important for the scalability of imaging methods to the big data setting of the next-generation telescopes. This article sheds new light on dimensionality reduction from the perspective of compressed sensing theory and studies its interplay with imaging algorithms designed in the context of convex optimization. We propose a post-gridding linear data embedding to the space spanned by the left singular vectors of the measurement operator, providing a dimensionality reduction below image size. This embedding preserves the null space of the measurement operator and hence also its sampling properties as per compressed sensing theory. We show that this can be approximated by first computing the dirty image and then applying a weighted subsampled discrete Fourier transform to obtain the final reduced data vector. This Fourier dimensionality reduction model ensures a fast implementation of the full measurement operator, essential for any iterative image reconstruction method. The proposed reduction also preserves the i.i.d. Gaussian properties of the original measurement noise. For convex optimization-based imaging algorithms, this is key to justify the use of the standard L2-norm as the data fidelity term. Our simulations confirm that this dimensionality reduction approach can be leveraged by convex optimization algorithms with no loss in imaging quality relative to reconstructing the image from the complete visibility data set. Reconstruction results in simulation settings with no direction dependent effects or calibration errors show promising performance of the proposed dimensionality reduction. MATLAB code implementing the proposed reduction method is available on GitHub.

20 pages, 10 figures, accepted in Monthly Notices of the Royal Astronomical Society

References in corpus (12)

• Sparsity and Incoherence in Compressive Sampling
• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• Multi-Scale CLEAN deconvolution of radio synthesis images
• Compressed sensing imaging techniques for radio interferometry
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• LOFAR Sparse Image Reconstruction
• Fast Holographic Deconvolution: a new technique for precision radio interferometry
• MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators. A sparse deconvolution algorithm for radio interferometric imaging
• Compressed sensing for wide-field radio interferometric imaging
• Radio Astronomical Image Formation using Constrained Least Squares and Krylov Subspaces
• Multi-frequency image reconstruction for radio interferometry. A regularized inverse problem approach

Cited by in corpus (5)

• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations
• A parallel & automatically tuned algorithm for multispectral image deconvolution
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): II. Code and real data proof of concept
• Fast algorithms to approximate the position-dependent point spread function responses in radio interferometric wide-field imaging