Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations

Abstract

Upcoming radio interferometers are aiming to image the sky at new levels of resolution and sensitivity, with wide-band image cubes reaching close to the Petabyte scale for SKA. Modern proximal optimization algorithms have shown a potential to significantly outperform CLEAN thanks to their ability to inject complex image models to regularize the inverse problem for image formation from visibility data. They were also shown to be parallelizable over large data volumes thanks to a splitting functionality enabling the decomposition of the data into blocks, for parallel processing of block-specific data-fidelity terms involved in the objective function. Focusing on intensity imaging, the splitting functionality is further exploited in this work to decompose the image cube into spatio-spectral facets, and enable parallel processing of facet-specific regularization terms in the objective function, leading to the "Faceted HyperSARA" algorithm. Reliable heuristics enabling an automatic setting of the regularization parameters involved in the objective are also introduced, based on estimates of the noise level, transferred from the visibility domain to the domains where the regularization is applied. Simulation results based on a MATLAB implementation and involving synthetic image cubes and data close to Gigabyte size confirm that faceting can provide a major increase in parallelization capability when compared to the non-faceted approach (HyperSARA).

arXiv commentPublished in MNRAS

References in corpus (35)

• LOFAR: The LOw-Frequency ARray
• 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
• Faceting for direction-dependent spectral deconvolution
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• PURIFY: a new approach to radio-interferometric imaging
• The Origin and Evolution of Cosmic Magnetism
• 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
• Sparsity Averaging for Compressive Imaging
• Probing the Dark Ages with the Square Kilometer Array
• 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
• Spread spectrum for imaging techniques in radio interferometry
• Robust sparse image reconstruction of radio interferometric observations with purify
• Unified Radio Interferometric Calibration and Imaging with Joint Uncertainty Quantification
• 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
• Probabilistic image reconstruction for radio interferometers
• Optimal Image Reconstruction in Radio Interferometry
• Sparse representations and convex optimization as tools for LOFAR radio interferometric imaging
• A new approach to multi-frequency synthesis in radio interferometry
• Galaxy Evolution and Cosmology with the Square Kilometre Array
• 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
• A Fourier dimensionality reduction model for big data interferometric imaging
• Joint Multichannel Deconvolution and Blind Source Separation
• The w-effect in interferometric imaging: from a fast sparse measurement operator to super-resolution

Cited by in corpus (8)

• Scalable precision wide-field imaging in radio interferometry: II. AIRI validated on ASKAP data
• 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
• Deep network series for large-scale high-dynamic range imaging
• A Distributed Block-Split Gibbs Sampler with Hypergraph Structure for High-Dimensional Inverse Problems
• A Decentralized Framework for Radio-interferometric Image Reconstruction