Wideband Super-resolution Imaging in Radio Interferometry via Low Rankness and Joint Average Sparsity Models (HyperSARA)

Abstract

We propose a new approach within the versatile framework of convex optimization to solve the radio-interferometric wideband imaging problem. Our approach, dubbed HyperSARA, solves a sequence of weighted nuclear norm and l21 minimization problems promoting low rankness and joint average sparsity of the wideband model cube. On the one hand, enforcing low rankness enhances the overall resolution of the reconstructed model cube by exploiting the correlation between the different channels. On the other hand, promoting joint average sparsity improves the overall sensitivity by rejecting artefacts present on the different channels. An adaptive Preconditioned Primal-Dual algorithm is adopted to solve the minimization problem. The algorithmic structure is highly scalable to large data sets and allows for imaging in the presence of unknown noise levels and calibration errors. We showcase the superior performance of the proposed approach, reflected in high-resolution images on simulations and real VLA observations with respect to single channel imaging and the CLEAN-based wideband imaging algorithm in the WSCLEAN software. Our MATLAB code is available online on GITHUB.

References in corpus (25)

• LOFAR: The LOw-Frequency ARray
• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• A multi-scale multi-frequency deconvolution algorithm for synthesis imaging in radio interferometry
• Multi-Scale CLEAN deconvolution of radio synthesis images
• The Expanded Very Large Array -- a New Telescope for New Science
• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical 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
• PURIFY: a new approach to radio-interferometric imaging
• The application of compressive sampling to radio astronomy I: Deconvolution
• LOFAR Sparse Image Reconstruction
• Scale sensitive deconvolution of interferometric images I. Adaptive Scale Pixel (Asp) decomposition
• Scalable splitting algorithms for big-data interferometric imaging in the SKA era
• Dynamical Imaging with Interferometry
• MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators. A sparse deconvolution algorithm for radio interferometric imaging
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• Robust sparse image reconstruction of radio interferometric observations with purify
• An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet
• Cygnus A super-resolved via convex optimisation from VLA data
• Sparse representations and convex optimization as tools for LOFAR radio interferometric imaging
• A new approach to multi-frequency synthesis in radio interferometry
• Joint Multichannel Deconvolution and Blind Source Separation
• Using SKA Rotation Measures to Reveal the Mysteries of the Magnetised Universe
• Multi-frequency image reconstruction for radio interferometry. A regularized inverse problem approach

Cited by in corpus (11)

• Image reconstruction algorithms in radio interferometry: from handcrafted to learned regularization denoisers
• Comparison of classical and Bayesian imaging in radio interferometry
• Deep Radio Interferometric Imaging with POLISH: DSA-2000 and weak lensing
• First AI for deep super-resolution wide-field imaging in radio astronomy: unveiling structure in ESO 137--006
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations
• Scalable precision wide-field imaging in radio interferometry: I. uSARA validated on ASKAP data
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): II. Code and real data proof of concept
• fast-resolve: Fast Bayesian Radio Interferometric Imaging
• PolyCLEAN: Atomic Optimization for Super-Resolution Imaging and Uncertainty Estimation in Radio Interferometry
• A novel radio imaging method for physical spectral index modelling
• High-Performance Gridding For Radio Interferometric Image Synthesis