57 makeDirect=True, makePsfMatched=False):
58 """Warp and optionally PSF-match exposure
59
60 Parameters
61 ----------
62 exposure : :cpp:class: `lsst::afw::image::Exposure`
63 Exposure to preprocess.
64 wcs : :cpp:class:`lsst::afw::image::Wcs`
65 Desired WCS of temporary images.
66 modelPsf : :cpp:class: `lsst::meas::algorithms::KernelPsf` or None
67 Target PSF to which to match.
68 maxBBox : :cpp:class:`lsst::afw::geom::Box2I` or None
69 Maximum allowed parent bbox of warped exposure.
70 If None then the warped exposure will be just big enough to contain all warped pixels;
71 if provided then the warped exposure may be smaller, and so missing some warped pixels;
72 ignored if destBBox is not None.
73 destBBox: :cpp:class: `lsst::afw::geom::Box2I` or None
74 Exact parent bbox of warped exposure.
75 If None then maxBBox is used to determine the bbox, otherwise maxBBox is ignored.
76 makeDirect : bool
77 Return an exposure that has been only warped?
78 makePsfMatched : bool
79 Return an exposure that has been warped and PSF-matched?
80
81 Returns
82 -------
83 An lsst.pipe.base.Struct with the following fields:
84
85 direct : :cpp:class:`lsst::afw::image::Exposure`
86 warped exposure
87 psfMatched : :cpp:class: `lsst::afw::image::Exposure`
88 warped and psf-Matched temporary exposure
89 """
90 if makePsfMatched and modelPsf is None:
91 raise RuntimeError("makePsfMatched=True, but no model PSF was provided")
92
93 if not makePsfMatched and not makeDirect:
94 self.log.warning("Neither makeDirect nor makePsfMatched requested")
95
96
98 psfWarped = WarpedPsf(exposure.getPsf(), xyTransform)
99
100 if makePsfMatched and maxBBox is not None:
101
102 pixToGrow = 2 *
max(self.psfMatch.kConfig.sizeCellX,
103 self.psfMatch.kConfig.sizeCellY)
104
106 maxBBox.grow(pixToGrow)
107
108 with self.timer("warp"):
109 exposure = self.warper.warpExposure(wcs, exposure, maxBBox=maxBBox, destBBox=destBBox)
110 exposure.setPsf(psfWarped)
111
112 if makePsfMatched:
113 try:
114 exposurePsfMatched = self.psfMatch.run(exposure, modelPsf).psfMatchedExposure
115 except Exception as e:
116 exposurePsfMatched = None
117 self.log.info("Cannot PSF-Match: %s", e)
118
119 return pipeBase.Struct(
120 direct=exposure if makeDirect else None,
121 psfMatched=exposurePsfMatched if makePsfMatched else None
122 )
An integer coordinate rectangle.
std::shared_ptr< TransformPoint2ToPoint2 > makeWcsPairTransform(SkyWcs const &src, SkyWcs const &dst)
A Transform obtained by putting two SkyWcs objects "back to back".