How to COAAD Images. I. Optimal Source Detection and Photometry of Point Sources Using Ensembles of Images

• Published in: The Astrophysical journal Vol. 836; no. 2; pp. 187 - 200
• Main Authors: Zackay, Barak, Ofek, Eran O.
• Format: Journal Article
• Language: English
• Published: Philadelphia The American Astronomical Society 20.02.2017; IOP Publishing
• Subjects: Algorithms; Apertures; Astrophysics; Background noise; Celestial bodies; Computer simulation; cosmology: observations; Digital imaging; Galaxies; gravitational lensing: weak; Ground-based observation; Image detection; Image filters; Matched filters; methods: data analysis; Photometry; Point sources; Point spread functions; Signal processing; Signal to noise ratio; Sky surveys (astronomy); Statistical tests; surveys; techniques: image processing; techniques: photometric
• ISSN: 0004-637X; 1538-4357
• DOI: 10.3847/1538-4357/836/2/187

Stacks of digital astronomical images are combined in order to increase image depth. The variable seeing conditions, sky background, and transparency of ground-based observations make the coaddition process nontrivial. We present image coaddition methods that maximize the signal-to-noise ratio (S/N) and optimized for source detection and flux measurement. We show that for these purposes, the best way to combine images is to apply a matched filter to each image using its own point-spread function (PSF) and only then to sum the images with the appropriate weights. Methods that either match the filter after coaddition or perform PSF homogenization prior to coaddition will result in loss of sensitivity. We argue that our method provides an increase of between a few and 25% in the survey speed of deep ground-based imaging surveys compared with weighted coaddition techniques. We demonstrate this claim using simulated data as well as data from the Palomar Transient Factory data release 2. We present a variant of this coaddition method, which is optimal for PSF or aperture photometry. We also provide an analytic formula for calculating the S/N for PSF photometry on single or multiple observations. In the next paper in this series, we present a method for image coaddition in the limit of background-dominated noise, which is optimal for any statistical test or measurement on the constant-in-time image (e.g., source detection, shape or flux measurement, or star-galaxy separation), making the original data redundant. We provide an implementation of these algorithms in MATLAB.