How to make CLEAN variants faster? Using clustered components informed by the autocorrelation function
arXiv:2504.16058 · doi:10.1051/0004-6361/202553990
Abstract
Deconvolution, imaging and calibration of data from radio interferometers is a challenging computational (inverse) problem. The upcoming generation of radio telescopes poses significant challenges to existing, and well proven data reduction pipelines due to the large data sizes expected from these experiments, and the high resolution and dynamic range. In this manuscript, we deal with the deconvolution problem. A variety of multiscalar variants to the classical CLEAN algorithm (the de-facto standard) have been proposed in the past, often outperforming CLEAN at the cost of significantly increasing numerical resources. In this work, we aim to combine some of these ideas for a new algorithm, Autocorr-CLEAN, to accelerate the deconvolution and prepare the data reduction pipelines for the data sizes expected by the upcoming generation of instruments. To this end, we propose to use a cluster of CLEAN components fitted to the autocorrelation function of the residual in a subminor loop, to derive continuously changing, and potentially non-radially symmetric, basis functions for CLEANing the residual. Autocorr-CLEAN allows for the superior reconstruction fidelity achieved by modern multiscalar approaches, and their superior convergence speed. It achieves this without utilizing any substep of super-linear complexity in the minor loops, keeping the single minor loop and subminor loop iterations at an execution time comparable to CLEAN. Combining these advantages, Autocorr-CLEAN is found to be up to a magnitude faster than the classical CLEAN procedure. Autocorr-CLEAN fits well in the algorithmic framework common for radio interferometry, making it relatively straightforward to include in future data reduction pipelines. With its accelerated convergence speed, and smaller residual, Autocorr-CLEAN may be an important asset for the data analysis in the future.
accepted for A&A
References in corpus (30)
- 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
- An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
- 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
- Scale sensitive deconvolution of interferometric images I. Adaptive Scale Pixel (Asp) decomposition
- Applying full polarization A-Projection to very wide field of view instruments: An imager for LOFAR
- AUTO-MULTITHRESH: A General Purpose Automasking Algorithm
- Scalable splitting algorithms for big-data interferometric imaging in the SKA era
- RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
- Comparison of classical and Bayesian imaging in radio interferometry
- An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet
- Deep Learning-based Imaging in Radio Interferometry
- Efficient implementation of the adaptive scale pixel decomposition algorithm
- Bayesian radio interferometric imaging with direction-dependent calibration
- Using multiobjective optimization to reconstruct interferometric data (I)
- Multi-scale and Multi-directional VLBI Imaging with CLEAN
- Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data
- Evidence for a toroidal magnetic field in the core of 3C 84
- Dynamic and polarimetric VLBI imaging with a multiscalar approach
- Imaging a ring-like structure and the extended jet of M87 at 86 GHz
- Swarm intelligence for full Stokes dynamic imaging reconstruction of interferometric data
- The Pointing Self Calibration algorithm for aperture synthesis radio telescopes
- PolyCLEAN: Atomic Optimization for Super-Resolution Imaging and Uncertainty Estimation in Radio Interferometry
- Prospects of using closure traces directly for imaging in Very Long Baseline Interferometry
- Probing circular polarization and magnetic field structure in AGN