lsst.pipe.tasks.measurePsf Namespace Reference

Classes
class	MeasurePsfConfig
class	MeasurePsfTask
	Measure the PSF. More...

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

Function Documentation

plotPsfCandidates()

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

Definition at line 356 of file measurePsf.py.

def plotPsfCandidates(cellSet, showBadCandidates=False, frame=1):
    stamps = []
    for cell in cellSet.getCellList():
        for cand in cell.begin(not showBadCandidates):  # maybe include bad candidates
            try:
                im = cand.getMaskedImage()
 
                chi2 = cand.getChi2()
                if chi2 < 1e100:
                    chi2 = "%.1f" % chi2
                else:
                    chi2 = float("nan")
 
                stamps.append((im, "%d%s" %
                               (maUtils.splitId(cand.getSource().getId(), True)["objId"], chi2),
                               cand.getStatus()))
            except Exception:
                continue
 
    mos = afwDisplay.utils.Mosaic()
    disp = afwDisplay.Display(frame=frame)
    for im, label, status in stamps:
        im = type(im)(im, True)
        try:
            im /= afwMath.makeStatistics(im, afwMath.MAX).getValue()
        except NotImplementedError:
            pass
 
        mos.append(im, label,
                   afwDisplay.GREEN if status == afwMath.SpatialCellCandidate.GOOD else
                   afwDisplay.YELLOW if status == afwMath.SpatialCellCandidate.UNKNOWN else afwDisplay.RED)
 
    if mos.images:
        disp.mtv(mos.makeMosaic(), title="Psf Candidates")
 
 

plotResiduals()

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

Definition at line 392 of file measurePsf.py.

def plotResiduals(exposure, cellSet, showBadCandidates=False, normalizeResiduals=True, frame=2):
    psf = exposure.getPsf()
    disp = afwDisplay.Display(frame=frame)
    while True:
        try:
            maUtils.showPsfCandidates(exposure, cellSet, psf=psf, display=disp,
                                      normalize=normalizeResiduals,
                                      showBadCandidates=showBadCandidates)
            frame += 1
            maUtils.showPsfCandidates(exposure, cellSet, psf=psf, display=disp,
                                      normalize=normalizeResiduals,
                                      showBadCandidates=showBadCandidates,
                                      variance=True)
            frame += 1
        except Exception:
            if not showBadCandidates:
                showBadCandidates = True
                continue
        break
 
    return frame

showPsfSpatialCells()

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

Definition at line 343 of file measurePsf.py.

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