A Novel CMB Component Separation Method: Hierarchical Generalized Morphological Component Analysis
arXiv:1910.08077 · doi:10.1093/mnras/staa744
Abstract
We present a novel technique for Cosmic Microwave Background (CMB) foreground subtraction based on the framework of blind source separation. Inspired by previous work incorporating local variation to Generalized Morphological Component Analysis (GMCA), we introduce Hierarchical GMCA (HGMCA), a Bayesian hierarchical graphical model for source separation. We test our method on simulated sky maps that include dust, synchrotron, free-free and anomalous microwave emission, and show that HGMCA reduces foreground contamination by over GMCA in both the regions included and excluded by the Planck UT78 mask, decreases the error in the measurement of the CMB temperature power spectrum to the level at (and for all ), and reduces correlation to all the foregrounds. We find equivalent or improved performance when compared to state-of-the-art Internal Linear Combination (ILC)-type algorithms on these simulations, suggesting that HGMCA may be a competitive alternative to foreground separation techniques previously applied to observed CMB data. Additionally, we show that our performance does not suffer when we perturb model parameters or alter the CMB realization, which suggests that our algorithm generalizes well beyond our simplified simulations. Our results open a new avenue for constructing CMB maps through Bayesian hierarchical analysis.
arXiv commentUpdated to reflect accepted MNRAS version
References in corpus (48)
- Planck 2015 results. XIII. Cosmological parameters
- Planck 2013 results. XVI. Cosmological parameters
- HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
- Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results
- Planck 2013 results. I. Overview of products and scientific results
- Planck 2018 results. I. Overview and the cosmological legacy of Planck
- Planck 2015 results. XI. CMB power spectra, likelihoods, and robustness of parameters
- A high resolution foreground cleaned CMB map from WMAP
- Planck 2013 results. XXIII. Isotropy and statistics of the CMB
- Planck 2018 results. IV. Diffuse component separation
- A full sky, low foreground, high resolution CMB map from WMAP
- Planck 2013 results. XII. Component separation
- Planck 2015 results. X. Diffuse component separation: Foreground maps
- Component separation methods for the Planck mission
- The Python Sky Model: software for simulating the Galactic microwave sky
- Power Spectrum of Primordial Inhomogeneity Determined from the 4-Year COBE DMR Sky Maps
- Power spectrum estimation from high-resolution maps by Gibbs sampling
- Multi-Detector Multi-Component spectral matching and applications for CMB data analysis
- How isotropic is the Universe?
- Foreground separation methods for satellite observations of the cosmic microwave background
- A novel sampling theorem on the sphere
- Planck 2013 results. XIX. The integrated Sachs-Wolfe effect
- Planck 2018 results. IX. Constraints on primordial non-Gaussianity
- Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Bayesian Estimation of CMB Polarization Maps
- Exact reconstruction with directional wavelets on the sphere
- Asymptotics for spherical needlets
- A needlet ILC analysis of WMAP 7-year data: estimation of CMB temperature map and power spectrum
- Fast and Efficient Template Fitting of Deterministic Anisotropic Cosmological Models Applied to WMAP Data
- Planck 2015 results. XII. Full Focal Plane simulations
- Spatial variations in the spectral index of polarized synchrotron emission in the 9 yr WMAP sky maps
- Joint Planck and WMAP CMB Map Reconstruction
- A needlet ILC analysis of WMAP 9-year polarisation data: CMB polarisation power spectra
- S2LET: A code to perform fast wavelet analysis on the sphere
- Template fitting and the large-angle CMB anomalies
- Sparse component separation for accurate CMB map estimation
- CMB power spectrum estimation and map reconstruction with the Expectation - Maximization algorithm
- Separation of correlated astrophysical sources using multiple-lag data covariance matrices
- Localisation of directional scale-discretised wavelets on the sphere
- Blind component separation in wavelet space. Application to CMB analysis
- Sparse and Non-Negative BSS for Noisy Data
- Second-Generation Curvelets on the Sphere
- Directional spin wavelets on the sphere
- Wavelet-Bayesian inference of cosmic strings embedded in the cosmic microwave background
- A novel sampling theorem on the rotation group
- SILC: a new Planck Internal Linear Combination CMB temperature map using directional wavelets
- Spin-SILC: CMB polarisation component separation with spin wavelets
- A multi-level solver for Gaussian constrained CMB realizations
- "Internal Linear Combination" method for the separation of CMB from Galactic foregrounds in the harmonic domain