MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators. A sparse deconvolution algorithm for radio interferometric imaging

Abstract

(arXiv abridged abstract) The current years are seeing huge developments of radio telescopes and a tremendous increase of their capabilities. Such systems make mandatory the design of more sophisticated techniques not only for transporting, storing and processing this new generation of radio interferometric data, but also for restoring the astrophysical information contained in such data. In this paper we present a new radio deconvolution algorithm named MORESANE and its application to fully realistic simulated data of MeerKAT, one of the SKA precursors. This method has been designed for the difficult case of restoring diffuse astronomical sources which are faint in brightness, complex in morphology and possibly buried in the dirty beam's side lobes of bright radio sources in the field. MORESANE is a greedy algorithm which combines complementary types of sparse recovery methods in order to reconstruct the most appropriate sky model from observed radio visibilities. A synthesis approach is used for the reconstruction of images, in which the synthesis atoms representing the unknown sources are learned using analysis priors. We apply this new deconvolution method to fully realistic simulations of radio observations of a galaxy cluster and of an HII region in M31. We show that MORESANE is able to efficiently reconstruct images composed from a wide variety of sources from radio interferometric data. Comparisons with other available algorithms, which include multi-scale CLEAN and the recently proposed methods by Li et al. (2011) and Carrillo et al. (2012), indicate that MORESANE provides competitive results in terms of both total flux/surface brightness conservation and fidelity of the reconstructed model. MORESANE seems particularly well suited for the recovery of diffuse and extended sources, as well as bright and compact radio sources known to be hosted in galaxy clusters.

17 pages, 11 figures, accepted for publication on A&A

References in corpus (23)

• LOFAR: The LOw-Frequency ARray
• Clusters of galaxies : observational properties of the diffuse radio emission
• Cosmic rays in galaxy clusters and their non-thermal emission
• Multi-Scale CLEAN deconvolution of radio synthesis images
• Shapelets: I. A Method for Image Analysis
• Observations of extended radio emission in clusters
• Magnetic Fields and Faraday Rotation in Clusters of Galaxies
• Compressed sensing imaging techniques for radio interferometry
• Global testing under sparse alternatives: ANOVA, multiple comparisons and the higher criticism
• Deconvolution with correct sampling
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• The intracluster magnetic field power spectrum in Abell 2255
• M87 at metre wavelengths: the LOFAR picture
• The application of compressive sampling to radio astronomy I: Deconvolution
• The MeqTrees software system and its use for third-generation calibration of radio interferometers
• LOFAR Sparse Image Reconstruction
• Applying full polarization A-Projection to very wide field of view instruments: An imager for LOFAR
• Radio Continuum Surveys with Square Kilometre Array Pathfinders
• Compressed sensing for wide-field radio interferometric imaging
• Positive deconvolution for superimposed extended source and point sources
• Spread spectrum for imaging techniques in radio interferometry
• Radio-continuum study of Large Magellanic Cloud Supernova Remnant J0509-6731
• Thermal and non-thermal components of the interstellar medium at sub-kiloparsec scales in galaxies

Cited by in corpus (39)

• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• Faceting for direction-dependent spectral deconvolution
• Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope using Sparse Modeling
• Super-resolution Full Polarimetric Imaging for Radio Interferometry with Sparse Modeling
• Scalable splitting algorithms for big-data interferometric imaging in the SKA era
• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• The MUSE Extremely Deep Field: the Cosmic Web in Emission at High Redshift
• Extragalactic Radio Continuum Surveys and the Transformation of Radio Astronomy
• Forming intracluster gas in a galaxy protocluster at a redshift of 2.16
• Superresolution Interferometric Imaging with Sparse Modeling Using Total Squared Variation --- Application to Imaging the Black Hole Shadow
• Robust sparse image reconstruction of radio interferometric observations with purify
• Image reconstruction algorithms in radio interferometry: from handcrafted to learned regularization denoisers
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet
• First AI for deep super-resolution wide-field imaging in radio astronomy: unveiling structure in ESO 137--006
• Cygnus A super-resolved via convex optimisation from VLA data
• Sparse representations and convex optimization as tools for LOFAR radio interferometric imaging
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations
• Automated detection of extended sources in radio maps: progress from the SCORPIO survey
• LOFAR 150-MHz observations of SS 433 and W 50
• DAWIS, a Detection Algorithm with Wavelets for Intracluster light Studies
• Wideband Super-resolution Imaging in Radio Interferometry via Low Rankness and Joint Average Sparsity Models (HyperSARA)
• A Fourier dimensionality reduction model for big data interferometric imaging
• 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
• The w-effect in interferometric imaging: from a fast sparse measurement operator to super-resolution
• Regularized Maximum Likelihood Image Synthesis and Validation for ALMA Continuum Observations of Protoplanetary Disks
• High resolution VLBI polarisation imaging of AGN with the Maximum Entropy Method
• A parallel & automatically tuned algorithm for multispectral image deconvolution
• Wide-band Rotation Measure Synthesis
• Deep learning-based deconvolution for interferometric radio transient reconstruction
• A novel radio imaging method for physical spectral index modelling
• Multi-Step Reconstruction of Radio-Interferometric Images
• PCAT-DE: Reconstructing point-like and diffuse signals in astronomical images using spatial and spectral information
• Convolutional Deep Denoising Autoencoders for Radio Astronomical Images
• A Decentralized Framework for Radio-interferometric Image Reconstruction
• Fast algorithms to approximate the position-dependent point spread function responses in radio interferometric wide-field imaging
• Prospects for radio weak lensing: studies using LOFAR observations in the ELAIS-N1 field
• An Unsupervised Learning Method for Radio Interferometry Deconvolution