Make a DiscreteSkyMap in a repository, using the bounding box of a set of calexps. More...

Inheritance diagram for lsst.pipe.tasks.makeDiscreteSkyMap.MakeDiscreteSkyMapTask:

Public Member Functions
def	__init__ (self, **kwargs)

def	runDataRef (self, butler, dataRefList)

def	run (self, wcs_bbox_tuple_list, oldSkyMap=None)

Static Public Attributes
	ConfigClass = MakeDiscreteSkyMapConfig

	RunnerClass = MakeDiscreteSkyMapRunner

Detailed Description

Make a DiscreteSkyMap in a repository, using the bounding box of a set of calexps.

The command-line and run signatures and config are sufficiently different from MakeSkyMapTask
that we don't inherit from it, but it is a replacement, so we use the same config/metadata names.

Definition at line 118 of file makeDiscreteSkyMap.py.

Constructor & Destructor Documentation

◆ init()

def lsst.pipe.tasks.makeDiscreteSkyMap.MakeDiscreteSkyMapTask.__init__	(		self,
		**	kwargs
	)

Definition at line 128 of file makeDiscreteSkyMap.py.

     def __init__(self, **kwargs):
         super().__init__(**kwargs)
  

Member Function Documentation

◆ run()

def lsst.pipe.tasks.makeDiscreteSkyMap.MakeDiscreteSkyMapTask.run	(	self,
		wcs_bbox_tuple_list,
		oldSkyMap = `None`
	)

Make a SkyMap from the bounds of the given set of calexp metadata.

Parameters
----------
wcs_bbox_tuple_list : iterable
   A list of tuples with each element expected to be a (Wcs, Box2I) pair
oldSkyMap : `lsst.skymap.DiscreteSkyMap`, option
   The SkyMap to extend if appending
Returns
-------
struct : `lsst.pipe.base.Struct
   The returned struct has one attribute, ``skyMap``, which holds the returned SkyMap

Definition at line 167 of file makeDiscreteSkyMap.py.

     def run(self, wcs_bbox_tuple_list, oldSkyMap=None):
         """Make a SkyMap from the bounds of the given set of calexp metadata.
  
         Parameters
         ----------
         wcs_bbox_tuple_list : iterable
            A list of tuples with each element expected to be a (Wcs, Box2I) pair
         oldSkyMap : `lsst.skymap.DiscreteSkyMap`, option
            The SkyMap to extend if appending
         Returns
         -------
         struct : `lsst.pipe.base.Struct
            The returned struct has one attribute, ``skyMap``, which holds the returned SkyMap
         """
         self.log.info("Extracting bounding boxes of %d images", len(wcs_bbox_tuple_list))
         points = []
         for wcs, boxI in wcs_bbox_tuple_list:
             boxD = geom.Box2D(boxI)
             points.extend(wcs.pixelToSky(corner).getVector() for corner in boxD.getCorners())
         if len(points) == 0:
             raise RuntimeError("No data found from which to compute convex hull")
         self.log.info("Computing spherical convex hull")
         polygon = lsst.sphgeom.ConvexPolygon.convexHull(points)
         if polygon is None:
             raise RuntimeError(
                 "Failed to compute convex hull of the vertices of all calexp bounding boxes; "
                 "they may not be hemispherical."
             )
         circle = polygon.getBoundingCircle()
  
         skyMapConfig = DiscreteSkyMap.ConfigClass()
         if oldSkyMap:
             skyMapConfig.raList.extend(oldSkyMap.config.raList)
             skyMapConfig.decList.extend(oldSkyMap.config.decList)
             skyMapConfig.radiusList.extend(oldSkyMap.config.radiusList)
         configIntersection = {k: getattr(self.config.skyMap, k)
                               for k in self.config.skyMap.toDict()
                               if k in skyMapConfig}
         skyMapConfig.update(**configIntersection)
         circleCenter = lsst.sphgeom.LonLat(circle.getCenter())
         skyMapConfig.raList.append(circleCenter[0].asDegrees())
         skyMapConfig.decList.append(circleCenter[1].asDegrees())
         circleRadiusDeg = circle.getOpeningAngle().asDegrees()
         skyMapConfig.radiusList.append(circleRadiusDeg + self.config.borderSize)
         skyMap = DiscreteSkyMap(skyMapConfig)
  
         for tractInfo in skyMap:
             wcs = tractInfo.getWcs()
             posBox = geom.Box2D(tractInfo.getBBox())
             pixelPosList = (
                 posBox.getMin(),
                 geom.Point2D(posBox.getMaxX(), posBox.getMinY()),
                 posBox.getMax(),
                 geom.Point2D(posBox.getMinX(), posBox.getMaxY()),
             )
             skyPosList = [wcs.pixelToSky(pos).getPosition(geom.degrees) for pos in pixelPosList]
             posStrList = ["(%0.3f, %0.3f)" % tuple(skyPos) for skyPos in skyPosList]
             self.log.info("tract %s has corners %s (RA, Dec deg) and %s x %s patches",
                           tractInfo.getId(), ", ".join(posStrList),
                           tractInfo.getNumPatches()[0], tractInfo.getNumPatches()[1])
         return pipeBase.Struct(
             skyMap=skyMap
         )
  

◆ runDataRef()

def lsst.pipe.tasks.makeDiscreteSkyMap.MakeDiscreteSkyMapTask.runDataRef	(	self,
		butler,
		dataRefList
	)

Make a skymap from the bounds of the given set of calexps using the butler.

Parameters
----------
butler : `lsst.daf.persistence.Butler`
   Gen2 data butler used to save the SkyMap
dataRefList : iterable
   A list of Gen2 data refs of calexps used to determin the size and pointing of the SkyMap
Returns
-------
struct : `lsst.pipe.base.Struct`
   The returned struct has one attribute, ``skyMap``, which holds the returned SkyMap

Definition at line 131 of file makeDiscreteSkyMap.py.

     def runDataRef(self, butler, dataRefList):
         """Make a skymap from the bounds of the given set of calexps using the butler.
  
         Parameters
         ----------
         butler : `lsst.daf.persistence.Butler`
            Gen2 data butler used to save the SkyMap
         dataRefList : iterable
            A list of Gen2 data refs of calexps used to determin the size and pointing of the SkyMap
         Returns
         -------
         struct : `lsst.pipe.base.Struct`
            The returned struct has one attribute, ``skyMap``, which holds the returned SkyMap
         """
         wcs_bbox_tuple_list = []
         oldSkyMap = None
         datasetName = self.config.coaddName + "Coadd_skyMap"
         for dataRef in dataRefList:
             if not dataRef.datasetExists("calexp"):
                 self.log.warning("CalExp for %s does not exist: ignoring", dataRef.dataId)
                 continue
             wcs_bbox_tuple_list.append((dataRef.get("calexp_wcs", immediate=True),
                                         dataRef.get("calexp_bbox", immediate=True)))
         if self.config.doAppend and butler.datasetExists(datasetName):
             oldSkyMap = butler.get(datasetName, immediate=True)
             if not isinstance(oldSkyMap.config, DiscreteSkyMap.ConfigClass):
                 raise TypeError("Cannot append to existing non-discrete skymap")
             compareLog = []
             if not self.config.skyMap.compare(oldSkyMap.config, output=compareLog.append):
                 raise ValueError("Cannot append to existing skymap - configurations differ:", *compareLog)
         result = self.run(wcs_bbox_tuple_list, oldSkyMap)
         if self.config.doWrite:
             butler.put(result.skyMap, datasetName)
         return result
  

Member Data Documentation

◆ ConfigClass

lsst.pipe.tasks.makeDiscreteSkyMap.MakeDiscreteSkyMapTask.ConfigClass = MakeDiscreteSkyMapConfig

static

Definition at line 124 of file makeDiscreteSkyMap.py.

◆ RunnerClass

lsst.pipe.tasks.makeDiscreteSkyMap.MakeDiscreteSkyMapTask.RunnerClass = MakeDiscreteSkyMapRunner

static

Definition at line 126 of file makeDiscreteSkyMap.py.

The documentation for this class was generated from the following file:

/j/snowflake/release/lsstsw/stack/lsst-scipipe-0.7.0/Linux64/pipe_tasks/21.0.0-172-gfb10e10a+18fedfabac/python/lsst/pipe/tasks/makeDiscreteSkyMap.py

Public Member Functions

Static Public Attributes