Best way to warp a PSF to match a warped image?

[ameisner (Aaron Meisner)]

Hi @charlotteaward – thanks for your message. Another community member recently posted an example of using WarpedPsf :

Issues with Image Subtraction

Hi everyone, So I think I’ve figured it out. I wanted to share my solution in case anyone else wants to do this. The first step involves obtaining the sources in the visit image and making sure that the final SourceCatalog has the right columns. from lsst.meas.algorithms.detection import SourceDetectionTask from lsst.meas.deblender import SourceDeblendTask from lsst.meas.base import SingleFrameMeasurementTask import lsst.daf.base as dafBase calexp.mask.removeAndClearMaskPlane('DETECTED') sch…

Specifically the lines:

xyTransform = afwGeom.makeWcsPairTransform(template.wcs,calexp.wcs)
warpedPsf = measAlgorithms.WarpedPsf(template.psf, xyTransform)

Does that approach (passing an output from afwGeom.makeWcsPairTransform as an input to WarpedPsf) work for you? Thanks again.