Optimal matched filter in the low-number count Poisson noise regime and implications for X-ray source detection

Abstract

Detection of templates (e.g., sources) embedded in low-number count Poisson noise is a common problem in astrophysics. Examples include source detection in X-ray images, gamma-rays, UV, neutrinos, and search for clusters of galaxies and stellar streams. However, the solutions in the X-ray-related literature are sub-optimal -- in some cases by considerable factors. Using the lemma of Neyman-Pearson we derive the optimal statistics for template detection in the presence of Poisson noise. We demonstrate that this method provides higher completeness, for a fixed false-alarm probability value, compared with filtering the image with the point-spread function (PSF). In turn, we find that filtering by the PSF is better than filtering the image using the Mexican-hat wavelet (used by wavedetect). For some background levels, our method improves the sensitivity of source detection by more than a factor of two over the popular Mexican-hat wavelet filtering. This filtering technique can also be used also for fast PSF photometry and flare detection, and it is efficient, as well as straight forward to implement. We provide an implementation in MATLAB.

Submitted to ApJ

Cited by in corpus (10)

• Optimal and Efficient Streak Detection in Astronomical Images
• A semi-supervised Machine Learning search for never-seen Gravitational-Wave sources
• Sifting Through the Static: Moving Object Detection in Difference Images
• AT2018lqh and the nature of the emerging population of day-scale duration optical transients
• Introducing constrained matched filters for improved separation of point sources from galaxy clusters
• PCAT-DE: Reconstructing point-like and diffuse signals in astronomical images using spatial and spectral information
• Matched filter in the low-number count Poisson noise regime: an efficient and effective implementation
• SRG/ART-XC all-sky X-ray survey: Catalog of sources detected during the first year
• Signal-to-noise-ratio and SNR-max detection statistics in template bank searches for exotic physics transients with networks of quantum sensors
• A New Search Pipeline for Short Gamma Ray Bursts in Fermi/GBM Data -- A 50% Increase in the Number of Detections