A new approach for the statistical denoising of Planck interstellar dust polarization data

Abstract

Dust emission is the main foreground for cosmic microwave background (CMB) polarization. Its statistical characterization must be derived from the analysis of observational data because the precision required for a reliable component separation is far greater than what is currently achievable with physical models of the turbulent magnetized interstellar medium. This letter takes a significant step toward this goal by proposing a method that retrieves non-Gaussian statistical characteristics of dust emission from noisy Planck polarization observations at 353 GHz. We devised a statistical denoising method based on wavelet phase harmonics (WPH) statistics, which characterize the coherent structures in non-Gaussian random fields and define a generative model of the data. The method was validated on mock data combining a dust map from a magnetohydrodynamic simulation and Planck noise maps. The denoised map reproduces the true power spectrum down to scales where the noise power is an order of magnitude larger than that of the signal. It remains highly correlated to the true emission and retrieves some of its non-Gaussian properties. Applied to Planck data, the method provides a new approach to building a generative model of dust polarization that will characterize the full complexity of the dust emission. We also release PyWPH, a public Python package, to perform GPU-accelerated WPH analyses on images.

arXiv comment10 pages, 9 figures, accepted by A&A

References in corpus (31)

• Beyond a Gaussian Denoiser: Residual Learning of Deep CNN for Image Denoising
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• A Joint Analysis of BICEP2/Keck Array and Planck Data
• Planck 2018 results. IV. Diffuse component separation
• Planck 2018 results. III. High Frequency Instrument data processing and frequency maps
• Wiener Reconstruction of The Large Scale Structure
• Global, Exact Cosmic Microwave Background Data Analysis Using Gibbs Sampling
• Planck 2018 results. XI. Polarized dust foregrounds
• Density Studies of MHD Interstellar Turbulence: Statistical Moments, Correlations and Bispectrum
• Planck 2018 results. XII. Galactic astrophysics using polarized dust emission
• A new approach to observational cosmology using the scattering transform
• Mapping the Magnetic Interstellar Medium in Three Dimensions Over the Full Sky with Neutral Hydrogen
• New Interpretable Statistics for Large Scale Structure Analysis and Generation
• Modeling the magnetized Local Bubble from dust data
• Cosmological parameters, shear maps and power spectra from CFHTLenS using Bayesian hierarchical inference
• The RWST, a comprehensive statistical description of the non-Gaussian structures in the ISM
• Efficient Wiener filtering without preconditioning
• Improving Weak Lensing Mass Map Reconstructions using Gaussian and Sparsity Priors: Application to DES SV
• Full-sky Cosmic Microwave Background Foreground Cleaning Using Machine Learning
• Mitigating Complex Dust Foregrounds in Future CMB Polarization Experiments
• Classification of Magnetohydrodynamic Simulations using Wavelet Scattering Transforms
• Statistical description of dust polarized emission from the diffuse interstellar medium -- A RWST approach
• Intermittency of interstellar turbulence: Parsec-scale coherent structure of intense velocity-shear
• Impact of Polarized Galactic Foreground Emission on CMB Lensing Reconstruction and Delensing of B-Modes
• Matrix-free Large Scale Bayesian inference in cosmology
• ForSE: a GAN based algorithm for extending CMB foreground models to sub-degree angular scales
• Cleaning our own Dust: Simulating and Separating Galactic Dust Foregrounds with Neural Networks
• Maximum Entropy Models from Phase Harmonic Covariances
• Wiener filtering and pure E/B decomposition of CMB maps with anisotropic correlated noise
• A Generative Model of Galactic Dust Emission Using Variational Inference
• Fast Sampling from Wiener Posteriors for Image Data with Dataflow Engines

Cited by in corpus (12)

• Non-Gaussian modelling and statistical denoising of Planck dust polarization full-sky maps using scattering transforms
• The Atacama Cosmology Telescope: Map-Based Noise Simulations for DR6
• Single frequency CMB B-mode inference with realistic foregrounds from a single training image
• Full-sky, arcminute-scale, 3D models of Galactic microwave foreground dust emission based on filaments
• Scattering Spectra Models for Physics
• Generative models of astrophysical fields with scattering transforms on the sphere
• Generative Models of Multi-channel Data from a Single Example -- Application to Dust Emission
• Self-supervised component separation for the extragalactic submillimeter sky
• From few to many maps: A fast map-level emulator for extreme augmentation of CMB systematics datasets
• Quantitative Morphology of Galactic Cirrus in Deep Optical Imaging
• Constraining Power of Wavelet vs. Power Spectrum Statistics for CMB Lensing and Weak Lensing with Learned Binning
• Bayesian imaging inverse problem with scattering transform