Accelerating the CLEAN algorithm of radio interferometry with convex optimization

Abstract

In radio-interferometry, we recover an image from an incompletely sampled Fourier data. The de-facto standard algorithm, the Cotton-Schwab CLEAN, is iteratively switching between computing a deconvolution (minor loop) and subtracting the model from the visibilities (major loop). The next generation of radio interferometers is expected to deal with much higher data rates, image sizes and sensitivity, making an acceleration of current data processing algorithms necessary. We aim to achieve this by evaluating the potential of various well-known acceleration techniques in convex optimization to the major loop. For the present manuscript, we limit the scope to study these techniques only in the CLEAN framework. To this end, we identify CLEAN with a Newton scheme, and use this chain of arguments backwards to express Nesterov acceleration and conjugate gradient orthogonalization in the major and minor loop framework. The resulting algorithms are simple extensions of the traditional framework, but converge multiple times faster than traditional techniques, and reduce the residual significantly deeper. These improvements achieved by accelerating the major loop are competitive to well-known improvements by replacing the minor loop with more advanced algorithms, but at lower numerical cost. The best performance is achieved by combining these two developments.CLEAN remains among the fastest and most robust algorithms for imaging in radio interferometry, and can be easily extended to an almost an order of magnitude faster convergence speed and dynamic range. The procedure outlined in this manuscript is relatively straightforward and could be easily extended.

accepted for publication A&A

References in corpus (33)

• First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole
• CASA, the Common Astronomy Software Applications for Radio Astronomy
• 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
• First Sagittarius A* Event Horizon Telescope Results. III: Imaging of the Galactic Center Supermassive Black Hole
• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• MOJAVE: XV. VLBA 15 GHz Total Intensity and Polarization Maps of 437 Parsec-Scale AGN Jets From 1996-2017
• Compressed sensing imaging techniques for radio interferometry
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• Interferometric Imaging Directly with Closure Phases and Closure Amplitudes
• Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope using Sparse Modeling
• A ring-like accretion structure in M87 connecting its black hole and jet
• LOFAR Sparse Image Reconstruction
• Scale sensitive deconvolution of interferometric images I. Adaptive Scale Pixel (Asp) decomposition
• The ALMA Interferometric Pipeline Heuristics
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• PHANGS-JWST First Results: Mapping the 3.3 micron Polycyclic Aromatic Hydrocarbon Vibrational Band in Nearby Galaxies with NIRCam Medium Bands
• PHANGS-JWST First Results: Tracing the Diffuse ISM with JWST Imaging of Polycyclic Aromatic Hydrocarbon Emission in Nearby Galaxies
• Comparison of classical and Bayesian imaging in radio interferometry
• An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet
• Bayesian radio interferometric imaging with direction-dependent calibration
• An Updated Formalism For Line-Driven Radiative Acceleration and Implications for Stellar Mass Loss
• Using multiobjective optimization to reconstruct interferometric data (I)
• Multi-scale and Multi-directional VLBI Imaging with CLEAN
• alpha-Deep Probabilistic Inference (alpha-DPI): efficient uncertainty quantification from exoplanet astrometry to black hole feature extraction
• Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data
• Dynamic and polarimetric VLBI imaging with a multiscalar approach
• Swarm intelligence for full Stokes dynamic imaging reconstruction of interferometric data
• fast-resolve: Fast Bayesian Radio Interferometric Imaging
• 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
• An Algorithm Architecture for Radio Interferometric Data Processing