CMB power spectrum estimation and map reconstruction with the Expectation - Maximization algorithm
arXiv:astro-ph/0302094 · doi:10.1046/j.1365-2966.2003.06885.x
Abstract
We apply the iterative Expectation-Maximization algorithm (EM) to estimate the power spectrum of the CMB from multifrequency microwave maps. In addition, we are also able to provide a reconstruction of the CMB map. By assuming that the combined emission of the foregrounds plus the instrumental noise is Gaussian distributed in Fourier space, we have simplified the EM procedure finding an analytical expression for the maximization step. By using the simplified expression the CPU time can be greatly reduced. We test the stability of our power spectrum estimator with realistic simulations of Planck data, including point sources and allowing for spatial variation of the frequency dependence of the Galactic emissions. Without prior information about any of the components, our new estimator can recover the CMB power spectrum up to scales l~1500 with less than 10% error. This result is significantly improved if the brightest point sources are removed before applying our estimator. In this way, the CMB power spectrum can be recovered up to l~1700 with 10% error and up to l~2100 with 50% error. This result is very close to the one that would be obtained in the ideal case of only CMB plus white noise, for which all our assumptions are satisfied. Moreover, the EM algorithm also provides an straightforward mechanism to reconstruct the CMB map. The recovered cosmological signal shows a high degree of correlation (r = 0.98) with the input map and low residuals.
Final version (minor changes). 10 pages, 6 figures, MNRAS accepted
References in corpus (30)
- A Line of Sight Approach to Cosmic Microwave Background Anisotropies
- MAXIMA-1: A Measurement of the Cosmic Microwave Background Anisotropy on angular scales of 10 arcminutes to 5 degrees
- A measurement by BOOMERANG of multiple peaks in the angular power spectrum of the cosmic microwave background
- DASI First Results: A Measurement of the Cosmic Microwave Background Angular Power Spectrum
- Detection of Polarization in the Cosmic Microwave Background using DASI
- Extrapolation of Galactic Dust Emission at 100 Microns to CMBR Frequencies Using FIRAS
- A Robotic Wide-Angle H-Alpha Survey of the Southern Sky
- Diffuse Galactic Emission from Spinning Dust Grains
- The Anisotropy of the Microwave Background to l = 3500: Deep Field Observations with the Cosmic Background Imager
- A method for subtracting foregrounds from multi-frequency CMB sky maps
- Extragalactic Source Counts and Contributions to the Anisotropies of the Cosmic Microwave Background. Predictions for the Planck Surveyor mission
- Foreground separation methods for satellite observations of the cosmic microwave background
- Foregrounds and CMB Experiments: I. Semi-analytical estimates of contamination
- Galactic microwave emission at degree angular scales
- Towards a model of full-sky Galactic synchrotron intensity and linear polarisation: a re-analysis of the Parkes data
- Filtering techniques for the detection of Sunyaev-Zel'dovich clusters in multifrequency CMB maps
- Removing point sources from CMB maps
- Measuring Polarization with DASI
- The CMB power spectrum out to l=1400 measured by the VSA
- Neural networks and separation of Cosmic Microwave Background and astrophysical signals in sky maps
- Isotropic Wavelets: a Powerful Tool to Extract Point Sources from CMB Maps
- Tentative Detection of Electric Dipole Emission from Rapidly Rotating Dust Grains
- Point Source Detection using the Spherical Mexican Hat Wavelet on simulated all-sky Planck maps
- All-sky component separation for the Planck mission
- Optimal detection of sources on a homogeneous and isotropic background
- Combining maximum-entropy and the mexican hat wavelet to reconstruct the microwave sky
- Scale-Adaptive filters for the detection/separation of compact sources
- Predicted Planck Extragalactic Point Source Catalogue
- Constraining our Universe with X-ray & Optical Cluster Data
- An adaptive filter for the PLANCK Compact Source Catalogue construction
Cited by in corpus (36)
- Planck 2013 results. XII. Component separation
- Predictions for high-frequency radio surveys of extragalactic sources
- Component separation methods for the Planck mission
- The pre-launch Planck Sky Model: a model of sky emission at submillimetre to centimetre wavelengths
- Multi-resolution internal template cleaning: An application to the Wilkinson Microwave Anisotropy Probe 7-yr polarization data
- Full-sky Cosmic Microwave Background Foreground Cleaning Using Machine Learning
- Mining the UKIDSS GPS: star formation and embedded clusters
- Direct detection of the kinetic Sunyaev-Zel'dovich effect in galaxy clusters
- Component separation with flexible models. Application to the separation of astrophysical emissions
- Foreground separation using a flexible maximum-entropy algorithm: an application to COBE data
- The hybrid SZ power spectrum: Combining cluster counts and SZ fluctuations to probe gas physics
- Multi-Clustering Needlet-ILC for CMB B-modes component separation
- SILC: a new Planck Internal Linear Combination CMB temperature map using directional wavelets
- Minkowski Functionals of CMB polarisation intensity with Pynkowski: theory and application to Planck and future data
- Foreground analysis from the 1-year Wilkinson Microwave Anisotropy Probe (WMAP) data
- CENN: A fully convolutional neural network for CMB recovery in realistic microwave sky simulations
- Spin-SILC: CMB polarisation component separation with spin wavelets
- Correlated Component Analysis for diffuse component separation with error estimation on simulated Planck polarization data
- Hierarchical Bayesian CMB Component Separation with the No-U-Turn Sampler
- Analysis of NILC performance on B-modes data of sub-orbital experiments
- Effect of Foregrounds on the CMBR Multipole Alignment
- Efficient ILC analysis on polarization maps after EB leakage correction
- Linearized iterative least-squares (LIL): A parameter fitting algorithm for component separation in multifrequency CMB experiments such as Planck
- A Novel CMB Component Separation Method: Hierarchical Generalized Morphological Component Analysis
- Effect of component separation on the temperature distribution of the CMB
- Testing the ABS method with the simulated Planck temperature maps
- Testing for foreground residuals in the Planck foreground cleaned maps: A new method for designing confidence masks
- Sparse point-source removal for full-sky CMB experiments: application to WMAP 9-year data
- Diffuse source separation in CMB observations
- A Constrained NILC method for CMB B mode observations
- The Cosmic Microwave Background anisotropies: open problems
- Limits on Polarized Dust Spectral Index Variations for CMB Foreground Analysis
- Dark matter search with CMB: a study of foregrounds
- Constraining the polarisation flux density and angle of point sources by training a convolutional neural network
- Unbiased Primordial Gravitational Wave Inference from the CMB with SMICA
- Multi-frequency point source detection with fully convolutional networks: Performance in realistic microwave sky simulations