A novel radio imaging method for physical spectral index modelling
arXiv:2308.07306 · doi:10.1093/mnras/stad2465
Abstract
We present a new method, called "forced-spectrum fitting", for physically-based spectral modelling of radio sources during deconvolution. This improves upon current common deconvolution fitting methods, which often produce inaccurate spectra. Our method uses any pre-existing spectral index map to assign spectral indices to each model component cleaned during the multi-frequency deconvolution of WSClean, where the pre-determined spectrum is fitted. The component magnitude is evaluated by performing a modified weighted linear least-squares fit. We test this method on a simulated LOFAR-HBA observation of the 3C196 QSO and a real LOFAR-HBA observation of the 4C+55.16 FRI galaxy. We compare the results from the forced-spectrum fitting with traditional joined-channel deconvolution using polynomial fitting. Because no prior spectral information was available for 4C+55.16, we demonstrate a method for extracting spectral indices in the observed frequency band using "clustering". The models generated by the forced-spectrum fitting are used to improve the calibration of the datasets. The final residuals are comparable to existing multi-frequency deconvolution methods, but the output model agrees with the provided spectral index map, embedding correct spectral information. While forced-spectrum fitting does not solve the determination of the spectral information itself, it enables the construction of accurate multi-frequency models that can be used for wide-band calibration and subtraction.
arXiv comment17 pages, 9 figures, 5 tables. Accepted for publication in MNRAS
References in corpus (64)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Array Programming with NumPy
- The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
- LOFAR: The LOw-Frequency ARray
- WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
- The LOFAR Two-metre Sky Survey - I. Survey Description and Preliminary Data Release
- The LOFAR Two-metre Sky Survey - II. First data release
- A multi-scale multi-frequency deconvolution algorithm for synthesis imaging in radio interferometry
- GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey I: A low-frequency extragalactic catalogue
- Multi-Scale CLEAN deconvolution of radio synthesis images
- The Expanded Very Large Array -- a New Telescope for New Science
- The LOFAR Two-metre Sky Survey -- V. Second data release
- Shapelets: I. A Method for Image Analysis
- An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
- The third realization of the International Celestial Reference Frame by very long baseline interferometry
- Foreground simulations for the LOFAR - Epoch of Reionization Experiment
- Compressed sensing imaging techniques for radio interferometry
- Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
- The Sloan Digital Sky Survey quasar catalog: tenth data release
- Data Analysis for Precision 21 cm Cosmology
- Real-Time Calibration of the Murchison Widefield Array
- Foreground Removal using FastICA: A Showcase of LOFAR-EoR
- Extragalactic Peaked-Spectrum Radio Sources at Low Frequencies
- Correcting direction-dependent gains in the deconvolution of radio interferometric images
- Initial deep LOFAR observations of Epoch of Reionization windows: I. The North Celestial Pole
- Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
- Radio Interferometric Calibration Using The SAGE Algorithm
- A radio spectral index map and catalogue at 147-1400 MHz covering 80 per cent of the sky
- Statistical 21-cm Signal Separation via Gaussian Process Regression Analysis
- The application of compressive sampling to radio astronomy I: Deconvolution
- LOFAR Sparse Image Reconstruction
- Applying full polarization A-Projection to very wide field of view instruments: An imager for LOFAR
- Anticipated Performance of the Square Kilometre Array -- Phase 1 (SKA1)
- Calibration of Low-Frequency, Wide-Field Radio Interferometers Using Delay/Delay-Rate Filtering
- Distributed Radio Interferometric Calibration
- Sub-arcsecond imaging with the International LOFAR Telescope I. Foundational calibration strategy and pipeline
- Parametrising Epoch of Reionization foregrounds: A deep survey of low-frequency point-source spectra with the MWA
- Can the reionization epoch be detected as a global signature in the cosmic background?
- Multi-Scale CLEAN: A comparison of its performance against classical CLEAN in galaxies using THINGS
- LOFAR imaging of Cygnus A -- Direct detection of a turnover in the hotspot radio spectra
- MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators. A sparse deconvolution algorithm for radio interferometric imaging
- Efficient wide-field radio interferometry response
- DAYENU: A Simple Filter of Smooth Foregrounds for Intensity Mapping Power Spectra
- Low-Frequency Radio Absorption in Cassiopeia A
- Foreground removal for Square Kilometre Array observations of the Epoch of Reionization with the Correlated Component Analysis
- A deep / wide 1-2 GHz snapshot survey of SDSS Stripe 82 using the Karl G. Jansky Very Large Array in a compact hybrid configuration
- 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
- GASP XXXVIII: The LOFAR-MeerKAT-VLA view on the non-thermal side of a jellyfish galaxy
- Measuring cavity powers of active galactic nuclei in clusters using a hybrid X-ray-radio method -- A new window on feedback opened by subarcsecond LOFAR-VLBI observations
- Origin of the ring structures in Hercules A -- Sub-arcsecond 144 MHz to 7 GHz observations
- A new approach to multi-frequency synthesis in radio interferometry
- Wideband Super-resolution Imaging in Radio Interferometry via Low Rankness and Joint Average Sparsity Models (HyperSARA)
- Deep wideband single pointings and mosaics in radio interferometry - How accurately do we reconstruct intensities and spectral indices of faint sources?
- A method for determining radio continuum spectra, and its application to large surveys
- Modelling and peeling extended sources with shapelets: a Fornax A case study
- Polca SARA - Full polarization, direction-dependent calibration and sparse imaging for radio interferometry
- Observing the Reionization : Effect of Calibration and Position Errors on Realistic Observation Conditions
- A new look at old friends. I. Imaging classical radio galaxies with uGMRT and MeerKAT
- A Polarimetric Study of 9 PG Quasars with the VLA
- Degree-Scale Galactic Radio Emission at 122 MHz around the North Celestial Pole with LOFAR-AARTFAAC
- Multi-frequency image reconstruction for radio interferometry. A regularized inverse problem approach
- Multi-frequency synthesis algorithm based on Generalized Maximum Entropy Method. Application to 0954+658
- Correction of Radio Interferometric Imaging for Antenna Patterns