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
lsst.obs.base.utils Namespace Reference

Classes

class  InitialSkyWcsError
 

Functions

def createInitialSkyWcs (visitInfo, detector, flipX=False)
 
def bboxFromIraf (irafBBoxStr)
 

Function Documentation

◆ bboxFromIraf()

def lsst.obs.base.utils.bboxFromIraf (   irafBBoxStr)
Return a Box2I corresponding to an IRAF-style BBOX

[x0:x1,y0:y1] where x0 and x1 are the one-indexed start and end columns, and correspondingly
y0 and y1 are the start and end rows.

Definition at line 84 of file utils.py.

84 def bboxFromIraf(irafBBoxStr):
85  """Return a Box2I corresponding to an IRAF-style BBOX
86 
87  [x0:x1,y0:y1] where x0 and x1 are the one-indexed start and end columns, and correspondingly
88  y0 and y1 are the start and end rows.
89  """
90 
91  mat = re.search(r"^\[([-\d]+):([-\d]+),([-\d]+):([-\d]+)\]$", irafBBoxStr)
92  if not mat:
93  raise RuntimeError("Unable to parse IRAF-style bbox \"%s\"" % irafBBoxStr)
94  x0, x1, y0, y1 = [int(_) for _ in mat.groups()]
95 
96  return geom.BoxI(geom.PointI(x0 - 1, y0 - 1), geom.PointI(x1 - 1, y1 - 1))
97 
def bboxFromIraf(irafBBoxStr)
Definition: utils.py:84
An integer coordinate rectangle.
Definition: Box.h:55

◆ createInitialSkyWcs()

def lsst.obs.base.utils.createInitialSkyWcs (   visitInfo,
  detector,
  flipX = False 
)
Create a SkyWcs from the telescope boresight and detector geometry.

A typical usecase for this is to create the initial WCS for a newly-read
raw exposure.


Parameters
----------
visitInfo : `lsst.afw.image.VisitInfo`
    Where to get the telescope boresight and rotator angle from.
detector : `lsst.afw.cameraGeom.Detector`
    Where to get the camera geomtry from.
flipX : `bool`, optional
    If False, +X is along W, if True +X is along E.

Returns
-------
skyWcs : `lsst.afw.geom.SkyWcs`
    The new composed WCS.

Raises
------
InitialSkyWcsError
    Raised if there is an error generating the SkyWcs, chained from the
    lower-level exception if available.

Definition at line 43 of file utils.py.

43 def createInitialSkyWcs(visitInfo, detector, flipX=False):
44  """Create a SkyWcs from the telescope boresight and detector geometry.
45 
46  A typical usecase for this is to create the initial WCS for a newly-read
47  raw exposure.
48 
49 
50  Parameters
51  ----------
52  visitInfo : `lsst.afw.image.VisitInfo`
53  Where to get the telescope boresight and rotator angle from.
54  detector : `lsst.afw.cameraGeom.Detector`
55  Where to get the camera geomtry from.
56  flipX : `bool`, optional
57  If False, +X is along W, if True +X is along E.
58 
59  Returns
60  -------
61  skyWcs : `lsst.afw.geom.SkyWcs`
62  The new composed WCS.
63 
64  Raises
65  ------
66  InitialSkyWcsError
67  Raised if there is an error generating the SkyWcs, chained from the
68  lower-level exception if available.
69  """
70  if visitInfo.getRotType() != RotType.SKY:
71  msg = (f"Cannot create SkyWcs from camera geometry: rotator angle defined using "
72  f"RotType={visitInfo.getRotType()} instead of SKY.")
73  raise InitialSkyWcsError(msg)
74  orientation = visitInfo.getBoresightRotAngle()
75  boresight = visitInfo.getBoresightRaDec()
76  try:
77  pixelsToFieldAngle = detector.getTransform(detector.makeCameraSys(PIXELS),
78  detector.makeCameraSys(FIELD_ANGLE))
80  raise InitialSkyWcsError("Cannot compute PIXELS to FIELD_ANGLE Transform.") from e
81  return makeSkyWcs(pixelsToFieldAngle, orientation, flipX, boresight)
82 
83 
def createInitialSkyWcs(visitInfo, detector, flipX=False)
Definition: utils.py:43
std::shared_ptr< SkyWcs > makeSkyWcs(daf::base::PropertySet &metadata, bool strip=false)
Construct a SkyWcs from FITS keywords.
Definition: SkyWcs.cc:506
Reports invalid arguments.
Definition: Runtime.h:66