skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Wavelet-based techniques for the gamma-ray sky

Abstract

Here, we demonstrate how the image analysis technique of wavelet decomposition can be applied to the gamma-ray sky to separate emission on different angular scales. New structures on scales that differ from the scales of the conventional astrophysical foreground and background uncertainties can be robustly extracted, allowing a model-independent characterization with no presumption of exact signal morphology. As a test case, we generate mock gamma-ray data to demonstrate our ability to extract extended signals without assuming a fixed spatial template. For some point source luminosity functions, our technique also allows us to differentiate a diffuse signal in gamma-rays from dark matter annihilation and extended gamma-ray point source populations in a data-driven way.

Authors:
 [1];  [2];  [3];  [4]
  1. SUNY Stony Brook, Stony Brook, NY (United States)
  2. Fermi National Accelerator Lab. (FNAL), Batavia, IL (United States)
  3. The Johns Hopkins Univ., Baltimore, MD (United States); Fermi National Accelerator Lab. (FNAL), Batavia, IL (United States)
  4. Broad Institute, Cambridge, MA (United States); Princeton Univ., Princeton, NJ (United States)
Publication Date:
Research Org.:
Fermi National Accelerator Lab. (FNAL), Batavia, IL (United States)
Sponsoring Org.:
USDOE Office of Science (SC), High Energy Physics (HEP) (SC-25)
OSTI Identifier:
1341885
Report Number(s):
FERMILAB-PUB-15-504-T; YITP-SB-15-43; arXiv:1512.00012
Journal ID: ISSN 1475-7516; 1407425; TRN: US1701809
Grant/Contract Number:
AC02-07CH11359
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Cosmology and Astroparticle Physics
Additional Journal Information:
Journal Volume: 2016; Journal Issue: 07; Journal ID: ISSN 1475-7516
Publisher:
Institute of Physics (IOP)
Country of Publication:
United States
Language:
English
Subject:
79 ASTRONOMY AND ASTROPHYSICS; 46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND TECHNOLOGY; 72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; dark matter theory; gamma ray detectors

Citation Formats

McDermott, Samuel D., Fox, Patrick J., Cholis, Ilias, and Lee, Samuel K. Wavelet-based techniques for the gamma-ray sky. United States: N. p., 2016. Web. doi:10.1088/1475-7516/2016/07/045.
McDermott, Samuel D., Fox, Patrick J., Cholis, Ilias, & Lee, Samuel K. Wavelet-based techniques for the gamma-ray sky. United States. doi:10.1088/1475-7516/2016/07/045.
McDermott, Samuel D., Fox, Patrick J., Cholis, Ilias, and Lee, Samuel K. 2016. "Wavelet-based techniques for the gamma-ray sky". United States. doi:10.1088/1475-7516/2016/07/045. https://www.osti.gov/servlets/purl/1341885.
@article{osti_1341885,
title = {Wavelet-based techniques for the gamma-ray sky},
author = {McDermott, Samuel D. and Fox, Patrick J. and Cholis, Ilias and Lee, Samuel K.},
abstractNote = {Here, we demonstrate how the image analysis technique of wavelet decomposition can be applied to the gamma-ray sky to separate emission on different angular scales. New structures on scales that differ from the scales of the conventional astrophysical foreground and background uncertainties can be robustly extracted, allowing a model-independent characterization with no presumption of exact signal morphology. As a test case, we generate mock gamma-ray data to demonstrate our ability to extract extended signals without assuming a fixed spatial template. For some point source luminosity functions, our technique also allows us to differentiate a diffuse signal in gamma-rays from dark matter annihilation and extended gamma-ray point source populations in a data-driven way.},
doi = {10.1088/1475-7516/2016/07/045},
journal = {Journal of Cosmology and Astroparticle Physics},
number = 07,
volume = 2016,
place = {United States},
year = 2016,
month = 7
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 2works
Citation information provided by
Web of Science

Save / Share:
  • Low-resolution isotope identifiers are widely deployed for nuclear security purposes, but these detectors currently demonstrate problems in making correct identifications in many typical usage scenarios. While there are many hardware alternatives and improvements that can be made, performance on existing low resolution isotope identifiers should be able to be improved by developing new identification algorithms. We have developed a wavelet-based peak extraction algorithm and an implementation of a Bayesian classifier for automated peak-based identification. The peak extraction algorithm has been extended to compute uncertainties in the peak area calculations. To build empirical joint probability distributions of the peak areas andmore » uncertainties, a large set of spectra were simulated in MCNP6 and processed with the wavelet-based feature extraction algorithm. Kernel density estimation was then used to create a new component of the likelihood function in the Bayesian classifier. Furthermore, identification performance is demonstrated on a variety of real low-resolution spectra, including Category I quantities of special nuclear material.« less
  • Datasets of the geometry of fracture systems were analysed at various scales in the western Arabian sedimentary platform by means of geostatistical, multifractal, and anisotropic-wavelet techniques. The investigations covered a wide range of scales, from regional to outcrops in a well-exposed area, and were based on field mapping of fractures, and the interpretation and digitizing of fracture patterns on aerial photographs and satellite images. As a first step, fracture data sets were used to examine the direction, size, spacing and density systematics, and the variability in these quantities with space and scale. Secondly, a multifractal analysis was carried out, whichmore » consists in estimating the moments of the spatial distribution of fractures at different resolutions. This global multifractal method was complemented by a local wavelet analysis, using a new anisotropic technique tailored to linear structures. For a map with a given scale of detail, this procedure permits to define integrated fracture patterns and their associated directions at a more regional scale. The main result of this combined approach is that fracturing is not a self-similar process from the centimeter scale up to the one-million-kilometer scale. Spatial distribution of faults appears as being highly controlled by the thickness of the different rheological layers that constitute the crust. A proceeding for upscaling fracture systems in sedimentary reservoirs can be proposed, based on (i) a power law for joint-length distribution, (ii) characteristic joint spacing depending on the critical sedimentary units, and (iii) fractal fault geometry for faults larger than the whole thickness of the sedimentary basin.« less
  • A method for applying goal-based adaptive methods to the angular resolution of the neutral particle transport equation is presented. The methods are applied to an octahedral wavelet discretisation of the spherical angular domain which allows for anisotropic resolution. The angular resolution is adapted across both the spatial and energy dimensions. The spatial domain is discretised using an inner-element sub-grid scale finite element method. The goal-based adaptive methods optimise the angular discretisation to minimise the error in a specific functional of the solution. The goal-based error estimators require the solution of an adjoint system to determine the importance to the specifiedmore » functional. The error estimators and the novel methods to calculate them are described. Several examples are presented to demonstrate the effectiveness of the methods. It is shown that the methods can significantly reduce the number of unknowns and computational time required to obtain a given error. The novelty of the work is the use of goal-based adaptive methods to obtain anisotropic resolution in the angular domain for solving the transport equation. -- Highlights: •Wavelet angular discretisation used to solve transport equation. •Adaptive method developed for the wavelet discretisation. •Anisotropic angular resolution demonstrated through the adaptive method. •Adaptive method provides improvements in computational efficiency.« less
  • The authors represent the standard ramp filter operator of the filtered-back-projection (FBP) reconstruction in different bases composed of Haar and Daubechies compactly supported wavelets. The resulting multiscale representation of the ramp-filter matrix operator is approximately diagonal. The accuracy of this diagonal approximation becomes better as wavelets with larger numbers of vanishing moments are used. This wavelet-based representation enables them to formulate a multiscale tomographic reconstruction technique in which the object is reconstructed at multiple scales or resolutions. A complete reconstruction is obtained by combining the reconstructions at different scales. The multiscale reconstruction technique has the same computational complexity as themore » FBP reconstruction method. It differs from other multiscale reconstruction techniques in that (1) the object is defined through a one-dimensional multiscale transformation of the projection domain, and (2) they explicitly account for noise in the projections data by calculating maximum a posteriori probability (MAP) multiscale reconstruction estimates based on a chosen fractal prior on the multiscale object coefficients. The computational complexity of this maximum a posteriori probability (MAP) solution is also the same as that of the FBP reconstruction. This result is in contrast to commonly used methods of statistical regularization, which result in computationally intensive optimization algorithms.« less
  • We show that the statistical properties of the large scales of the Kuramoto-Sivashinsky equation in the extended system limit can be understood in terms of the dynamical behavior of the same equation in a small finite domain. Our method relies on the description of the solutions of this equation in terms of wavelets, and allows us to model the energy transfer between small and large scales. We show that the effective equation obtained in this way can be consistently approximated by a forced Burgers equation only for scales far from the cutoff between small and large wavelengths. {copyright} {ital 1996more » The American Physical Society.}« less