When extracting the psf (
lsst.afw.detection.detectionLib.Psf) as an image (
Psf.computeImage().getArray()) in each band for a set of calexps (
lsst.afw.image.imageLib.ExposureF), while the calexps are spatially aligned, the psf images are not (and in fact they are often different sizes). This is an issue when trying to use the new deblender that operates on de-convolved templates, resulting in color gradients that do not exist.
The dirty way to fix this is to calculate the maximum pixel in each psf image and project them onto a common image, trimming the pixels in each psf outside of the common psf image. But I don’t know enough about how the psfs are generated in the first place, and what invalid assumptions I might be making while implementing this procedure. So I was wondering if there is a more mathematically rigorous solution (that ideally is implemented in the code already!). Does anyone know of such a thing?