LSST Applications  21.0.0+04719a4bac,21.0.0-1-ga51b5d4+f5e6047307,21.0.0-11-g2b59f77+a9c1acf22d,21.0.0-11-ga42c5b2+86977b0b17,21.0.0-12-gf4ce030+76814010d2,21.0.0-13-g1721dae+760e7a6536,21.0.0-13-g3a573fe+768d78a30a,21.0.0-15-g5a7caf0+f21cbc5713,21.0.0-16-g0fb55c1+b60e2d390c,21.0.0-19-g4cded4ca+71a93a33c0,21.0.0-2-g103fe59+bb20972958,21.0.0-2-g45278ab+04719a4bac,21.0.0-2-g5242d73+3ad5d60fb1,21.0.0-2-g7f82c8f+8babb168e8,21.0.0-2-g8f08a60+06509c8b61,21.0.0-2-g8faa9b5+616205b9df,21.0.0-2-ga326454+8babb168e8,21.0.0-2-gde069b7+5e4aea9c2f,21.0.0-2-gecfae73+1d3a86e577,21.0.0-2-gfc62afb+3ad5d60fb1,21.0.0-25-g1d57be3cd+e73869a214,21.0.0-3-g357aad2+ed88757d29,21.0.0-3-g4a4ce7f+3ad5d60fb1,21.0.0-3-g4be5c26+3ad5d60fb1,21.0.0-3-g65f322c+e0b24896a3,21.0.0-3-g7d9da8d+616205b9df,21.0.0-3-ge02ed75+a9c1acf22d,21.0.0-4-g591bb35+a9c1acf22d,21.0.0-4-g65b4814+b60e2d390c,21.0.0-4-gccdca77+0de219a2bc,21.0.0-4-ge8a399c+6c55c39e83,21.0.0-5-gd00fb1e+05fce91b99,21.0.0-6-gc675373+3ad5d60fb1,21.0.0-64-g1122c245+4fb2b8f86e,21.0.0-7-g04766d7+cd19d05db2,21.0.0-7-gdf92d54+04719a4bac,21.0.0-8-g5674e7b+d1bd76f71f,master-gac4afde19b+a9c1acf22d,w.2021.13
LSST Data Management Base Package
Functions
lsst.jointcal.testUtils Namespace Reference

Functions

def canRunTests ()
 
def createTwoFakeCcdImages (num1=4, num2=4, seed=100, fakeCcdId=12, photoCalibMean1=1e-2, photoCalibMean2=1.2e-2, fakeWcses=(None, None), fakeVisitInfos=(None, None))
 
def createFakeCcdImage (butler, visit, num, fluxFieldName, photoCalibMean=1e-2, photoCalibErr=1.0, fakeCcdId=12, fakeWcs=None, fakeVisitInfo=None)
 
def createFakeCatalog (num, bbox, fluxFieldName, skyWcs=None, refCat=False)
 
def fillCatalog (schema, num, bbox, centroidKey, xErrKey, yErrKey, shapeKey, fluxFieldName, skyWcs=None, fluxErrFraction=0.05, refCat=False)
 
def getMeasuredStarsFromCatalog (catalog, pixToFocal)
 

Function Documentation

◆ canRunTests()

def lsst.jointcal.testUtils.canRunTests ( )
Returns True if the necessary packages and files are available.

We need ``obs_cfht`` to load the test/data/cfht_minimal dataset, which
includes the metadata that is used to build the fake catalogs.

Definition at line 39 of file testUtils.py.

39 def canRunTests():
40  """Returns True if the necessary packages and files are available.
41 
42  We need ``obs_cfht`` to load the test/data/cfht_minimal dataset, which
43  includes the metadata that is used to build the fake catalogs.
44  """
45  try:
46  import lsst.obs.cfht # noqa: F401
47  return True
48  except ImportError:
49  return False
50 
51 

◆ createFakeCatalog()

def lsst.jointcal.testUtils.createFakeCatalog (   num,
  bbox,
  fluxFieldName,
  skyWcs = None,
  refCat = False 
)
Return a fake minimally-useful catalog for jointcal.

Parameters
----------
num : `int`
    Number of sources to put in the catalogs. Should be
    a square, to have sqrt(num) centroids on a grid.
bbox : `lsst.geom.Box2I`
    Bounding Box of the detector to populate.
fluxFieldName : `str`
    Name of the flux field to populate in the catalog, without `_instFlux`
    (e.g. "slot_CalibFlux").
skyWcs : `lsst.afw.geom.SkyWcs` or None, optional
    If supplied, use this to fill in coordinates from centroids.
refCat : `bool`, optional
    Return a ``SimpleCatalog`` so that it behaves like a reference catalog?

Returns
-------
catalog : `lsst.afw.table.SourceCatalog`
    A populated source catalog.

Definition at line 188 of file testUtils.py.

188 def createFakeCatalog(num, bbox, fluxFieldName, skyWcs=None, refCat=False):
189  """Return a fake minimally-useful catalog for jointcal.
190 
191  Parameters
192  ----------
193  num : `int`
194  Number of sources to put in the catalogs. Should be
195  a square, to have sqrt(num) centroids on a grid.
196  bbox : `lsst.geom.Box2I`
197  Bounding Box of the detector to populate.
198  fluxFieldName : `str`
199  Name of the flux field to populate in the catalog, without `_instFlux`
200  (e.g. "slot_CalibFlux").
201  skyWcs : `lsst.afw.geom.SkyWcs` or None, optional
202  If supplied, use this to fill in coordinates from centroids.
203  refCat : `bool`, optional
204  Return a ``SimpleCatalog`` so that it behaves like a reference catalog?
205 
206  Returns
207  -------
208  catalog : `lsst.afw.table.SourceCatalog`
209  A populated source catalog.
210  """
212  # centroid
213  centroidKey = lsst.afw.table.Point2DKey.addFields(schema, "centroid", "centroid", "pixels")
214  xErrKey = schema.addField("centroid_xErr", type="F")
215  yErrKey = schema.addField("centroid_yErr", type="F")
216  # shape
217  shapeKey = lsst.afw.table.QuadrupoleKey.addFields(schema, "shape", "",
218  lsst.afw.table.CoordinateType.PIXEL)
219  # Put the fake sources in the minimal catalog.
220  schema.addField(fluxFieldName+"_instFlux", type="D", doc="post-ISR instFlux")
221  schema.addField(fluxFieldName+"_instFluxErr", type="D", doc="post-ISR instFlux stddev")
222  schema.addField(fluxFieldName+"_flux", type="D", doc="source flux (nJy)")
223  schema.addField(fluxFieldName+"_fluxErr", type="D", doc="flux stddev (nJy)")
224  schema.addField(fluxFieldName+"_mag", type="D", doc="magnitude")
225  schema.addField(fluxFieldName+"_magErr", type="D", doc="magnitude stddev")
226  return fillCatalog(schema, num, bbox,
227  centroidKey, xErrKey, yErrKey, shapeKey, fluxFieldName,
228  skyWcs=skyWcs, refCat=refCat)
229 
230 
static QuadrupoleKey addFields(Schema &schema, std::string const &name, std::string const &doc, CoordinateType coordType=CoordinateType::PIXEL)
Add a set of quadrupole subfields to a schema and return a QuadrupoleKey that points to them.
Definition: aggregates.cc:100
static Schema makeMinimalSchema()
Return a minimal schema for Source tables and records.
Definition: Source.h:260
def fillCatalog(schema, num, bbox, centroidKey, xErrKey, yErrKey, shapeKey, fluxFieldName, skyWcs=None, fluxErrFraction=0.05, refCat=False)
Definition: testUtils.py:233
def createFakeCatalog(num, bbox, fluxFieldName, skyWcs=None, refCat=False)
Definition: testUtils.py:188

◆ createFakeCcdImage()

def lsst.jointcal.testUtils.createFakeCcdImage (   butler,
  visit,
  num,
  fluxFieldName,
  photoCalibMean = 1e-2,
  photoCalibErr = 1.0,
  fakeCcdId = 12,
  fakeWcs = None,
  fakeVisitInfo = None 
)
Create a fake CcdImage by making a fake catalog.

Parameters
----------
butler : `lsst.daf.persistence.Butler`
    Butler to load metadata from.
visit : `int`
    Visit identifier to build a butler dataId.
num : `int`
    Number of sources to put in the catalogs. Should be
    a square, to have sqrt(num) centroids on a grid.
fluxFieldName : `str`
    Name of the flux field to populate in the catalog, without `_instFlux`
    (e.g. "slot_CalibFlux").
photoCalibMean : `float`, optional
    Value to set for calibrationMean in the created PhotoCalib.
    Note: this value is 1/instFluxMag0, so it should be less than 1.
photoCalibErr : `float`, optional
    Value to set for calibrationErr in the created PhotoCalib.
fakeCcdId : `int`, optional
    Use this as the ccdId in the returned CcdImage.
fakeWcs : `lsst.afw.geom.SkyWcs`, optional
    A SkyWcs to use instead of one read from disk.
fakeVisitInfo : `lsst.afw.image.VisitInfo`, optional
    A VisitInfo to use instead of one read from disk.

Returns
-------
struct : `lsst.pipe.base.Struct`
   Result struct with components:

   - `catalog` : Catalogs containing fake sources
       (`lsst.afw.table.SourceCatalog`).
   - `ccdImage` : CcdImage containing the metadata and fake sources
       (`lsst.jointcal.CcdImage`).
   - `bbox` : Bounding Box of the image (`lsst.geom.Box2I`).
   - `skyWcs` : SkyWcs of the image (`lsst.afw.geom.SkyWcs`).

Definition at line 130 of file testUtils.py.

132  fakeWcs=None, fakeVisitInfo=None):
133  """Create a fake CcdImage by making a fake catalog.
134 
135  Parameters
136  ----------
137  butler : `lsst.daf.persistence.Butler`
138  Butler to load metadata from.
139  visit : `int`
140  Visit identifier to build a butler dataId.
141  num : `int`
142  Number of sources to put in the catalogs. Should be
143  a square, to have sqrt(num) centroids on a grid.
144  fluxFieldName : `str`
145  Name of the flux field to populate in the catalog, without `_instFlux`
146  (e.g. "slot_CalibFlux").
147  photoCalibMean : `float`, optional
148  Value to set for calibrationMean in the created PhotoCalib.
149  Note: this value is 1/instFluxMag0, so it should be less than 1.
150  photoCalibErr : `float`, optional
151  Value to set for calibrationErr in the created PhotoCalib.
152  fakeCcdId : `int`, optional
153  Use this as the ccdId in the returned CcdImage.
154  fakeWcs : `lsst.afw.geom.SkyWcs`, optional
155  A SkyWcs to use instead of one read from disk.
156  fakeVisitInfo : `lsst.afw.image.VisitInfo`, optional
157  A VisitInfo to use instead of one read from disk.
158 
159  Returns
160  -------
161  struct : `lsst.pipe.base.Struct`
162  Result struct with components:
163 
164  - `catalog` : Catalogs containing fake sources
165  (`lsst.afw.table.SourceCatalog`).
166  - `ccdImage` : CcdImage containing the metadata and fake sources
167  (`lsst.jointcal.CcdImage`).
168  - `bbox` : Bounding Box of the image (`lsst.geom.Box2I`).
169  - `skyWcs` : SkyWcs of the image (`lsst.afw.geom.SkyWcs`).
170  """
171  ccdId = 12 # we only have data for ccd=12
172 
173  dataId = dict(visit=visit, ccd=ccdId)
174  skyWcs = fakeWcs if fakeWcs is not None else butler.get('calexp_wcs', dataId=dataId)
175  visitInfo = fakeVisitInfo if fakeVisitInfo is not None else butler.get('calexp_visitInfo', dataId=dataId)
176  bbox = butler.get('calexp_bbox', dataId=dataId)
177  detector = butler.get('calexp_detector', dataId=dataId)
178  filt = butler.get("calexp_filter", dataId=dataId).getName()
179  photoCalib = lsst.afw.image.PhotoCalib(photoCalibMean, photoCalibErr)
180 
181  catalog = createFakeCatalog(num, bbox, fluxFieldName, skyWcs=skyWcs)
182  ccdImage = lsst.jointcal.ccdImage.CcdImage(catalog, skyWcs, visitInfo, bbox, filt, photoCalib,
183  detector, visit, fakeCcdId, fluxFieldName)
184 
185  return lsst.pipe.base.Struct(catalog=catalog, ccdImage=ccdImage, bbox=bbox, skyWcs=skyWcs)
186 
187 
The photometric calibration of an exposure.
Definition: PhotoCalib.h:114
std::string const & getName() const noexcept
Return a filter's name.
Definition: Filter.h:78

◆ createTwoFakeCcdImages()

def lsst.jointcal.testUtils.createTwoFakeCcdImages (   num1 = 4,
  num2 = 4,
  seed = 100,
  fakeCcdId = 12,
  photoCalibMean1 = 1e-2,
  photoCalibMean2 = 1.2e-2,
  fakeWcses = (None, None),
  fakeVisitInfos = (None, None) 
)
Return two fake ccdImages built on CFHT Megacam metadata.

If ``num1 == num2``, the catalogs will align on-sky so each source will
have a match in the other catalog.

This uses the butler dataset stored in `tests/data/cfht_minimal` to
bootstrap the metadata.

Parameters
----------
num1, num2 : `int`, optional
    Number of sources to put in the first and second catalogs. Should be
    a square, to have sqrt(num) centroids on a grid.
seed : `int`, optional
    Seed value for np.random.
fakeCcdId : `int`, optional
    Sensor identifier to use for both CcdImages. The wcs, bbox, photoCalib, etc.
    will still be drawn from the CFHT ccd=12 files, as that is the only
    testdata that is included in this simple test dataset.
photoCalibMean1, photoCalibMean2: `float`, optional
    The mean photometric calibration to pass to each ccdImage construction.
    Note: this value is 1/instFluxMag0, so it should be less than 1.
fakeWcses : `list` [`lsst.afw.geom.SkyWcs`], optional
    The SkyWcses to use instead of the ones read from disk.
fakeWcses : `list` [`lsst.afw.image.VisitInfo`], optional
    The VisitInfos to use instead of the ones read from disk.

Returns
-------
struct : `lsst.pipe.base.Struct`
   Result struct with components:

   - `camera` : Camera representing these catalogs
       (`lsst.afw.cameraGeom.Camera`).
   - `catalogs` : Catalogs containing fake sources
       (`list` of `lsst.afw.table.SourceCatalog`).
   - `ccdImageList` : CcdImages containing the metadata and fake sources
       (`list` of `lsst.jointcal.CcdImage`).
   - `bbox` : Bounding Box of the image (`lsst.geom.Box2I`).
   - 'fluxFieldName' : name of the instFlux field in the catalogs ('str').

Definition at line 52 of file testUtils.py.

55  fakeVisitInfos=(None, None)):
56  """Return two fake ccdImages built on CFHT Megacam metadata.
57 
58  If ``num1 == num2``, the catalogs will align on-sky so each source will
59  have a match in the other catalog.
60 
61  This uses the butler dataset stored in `tests/data/cfht_minimal` to
62  bootstrap the metadata.
63 
64  Parameters
65  ----------
66  num1, num2 : `int`, optional
67  Number of sources to put in the first and second catalogs. Should be
68  a square, to have sqrt(num) centroids on a grid.
69  seed : `int`, optional
70  Seed value for np.random.
71  fakeCcdId : `int`, optional
72  Sensor identifier to use for both CcdImages. The wcs, bbox, photoCalib, etc.
73  will still be drawn from the CFHT ccd=12 files, as that is the only
74  testdata that is included in this simple test dataset.
75  photoCalibMean1, photoCalibMean2: `float`, optional
76  The mean photometric calibration to pass to each ccdImage construction.
77  Note: this value is 1/instFluxMag0, so it should be less than 1.
78  fakeWcses : `list` [`lsst.afw.geom.SkyWcs`], optional
79  The SkyWcses to use instead of the ones read from disk.
80  fakeWcses : `list` [`lsst.afw.image.VisitInfo`], optional
81  The VisitInfos to use instead of the ones read from disk.
82 
83  Returns
84  -------
85  struct : `lsst.pipe.base.Struct`
86  Result struct with components:
87 
88  - `camera` : Camera representing these catalogs
89  (`lsst.afw.cameraGeom.Camera`).
90  - `catalogs` : Catalogs containing fake sources
91  (`list` of `lsst.afw.table.SourceCatalog`).
92  - `ccdImageList` : CcdImages containing the metadata and fake sources
93  (`list` of `lsst.jointcal.CcdImage`).
94  - `bbox` : Bounding Box of the image (`lsst.geom.Box2I`).
95  - 'fluxFieldName' : name of the instFlux field in the catalogs ('str').
96  """
97  if not canRunTests():
98  msg = "Necessary packages not available to run tests that use the cfht_minimal dataset."
99  raise unittest.SkipTest(msg)
100 
101  np.random.seed(seed)
102 
103  visit1 = 849375
104  visit2 = 850587
105  fluxFieldName = "SomeFlux"
106 
107  # Load or fake the necessary metadata for each CcdImage
108  dataDir = lsst.utils.getPackageDir('jointcal')
109  inputDir = os.path.join(dataDir, 'tests/data/cfht_minimal')
110  butler = lsst.daf.persistence.Butler(inputDir)
111 
112  # so we can access parts of the camera later (e.g. focal plane)
113  camera = butler.get('camera', visit=visit1)
114 
115  struct1 = createFakeCcdImage(butler, visit1, num1, fluxFieldName,
116  photoCalibMean=photoCalibMean1, photoCalibErr=1.0, fakeCcdId=fakeCcdId,
117  fakeWcs=fakeWcses[0], fakeVisitInfo=fakeVisitInfos[0])
118  struct2 = createFakeCcdImage(butler, visit2, num2, fluxFieldName,
119  photoCalibMean=photoCalibMean2, photoCalibErr=5.0, fakeCcdId=fakeCcdId,
120  fakeWcs=fakeWcses[1], fakeVisitInfo=fakeVisitInfos[1])
121 
122  return lsst.pipe.base.Struct(camera=camera,
123  catalogs=[struct1.catalog, struct2.catalog],
124  ccdImageList=[struct1.ccdImage, struct2.ccdImage],
125  bbox=struct1.bbox,
126  skyWcs=[struct1.skyWcs, struct2.skyWcs],
127  fluxFieldName=fluxFieldName)
128 
129 
def createFakeCcdImage(butler, visit, num, fluxFieldName, photoCalibMean=1e-2, photoCalibErr=1.0, fakeCcdId=12, fakeWcs=None, fakeVisitInfo=None)
Definition: testUtils.py:132
std::string getPackageDir(std::string const &packageName)
return the root directory of a setup package
Definition: packaging.cc:33

◆ fillCatalog()

def lsst.jointcal.testUtils.fillCatalog (   schema,
  num,
  bbox,
  centroidKey,
  xErrKey,
  yErrKey,
  shapeKey,
  fluxFieldName,
  skyWcs = None,
  fluxErrFraction = 0.05,
  refCat = False 
)
Return a catalog populated with fake, but reasonable, sources.

Centroids are placed on a uniform grid, errors are normally distributed.

Parameters
----------
schema : `lsst.afw.table.Schema`
    Pre-built schema to make the catalog from.
num : `int`
    Number of sources to put in the catalog.
bbox : `lsst.geom.Box2I`
    Bounding box of the ccd to put sources in.
centroidKey : `lsst.afw.table.Key`
    Key for the centroid field to populate.
xErrKey : `lsst.afw.table.Key`
    Key for the xErr field to populate.
yErrKey : `lsst.afw.table.Key`
    Key for the yErr field to populate.
shapeKey : `lsst.afw.table.Key`
    Key for the shape field to populate.
fluxFieldName : `str`
    Name of the flux field to populate in the catalog, without `_instFlux`
    (e.g. "slot_CalibFlux").
skyWcs : `lsst.afw.geom.SkyWcs` or None, optional
    If supplied, use this to fill in coordinates from centroids.
fluxErrFraction : `float`, optional
    Fraction of instFlux to use for the instFluxErr.
refCat : `bool`, optional
    Return a ``SimpleCatalog`` so that it behaves like a reference catalog?

Returns
-------
catalog : `lsst.afw.table.SourceCatalog`
    The filled catalog.

Definition at line 231 of file testUtils.py.

233  skyWcs=None, fluxErrFraction=0.05, refCat=False):
234  """Return a catalog populated with fake, but reasonable, sources.
235 
236  Centroids are placed on a uniform grid, errors are normally distributed.
237 
238  Parameters
239  ----------
240  schema : `lsst.afw.table.Schema`
241  Pre-built schema to make the catalog from.
242  num : `int`
243  Number of sources to put in the catalog.
244  bbox : `lsst.geom.Box2I`
245  Bounding box of the ccd to put sources in.
246  centroidKey : `lsst.afw.table.Key`
247  Key for the centroid field to populate.
248  xErrKey : `lsst.afw.table.Key`
249  Key for the xErr field to populate.
250  yErrKey : `lsst.afw.table.Key`
251  Key for the yErr field to populate.
252  shapeKey : `lsst.afw.table.Key`
253  Key for the shape field to populate.
254  fluxFieldName : `str`
255  Name of the flux field to populate in the catalog, without `_instFlux`
256  (e.g. "slot_CalibFlux").
257  skyWcs : `lsst.afw.geom.SkyWcs` or None, optional
258  If supplied, use this to fill in coordinates from centroids.
259  fluxErrFraction : `float`, optional
260  Fraction of instFlux to use for the instFluxErr.
261  refCat : `bool`, optional
262  Return a ``SimpleCatalog`` so that it behaves like a reference catalog?
263 
264  Returns
265  -------
266  catalog : `lsst.afw.table.SourceCatalog`
267  The filled catalog.
268  """
269  table = lsst.afw.table.SourceTable.make(schema)
270  table.defineCentroid('centroid')
271  table.defineShape('shape')
272  table.defineCalibFlux(fluxFieldName)
273  if refCat:
274  catalog = lsst.afw.table.SimpleCatalog(table)
275  else:
276  catalog = lsst.afw.table.SourceCatalog(table)
277 
278  instFlux = np.random.random(num)*10000
279  instFluxErr = np.abs(instFlux * np.random.normal(fluxErrFraction, scale=0.1, size=num))
280  xx = np.linspace(bbox.getMinX(), bbox.getMaxX(), int(np.sqrt(num)))
281  yy = np.linspace(bbox.getMinY(), bbox.getMaxY(), int(np.sqrt(num)))
282  xv, yv = np.meshgrid(xx, yy)
283  vx = np.random.normal(scale=0.1, size=num)
284  vy = np.random.normal(scale=0.1, size=num)
285 
286  # make all the sources perfectly spherical, for simplicity.
287  mxx = 1
288  myy = 1
289  mxy = 0
290 
291  for i, (x, y) in enumerate(zip(xv.ravel(), yv.ravel())):
292  record = catalog.addNew()
293  record.set('id', i)
294  record.set(centroidKey, lsst.geom.Point2D(x, y))
295  record.set(shapeKey, lsst.afw.geom.ellipses.Quadrupole(mxx, myy, mxy))
296 
297  if skyWcs is not None:
298  lsst.afw.table.updateSourceCoords(skyWcs, catalog)
299 
300  catalog[xErrKey] = vx
301  catalog[yErrKey] = vy
302  catalog[fluxFieldName + '_instFlux'] = instFlux
303  catalog[fluxFieldName + '_instFluxErr'] = instFluxErr
304 
305  return catalog
306 
307 
An ellipse core with quadrupole moments as parameters.
Definition: Quadrupole.h:47
Custom catalog class for record/table subclasses that are guaranteed to have an ID,...
Definition: SortedCatalog.h:42
static std::shared_ptr< SourceTable > make(Schema const &schema, std::shared_ptr< IdFactory > const &idFactory)
Construct a new table.
Definition: Source.cc:382
void updateSourceCoords(geom::SkyWcs const &wcs, SourceCollection &sourceList)
Update sky coordinates in a collection of source objects.
Definition: wcsUtils.cc:96

◆ getMeasuredStarsFromCatalog()

def lsst.jointcal.testUtils.getMeasuredStarsFromCatalog (   catalog,
  pixToFocal 
)
Return a list of measuredStars built from a catalog.

Parameters
----------
catalog : `lsst.afw.table.SourceCatalog`
    The table to get sources from.
pixToFocal : `lsst.afw.geom.TransformPoint2ToPoint2`
    Transform that goes from pixel to focal plane coordinates, to set the
    MeasuredStar x/y focal points.

Returns
-------
stars : `list` of `lsst.jointcal.MeasuredStar`
    MeasuredStars built from the catalog sources.

Definition at line 308 of file testUtils.py.

308 def getMeasuredStarsFromCatalog(catalog, pixToFocal):
309  """Return a list of measuredStars built from a catalog.
310 
311  Parameters
312  ----------
313  catalog : `lsst.afw.table.SourceCatalog`
314  The table to get sources from.
315  pixToFocal : `lsst.afw.geom.TransformPoint2ToPoint2`
316  Transform that goes from pixel to focal plane coordinates, to set the
317  MeasuredStar x/y focal points.
318 
319  Returns
320  -------
321  stars : `list` of `lsst.jointcal.MeasuredStar`
322  MeasuredStars built from the catalog sources.
323  """
324  stars = []
325  for record in catalog:
326  star = lsst.jointcal.star.MeasuredStar()
327  star.x = record.getX()
328  star.y = record.getY()
329  star.setInstFluxAndErr(record.getCalibInstFlux(), record.getCalibInstFluxErr())
330  # TODO: cleanup after DM-4044
331  point = lsst.geom.Point2D(star.x, star.y)
332  pointFocal = pixToFocal.applyForward(point)
333  star.setXFocal(pointFocal.getX())
334  star.setYFocal(pointFocal.getY())
335  stars.append(star)
336 
337  return stars
def getMeasuredStarsFromCatalog(catalog, pixToFocal)
Definition: testUtils.py:308