LSST Applications g013ef56533+d2224463a4,g199a45376c+0ba108daf9,g19c4beb06c+9f335b2115,g1fd858c14a+2459ca3e43,g210f2d0738+2d3d333a78,g262e1987ae+abbb004f04,g2825c19fe3+eedc38578d,g29ae962dfc+0cb55f06ef,g2cef7863aa+aef1011c0b,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+19c3a54948,g47891489e3+501a489530,g4cdb532a89+a047e97985,g511e8cfd20+ce1f47b6d6,g53246c7159+8c5ae1fdc5,g54cd7ddccb+890c8e1e5d,g5fd55ab2c7+951cc3f256,g64539dfbff+2d3d333a78,g67b6fd64d1+501a489530,g67fd3c3899+2d3d333a78,g74acd417e5+0ea5dee12c,g786e29fd12+668abc6043,g87389fa792+8856018cbb,g89139ef638+501a489530,g8d7436a09f+5ea4c44d25,g8ea07a8fe4+81eaaadc04,g90f42f885a+34c0557caf,g9486f8a5af+165c016931,g97be763408+d5e351dcc8,gbf99507273+8c5ae1fdc5,gc2a301910b+2d3d333a78,gca7fc764a6+501a489530,gce8aa8abaa+8c5ae1fdc5,gd7ef33dd92+501a489530,gdab6d2f7ff+0ea5dee12c,ge410e46f29+501a489530,geaed405ab2+e3b4b2a692,gf9a733ac38+8c5ae1fdc5,w.2025.41
LSST Data Management Base Package
Loading...
Searching...
No Matches
lsst.pipe.tasks.measurePsf Namespace Reference

Classes

class  MeasurePsfConfig
 
class  MeasurePsfTask
 
class  NonfinitePsfShapeError
 

Functions

 showPsfSpatialCells (exposure, cellSet, showBadCandidates, frame=1)
 
 plotPsfCandidates (cellSet, showBadCandidates=False, frame=1)
 
 plotResiduals (exposure, cellSet, showBadCandidates=False, normalizeResiduals=True, frame=2)
 

Function Documentation

◆ plotPsfCandidates()

lsst.pipe.tasks.measurePsf.plotPsfCandidates ( cellSet,
showBadCandidates = False,
frame = 1 )

Definition at line 330 of file measurePsf.py.

330def plotPsfCandidates(cellSet, showBadCandidates=False, frame=1):
331 stamps = []
332 for cell in cellSet.getCellList():
333 for cand in cell.begin(not showBadCandidates): # maybe include bad candidates
334 try:
335 im = cand.getMaskedImage()
336
337 chi2 = cand.getChi2()
338 if chi2 < 1e100:
339 chi2 = "%.1f" % chi2
340 else:
341 chi2 = float("nan")
342
343 stamps.append((im, "%d%s" %
344 (maUtils.splitId(cand.getSource().getId(), True)["objId"], chi2),
345 cand.getStatus()))
346 except Exception:
347 continue
348
349 mos = afwDisplay.utils.Mosaic()
350 disp = afwDisplay.Display(frame=frame)
351 for im, label, status in stamps:
352 im = type(im)(im, True)
353 try:
354 im /= afwMath.makeStatistics(im, afwMath.MAX).getValue()
355 except NotImplementedError:
356 pass
357
358 mos.append(im, label,
359 afwDisplay.GREEN if status == afwMath.SpatialCellCandidate.GOOD else
360 afwDisplay.YELLOW if status == afwMath.SpatialCellCandidate.UNKNOWN else afwDisplay.RED)
361
362 if mos.images:
363 disp.mtv(mos.makeMosaic(), title="Psf Candidates")
364
365
Statistics makeStatistics(lsst::afw::image::Image< Pixel > const &img, lsst::afw::image::Mask< image::MaskPixel > const &msk, int const flags, StatisticsControl const &sctrl=StatisticsControl())
Handle a watered-down front-end to the constructor (no variance)
Definition Statistics.h:361

◆ plotResiduals()

lsst.pipe.tasks.measurePsf.plotResiduals ( exposure,
cellSet,
showBadCandidates = False,
normalizeResiduals = True,
frame = 2 )

Definition at line 366 of file measurePsf.py.

366def plotResiduals(exposure, cellSet, showBadCandidates=False, normalizeResiduals=True, frame=2):
367 psf = exposure.getPsf()
368 disp = afwDisplay.Display(frame=frame)
369 while True:
370 try:
371 maUtils.showPsfCandidates(exposure, cellSet, psf=psf, display=disp,
372 normalize=normalizeResiduals,
373 showBadCandidates=showBadCandidates)
374 frame += 1
375 maUtils.showPsfCandidates(exposure, cellSet, psf=psf, display=disp,
376 normalize=normalizeResiduals,
377 showBadCandidates=showBadCandidates,
378 variance=True)
379 frame += 1
380 except Exception:
381 if not showBadCandidates:
382 showBadCandidates = True
383 continue
384 break
385
386 return frame

◆ showPsfSpatialCells()

lsst.pipe.tasks.measurePsf.showPsfSpatialCells ( exposure,
cellSet,
showBadCandidates,
frame = 1 )

Definition at line 317 of file measurePsf.py.

317def showPsfSpatialCells(exposure, cellSet, showBadCandidates, frame=1):
318 disp = afwDisplay.Display(frame=frame)
319 maUtils.showPsfSpatialCells(exposure, cellSet,
320 symb="o", ctype=afwDisplay.CYAN, ctypeUnused=afwDisplay.YELLOW,
321 size=4, display=disp)
322 for cell in cellSet.getCellList():
323 for cand in cell.begin(not showBadCandidates): # maybe include bad candidates
324 status = cand.getStatus()
325 disp.dot('+', *cand.getSource().getCentroid(),
326 ctype=afwDisplay.GREEN if status == afwMath.SpatialCellCandidate.GOOD else
327 afwDisplay.YELLOW if status == afwMath.SpatialCellCandidate.UNKNOWN else afwDisplay.RED)
328
329