Multi-scale and Multi-directional VLBI Imaging with CLEAN

Abstract

Very long baseline interferometry (VLBI) is a radio-astronomical technique in which the correlated signal from various baselines is combined into an image of highest angular resolution. Due to sparsity of the measurements, this imaging procedure constitutes an ill-posed inverse problem. For decades the CLEAN algorithm was the standard choice in VLBI studies, although having some serious disadvantages and pathologies that are challenged by the requirements of modern frontline VLBI applications. We develop a novel multi-scale CLEAN deconvolution method (DoB-CLEAN) based on continuous wavelet transforms that address several pathologies in CLEAN imaging. We benchmark this novel algorithm against CLEAN reconstructions on synthetic data and reanalyze BL Lac observations of RadioAstron with DoB-CLEAN. DoB-CLEAN approaches the image by multi-scalar and multi-directional wavelet dictionaries. Two different dictionaries are used. Firstly, a difference of elliptical spherical Bessel functions dictionary fitted to the uv-coverage of the observation that is used to sparsely represent the features in the dirty image. Secondly, a difference of elliptical Gaussian wavelet dictionary that is well suited to represent relevant image features cleanly. The deconvolution is performed by switching between the dictionaries. DoB-CLEAN achieves super-resolution compared to CLEAN and remedies the spurious regularization properties of CLEAN. In contrast to CLEAN, the representation by basis functions has a physical meaning. Hence, the computed deconvolved image still fits the observed visibilities, opposed to CLEAN. State-of-the-art multi-scalar imaging approaches seem to outperform single-scalar standard approaches in VLBI and are well suited to maximize the extraction of information in ongoing frontline VLBI applications.

arXiv commentAccepted for publication in A&A

References in corpus (19)

• First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole
• A multi-scale multi-frequency deconvolution algorithm for synthesis imaging in radio interferometry
• Multi-Scale CLEAN deconvolution of radio synthesis images
• Interferometric Imaging Directly with Closure Phases and Closure Amplitudes
• Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope using Sparse Modeling
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• Probing the innermost regions of AGN jets and their magnetic fields with RadioAstron. I. Imaging BL Lacertae at 21 microarcsecond resolution
• LOFAR Sparse Image Reconstruction
• Super-resolution Full Polarimetric Imaging for Radio Interferometry with Sparse Modeling
• Scale sensitive deconvolution of interferometric images I. Adaptive Scale Pixel (Asp) decomposition
• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• Hybrid Very Long Baseline Interferometry Imaging and Modeling with Themis
• Comparison of classical and Bayesian imaging in radio interferometry
• Wavelet-based decomposition and analysis of structural patterns in astronomical images
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• Unified Radio Interferometric Calibration and Imaging with Joint Uncertainty Quantification
• VLBI imaging of a flare in the Crab Nebula: More than just a spot
• Modelling and peeling extended sources with shapelets: a Fornax A case study

Cited by in corpus (13)

• Using multiobjective optimization to reconstruct interferometric data (I)
• Dynamic and polarimetric VLBI imaging with a multiscalar approach
• CLEAN imaging systematics of M87 radio jet
• 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
• Event-horizon-scale Imaging of M87* under Different Assumptions via Deep Generative Image Priors
• PolyCLEAN: Atomic Optimization for Super-Resolution Imaging and Uncertainty Estimation in Radio Interferometry
• How to make CLEAN variants faster? Using clustered components informed by the autocorrelation function
• Prospects of using closure traces directly for imaging in Very Long Baseline Interferometry
• Probing circular polarization and magnetic field structure in AGN
• Accelerating the CLEAN algorithm of radio interferometry with convex optimization
• Radio properties of the quasi-periodic eruption source RXJ1301.9+2747 at parsec scales
• Multiobjective optimization for scattering mitigation and scattering screen reconstruction in VLBI observations of the Galactic Center