Repository logo

Infoscience

  • English
  • French
Log In
Logo EPFL, École polytechnique fédérale de Lausanne

Infoscience

  • English
  • French
Log In
  1. Home
  2. Academic and Research Output
  3. Journal articles
  4. PolyCLEAN: Atomic optimization for super-resolution imaging and uncertainty estimation in radio interferometry
 
research article

PolyCLEAN: Atomic optimization for super-resolution imaging and uncertainty estimation in radio interferometry

Jarret, Adrian  
•
Kashani, Sepand  
•
Rue-Queralt, Joan  
Show more
January 21, 2025
Astronomy & Astrophysics

Context. Imaging in radio interferometry requires 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 demonstrate excellent reconstruction quality thanks to their ability to enforce Bayesian priors. Nowadays, the main limitation to their adoption is 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, we propose a method for a fine resolution setup, which scales naturally in terms of memory usage and reconstruction speed. Second, we develop a new tool to localize a region of uncertainty, paving the way for quantitative imaging in radio interferometry. Methods. The classical & ell;(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 concept of 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.

  • Files
  • Details
  • Metrics
Loading...
Thumbnail Image
Name

10.1051_0004-6361_202450970.pdf

Type

Main Document

Version

Published version

Access type

openaccess

License Condition

CC BY

Size

22.86 MB

Format

Adobe PDF

Checksum (MD5)

7f5e19eb73ab7cbde299cfdf54cd96cc

Logo EPFL, École polytechnique fédérale de Lausanne
  • Contact
  • infoscience@epfl.ch

  • Follow us on Facebook
  • Follow us on Instagram
  • Follow us on LinkedIn
  • Follow us on X
  • Follow us on Youtube
AccessibilityLegal noticePrivacy policyCookie settingsEnd User AgreementGet helpFeedback

Infoscience is a service managed and provided by the Library and IT Services of EPFL. © EPFL, tous droits réservés