Non-convex optimization for self-calibration of direction-dependent effects in radio interferometric imaging

Abstract

Radio interferometric imaging aims to estimate an unknown sky intensity image from degraded observations, acquired through an antenna array. In the theoretical case of a perfectly calibrated array, it has been shown that solving the corresponding imaging problem by iterative algorithms based on convex optimization and compressive sensing theory can be competitive with classical algorithms such as CLEAN. However, in practice, antenna-based gains are unknown and have to be calibrated. Future radio telescopes, such as the SKA, aim at improving imaging resolution and sensitivity by orders of magnitude. At this precision level, the direction-dependency of the gains must be accounted for, and radio interferometric imaging can be understood as a blind deconvolution problem. In this context, the underlying minimization problem is non-convex, and adapted techniques have to be designed. In this work, leveraging recent developments in non-convex optimization, we propose the first joint calibration and imaging method in radio interferometry, with proven convergence guarantees. Our approach, based on a block-coordinate forward-backward algorithm, jointly accounts for visibilities and suitable priors on both the image and the direction-dependent effects (DDEs). As demonstrated in recent works, sparsity remains the prior of choice for the image, while DDEs are modelled as smooth functions of the sky, i.e. spatially band-limited. Finally, we show through simulations the efficiency of our method, for the reconstruction of both images of point sources and complex extended sources. MATLAB code is available on GitHub.

References in corpus (20)

• Multi-Scale CLEAN deconvolution of radio synthesis images
• Ionospheric Calibration of Low Frequency Radio Interferometric Observations using the Peeling Scheme: I. Method Description and First Results
• Compressed sensing imaging techniques for radio interferometry
• Compressed Sensing in Astronomy
• Radio interferometric gain calibration as a complex optimization problem
• Splitting methods with variable metric for KL functions
• Real-Time Calibration of the Murchison Widefield Array
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• PURIFY: a new approach to radio-interferometric imaging
• Radio Interferometric Calibration Using The SAGE Algorithm
• LOFAR Sparse Image Reconstruction
• Scale sensitive deconvolution of interferometric images I. Adaptive Scale Pixel (Asp) decomposition
• Advances in Calibration and Imaging Techniques in Radio Interferometry
• Scalable splitting algorithms for big-data interferometric imaging in the SKA era
• Fast gain calibration in radio astronomy using alternating direction implicit methods: Analysis and applications
• Compressed sensing for wide-field radio interferometric imaging
• Spread spectrum for imaging techniques in radio interferometry
• Calibration artefacts in radio interferometry. I. Ghost sources in WSRT data
• Compressed sensing reconstruction of a string signal from interferometric observations of the cosmic microwave background
• Tensor optimisation for optical-interferometric imaging

Cited by in corpus (13)

• The LOFAR Two Meter Sky Survey: Deep Fields, I -- Direction-dependent calibration and imaging
• Image reconstruction algorithms in radio interferometry: from handcrafted to learned regularization denoisers
• Cygnus A super-resolved via convex optimisation from VLA data
• First AI for deep super-resolution wide-field imaging in radio astronomy: unveiling structure in ESO 137--006
• Bayesian radio interferometric imaging with direction-dependent calibration
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations
• Sparse interferometric Stokes imaging under polarization constraint (Polarized SARA)
• Polca SARA - Full polarization, direction-dependent calibration and sparse imaging for radio interferometry
• Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data
• Scalable precision wide-field imaging in radio interferometry: I. uSARA validated on ASKAP data
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): II. Code and real data proof of concept
• Robust direction-dependent gain-calibration of beam-modelling errors far from the target field
• Primal-dual splitting scheme with backtracking for handling with epigraphic constraint and sparse analysis regularization