Efficient Wiener filtering without preconditioning

Abstract

We present a new approach to calculate the Wiener filter solution of general data sets. It is trivial to implement, flexible, numerically absolutely stable, and guaranteed to converge. Most importantly, it does not require an ingenious choice of preconditioner to work well. The method is capable of taking into account inhomogeneous noise distributions and arbitrary mask geometries. It iteratively builds up the signal reconstruction by means of a messenger field, introduced to mediate between the different preferred bases in which signal and noise properties can be specified most conveniently. Using cosmic microwave background (CMB) radiation data as a showcase, we demonstrate the capabilities of our scheme by computing Wiener filtered WMAP7 temperature and polarization maps at full resolution for the first time. We show how the algorithm can be modified to synthesize fluctuation maps, which, combined with the Wiener filter solution, result in unbiased constrained signal realizations, consistent with the observations. The algorithm performs well even on simulated CMB maps with Planck resolution and dynamic range.

5 pages, 2 figures. Submitted to Astronomy and Astrophysics. Replaced to match published version

References in corpus (27)

• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Likelihoods and Parameters from the WMAP data
• Three Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Polarization Analysis
• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Sky Maps, Systematic Errors, and Basic Results
• Three-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Temperature Analysis
• Planck Early Results: The Planck mission
• Estimating the Power Spectrum of the Cosmic Microwave Background
• How to measure CMB power spectra without losing information
• Detection of Gravitational Lensing in the Cosmic Microwave Background
• Measuring primordial non-Gaussianity in the cosmic microwave background
• How to make maps from CMB data without losing information
• Global, Exact Cosmic Microwave Background Data Analysis Using Gibbs Sampling
• Optimal limits on f_{NL}^{local} from WMAP 5-year data
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Data Processing Methods and Systematic Errors Limits
• An Efficient Technique to Determine the Power Spectrum from Cosmic Microwave Background Sky Maps
• Bayesian power-spectrum inference for Large Scale Structure data
• Cross-correlation of CMB with large-scale structure: weak gravitational lensing
• Estimation of Polarized Power Spectra by Gibbs sampling
• Reconstruction of signals with unknown spectra in information field theory with parameter uncertainty
• Non-Gaussianity and the CMB Bispectrum: confusion between Primordial and Lensing-Rees Sciama contribution?
• A high-resolution map of the cosmic microwave background around the north celestial pole
• Reconstructing signals from noisy data with unknown signal and noise covariance
• Local non-Gaussianity in the Cosmic Microwave Background the Bayesian way
• Likelihood, Fisher information, and systematics of cosmic microwave background experiments
• Fast calculation of the Fisher matrix for cosmic microwave background experiments
• ARKCoS: Artifact-Suppressed Accelerated Radial Kernel Convolution on the Sphere

Cited by in corpus (55)

• Planck 2013 results. XXIII. Isotropy and statistics of the CMB
• Planck 2013 Results. XXIV. Constraints on primordial non-Gaussianity
• Detection of B-mode Polarization in the Cosmic Microwave Background with Data from the South Pole Telescope
• Physical Bayesian modelling of the non-linear matter distribution: new insights into the Nearby Universe
• Towards optimal extraction of cosmological information from nonlinear data
• Review: Deep Learning in Electron Microscopy
• Dark Energy Survey Year 3 results: curved-sky weak lensing mass map reconstruction
• Unmasking the Masked Universe: the 2M++ catalogue through Bayesian eyes
• Bayesian optimisation for likelihood-free cosmological inference
• Hierarchical Cosmic Shear Power Spectrum Inference
• Cosmological parameters, shear maps and power spectra from CFHTLenS using Bayesian hierarchical inference
• Bayesian delensing delight: sampling-based inference of the primordial CMB and gravitational lensing
• Bayesian delensing of CMB temperature and polarization
• Reconstructing cosmic growth with kSZ observations in the era of Stage IV experiments
• Improving Weak Lensing Mass Map Reconstructions using Gaussian and Sparsity Priors: Application to DES SV
• The skewed weak lensing likelihood: why biases arise, despite data and theory being sound
• Unbiased methods for removing systematics from galaxy clustering measurements
• Unbiased pseudo-Cl power spectrum estimation with mode projection
• Probabilistic Mass Mapping with Neural Score Estimation
• Bayesian inference of CMB gravitational lensing
• Using inpainting to construct accurate cut-sky CMB estimators
• Matrix-free Large Scale Bayesian inference in cosmology
• A new approach for the statistical denoising of Planck interstellar dust polarization data
• Efficient Optimal Reconstruction of Linear Fields and Band-powers from Cosmological Data
• Weak lensing mass reconstruction using sparsity and a Gaussian random field
• Multi-resolution Bayesian CMB component separation through Wiener-filtering with a pseudo-inverse preconditioner
• A new probe of Axion-Like Particles: CMB polarization distortions due to cluster magnetic fields
• Wiener filter reloaded: fast signal reconstruction without preconditioning
• Non-Gaussianity of secondary anisotropies from ACTPol and Planck
• Wiener filtering and pure E/B decomposition of CMB maps with anisotropic correlated noise
• Pure E and B polarization maps via Wiener filtering
• Optimal bispectrum estimator and simulations of the the CMB Lensing-ISW non-Gaussian signal
• Bayesian power-spectrum inference with foreground and target contamination treatment
• A multi-level solver for Gaussian constrained CMB realizations
• Semi-blind Bayesian inference of CMB map and power spectrum
• Extracting the Signal of Cosmic String Wakes from 21-cm Observations
• CMB bounds on tensor-scalar-scalar inflationary correlations
• Fast Sampling from Wiener Posteriors for Image Data with Dataflow Engines
• The causal effect of environment on halo mass and concentration
• Optimal and fast E/B separation with a dual messenger field
• Cosmic Microwave Background Mapmaking with a Messenger Field
• Sparse Inpainting and Isotropy
• Solving linear equations with messenger-field and conjugate gradients techniques - an application to CMB data analysis
• A close examination of cosmic microwave background mirror-parity after Planck
• Preconditioner-free Wiener filtering with a dense noise matrix
• Towards efficient and optimal analysis of CMB anisotropies on a masked sky
• Probabilistic Cosmological Mass Mapping from Weak Lensing Shear
• Enhancing CMB map reconstruction and power spectrum estimation with convolutional neural networks
• Testing Quadratic Maximum Likelihood estimators for forthcoming Stage-IV weak lensing surveys
• Optimal estimator for the amplitude of the bispectrum from IR sources
• Study on the filters of atmospheric contamination in ground based CMB observation
• Improved Gibbs samplers for Cosmic Microwave Background power spectrum estimation
• Compressed convolution
• Pixel domain multi-resolution minimum variance painting of CMB maps
• Cooling Improves Cosmic Microwave Background Map-Making When Low-Frequency Noise is Large