PURIFY: a new approach to radio-interferometric imaging

Abstract

In a recent article series, the authors have promoted convex optimization algorithms for radio-interferometric imaging in the framework of compressed sensing, which leverages sparsity regularization priors for the associated inverse problem and defines a minimization problem for image reconstruction. This approach was shown, in theory and through simulations in a simple discrete visibility setting, to have the potential to outperform significantly CLEAN and its evolutions. In this work, we leverage the versatility of convex optimization in solving minimization problems to both handle realistic continuous visibilities and offer a highly parallelizable structure paving the way to significant acceleration of the reconstruction and high-dimensional data scalability. The new algorithmic structure promoted relies on the simultaneous-direction method of multipliers (SDMM), and contrasts with the current major-minor cycle structure of CLEAN and its evolutions, which in particular cannot handle the state-of-the-art minimization problems under consideration where neither the regularization term nor the data term are differentiable functions. We release a beta version of an SDMM-based imaging software written in C and dubbed PURIFY (http://basp-group.github.io/purify/) that handles various sparsity priors, including our recent average sparsity approach SARA. We evaluate the performance of different priors through simulations in the continuous visibility setting, confirming the superiority of SARA.

14 pages, 5 figures. Accepted in MNRAS

References in corpus (10)

• 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
• Correcting direction-dependent gains in the deconvolution of radio interferometric images
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• Wide-field wide-band interferometric imaging:The WB A-Projection and hybrid algorithms
• Fast Holographic Deconvolution: a new technique for precision radio interferometry
• Sparsity Averaging for Compressive Imaging
• Compressed sensing for wide-field radio interferometric imaging
• Revisiting the spread spectrum effect in radio interferometric imaging: a sparse variant of the w-projection algorithm

Cited by in corpus (56)

• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope using Sparse Modeling
• LOFAR Sparse Image Reconstruction
• Super-resolution Full Polarimetric Imaging for Radio Interferometry with Sparse Modeling
• Scalable splitting algorithms for big-data interferometric imaging in the SKA era
• Distributed Radio Interferometric Calibration
• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• Robust sparse image reconstruction of radio interferometric observations with purify
• Superresolution Interferometric Imaging with Sparse Modeling Using Total Squared Variation --- Application to Imaging the Black Hole Shadow
• Resolving the extended atmosphere and the inner wind of Mira ($o$ Ceti) with long ALMA baselines
• Image reconstruction algorithms in radio interferometry: from handcrafted to learned regularization denoisers
• Comparison of classical and Bayesian imaging in radio interferometry
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet
• Non-convex optimization for self-calibration of direction-dependent effects in radio interferometric imaging
• Cygnus A super-resolved via convex optimisation from VLA data
• Data Combination: Interferometry and Single-dish Imaging in Radio Astronomy
• Sparse representations and convex optimization as tools for LOFAR radio interferometric imaging
• An improved method for polarimetric image restoration in interferometry
• 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)
• Weak gravitational lensing with the Square Kilometre Array
• Sparse interferometric Stokes imaging under polarization constraint (Polarized SARA)
• Polca SARA - Full polarization, direction-dependent calibration and sparse imaging for radio interferometry
• A Fourier dimensionality reduction model for big data interferometric imaging
• 21cm Signal Recovery via the Robust Principle Component Analysis
• Efficient Computation of Slepian Functions for Arbitrary Regions on the Sphere
• Revisiting the spread spectrum effect in radio interferometric imaging: a sparse variant of the w-projection algorithm
• 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
• Slepian Spatial-Spectral Concentration Problem on the Sphere: Analytical Formulation for Limited Colatitude-Longitude Spatial Region
• Regularized Maximum Likelihood Image Synthesis and Validation for ALMA Continuum Observations of Protoplanetary Disks
• Radio Astronomical Image Formation using Constrained Least Squares and Krylov Subspaces
• High resolution VLBI polarisation imaging of AGN with the Maximum Entropy Method
• A parallel & automatically tuned algorithm for multispectral image deconvolution
• Posterior sampling for inverse imaging problems on the sphere in seismology and cosmology
• Multi-frequency image reconstruction for radio interferometry. A regularized inverse problem approach
• Learned Interferometric Imaging for the SPIDER Instrument
• $w$-stacking $w$-projection hybrid algorithm for wide-field interferometric imaging: implementation details and improvements
• Parallel Calibration for Sensor Array Radio Interferometers
• Multi-Step Reconstruction of Radio-Interferometric Images
• Distributed and parallel sparse convex optimization for radio interferometry with PURIFY
• Convolutional Deep Denoising Autoencoders for Radio Astronomical Images
• A randomised primal-dual algorithm for distributed radio-interferometric imaging
• A Decentralized Framework for Radio-interferometric Image Reconstruction
• PURIFY: a new algorithmic framework for next-generation radio-interferometric imaging
• Multi-frequency image reconstruction for radio-interferometry with self-tuned regularization parameters
• Fast algorithms to approximate the position-dependent point spread function responses in radio interferometric wide-field imaging
• Why CLEAN when you can PURIFY? A new approach for next-generation radio-interferometric imaging
• CS-ROMER: A novel compressed sensing framework for Faraday depth reconstruction
• Quantifying Uncertainty in High Dimensional Inverse Problems by Convex Optimisation
• DoG-HiT: A novel VLBI Multiscale Imaging Approach
• High-Performance Gridding For Radio Interferometric Image Synthesis
• Large Scale 3D Image Reconstruction in Optical Interferometry
• Prospects for radio weak lensing: studies using LOFAR observations in the ELAIS-N1 field