Image subtraction in astronomy is a tool for transient object discovery and characterization, particularly useful in wide fields, and is well suited for moving or photometrically varying objects such as asteroids, extra-solar planets and supernovae. A convolution technique is used to match point spread functions (PSFs) between images of the same field taken at different times prior to pixel-by-pixel subtraction. Particularly suitable for large-scale images is a spatially-varying kernel, where the convolution is allowed to adapt to PSF changes as a function of position within the images. The most versatile basis for fitting the spatially-varying kernel is the Dirac delta function. However, the convolution kernel based on the delta function does not discriminate between pixel scale noise variations and the intended stellar point spread function signals. The situation can frequently lead to reduced signal to noise ratios for variable objects detectable in the resulting subtraction. This work presents a cross-correlation method for reducing noise effects on the delta function derived convolution kernels, thus yielding significantly improved signal to noise in the resulting subtraction.

Image subtraction in astronomy is a tool for transient object discovery such as asteroids, extra-solar planets and supernovae. To match point spread functions (PSFs) between images of the same field taken at different times a convolution technique is used. Particularly suitable for large-scale images is a computationally intensive spatially-varying kernel. The underlying algorithm is inherently massively parallel due to unique kernel generation at every pixel location. The spatially-varying kernel cannot be efficiently computed through the Convolution Theorem, and thus does not lend itself to acceleration by Fast Fourier Transform (FFT). This work presents results of accelerated implementation of the spatially-varying kernel image convolution in multi-cores with OpenMP and graphic processing units (GPUs). Typical speedups over ANSI-C were a factor of 50 and a factor of 1000 over the initial IDL implementation, demonstrating that the techniques are a practical and high impact path to terabyte-per-night image pipelines and petascale processing.