LSSTApplications  17.0+50,17.0+84,17.0+9,18.0.0+14,18.0.0+2,18.0.0+30,18.0.0+4,18.0.0+9,18.0.0-2-ge43143a+4,18.1.0-1-g0001055,18.1.0-1-g0896a44+6,18.1.0-1-g1349e88+4,18.1.0-1-g2505f39+3,18.1.0-1-g380d4d4+4,18.1.0-1-g5e4b7ea,18.1.0-1-g85f8cd4+3,18.1.0-1-g9a6769a+2,18.1.0-1-ga1a4c1a+2,18.1.0-1-gc037db8,18.1.0-1-gd55f500+1,18.1.0-1-ge10677a+3,18.1.0-10-g73b8679e+7,18.1.0-11-g311e899+3,18.1.0-12-g0d73a3591,18.1.0-12-gc95f69a+3,18.1.0-2-g000ad9a+3,18.1.0-2-g31c43f9+3,18.1.0-2-g9c63283+4,18.1.0-2-gdf0b915+4,18.1.0-2-gf03bb23,18.1.0-3-g2e29e3d+6,18.1.0-3-g52aa583+2,18.1.0-3-g9cb968e+3,18.1.0-4-gd2e8982+6,18.1.0-5-g510c42a+3,18.1.0-5-gaeab27e+4,18.1.0-6-gdda7f3e+6,18.1.0-7-g89824ecc+4,w.2019.32
LSSTDataManagementBasePackage
Classes | Functions
lsst.meas.astrom.approximateWcs Namespace Reference

Classes

class  _MockTestCase
 

Functions

def approximateWcs (wcs, bbox, order=3, nx=20, ny=20, iterations=3, skyTolerance=0.001 *lsst.geom.arcseconds, pixelTolerance=0.02, useTanWcs=False)
 

Function Documentation

◆ approximateWcs()

def lsst.meas.astrom.approximateWcs.approximateWcs (   wcs,
  bbox,
  order = 3,
  nx = 20,
  ny = 20,
  iterations = 3,
  skyTolerance = 0.001*lsst.geom.arcseconds,
  pixelTolerance = 0.02,
  useTanWcs = False 
)
Approximate an existing WCS as a TAN-SIP WCS

The fit is performed by evaluating the WCS at a uniform grid of points
within a bounding box.

Parameters
----------
wcs : `lsst.afw.geom.SkyWcs`
    wcs to approximate
bbox : `lsst.geom.Box2I`
    the region over which the WCS will be fit
order : `int`
    order of SIP fit
nx : `int`
    number of grid points along x
ny : `int`
    number of grid points along y
iterations : `int`
    number of times to iterate over fitting
skyTolerance : `lsst.geom.Angle`
    maximum allowed difference in world coordinates between
    input wcs and approximate wcs (default is 0.001 arcsec)
pixelTolerance : `float`
    maximum allowed difference in pixel coordinates between
    input wcs and approximate wcs (default is 0.02 pixels)
useTanWcs : `bool`
    send a TAN version of wcs to the fitter? It is documented to require that,
    but I don't think the fitter actually cares

Returns
-------
fitWcs : `lsst.afw.geom.SkyWcs`
    the fit TAN-SIP WCS

Definition at line 44 of file approximateWcs.py.

44  skyTolerance=0.001*lsst.geom.arcseconds, pixelTolerance=0.02, useTanWcs=False):
45  """Approximate an existing WCS as a TAN-SIP WCS
46 
47  The fit is performed by evaluating the WCS at a uniform grid of points
48  within a bounding box.
49 
50  Parameters
51  ----------
52  wcs : `lsst.afw.geom.SkyWcs`
53  wcs to approximate
54  bbox : `lsst.geom.Box2I`
55  the region over which the WCS will be fit
56  order : `int`
57  order of SIP fit
58  nx : `int`
59  number of grid points along x
60  ny : `int`
61  number of grid points along y
62  iterations : `int`
63  number of times to iterate over fitting
64  skyTolerance : `lsst.geom.Angle`
65  maximum allowed difference in world coordinates between
66  input wcs and approximate wcs (default is 0.001 arcsec)
67  pixelTolerance : `float`
68  maximum allowed difference in pixel coordinates between
69  input wcs and approximate wcs (default is 0.02 pixels)
70  useTanWcs : `bool`
71  send a TAN version of wcs to the fitter? It is documented to require that,
72  but I don't think the fitter actually cares
73 
74  Returns
75  -------
76  fitWcs : `lsst.afw.geom.SkyWcs`
77  the fit TAN-SIP WCS
78  """
79  if useTanWcs:
80  crpix = wcs.getPixelOrigin()
81  crval = wcs.getSkyOrigin()
82  cdMatrix = wcs.getCdMatrix(crpix)
83  tanWcs = afwGeom.makeSkyWcs(crpix=crpix, crval=crval, cdMatrix=cdMatrix)
84  else:
85  tanWcs = wcs
86 
87  # create a matchList consisting of a grid of points covering the bbox
88  refSchema = afwTable.SimpleTable.makeMinimalSchema()
89  refCoordKey = afwTable.CoordKey(refSchema["coord"])
90  refCat = afwTable.SimpleCatalog(refSchema)
91 
92  sourceSchema = afwTable.SourceTable.makeMinimalSchema()
93  SingleFrameMeasurementTask(schema=sourceSchema) # expand the schema
94  sourceCentroidKey = afwTable.Point2DKey(sourceSchema["slot_Centroid"])
95 
96  sourceCat = afwTable.SourceCatalog(sourceSchema)
97 
98  matchList = []
99 
100  bboxd = lsst.geom.Box2D(bbox)
101  for x in np.linspace(bboxd.getMinX(), bboxd.getMaxX(), nx):
102  for y in np.linspace(bboxd.getMinY(), bboxd.getMaxY(), ny):
103  pixelPos = lsst.geom.Point2D(x, y)
104  skyCoord = wcs.pixelToSky(pixelPos)
105 
106  refObj = refCat.addNew()
107  refObj.set(refCoordKey, skyCoord)
108 
109  source = sourceCat.addNew()
110  source.set(sourceCentroidKey, pixelPos)
111 
112  matchList.append(afwTable.ReferenceMatch(refObj, source, 0.0))
113 
114  # The TAN-SIP fitter is fitting x and y separately, so we have to iterate to make it converge
115  for indx in range(iterations):
116  sipObject = makeCreateWcsWithSip(matchList, tanWcs, order, bbox)
117  tanWcs = sipObject.getNewWcs()
118  fitWcs = sipObject.getNewWcs()
119 
120  mockTest = _MockTestCase()
121  assertWcsAlmostEqualOverBBox(mockTest, wcs, fitWcs, bbox, maxDiffSky=skyTolerance,
122  maxDiffPix=pixelTolerance)
123 
124  return fitWcs
125 
A floating-point coordinate rectangle geometry.
Definition: Box.h:305
def assertWcsAlmostEqualOverBBox(testCase, wcs0, wcs1, bbox, maxDiffSky=0.01 *lsst.geom.arcseconds, maxDiffPix=0.01, nx=5, ny=5, msg="WCSs differ")
Definition: utils.py:155
Custom catalog class for record/table subclasses that are guaranteed to have an ID, and should generally be sorted by that ID.
Definition: fwd.h:63
Lightweight representation of a geometric match between two records.
Definition: fwd.h:101
std::shared_ptr< SkyWcs > makeSkyWcs(TransformPoint2ToPoint2 const &pixelsToFieldAngle, lsst::geom::Angle const &orientation, bool flipX, lsst::geom::SpherePoint const &boresight, std::string const &projection="TAN")
Construct a FITS SkyWcs from camera geometry.
Definition: SkyWcs.cc:516
CreateWcsWithSip< MatchT > makeCreateWcsWithSip(std::vector< MatchT > const &matches, afw::geom::SkyWcs const &linearWcs, int const order, geom::Box2I const &bbox=geom::Box2I(), int const ngrid=0)
Factory function for CreateWcsWithSip.
A FunctorKey used to get or set celestial coordinates from a pair of lsst::geom::Angle keys...
Definition: aggregates.h:210