An official website of the United States government
Data from high-energy observations are usually obtained as lists of photon events. A common analysis task for such data is to identify whether diffuse emission exists, and to estimate its surface brightness, even in the presence of point sources that may be superposed. We have developed a novel nonparametric event list segmentation algorithm to divide up the field of view into distinct emission components. We use photon location data directly, without binning them into an image. We first construct a graph from the Voronoi tessellation of the observed photon locations and then grow segments using a new adaptation of seeded region growing that we call Seeded Region Growing on Graph, after which the overall method is named SRGonG. Starting with a set of seed locations, this results in an oversegmented data set, which SRGonG then coalesces using a greedy algorithm where adjacent segments are merged to minimize a model comparison statistic; we use the Bayesian Information Criterion. Using SRGonG we are able to identify point-like and diffuse extended sources in the data with equal facility. We validate SRGonG using simulations, demonstrating that it is capable of discerning irregularly shaped low-surface-brightness emission structures as well as point-like sources with strengths comparable to that seen in typical X-ray data. We demonstrate SRGonG's use on the Chandra data of the Antennae galaxies and show that it segments the complex structures appropriately.
more »
« less
Identifying Diffuse Spatial Structures in High-energy Photon Lists
Abstract Data from high-energy observations are usually obtained as lists of photon events. A common analysis task for such data is to identify whether diffuse emission exists, and to estimate its surface brightness, even in the presence of point sources that may be superposed. We have developed a novel nonparametric event list segmentation algorithm to divide up the field of view into distinct emission components. We use photon location data directly, without binning them into an image. We first construct a graph from the Voronoi tessellation of the observed photon locations and then grow segments using a new adaptation of seeded region growing that we callSeeded Region Growing on Graph, after which the overall method is namedSRGonG. Starting with a set of seed locations, this results in an oversegmented data set, whichSRGonGthen coalesces using a greedy algorithm where adjacent segments are merged to minimize a model comparison statistic; we use the Bayesian Information Criterion. UsingSRGonGwe are able to identify point-like and diffuse extended sources in the data with equal facility. We validateSRGonGusing simulations, demonstrating that it is capable of discerning irregularly shaped low-surface-brightness emission structures as well as point-like sources with strengths comparable to that seen in typical X-ray data. We demonstrateSRGonG’s use on the Chandra data of the Antennae galaxies and show that it segments the complex structures appropriately.
more »
« less
- PAR ID:
- 10392439
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astronomical Journal
- Volume:
- 165
- Issue:
- 2
- ISSN:
- 0004-6256
- Format(s):
- Medium: X Size: Article No. 66
- Size(s):
- Article No. 66
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Galactic gamma-ray diffuse emission (GDE) is emitted by cosmic rays (CRs), ultra-relativistic protons, and electrons, interacting with gas and electromagnetic radiation fields in the interstellar medium. Here we present the analysis of teraelectronvolt diffuse emission from a region of the Galactic plane over the range in longitude ofl∈ [43°, 73°], using data collected with the High Altitude Water Cherenkov (HAWC) detector. Spectral, longitudinal, and latitudinal distributions of the teraelectronvolt diffuse emission are shown. The radiation spectrum is compatible with the spectrum of the emission arising from a CR population with anindexsimilar to that of the observed CRs. When comparing with theDRAGONbase model, the HAWC GDE flux is higher by about a factor of 2. Unresolved sources such as pulsar wind nebulae and teraelectronvolt halos could explain the excess emission. Finally, deviations of the Galactic CR flux from the locally measured CR flux may additionally explain the difference between the predicted and measured diffuse fluxes.more » « less
-
We present a new method to estimate the boundary of extended sources in high-energy photon lists and to quantify the uncertainty in the boundary. This method extends the graphed seeded region growing method developed by M. Fan et al. Here, we describe how an unambiguous boundary of a centrally concentrated astronomical source may be defined by first spatially segmenting the photon list, then forcibly merging the segments until only two segments—an extended source and its background—remain, and finally constructing a boundary as the connected outer edges of the Voronoi tessellation of the photons included in the source segment. The resulting boundary is then modeled using Fourier descriptors to generate a smooth curve, and this curve is bootstrapped to generate uncertainties. We apply the method to photon event lists obtained during the observations of galaxies NGC 2300 and Arp 299. We demonstrate how the derived extent and enclosed flux of NGC 2300 obtained with Chandra and XMM-Newton are comparable. We also show how complex internal structure, as in the case of Arp 299, may be subsumed to construct a compact boundary of the object.more » « less
-
The Green Bank Telescope (GBT) Diffuse Ionized Gas Survey (GDIGS) traces ionized gas in the Galactic midplane by observing radio recombination line (RRL) emission from 4–8 GHz. The nominal survey zone is 32.3° > l > -5°, |b| < 0.5°. Here, we analyze GDIGS Hnα ionized gas emission toward discrete sources with sizes comparable to the 2.065' GDIGS Hnα beam. We use GDIGS data to identify the correct velocity of 39 H II regions that have multiple RRL velocity components. We identify and characterize RRL emission from 88 H II regions that previously lacked measured ionized gas velocities. We additionally identify and characterize RRL emission from eight locations that appear to be previously-unidentified H II regions and 41 locations of RRL emission that do not appear to be H II regions based on their lack of mid-infrared emission. We identify 10 discrete sources that have anomalously high RRL velocities for their locations in the Galactic plane and we compare the objects’ RRL data to 13CO, H I and mid-infrared data. These sources do not have the expected 24 μm emission characteristic of H II regions. Based on this comparison we do not think these objects are H II regions, but we are unable to classify them as a known type of object.more » « less
-
Abstract The Green Bank Telescope Diffuse Ionized Gas Survey (GDIGS) traces ionized gas in the Galactic midplane by observing radio recombination line (RRL) emission from 4 to 8 GHz. The nominal survey zone is 32.°3 >ℓ> −5°, ∣b∣ < 0.°5. Here, we analyze GDIGS Hnαionized gas emission toward discrete sources. Using GDIGS data, we identify the velocity of 35 Hiiregions that have multiple detected RRL velocity components. We identify and characterize RRL emission from 88 Hiiregions that previously lacked measured ionized gas velocities. We also identify and characterize RRL emission from eight locations that appear to be previously unidentified Hiiregions and 30 locations of RRL emission that do not appear to be Hiiregions based on their lack of mid-infrared emission. This latter group may be a compact component of the Galactic Diffuse Ionized Gas. There are an additional 10 discrete sources that have anomalously high RRL velocities for their locations in the Galactic plane. We compare these objects’ RRL data to13CO, Hi,and mid-infrared data, and find that these sources do not have the expected 24μm emission characteristic of Hiiregions. Based on this comparison we do not think these objects are Hiiregions, but we are unable to classify them as a known type of object.more » « less
