PolyCLEAN: When Högbom meets Bayes ⋆ Fast Super-Resolution Imaging with Bayesian MAP Estimation
Context. Imaging in radio interferometry amounts to solving an ill-posed noisy inverse problem, for which the most adopted algorithm is the original CLEAN algorithm and its variants. Alternative explicit optimization methods have gained increasing attention, as they have demonstrated excellent reconstruction quality due to their ability to enforce Bayesian priors. Nowadays, the main limitation to their adoption lies in run-time speed. Additionally, uncertainty quantification is difficult for both CLEAN and convex optimization techniques. Aims. We address two issues for the adoption of convex optimization in radio interferometric imaging. First, a method for a fine resolution setup is proposed which scales naturally in terms of memory usage and reconstruction speed. Second, a new tool to localize a region of uncertainty is developed, paving the way for quantitative imaging in radio interferometry. Methods. The classical ℓ 1 penalty is used to turn the inverse problem into a sparsity-promoting optimization. For efficient implementation, the so-called Frank-Wolfe algorithm is used together with a polyatomic refinement. The algorithm naturally produces sparse images at each iteration, leveraged to reduce memory and computational requirements. In that regard, PolyCLEAN reproduces the numerical behavior of CLEAN while guaranteeing that it solves the minimization problem of interest. Additionally, we introduce the dual certificate image, which appears as a numerical byproduct of the Frank-Wolfe algorithm. This image is proposed as a tool for uncertainty quantification on the location of the recovered sources. Results. PolyCLEAN demonstrates good scalability performance, in particular for fine-resolution grids. On simulations, the Pythonbased implementation is competitive with the fast numerically-optimized CLEAN solver. This acceleration does not affect image reconstruction quality: PolyCLEAN images are consistent with CLEAN-obtained ones for both point sources and diffuse emission recovery. We also highlight PolyCLEAN reconstruction capabilities on observed radio measurements. Conclusions. PolyCLEAN can be considered as an alternative to CLEAN in the radio interferometric imaging pipeline, as it enables the use of Bayesian priors without impacting the scalability and numerical performance of the imaging method.
PolyCLEAN___A_A_Format.pdf
main document
openaccess
CC BY-NC-SA
5 MB
Adobe PDF
36e9754d3934e625284e949b1562e60a