The current warp code works nicely, but it assumes that the interface is two Exposures (or (Image, Wcs) tuple) and warps one to the other.
Yasuda-san here at IPMU is looking at dipoles in his SNe difference imaging which uses the LSST ip_diffim code. The astrometry isn’t quite good enough, so he’d like to
- Solve for the relative astrometry between two coadds
- Warp one onto the other
Is there an easy way to do this? It could clearly be done by transforming the polynomial to a TAN-SIP Wcs, but that’s painful.