A hybrid FourierReal Gaussian Mixture method for fast galaxyPSF convolution
Abstract
I present a method for the fast convolution of a model galaxy profile by a pointspread function (PSF) model represented as a pixel grid. The method relies upon three observations: First, most simple radial galaxy profiles of common interest (deVaucouleurs, exponential, Sersic) can be approximated as mixtures of Gaussians. Second, the Fourier transform of a Gaussian is a Gaussian, thus the Fourier transform of a mixtureofGausssian approximation of a galaxy can be directly evaluated as a mixture of Gaussians in Fourier space. Third, if a mixture component would result in Fourierspace aliasing, that component can be evaluated in real space. For mixture components to be evaluated in Fourier space, we use the FFT for the PSF model, direct evaluation of the Fourier transform for the galaxy, and the inverseFFT to return the result to pixel space. For mixture components to be evaluated in real spacewhich only happens when the mixture components is much larger than the PSFwe use a simple Gaussian approximation of the PSF, perform the convolution analytically, and evaluate in real pixel space. The method is fast and exact (to the extent that the mixtureofGaussians approximation of the galaxy profile is exact) as long as the pixelized PSF model is well sampled. This Fourier method can be seen as a way of applying a perfect lowpass filter to the (typically strongly undersampled) galaxy profile before convolution by the PSF, at exactly the Nyquist frequency of the PSF pixel model grid. In this way, it avoids the computational expense of a traditional superresolution approach. This method allows the efficient use of pixelized PSF models (ie, a PSF represented as a grid of pixel values) in galaxy forward modelfitting approaches such as the Tractor.
 Publication:

arXiv eprints
 Pub Date:
 December 2020
 arXiv:
 arXiv:2012.15797
 Bibcode:
 2020arXiv201215797L
 Keywords:

 Astrophysics  Instrumentation and Methods for Astrophysics