LSST Applications  21.0.0+04719a4bac,21.0.0-1-ga51b5d4+f5e6047307,21.0.0-11-g2b59f77+a9c1acf22d,21.0.0-11-ga42c5b2+86977b0b17,21.0.0-12-gf4ce030+76814010d2,21.0.0-13-g1721dae+760e7a6536,21.0.0-13-g3a573fe+768d78a30a,21.0.0-15-g5a7caf0+f21cbc5713,21.0.0-16-g0fb55c1+b60e2d390c,21.0.0-19-g4cded4ca+71a93a33c0,21.0.0-2-g103fe59+bb20972958,21.0.0-2-g45278ab+04719a4bac,21.0.0-2-g5242d73+3ad5d60fb1,21.0.0-2-g7f82c8f+8babb168e8,21.0.0-2-g8f08a60+06509c8b61,21.0.0-2-g8faa9b5+616205b9df,21.0.0-2-ga326454+8babb168e8,21.0.0-2-gde069b7+5e4aea9c2f,21.0.0-2-gecfae73+1d3a86e577,21.0.0-2-gfc62afb+3ad5d60fb1,21.0.0-25-g1d57be3cd+e73869a214,21.0.0-3-g357aad2+ed88757d29,21.0.0-3-g4a4ce7f+3ad5d60fb1,21.0.0-3-g4be5c26+3ad5d60fb1,21.0.0-3-g65f322c+e0b24896a3,21.0.0-3-g7d9da8d+616205b9df,21.0.0-3-ge02ed75+a9c1acf22d,21.0.0-4-g591bb35+a9c1acf22d,21.0.0-4-g65b4814+b60e2d390c,21.0.0-4-gccdca77+0de219a2bc,21.0.0-4-ge8a399c+6c55c39e83,21.0.0-5-gd00fb1e+05fce91b99,21.0.0-6-gc675373+3ad5d60fb1,21.0.0-64-g1122c245+4fb2b8f86e,21.0.0-7-g04766d7+cd19d05db2,21.0.0-7-gdf92d54+04719a4bac,21.0.0-8-g5674e7b+d1bd76f71f,master-gac4afde19b+a9c1acf22d,w.2021.13
LSST Data Management Base Package
warpAndPsfMatch.py
Go to the documentation of this file.
1 #
2 # LSST Data Management System
3 # Copyright 2008, 2009, 2010, 2011, 2012 LSST Corporation.
4 #
5 # This product includes software developed by the
6 # LSST Project (http://www.lsst.org/).
7 #
8 # This program is free software: you can redistribute it and/or modify
9 # it under the terms of the GNU General Public License as published by
10 # the Free Software Foundation, either version 3 of the License, or
11 # (at your option) any later version.
12 #
13 # This program is distributed in the hope that it will be useful,
14 # but WITHOUT ANY WARRANTY; without even the implied warranty of
15 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 # GNU General Public License for more details.
17 #
18 # You should have received a copy of the LSST License Statement and
19 # the GNU General Public License along with this program. If not,
20 # see <http://www.lsstcorp.org/LegalNotices/>.
21 #
22 import lsst.pex.config as pexConfig
23 import lsst.afw.math as afwMath
24 import lsst.geom as geom
25 import lsst.afw.geom as afwGeom
26 import lsst.pipe.base as pipeBase
27 from lsst.ip.diffim import ModelPsfMatchTask
28 from lsst.meas.algorithms import WarpedPsf
29 
30 __all__ = ["WarpAndPsfMatchTask"]
31 
32 
33 class WarpAndPsfMatchConfig(pexConfig.Config):
34  """Config for WarpAndPsfMatchTask
35  """
36  psfMatch = pexConfig.ConfigurableField(
37  target=ModelPsfMatchTask,
38  doc="PSF matching model to model task",
39  )
40  warp = pexConfig.ConfigField(
41  dtype=afwMath.Warper.ConfigClass,
42  doc="warper configuration",
43  )
44 
45 
46 class WarpAndPsfMatchTask(pipeBase.Task):
47  """A task to warp and PSF-match an exposure
48  """
49  ConfigClass = WarpAndPsfMatchConfig
50 
51  def __init__(self, *args, **kwargs):
52  pipeBase.Task.__init__(self, *args, **kwargs)
53  self.makeSubtask("psfMatch")
54  self.warperwarper = afwMath.Warper.fromConfig(self.config.warp)
55 
56  def run(self, exposure, wcs, modelPsf=None, maxBBox=None, destBBox=None,
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.warn("Neither makeDirect nor makePsfMatched requested")
95 
96  # Warp PSF before overwriting exposure
97  xyTransform = afwGeom.makeWcsPairTransform(exposure.getWcs(), wcs)
98  psfWarped = WarpedPsf(exposure.getPsf(), xyTransform)
99 
100  if makePsfMatched and maxBBox is not None:
101  # grow warped region to provide sufficient area for PSF-matching
102  pixToGrow = 2 * max(self.psfMatch.kConfig.sizeCellX,
103  self.psfMatch.kConfig.sizeCellY)
104  # replace with copy
105  maxBBox = geom.Box2I(maxBBox)
106  maxBBox.grow(pixToGrow)
107 
108  with self.timer("warp"):
109  exposure = self.warperwarper.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  )
int max
An integer coordinate rectangle.
Definition: Box.h:55
A Psf class that maps an arbitrary Psf through a coordinate transformation.
Definition: WarpedPsf.h:50
def run(self, exposure, wcs, modelPsf=None, maxBBox=None, destBBox=None, makeDirect=True, makePsfMatched=False)
std::shared_ptr< TransformPoint2ToPoint2 > makeWcsPairTransform(SkyWcs const &src, SkyWcs const &dst)
A Transform obtained by putting two SkyWcs objects "back to back".
Definition: SkyWcs.cc:151
int warpExposure(DestExposureT &destExposure, SrcExposureT const &srcExposure, WarpingControl const &control, typename DestExposureT::MaskedImageT::SinglePixel padValue=lsst::afw::math::edgePixel< typename DestExposureT::MaskedImageT >(typename lsst::afw::image::detail::image_traits< typename DestExposureT::MaskedImageT >::image_category()))
Warp (remap) one exposure to another.
Fit spatial kernel using approximate fluxes for candidates, and solving a linear system of equations.