D${}^3$PO - Denoising, Deconvolving, and Decomposing Photon Observations

Abstract

The analysis of astronomical images is a non-trivial task. The D3PO algorithm addresses the inference problem of denoising, deconvolving, and decomposing photon observations. Its primary goal is the simultaneous but individual reconstruction of the diffuse and point-like photon flux given a single photon count image, where the fluxes are superimposed. In order to discriminate between these morphologically different signal components, a probabilistic algorithm is derived in the language of information field theory based on a hierarchical Bayesian parameter model. The signal inference exploits prior information on the spatial correlation structure of the diffuse component and the brightness distribution of the spatially uncorrelated point-like sources. A maximum a posteriori solution and a solution minimizing the Gibbs free energy of the inference problem using variational Bayesian methods are discussed. Since the derivation of the solution is not dependent on the underlying position space, the implementation of the D3PO algorithm uses the NIFTY package to ensure applicability to various spatial grids and at any resolution. The fidelity of the algorithm is validated by the analysis of simulated data, including a realistic high energy photon count image showing a 32 x 32 arcmin^2 observation with a spatial resolution of 0.1 arcmin. In all tests the D3PO algorithm successfully denoised, deconvolved, and decomposed the data into a diffuse and a point-like signal estimate for the respective photon flux components.

22 pages, 8 figures, 2 tables, accepted by Astronomy & Astrophysics; refereed version, 1 figure added, results unchanged, software available at http://www.mpa-garching.mpg.de/ift/d3po/

References in corpus (13)

• Multi-Scale CLEAN deconvolution of radio synthesis images
• A proximal iteration for deconvolving Poisson noisy images using sparse representations
• The Scale of the Problem : Recovering Images of Reionization with GMCA
• 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
• Reconstruction of signals with unknown spectra in information field theory with parameter uncertainty
• 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
• Reply to "Comment on `Inference with minimal Gibbs free energy in information field theory'" by Iatsenko, Stefanovska and McClintock
• Rapid and deterministic estimation of probability densities using scale-free field theories
• Astrophysical data analysis with information field theory
• Simultaneous analysis of large INTEGRAL/SPI datasets: optimizing the computation of the solution and its variance using sparse matrix algorithms

Cited by in corpus (21)

• IMAGINE: A comprehensive view of the interstellar medium, Galactic magnetic fields and cosmic rays
• The Denoised, Deconvolved, and Decomposed Fermi $\gamma$-ray sky - An application of the D${}^3$PO algorithm
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• Distinguishing Dark Matter from Unresolved Point Sources in the Inner Galaxy with Photon Statistics
• Information theory for fields
• Firedec: a two-channel finite-resolution image deconvolution algorithm
• A new approach to multi-frequency synthesis in radio interferometry
• Image deconvolution and PSF reconstruction with STARRED: a wavelet-based two-channel method optimized for light-curve extraction
• Search for quasi-periodic signals in magnetar giant flares
• STARRED: a two-channel deconvolution method with Starlet regularization
• Generic inference of inflation models by non-Gaussianity and primordial power spectrum reconstruction
• Denoising, deconvolving and decomposing multi-domain photon observations- The D4PO algorithm
• Multi-Component Imaging of the Fermi Gamma-ray Sky in the Spatio-spectral Domain
• First spatio-spectral Bayesian imaging of SN1006 in X-ray
• Sharpening up Galactic all-sky maps with complementary data - A machine learning approach
• A study of Neural networks point source extraction on simulated Fermi/LAT Telescope images
• Exploring the MeV Sky with a Combined Coded Mask and Compton Telescope: The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO)
• GeV excess and phenomenological astrophysics modeling
• Latent-space Field Tension for Astrophysical Component Detection An application to X-ray imaging
• Bayesian Multi-wavelength Imaging of the LMC SN1987A with SRG/eROSITA
• Machine learning application to Fermi-LAT data: sharpening all-sky map and emphasizing variable sources