lsst.pipe.tasks.selectImages.WcsSelectImagesTask Class Reference

Inheritance diagram for lsst.pipe.tasks.selectImages.WcsSelectImagesTask:

Public Member Functions
	run (self, wcsList, bboxList, coordList, dataIds=None, **kwargs)
	getValidImageCorners (self, imageWcs, imageBox, patchPoly, dataId=None)

Detailed Description

Select images using their Wcs.

We use the "convexHull" method of lsst.sphgeom.ConvexPolygon to define
polygons on the celestial sphere, and test the polygon of the
patch for overlap with the polygon of the image.

We use "convexHull" instead of generating a ConvexPolygon
directly because the standard for the inputs to ConvexPolygon
are pretty high and we don't want to be responsible for reaching them.

Definition at line 162 of file selectImages.py.

Member Function Documentation

getValidImageCorners()

lsst.pipe.tasks.selectImages.WcsSelectImagesTask.getValidImageCorners	(		self,
			imageWcs,
			imageBox,
			patchPoly,
			dataId = None )

Return corners or `None` if bad.

Parameters
----------
imageWcs : `Unknown`
imageBox : `Unknown`
patchPoly : `Unknown`
dataId : `Unknown`

Definition at line 209 of file selectImages.py.

    def getValidImageCorners(self, imageWcs, imageBox, patchPoly, dataId=None):
        """Return corners or `None` if bad.
 
        Parameters
        ----------
        imageWcs : `Unknown`
        imageBox : `Unknown`
        patchPoly : `Unknown`
        dataId : `Unknown`
        """
        try:
            imageCorners = [imageWcs.pixelToSky(pix) for pix in geom.Box2D(imageBox).getCorners()]
        except (pexExceptions.DomainError, pexExceptions.RuntimeError) as e:
            # Protecting ourselves from awful Wcs solutions in input images
            self.log.debug("WCS error in testing calexp %s (%s): deselecting", dataId, e)
            return None
 
        imagePoly = lsst.sphgeom.ConvexPolygon.convexHull([coord.getVector() for coord in imageCorners])
        if imagePoly is None:
            self.log.debug("Unable to create polygon from image %s: deselecting", dataId)
            return None
 
        if patchPoly.intersects(imagePoly):
            # "intersects" also covers "contains" or "is contained by"
            self.log.info("Selecting calexp %s", dataId)
            return imageCorners
 
        return None
 
 

run()

lsst.pipe.tasks.selectImages.WcsSelectImagesTask.run	(		self,
			wcsList,
			bboxList,
			coordList,
			dataIds = None,
		**	kwargs )

Return indices of provided lists that meet the selection criteria.

Parameters
----------
wcsList : `list` [`lsst.afw.geom.SkyWcs`]
    Specifying the WCS's of the input ccds to be selected.
bboxList : `list` [`lsst.geom.Box2I`]
    Specifying the bounding boxes of the input ccds to be selected.
coordList : `list` [`lsst.geom.SpherePoint`]
    ICRS coordinates specifying boundary of the patch.
dataIds : iterable [`lsst.daf.butler.dataId`] or `None`, optional
    An iterable object of dataIds which point to reference catalogs.
**kwargs
    Additional keyword arguments.

Returns
-------
result : `list` [`int`]
    The indices of selected ccds.

Reimplemented from lsst.pipe.tasks.selectImages.BaseSelectImagesTask.

Definition at line 174 of file selectImages.py.

    def run(self, wcsList, bboxList, coordList, dataIds=None, **kwargs):
        """Return indices of provided lists that meet the selection criteria.
 
        Parameters
        ----------
        wcsList : `list` [`lsst.afw.geom.SkyWcs`]
            Specifying the WCS's of the input ccds to be selected.
        bboxList : `list` [`lsst.geom.Box2I`]
            Specifying the bounding boxes of the input ccds to be selected.
        coordList : `list` [`lsst.geom.SpherePoint`]
            ICRS coordinates specifying boundary of the patch.
        dataIds : iterable [`lsst.daf.butler.dataId`] or `None`, optional
            An iterable object of dataIds which point to reference catalogs.
        **kwargs
            Additional keyword arguments.
 
        Returns
        -------
        result : `list` [`int`]
            The indices of selected ccds.
        """
        if dataIds is None:
            dataIds = [None] * len(wcsList)
        patchVertices = [coord.getVector() for coord in coordList]
        patchPoly = lsst.sphgeom.ConvexPolygon.convexHull(patchVertices)
        result = []
        for i, (imageWcs, imageBox, dataId) in enumerate(zip(wcsList, bboxList, dataIds)):
            if imageWcs is None:
                self.log.info("De-selecting exposure %s:  Exposure has no WCS.", dataId)
            else:
                imageCorners = self.getValidImageCorners(imageWcs, imageBox, patchPoly, dataId)
                if imageCorners:
                    result.append(i)
        return result
 

---

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

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-8.0.0/Linux64/pipe_tasks/g2a9a014e59+6f9be1b9cd/python/lsst/pipe/tasks/selectImages.py