lsst.pipe.tasks.makeDiscreteSkyMap.MakeDiscreteSkyMapTask Class Reference

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

Public Member Functions
	run (self, wcs_bbox_tuple_list, oldSkyMap=None)

Static Public Attributes
	ConfigClass = MakeDiscreteSkyMapConfig

Static Protected Attributes
str	_DefaultName = "makeDiscreteSkyMap"

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 66 of file makeDiscreteSkyMap.py.

Member Function Documentation

run()

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`, optional
   The SkyMap to extend if appending.

Returns
-------
skyMap : `lsst.pipe.base.Struct`
    Sky map returned as a struct with attributes:

    ``skyMap``
        The returned SkyMap (`lsst.skyMap.SkyMap`).

Definition at line 77 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`, optional
           The SkyMap to extend if appending.
 
        Returns
        -------
        skyMap : `lsst.pipe.base.Struct`
            Sky map returned as a struct with attributes:
 
            ``skyMap``
                The returned SkyMap (`lsst.skyMap.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
        )

Member Data Documentation

_DefaultName

str lsst.pipe.tasks.makeDiscreteSkyMap.MakeDiscreteSkyMapTask._DefaultName = "makeDiscreteSkyMap"

staticprotected

Definition at line 74 of file makeDiscreteSkyMap.py.

ConfigClass

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

static

Definition at line 73 of file makeDiscreteSkyMap.py.

---

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

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-8.0.0/Linux64/pipe_tasks/g29320951ab+012e1474a1/python/lsst/pipe/tasks/makeDiscreteSkyMap.py