LSSTApplications  17.0+124,17.0+14,17.0+73,18.0.0+37,18.0.0+80,18.0.0-4-g68ffd23+4,18.1.0-1-g0001055+12,18.1.0-1-g03d53ef+5,18.1.0-1-g1349e88+55,18.1.0-1-g2505f39+44,18.1.0-1-g5315e5e+4,18.1.0-1-g5e4b7ea+14,18.1.0-1-g7e8fceb+4,18.1.0-1-g85f8cd4+48,18.1.0-1-g8ff0b9f+4,18.1.0-1-ga2c679d+1,18.1.0-1-gd55f500+35,18.1.0-10-gb58edde+2,18.1.0-11-g0997b02+4,18.1.0-13-gfe4edf0b+12,18.1.0-14-g259bd21+21,18.1.0-19-gdb69f3f+2,18.1.0-2-g5f9922c+24,18.1.0-2-gd3b74e5+11,18.1.0-2-gfbf3545+32,18.1.0-26-g728bddb4+5,18.1.0-27-g6ff7ca9+2,18.1.0-3-g52aa583+25,18.1.0-3-g8ea57af+9,18.1.0-3-gb69f684+42,18.1.0-3-gfcaddf3+6,18.1.0-32-gd8786685a,18.1.0-4-gf3f9b77+6,18.1.0-5-g1dd662b+2,18.1.0-5-g6dbcb01+41,18.1.0-6-gae77429+3,18.1.0-7-g9d75d83+9,18.1.0-7-gae09a6d+30,18.1.0-9-gc381ef5+4,w.2019.45
LSSTDataManagementBasePackage
Classes | Functions | Variables
lsst.meas.base.forcedPhotCcd Namespace Reference

Classes

class  ForcedPhotCcdConnections
 
class  PerTractCcdDataIdContainer
 

Functions

def imageOverlapsTract (tract, imageWcs, imageBox)
 

Variables

 applyMosaicResults = None
 

Function Documentation

◆ imageOverlapsTract()

def lsst.meas.base.forcedPhotCcd.imageOverlapsTract (   tract,
  imageWcs,
  imageBox 
) 
Return whether the given bounding box overlaps the tract given a WCS.

Parameters
----------
tract : `lsst.skymap.TractInfo`
    TractInfo specifying a tract.
imageWcs : `lsst.afw.geom.SkyWcs`
    World coordinate system for the image.
imageBox : `lsst.geom.Box2I`
    Bounding box for the image.

Returns
-------
overlap : `bool`
    `True` if the bounding box overlaps the tract; `False` otherwise.

Definition at line 110 of file forcedPhotCcd.py.

110 def imageOverlapsTract(tract, imageWcs, imageBox):
111  """Return whether the given bounding box overlaps the tract given a WCS.
112 
113  Parameters
114  ----------
115  tract : `lsst.skymap.TractInfo`
116  TractInfo specifying a tract.
117  imageWcs : `lsst.afw.geom.SkyWcs`
118  World coordinate system for the image.
119  imageBox : `lsst.geom.Box2I`
120  Bounding box for the image.
121 
122  Returns
123  -------
124  overlap : `bool`
125  `True` if the bounding box overlaps the tract; `False` otherwise.
126  """
127  tractPoly = tract.getOuterSkyPolygon()
128 
129  imagePixelCorners = lsst.geom.Box2D(imageBox).getCorners()
130  try:
131  imageSkyCorners = imageWcs.pixelToSky(imagePixelCorners)
132  except lsst.pex.exceptions.LsstCppException as e:
133  # Protecting ourselves from awful Wcs solutions in input images
134  if (not isinstance(e.message, lsst.pex.exceptions.DomainErrorException) and
135  not isinstance(e.message, lsst.pex.exceptions.RuntimeErrorException)):
136  raise
137  return False
138 
139  imagePoly = lsst.sphgeom.ConvexPolygon.convexHull([coord.getVector() for coord in imageSkyCorners])
140  return tractPoly.intersects(imagePoly) # "intersects" also covers "contains" or "is contained by"
141 
142 
lsst::geom::Box2D
A floating-point coordinate rectangle geometry.
Definition: Box.h:413
lsst::sphgeom::ConvexPolygon::convexHull
static ConvexPolygon convexHull(std::vector< UnitVector3d > const &points)
convexHull returns the convex hull of the given set of points if it exists and throws an exception ot...
Definition: ConvexPolygon.h:65
lsst::meas::base.forcedPhotCcd.imageOverlapsTract
def imageOverlapsTract(tract, imageWcs, imageBox)
Definition: forcedPhotCcd.py:110

Variable Documentation

◆ applyMosaicResults

lsst.meas.base.forcedPhotCcd.applyMosaicResults = None

Definition at line 41 of file forcedPhotCcd.py.

Generated on Sat Nov 16 2019 09:15:43 for LSSTApplications by   doxygen 1.8.13