-
Signal-background separation and energy reconstruction of gamma rays using pattern spectra and convolutional neural networks for the Small-Sized Telescopes of the Cherenkov Telescope Array
Authors:
J. Aschersleben,
T. T. H. Arnesen,
R. F. Peletier,
M. Vecchi,
C. Vlasakidis,
M. H. F. Wilkinson
Abstract:
Imaging Atmospheric Cherenkov Telescopes (IACTs) detect very-high-energy gamma rays from ground level by capturing the Cherenkov light of the induced particle showers. Convolutional neural networks (CNNs) can be trained on IACT camera images of such events to differentiate the signal from the background and to reconstruct the energy of the initial gamma ray. Pattern spectra provide a 2-dimensional…
▽ More
Imaging Atmospheric Cherenkov Telescopes (IACTs) detect very-high-energy gamma rays from ground level by capturing the Cherenkov light of the induced particle showers. Convolutional neural networks (CNNs) can be trained on IACT camera images of such events to differentiate the signal from the background and to reconstruct the energy of the initial gamma ray. Pattern spectra provide a 2-dimensional histogram of the sizes and shapes of features comprising an image and they can be used as an input for a CNN to significantly reduce the computational power required to train it. In this work, we generate pattern spectra from simulated gamma-ray and proton images to train a CNN for signal-background separation and energy reconstruction for the Small-Sized Telescopes (SSTs) of the Cherenkov Telescope Array (CTA). A comparison of our results with a CNN directly trained on CTA images shows that the pattern spectra-based analysis is about a factor of three less computationally expensive but not able to compete with the performance of an CTA image-based analysis. Thus, we conclude that the CTA images must be comprised of additional information not represented by the pattern spectra.
△ Less
Submitted 30 November, 2023; v1 submitted 11 July, 2023;
originally announced July 2023.
-
Event reconstruction using pattern spectra and convolutional neural networks for the Cherenkov Telescope Array
Authors:
J. Aschersleben,
M. Vecchi,
M. H. F. Wilkinson,
R. F. Peletier
Abstract:
The Cherenkov Telescope Array (CTA) is the future observatory for ground-based imaging atmospheric Cherenkov telescopes. Each telescope will provide a snapshot of gamma-ray induced particle showers by capturing the induced Cherenkov emission at ground level. The simulation of such events provides camera images that can be used as training data for convolutional neural networks (CNNs) to differenti…
▽ More
The Cherenkov Telescope Array (CTA) is the future observatory for ground-based imaging atmospheric Cherenkov telescopes. Each telescope will provide a snapshot of gamma-ray induced particle showers by capturing the induced Cherenkov emission at ground level. The simulation of such events provides camera images that can be used as training data for convolutional neural networks (CNNs) to differentiate signals from background events and to determine the energy of the initial gamma-ray events. Pattern spectra are commonly used tools for image classification and provide the distributions of the sizes and shapes of features comprising an image. The application of pattern spectra on a CNN allows the selection of relevant combinations of features within an image. In this work, we generate pattern spectra from simulated gamma-ray images to train a CNN for signal-background separation and energy reconstruction for CTA. We compare our results to a CNN trained with CTA images and find that the pattern spectra-based analysis is computationally less expensive but not competitive with the purely CTA images-based analysis. Thus, we conclude that the CNN must rely on additional features in the CTA images not captured by the pattern spectra.
△ Less
Submitted 23 February, 2023;
originally announced February 2023.
-
A comparative study of source-finding techniques in HI emission line cubes using SoFiA, MTObjects, and supervised deep learning
Authors:
J. A. Barkai,
M. A. W. Verheijen,
E. T. Martínez,
M. H. F. Wilkinson
Abstract:
The 21 cm spectral line emission of atomic neutral hydrogen (HI) is one of the primary wavelengths observed in radio astronomy. However, the signal is intrinsically faint and the HI content of galaxies depends on the cosmic environment, requiring large survey volumes and survey depth to investigate the HI Universe. As the amount of data coming from these surveys continues to increase with technolo…
▽ More
The 21 cm spectral line emission of atomic neutral hydrogen (HI) is one of the primary wavelengths observed in radio astronomy. However, the signal is intrinsically faint and the HI content of galaxies depends on the cosmic environment, requiring large survey volumes and survey depth to investigate the HI Universe. As the amount of data coming from these surveys continues to increase with technological improvements, so does the need for automatic techniques for identifying and characterising HI sources while considering the tradeoff between completeness and purity. This study aimed to find the optimal pipeline for finding and masking the most sources with the best mask quality and the fewest artefacts in 3D neutral hydrogen cubes. Various existing methods were explored in an attempt to create a pipeline to optimally identify and mask the sources in 3D neutral hydrogen 21 cm spectral line data cubes. Two traditional source-finding methods were tested, SoFiA and MTObjects, as well as a new supervised deep learning approach, in which a 3D convolutional neural network architecture, known as V-Net was used. These three source-finding methods were further improved by adding a classical machine learning classifier as a post-processing step to remove false positive detections. The pipelines were tested on HI data cubes from the Westerbork Synthesis Radio Telescope with additional inserted mock galaxies. SoFiA combined with a random forest classifier provided the best results, with the V-Net-random forest combination a close second. We suspect this is due to the fact that there are many more mock sources in the training set than real sources. There is, therefore, room to improve the quality of the V-Net network with better-labelled data such that it can potentially outperform SoFiA.
△ Less
Submitted 23 November, 2022;
originally announced November 2022.
-
The Fornax Deep Survey (FDS) with VST XII: Low surface brightness dwarf galaxies in the Fornax cluster
Authors:
Aku Venhola,
Reynier F. Peletier,
Heikki Salo,
Eija Laurikainen,
Joachim Janz,
Caroline Haigh,
Michael H. F. Wilkinson,
Enrichetta Iodice,
Michael Hilker,
Steffen Mieske,
Michele Cantiello,
Marilena Spavone
Abstract:
In this work we use Max-Tree Objects, (MTO) on the FDS data in order to detect previously undetected Low surface brightness (LSB) galaxies. After extending the existing Fornax dwarf galaxy catalogs with this sample, our goal is to understand the evolution of LSB dwarfs in the cluster. We also study the contribution of the newly detected galaxies to the faint end of the luminosity function. We test…
▽ More
In this work we use Max-Tree Objects, (MTO) on the FDS data in order to detect previously undetected Low surface brightness (LSB) galaxies. After extending the existing Fornax dwarf galaxy catalogs with this sample, our goal is to understand the evolution of LSB dwarfs in the cluster. We also study the contribution of the newly detected galaxies to the faint end of the luminosity function. We test the detection completeness and parameter extraction accuracy of MTO. We then apply MTO to the FDS images to identify LSB candidates. The identified objects are fitted with 2D Sérsic models using GALFIT and classified based on their morphological appearance, colors, and structure. With MTO, we are able to increase the completeness of our earlier FDS dwarf catalog (FDSDC) 0.5-1 mag deeper in terms of total magnitude and surface brightness. Due to the increased accuracy in measuring sizes of the detected objects, we also add many small galaxies to the catalog that were previously excluded as their outer parts had been missed in detection. We detect 265 new LSB dwarf galaxies in the Fornax cluster, which increases the total number of known dwarfs in Fornax to 821. Using the extended catalog, we show that the luminosity function has a faint-end slope of -1.38+/-0.02. We compare the obtained luminosity function with different environments studied earlier using deep data but do not find any significant differences. On the other hand, the Fornax-like simulated clusters in the IllustrisTNG cosmological simulation have shallower slopes than found in the observational data. We also find several trends in the galaxy colors, structure, and morphology that support the idea that the number of LSB galaxies is higher in the cluster center due to tidal forces and the age dimming of the stellar populations. The same result also holds for the subgroup of large LSB galaxies, so-called ultra-diffuse galaxies.
△ Less
Submitted 2 November, 2021;
originally announced November 2021.
-
Application of pattern spectra and convolutional neural networks to the analysis of simulated Cherenkov Telescope Array data
Authors:
J. Aschersleben,
R. F. Peletier,
M. Vecchi,
M. H. F. Wilkinson
Abstract:
The Cherenkov Telescope Array (CTA) will be the next generation gamma-ray observatory and will be the major global instrument for very-high-energy astronomy over the next decade, offering 5 - 10 x better flux sensitivity than current generation gamma-ray telescopes. Each telescope will provide a snapshot of gamma-ray induced particle showers by capturing the induced Cherenkov emission at ground le…
▽ More
The Cherenkov Telescope Array (CTA) will be the next generation gamma-ray observatory and will be the major global instrument for very-high-energy astronomy over the next decade, offering 5 - 10 x better flux sensitivity than current generation gamma-ray telescopes. Each telescope will provide a snapshot of gamma-ray induced particle showers by capturing the induced Cherenkov emission at ground level. The simulation of such events provides images that can be used as training data for convolutional neural networks (CNNs) to determine the energy of the initial gamma rays. Compared to other state-of-the-art algorithms, analyses based on CNNs promise to further enhance the performance to be achieved by CTA. Pattern spectra are commonly used tools for image classification and provide the distributions of the shapes and sizes of various objects comprising an image. The use of relatively shallow CNNs on pattern spectra would automatically select relevant combinations of features within an image, taking advantage of the 2D nature of pattern spectra. In this work, we generate pattern spectra from simulated gamma-ray events instead of using the raw images themselves in order to train our CNN for energy reconstruction. This is different from other relevant learning and feature selection methods that have been tried in the past. Thereby, we aim to obtain a significantly faster and less computationally intensive algorithm, with minimal loss of performance.
△ Less
Submitted 2 August, 2021;
originally announced August 2021.
-
Inferring the properties of the sources of reionization using the morphological spectra of the ionized regions
Authors:
Simon Gazagnes,
Léon V. E. Koopmans,
Michael H. F. Wilkinson
Abstract:
High-redshift 21-cm observations will provide crucial insights into the physical processes of the Epoch of Reionization. Next-generation interferometers such as the Square Kilometer Array will have enough sensitivity to directly image the 21-cm fluctuations and trace the evolution of the ionizing fronts. In this work, we develop an inferential approach to recover the sources and IGM properties of…
▽ More
High-redshift 21-cm observations will provide crucial insights into the physical processes of the Epoch of Reionization. Next-generation interferometers such as the Square Kilometer Array will have enough sensitivity to directly image the 21-cm fluctuations and trace the evolution of the ionizing fronts. In this work, we develop an inferential approach to recover the sources and IGM properties of the process of reionization using the number and, in particular, the morphological pattern spectra of the ionized regions extracted from realistic mock observations. To do so, we extend the Markov Chain Monte Carlo analysis tool 21CMMC by including these 21-cm tomographic statistics and compare this method to only using the power spectrum. We demonstrate that the evolution of the number-count and morphology of the ionized regions as a function of redshift provides independent information to disentangle multiple reionization scenarios because it probes the average ionizing budget per baryon. Although less precise, we find that constraints inferred using 21-cm tomographic statistics are more robust to the presence of contaminants such as foreground residuals. This work highlights that combining power spectrum and tomographic analyses more accurately recovers the astrophysics of reionization.
△ Less
Submitted 12 January, 2021; v1 submitted 16 November, 2020;
originally announced November 2020.
-
Optimising and comparing source extraction tools using objective segmentation quality criteria
Authors:
Caroline Haigh,
Nushkia Chamba,
Aku Venhola,
Reynier Peletier,
Lars Doorenbos,
Matthew Watkins,
Michael H. F. Wilkinson
Abstract:
With the growth of the scale, depth, and resolution of astronomical imaging surveys, there is an increased need for highly accurate automated detection and extraction of astronomical sources from images. This also means there is a need for objective quality criteria, and automated methods to optimise parameter settings for these software tools.
We present a comparison of several tools which have…
▽ More
With the growth of the scale, depth, and resolution of astronomical imaging surveys, there is an increased need for highly accurate automated detection and extraction of astronomical sources from images. This also means there is a need for objective quality criteria, and automated methods to optimise parameter settings for these software tools.
We present a comparison of several tools which have been developed to perform this task: namely SExtractor, ProFound, NoiseChisel, and MTObjects. In particular, we focus on evaluating performance in situations which present challenges for detection -- for example, faint and diffuse galaxies; extended structures, such as streams; and objects close to bright sources. Furthermore, we develop an automated method to optimise the parameters for the above tools.
We present four different objective segmentation quality measures, based on precision, recall, and a new measure for the correctly identified area of sources. Bayesian optimisation is used to find optimal parameter settings for each of the four tools on simulated data, for which a ground truth is known. After training, the tools are tested on similar simulated data, to provide a performance baseline. We then qualitatively assess tool performance on real astronomical images from two different surveys.
We determine that when area is disregarded, all four tools are capable of broadly similar levels of detection completeness, while only NoiseChisel and MTObjects are capable of locating the faint outskirts of objects. MTObjects produces the highest scores on all tests on all four quality measures, whilst SExtractor obtains the highest speeds. No tool has sufficient speed and accuracy to be well-suited to large-scale automated segmentation in its current form.
△ Less
Submitted 16 November, 2020; v1 submitted 16 September, 2020;
originally announced September 2020.