Denoising, deconvolving and decomposing multi-domain photon observations- The D4PO algorithm

Abstract

Astronomical imaging based on photon count data is a non-trivial task. In this context we show how to denoise, deconvolve, and decompose multi-domain photon observations. The primary objective is to incorporate accurate and well motivated likelihood and prior models in order to give reliable estimates about morphologically different but superimposed photon flux components present in the data set. Thereby we denoise and deconvolve photon counts, while simultaneously decomposing them into diffuse, point-like and uninteresting background radiation fluxes. The decomposition is based on a probabilistic hierarchical Bayesian parameter model within the framework of information field theory (IFT). In contrast to its predecessor D${}^3$PO, D${}^4$PO reconstructs multi-domain components. Thereby each component is defined over its own direct product of multiple independent domains, for example location and energy. D${}^4$PO has the capability to reconstruct correlation structures over each of the sub-domains of a component separately. Thereby the inferred correlations implicitly define the morphologically different source components, except for the spatial correlations of the point-like flux. Point-like source fluxes are spatially uncorrelated by definition. The capabilities of the algorithm are demonstrated by means of a synthetic, but realistic, mock data set, providing spectral and spatial information about each detected photon. D${}^4$PO successfully denoised, deconvolved, and decomposed a photon count image into diffuse, point-like and background flux, each being functions of location as well as energy. Moreover, uncertainty estimates of the reconstructed fields as well as of their correlation structure are provided employing their posterior density function and accounting for the manifolds the domains reside on.

References in corpus (36)

• The Large Area Telescope on the Fermi Gamma-ray Space Telescope Mission
• A multi-scale multi-frequency deconvolution algorithm for synthesis imaging in radio interferometry
• Restoration of Poissonian Images Using Alternating Direction Optimization
• Multiplicative Noise Removal Using Variable Splitting and Constrained Optimization
• Planck Early Results. VII. The Early Release Compact Source Catalog
• A proximal iteration for deconvolving Poisson noisy images using sparse representations
• Rejuvenating the matter power spectrum: restoring information with a logarithmic density mapping
• Information field theory for cosmological perturbation reconstruction and non-linear signal analysis
• Probability Distribution Function of Cosmological Density Fluctuations from Gaussian Initial Condition: Comparison of One- and Two-point Log-normal Model Predictions with N-body Simulations
• The Scale of the Problem : Recovering Images of Reionization with GMCA
• Fast Hamiltonian sampling for large scale structure inference
• Bayesian power-spectrum inference for Large Scale Structure data
• The Mexican Hat Wavelet Family. Application to point source detection in CMB maps
• A fast Bayesian approach to discrete object detection in astronomical datasets - PowellSnakes I
• Reconstruction of Gaussian and log-normal fields with spectral smoothness
• NIFTY - Numerical Information Field Theory - a versatile Python library for signal inference
• Cosmic Cartography of the Large-Scale Structure with Sloan Digital Sky Survey Data Release 6
• The Denoised, Deconvolved, and Decomposed Fermi $\gamma$-ray sky - An application of the D${}^3$PO algorithm
• SZ and CMB reconstruction using Generalized Morphological Component Analysis
• Reconstruction of signals with unknown spectra in information field theory with parameter uncertainty
• Positive deconvolution for superimposed extended source and point sources
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• Comparison of source detection procedures for XMM-Newton images
• Maximum-A-Posteriori Estimates in Linear Inverse Problems with Log-concave Priors are Proper Bayes Estimators
• Methods for Bayesian power spectrum inference with galaxy surveys
• Inference with minimal Gibbs free energy in information field theory
• Powellsnakes II: a fast Bayesian approach to discrete object detection in multi-frequency astronomical data sets
• Numerical Simulation of Non-Gaussian Random Fields with Prescribed Correlation Structure
• D${}^3$PO - Denoising, Deconvolving, and Decomposing Photon Observations
• Background-Source Separation in astronomical images with Bayesian probability theory (I): the method
• Bayesian weak lensing tomography: Reconstructing the 3D large-scale distribution of matter with a lognormal prior
• Constrained realizations and minimum variance reconstruction of non-Gaussian random fields
• Rapid and deterministic estimation of probability densities using scale-free field theories
• Search for quasi-periodic signals in magnetar giant flares
• Dynamic system classifier
• Simultaneous analysis of large INTEGRAL/SPI datasets: optimizing the computation of the solution and its variance using sparse matrix algorithms

Cited by in corpus (7)

• Information theory for fields
• Information Field Theory and Artificial Intelligence
• Multi-Component Imaging of the Fermi Gamma-ray Sky in the Spatio-spectral Domain
• First spatio-spectral Bayesian imaging of SN1006 in X-ray
• Radio Imaging With Information Field Theory
• Exploring the MeV Sky with a Combined Coded Mask and Compton Telescope: The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO)
• Bayesian Multi-wavelength Imaging of the LMC SN1987A with SRG/eROSITA