Filling in CMB map missing data using constrained Gaussian realizations

Abstract

For analyzing maps of the cosmic microwave background sky, it is necessary to mask out the region around the galactic equator where the parasitic foreground emission is strongest as well as the brightest compact sources. Since many of the analyses of the data, particularly those searching for non-Gaussianity of a primordial origin, are most straightforwardly carried out on full-sky maps, it is of great interest to develop efficient algorithms for filling in the missing information in a plausible way. We explore practical algorithms for filling in based on constrained Gaussian realizations. Although carrying out such realizations is in principle straightforward, for finely pixelized maps as will be required for the Planck analysis a direct brute force method is not numerically tractable. We present some concrete solutions to this problem, both on a spatially flat sky with periodic boundary conditions and on the pixelized sphere. One approach is to solve the linear system with an appropriately preconditioned conjugate gradient method. While this approach was successfully implemented on a rectangular domain with periodic boundary conditions and worked even for very wide masked regions, we found that the method failed on the pixelized sphere for reasons that we explain here. We present an approach that works for full-sky pixelized maps on the sphere involving a kernel-based multi-resolution Laplace solver followed by a series of conjugate gradient corrections near the boundary of the mask.

22 pages, 14 figures, minor changes, a few missing references added

References in corpus (4)

• A Measurement of the Damping Tail of the Cosmic Microwave Background Power Spectrum with the South Pole Telescope
• CMB data analysis and sparsity
• Avoiding bias in reconstructing the largest observable scales from partial-sky data
• Harmonic Inpainting of the Cosmic Microwave Background Sky: Formulation and Error Estimate

Cited by in corpus (32)

• Planck 2018 results. VII. Isotropy and Statistics of the CMB
• The Atacama Cosmology Telescope: A Measurement of the DR6 CMB Lensing Power Spectrum and its Implications for Structure Growth
• The Atacama Cosmology Telescope: Two-Season ACTPol Lensing Power Spectrum
• CMB Lensing Power Spectrum Biases from Galaxies and Clusters using High-angular Resolution Temperature Maps
• The Atacama Cosmology Telescope: A CMB lensing mass map over 2100 square degrees of sky and its cross-correlation with BOSS-CMASS galaxies
• The Atacama Cosmology Telescope: Component-separated maps of CMB temperature and the thermal Sunyaev-Zel'dovich effect
• Large-scale alignments from WMAP and Planck
• The Atacama Cosmology Telescope: Lensing of CMB Temperature and Polarization Derived from Cosmic Infrared Background Cross-Correlation
• Exploring Cosmic Origins with CORE: B-mode Component Separation
• SPT-3G D1: CMB temperature and polarization power spectra and cosmology from 2019 and 2020 observations of the SPT-3G Main field
• Reconstructing cosmic growth with kSZ observations in the era of Stage IV experiments
• Cross-correlating Astrophysical and Cosmological Gravitational Wave Backgrounds with the Cosmic Microwave Background
• Planck CMB Anomalies: Astrophysical and Cosmological Secondary Effects and the Curse of Masking
• Methods for pixel domain correction of EB leakage
• Sparsity and the Bayesian Perspective
• The Atacama Cosmology Telescope: Map-Based Noise Simulations for DR6
• Inpainting Galactic Foreground Intensity and Polarization maps using Convolutional Neural Network
• On Preferred Axes in WMAP Cosmic Microwave Background Data after Subtraction of the Integrated Sachs-Wolfe Effect
• A Generative Model of Galactic Dust Emission Using Variational Inference
• Inpainting CMB maps using Partial Convolutional Neural Networks
• Sparse point-source removal for full-sky CMB experiments: application to WMAP 9-year data
• Unraveling the CMB lack-of-correlation anomaly with the cosmological gravitational wave background
• Slepian models for Gaussian Random Landscapes
• Level correlations of the CMB temperature angular power spectrum
• Reconstruction of Missing Data using Iterative Harmonic Expansion
• Free-flight experiments in LISA Pathfinder
• The slingshot effect as a probe of transverse motions of galaxies
• CMB-PAInT: An inpainting tool for the cosmic microwave background
• Filtering in CMB data analysis with application to ACT DR4 and Planck observations
• Mitigating point-source contamination in CMB polarization: a Generalized Point Spread Function fitting approach
• Removing the mask -- reconstructing a scalar field on the sphere from a masked field
• Test of the Statistical Isotropy of the Universe using Gravitational Waves