Cleaning foregrounds from single-dish 21cm intensity maps with Kernel Principal Component Analysis
arXiv:2107.02267 · doi:10.1093/mnras/stab2855
Abstract
The high dynamic range between contaminating foreground emission and the fluctuating 21cm brightness temperature field is one of the most problematic characteristics of 21cm intensity mapping data. While these components would ordinarily have distinctive frequency spectra, making it relatively easy to separate them, instrumental effects and calibration errors further complicate matters by modulating and mixing them together. A popular class of foreground cleaning method are unsupervised techniques related to Principal Component Analysis (PCA), which exploit the different shapes and amplitudes of each component's contribution to the covariance of the data in order to segregate the signals. These methods have been shown to be effective at removing foregrounds, while also unavoidably filtering out some of the 21cm signal too. In this paper we examine, for the first time in the context of 21cm intensity mapping, a generalised method called Kernel PCA, which instead operates on the covariance of non-linear transformations of the data. This allows more flexible functional bases to be constructed, in principle allowing a cleaner separation between foregrounds and the 21cm signal to be found. We show that Kernel PCA is effective when applied to simulated single-dish (autocorrelation) 21cm data under a variety of assumptions about foregrounds models, instrumental effects etc. It presents a different set of behaviours to PCA, e.g. in terms of sensitivity to the data resolution and smoothing scale, outperforming it on intermediate to large scales in most scenarios.
arXiv comment14 pages + 3 new pages, 11 + 5 new figures. Updated to match version accepted by MNRAS
References in corpus (43)
- Scikit-learn: Machine Learning in Python
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- The NumPy array: a structure for efficient numerical computation
- HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
- Baryon Acoustic Oscillation Intensity Mapping as a Test of Dark Energy
- Measurement of 21 cm brightness fluctuations at z ~ 0.8 in cross-correlation
- Late-time cosmology with 21cm intensity mapping experiments
- Core Cosmology Library: Precision Cosmological Predictions for LSST
- nbodykit: an open-source, massively parallel toolkit for large-scale structure
- Determination of z~0.8 neutral hydrogen fluctuations using the 21 cm intensity mapping auto-correlation
- Separation of anomalous and synchrotron emissions using WMAP polarization data
- An Improved Model of Diffuse Galactic Radio Emission from 10 MHz to 5 THz
- Lack of clustering in low-redshift 21-cm intensity maps cross-correlated with 2dF galaxy densities
- Coaxing Cosmic 21cm Fluctuations from the Polarized Sky using m-mode Analysis
- Foreground component separation with generalised ILC
- Blind foreground subtraction for intensity mapping experiments
- The Scale of the Problem : Recovering Images of Reionization with GMCA
- The 1.28 GHz MeerKAT DEEP2 Image
- Statistical 21-cm Signal Separation via Gaussian Process Regression Analysis
- The effect of foreground subtraction on cosmological measurements from Intensity Mapping
- nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics
- Interpreting the unresolved intensity of cosmologically redshifted line radiation
- HI intensity mapping with MeerKAT: Calibration pipeline for multi-dish autocorrelation observations
- Fast simulations for intensity mapping experiments
- SZ and CMB reconstruction using Generalized Morphological Component Analysis
- Improved foreground removal in GMRT 610 MHz observations towards redshifted 21-cm tomography
- 21cm foregrounds and polarization leakage: a user's guide on cleaning and mitigation strategies
- Erasing the Milky Way: new cleaning technique applied to GBT intensity mapping data
- Extracting HI cosmological signal with Generalized Needlet Internal Linear Combination
- MeerKLASS: MeerKAT Large Area Synoptic Survey
- Recovery of 21 cm intensity maps with sparse component separation
- Generating Log-normal Mock Catalog of Galaxies in Redshift Space
- Primary beam effects of radio astronomy antennas -- II. Modelling the MeerKAT L-band beam
- Determining the HI content of galaxies via intensity mapping cross-correlations
- Impact of Foregrounds on HI Intensity Mapping Cross-Correlations with Optical Surveys
- Gaussian Process Foreground Subtraction and Power Spectrum Estimation for 21 cm Cosmology
- Hi intensity mapping with MeerKAT: Primary beam effects on foreground cleaning
- Prospects for clustering and lensing measurements with forthcoming intensity mapping and optical surveys
- Gaussian Process Regression for foreground removal in HI intensity mapping experiments
- Foreground modelling via Gaussian process regression: an application to HERA data
- Power spectrum multipole expansion for HI intensity mapping experiments: unbiased parameter estimation
- Cross-correlating 21 cm and galaxy surveys: implications for cosmology and astrophysics
- Statistical recovery of the BAO scale from multipoles of the beam-convolved 21cm correlation function
Cited by in corpus (11)
- SKAO HI Intensity Mapping: Blind Foreground Subtraction Challenge
- A review of unsupervised learning in astronomy
- The foreground transfer function for HI intensity mapping signal reconstruction: MeerKLASS and precision cosmology applications
- Detecting the power spectrum turnover with HI intensity mapping
- Cosmological Probes of Structure Growth and Tests of Gravity
- Multipole expansion for 21cm Intensity Mapping power spectrum: forecasted cosmological parameters estimation for the SKA Observatory
- Square Kilometre Array Science Data Challenge 3a: foreground removal for an EoR experiment
- Unveiling the Universe with Emerging Cosmological Probes
- Baryon acoustic oscillations from HI intensity mapping: the importance of cross-correlations in the monopole and quadrupole
- Spin-based removal of instrumental systematics in 21cm intensity mapping surveys
- Simulating a full-sky high resolution Galactic synchrotron spectral index map using neural networks