Poisson Denoising on the Sphere: Application to the Fermi Gamma Ray Space Telescope

Abstract

The Large Area Telescope (LAT), the main instrument of the Fermi Gamma-Ray Space Telescope, detects high energy gamma rays with energies from 20 MeV to more than 300 GeV. The two main scientific ob jectives, the study of the Milky Way diffuse background and the detection of point sources, are complicated by the lack of photons. That is why we need a powerful Poisson noise removal method on the sphere which is efficient on low count Poisson data. This paper presents a new multiscale decomposition on the sphere for data with Poisson noise, called Multi-Scale Variance Stabilizing Transform on the Sphere (MS-VSTS). This method is based on a Variance Stabilizing Transform (VST), a transform which aims to stabilize a Poisson data set such that each stabilized sample has a quasi constant variance. In addition, for the VST used in the method, the transformed data are asymptotically Gaussian. MS-VSTS consists of decomposing the data into a sparse multi-scale dictionary like wavelets or curvelets, and then applying a VST on the coefficients in order to get almost Gaussian stabilized coefficients. In this work, we use the Isotropic Undecimated Wavelet Transform (IUWT) and the Curvelet Transform as spherical multi-scale transforms. Then, binary hypothesis testing is carried out to detect significant coefficients, and the denoised image is reconstructed with an iterative algorithm based on Hybrid Steepest Descent (HSD). To detect point sources, we have to extract the Galactic diffuse background: an extension of the method to background separation is then proposed. In contrary, to study the Milky Way diffuse background, we remove point sources with a binary mask. The gaps have to be interpolated: an extension to inpainting is then proposed. The method, applied on simulated Fermi LAT data, proves to be adaptive, fast and easy to implement.

16 pages, 8 figures

References in corpus (3)

• Fermi Large Area Telescope Bright Gamma-ray Source List
• CMB data analysis and sparsity
• Source detection using a 3D sparse representation: application to the Fermi gamma-ray space telescope

Cited by in corpus (15)

• The Denoised, Deconvolved, and Decomposed Fermi $\gamma$-ray sky - An application of the D${}^3$PO algorithm
• Exact Wavelets on the Ball
• Distinguishing Dark Matter from Unresolved Point Sources in the Inner Galaxy with Photon Statistics
• S2LET: A code to perform fast wavelet analysis on the sphere
• Localisation of directional scale-discretised wavelets on the sphere
• Planck CMB Anomalies: Astrophysical and Cosmological Secondary Effects and the Curse of Masking
• Astronomical Image Denoising Using Dictionary Learning
• D${}^3$PO - Denoising, Deconvolving, and Decomposing Photon Observations
• Low-l CMB Analysis and Inpainting
• Spherical 3D Isotropic Wavelets
• Wavelet-Based Techniques for the Gamma-Ray Sky
• Fast directional spatially localized spherical harmonic transform
• Multichannel Poisson denoising and deconvolution on the sphere : Application to the Fermi Gamma Ray Space Telescope
• Prospects on high-energy source searches based on pattern recognition: Object detection in the H.E.S.S. Galactic Plane Survey and catalogue cross-matches
• Learning sparse representations on the sphere