Accelerating the CLEAN algorithm of radio interferometry with convex optimization

TL;DR

This paper enhances the CLEAN algorithm in radio interferometry by integrating convex optimization acceleration techniques, significantly improving convergence speed and residual reduction, thus enabling faster and more robust imaging.

Contribution

It introduces a novel approach to accelerate the CLEAN algorithm using convex optimization methods like Nesterov acceleration and conjugate gradient, achieving faster convergence with lower computational cost.

Findings

Algorithms converge multiple times faster than traditional methods.

Residuals are significantly reduced with the accelerated algorithms.

Combining acceleration techniques yields the best performance.

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.