DoG-HiT: A novel VLBI Multiscale Imaging Approach

Abstract

Reconstructing images from very long baseline interferometry (VLBI) data with sparse sampling of the Fourier domain (uv-coverage) constitutes an ill-posed deconvolution problem. It requires application of robust algorithms maximizing the information extraction from all of the sampled spatial scales and minimizing the influence of the unsampled scales on image quality. We develop a new multiscale wavelet deconvolution algorithm DoG-HiT for imaging sparsely sampled interferometric data which combines the difference of Gaussian (DoG) wavelets and hard image thresholding (HiT). Based on DoG-HiT, we propose a multi-step imaging pipeline for analysis of interferometric data. DoG-HiT applies the compressed sensing approach to imaging by employing a flexible DoG wavelet dictionary which is designed to adapt smoothly to the uv-coverage. It uses closure properties as data fidelity terms only initially and perform non-convex, non-smooth optimization by an amplitude conserving and total flux conserving hard thresholding splitting. DoG-HiT calculates a multiresolution support as a side product. The final reconstruction is refined through self-calibration loops and imaging with amplitude and phase information applied for the multiresolution support only. We demonstrate the stability of DoG-HiT and benchmark its performance against image reconstructions made with CLEAN and Regularized Maximum-Likelihood (RML) methods using synthetic data. The comparison shows that DoG-HiT matches the superresolution achieved by the RML reconstructions and surpasses the sensitivity to extended emission reached by CLEAN. Application of regularized maximum likelihood methods outfitted with flexible multiscale wavelet dictionaries to imaging of interferometric data matches the performance of state-of-the art convex optimization imaging algorithms and requires fewer prior and user defined constraints.

arXiv comment20 pages, 11 Figures, accepted for publication in A&A, the imaging software underlying this work will be made publicly available soon under the link https://github.com/hmuellergoe/mrbeam

References in corpus (27)

• 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
• Compressed sensing imaging techniques for radio interferometry
• 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
• PURIFY: a new approach to radio-interferometric imaging
• The application of compressive sampling to radio astronomy I: Deconvolution
• 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
• Evaluation of New Submillimeter VLBI Sites for the Event Horizon Telescope
• Scalable splitting algorithms for big-data interferometric imaging in the SKA era
• Variable structures in M87* from space, time and frequency resolved interferometry
• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• Closure statistics in interferometric data
• Comparison of the ion-to-electron temperature ratio prescription: GRMHD simulations with electron thermodynamics
• Robust sparse image reconstruction of radio interferometric observations with purify
• Principles of image reconstruction in optical interferometry: tutorial
• Wavelet-based decomposition and analysis of structural patterns in astronomical images
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet
• Sparse representations and convex optimization as tools for LOFAR radio interferometric imaging
• Modelling and peeling extended sources with shapelets: a Fornax A case study
• PRECL: A new method for interferometry imaging from closure phase
• Radio Astronomical Image Formation using Constrained Least Squares and Krylov Subspaces