CS-ROMER: A novel compressed sensing framework for Faraday depth reconstruction
arXiv:2205.01413 · doi:10.1093/mnras/stac3031
Abstract
The reconstruction of Faraday depth structure from incomplete spectral polarization radio measurements using the RM Synthesis technique is an under-constrained problem requiring additional regularisation. In this paper we present cs-romer: a novel object-oriented compressed sensing framework to reconstruct Faraday depth signals from spectro-polarization radio data. Unlike previous compressed sensing applications, this framework is designed to work directly with data that are irregularly sampled in wavelength-squared space and to incorporate multiple forms of compressed sensing regularisation. We demonstrate the framework using simulated data for the VLA telescope under a variety of observing conditions, and we introduce a methodology for identifying the optimal basis function for reconstruction of these data, using an approach that can also be applied to datasets from other telescopes and over different frequency ranges. In this work we show that the delta basis function provides optimal reconstruction for VLA L-band data and we use this basis with observations of the low-mass galaxy cluster Abell 1314 in order to reconstruct the Faraday depth of its constituent cluster galaxies. We use the cs-romer framework to de-rotate the Galactic Faraday depth contribution directly from the wavelength-squared data and to handle the spectral behaviour of different radio sources in a direction-dependent manner. The results of this analysis show that individual galaxies within Abell 1314 deviate from the behaviour expected for a Faraday-thin screen such as the intra-cluster medium and instead suggest that the Faraday rotation exhibited by these galaxies is dominated by their local environments.
Final version after peer reviewing
References in corpus (14)
- Pathwise coordinate optimization
- First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole
- An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
- Compressed sensing imaging techniques for radio interferometry
- Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
- The 1.28 GHz MeerKAT DEEP2 Image
- Super-resolution Full Polarimetric Imaging for Radio Interferometry with Sparse Modeling
- The intra-cluster magnetic field in the double relic galaxy cluster Abell 2345
- Direction Dependent Corrections in Polarimetric Radio Imaging I : Characterizing the effects of the primary beam on full Stokes imaging
- Giant Metrewave Radio Telescope Observations of Head-Tail Radio Galaxies
- Low-frequency observations of the Giant Radio Galaxy NGC 6251
- An Iterative Reconstruction Algorithm for Faraday Tomography
- Gaussian Process Modelling for Improved Resolution in Faraday Depth Reconstruction
- Wavelets and sparsity for Faraday tomography