Scalable splitting algorithms for big-data interferometric imaging in the SKA era

Abstract

In the context of next generation radio telescopes, like the Square Kilometre Array, the efficient processing of large-scale datasets is extremely important. Convex optimisation tasks under the compressive sensing framework have recently emerged and provide both enhanced image reconstruction quality and scalability to increasingly larger data sets. We focus herein mainly on scalability and propose two new convex optimisation algorithmic structures able to solve the convex optimisation tasks arising in radio-interferometric imaging. They rely on proximal splitting and forward-backward iterations and can be seen, by analogy with the CLEAN major-minor cycle, as running sophisticated CLEAN-like iterations in parallel in multiple data, prior, and image spaces. Both methods support any convex regularisation function, in particular the well studied l1 priors promoting image sparsity in an adequate domain. Tailored for big-data, they employ parallel and distributed computations to achieve scalability, in terms of memory and computational requirements. One of them also exploits randomisation, over data blocks at each iteration, offering further flexibility. We present simulation results showing the feasibility of the proposed methods as well as their advantages compared to state-of-the-art algorithmic solvers. Our Matlab code is available online on GitHub.

Monthly Notices of the Royal Astronomical Society (2016)

References in corpus (10)

• 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
• Discovery of 35 New Supernova Remnants in the Inner Galaxy
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• LOFAR Sparse Image Reconstruction
• Distributed Radio Interferometric Calibration
• MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators. A sparse deconvolution algorithm for radio interferometric imaging
• Compressed sensing for wide-field radio interferometric imaging
• Distributed image reconstruction for very large arrays in radio astronomy

Cited by in corpus (35)

• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• Superresolution Interferometric Imaging with Sparse Modeling Using Total Squared Variation --- Application to Imaging the Black Hole Shadow
• Robust sparse image reconstruction of radio interferometric observations with purify
• Image reconstruction algorithms in radio interferometry: from handcrafted to learned regularization denoisers
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet
• Cygnus A super-resolved via convex optimisation from VLA data
• Non-convex optimization for self-calibration of direction-dependent effects in radio interferometric imaging
• Spherical Bayesian mass-mapping with uncertainties: full sky observations on the celestial sphere
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations
• Wideband Super-resolution Imaging in Radio Interferometry via Low Rankness and Joint Average Sparsity Models (HyperSARA)
• Sparse interferometric Stokes imaging under polarization constraint (Polarized SARA)
• Polca SARA - Full polarization, direction-dependent calibration and sparse imaging for radio interferometry
• 21cm Signal Recovery via the Robust Principle Component Analysis
• Scalable precision wide-field imaging in radio interferometry: I. uSARA validated on ASKAP data
• Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data
• A fast and exact $w$-stacking and $w$-projection hybrid algorithm for wide-field interferometric imaging
• The w-effect in interferometric imaging: from a fast sparse measurement operator to super-resolution
• Regularized Maximum Likelihood Image Synthesis and Validation for ALMA Continuum Observations of Protoplanetary Disks
• 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
• Learned Interferometric Imaging for the SPIDER Instrument
• Statistical Performance of Radio Interferometric Calibration
• Deep network series for large-scale high-dynamic range imaging
• Sparse image reconstruction on the sphere: a general approach with uncertainty quantification
• $w$-stacking $w$-projection hybrid algorithm for wide-field interferometric imaging: implementation details and improvements
• How to make CLEAN variants faster? Using clustered components informed by the autocorrelation function
• Multi-Step Reconstruction of Radio-Interferometric Images
• Distributed and parallel sparse convex optimization for radio interferometry with PURIFY
• Toward a Robust R2D2 Paradigm for Radio-interferometric Imaging: Revisiting Deep Neural Network Training and Architecture
• A Decentralized Framework for Radio-interferometric Image Reconstruction
• Block Distributed Majorize-Minimize Memory Gradient Algorithm and its application to 3D image restoration
• Multi-frequency image reconstruction for radio-interferometry with self-tuned regularization parameters
• DoG-HiT: A novel VLBI Multiscale Imaging Approach
• Bayesian variational regularization on the ball