An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet

Abstract

Next generation radio-interferometers, like the Square Kilometre Array, will acquire tremendous amounts of data with the goal of improving the size and sensitivity of the reconstructed images by orders of magnitude. The efficient processing of large-scale data sets is of great importance. We propose an acceleration strategy for a recently proposed primal-dual distributed algorithm. A preconditioning approach can incorporate into the algorithmic structure both the sampling density of the measured visibilities and the noise statistics. Using the sampling density information greatly accelerates the convergence speed, especially for highly non-uniform sampling patterns, while relying on the correct noise statistics optimises the sensitivity of the reconstruction. In connection to CLEAN, our approach can be seen as including in the same algorithmic structure both natural and uniform weighting, thereby simultaneously optimising both the resolution and the sensitivity. The method relies on a new non-Euclidean proximity operator for the data fidelity term, that generalises the projection onto the ${\ell }_2$ ball where the noise lives for naturally weighted data, to the projection onto a generalised ellipsoid incorporating sampling density information through uniform weighting. Importantly, this non-Euclidean modification is only an acceleration strategy to solve the convex imaging problem with data fidelity dictated only by noise statistics. We showcase through simulations with realistic sampling patterns the acceleration obtained using the preconditioning. We also investigate the algorithm performance for the reconstruction of the 3C129 radio galaxy from real visibilities and compare with multi-scale CLEAN, showing better sensitivity and resolution. Our MATLAB code is available online on GitHub.

References in corpus (9)

• 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
• 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
• Adaptive Weighting in Radio Interferometric Imaging

Cited by in corpus (20)

• 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
• Image reconstruction algorithms in radio interferometry: from handcrafted to learned regularization denoisers
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• First AI for deep super-resolution wide-field imaging in radio astronomy: unveiling structure in ESO 137--006
• Cygnus A super-resolved via convex optimisation from VLA data
• 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
• Scalable precision wide-field imaging in radio interferometry: I. uSARA validated on ASKAP data
• The w-effect in interferometric imaging: from a fast sparse measurement operator to super-resolution
• Data Multiplexing in Radio Interferometric Calibration
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): II. Code and real data proof of concept
• Statistical Performance of Radio Interferometric Calibration
• PolyCLEAN: Atomic Optimization for Super-Resolution Imaging and Uncertainty Estimation in Radio Interferometry
• How to make CLEAN variants faster? Using clustered components informed by the autocorrelation function
• Toward a Robust R2D2 Paradigm for Radio-interferometric Imaging: Revisiting Deep Neural Network Training and Architecture
• Accelerating the CLEAN algorithm of radio interferometry with convex optimization
• DoG-HiT: A novel VLBI Multiscale Imaging Approach