lsst.meas.algorithms.psfSelectionFromMatchList Namespace Reference

Functions
def	selectPsfSources

Variables
list	args = [None, "MatchList", None]

Function Documentation

def lsst.meas.algorithms.psfSelectionFromMatchList.selectPsfSources	(		exposure,
			matches,
			psfPolicy
	)

Get a list of suitable stars to construct a PSF.

Definition at line 29 of file psfSelectionFromMatchList.py.

def selectPsfSources(exposure, matches, psfPolicy):
    """Get a list of suitable stars to construct a PSF."""

    import lsstDebug
    display = lsstDebug.Info(__name__).display
    displayExposure = lsstDebug.Info(__name__).displayExposure     # display the Exposure + spatialCells
    #
    # Unpack policy
    #
    kernelSize   = psfPolicy.get("kernelSize")
    borderWidth  = psfPolicy.get("borderWidth")
    sizePsfCellX = psfPolicy.get("sizeCellX")
    sizePsfCellY = psfPolicy.get("sizeCellY")
    #
    mi = exposure.getMaskedImage()

    if display and displayExposure:
        frame = 0
        ds9.mtv(mi, frame=frame, title="PSF candidates")

    psfCellSet = afwMath.SpatialCellSet(mi.getBBox(), sizePsfCellX, sizePsfCellY)
    psfStars = []

    for val in matches:
        ref, source = val[0:2]
        if not (ref.getFlagForDetection() & measAlg.Flags.STAR) or \
               (source.getFlagForDetection() & measAlg.Flags.BAD):
            continue

        try:
            cand = measAlg.makePsfCandidate(source, mi)
            #
            # The setXXX methods are class static, but it's convenient to call them on
            # an instance as we don't know Exposure's pixel type (and hence cand's exact type)
            if cand.getWidth() == 0:
                cand.setBorderWidth(borderWidth)
                cand.setWidth(kernelSize + 2*borderWidth)
                cand.setHeight(kernelSize + 2*borderWidth)

            im = cand.getMaskedImage().getImage()
            max = afwMath.makeStatistics(im, afwMath.MAX).getValue()
            if not numpy.isfinite(max):
                continue

            psfCellSet.insertCandidate(cand)

            if display and displayExposure:
                ds9.dot("+", source.getXAstrom() - mi.getX0(), source.getYAstrom() - mi.getY0(),
                        size=4, frame=frame, ctype=ds9.CYAN)
                ds9.dot("o", source.getXAstrom() - mi.getX0(), source.getYAstrom() - mi.getY0(),
                        size=4, frame=frame, ctype=ds9.CYAN)
        except Exception:
            continue

        source.setFlagForDetection(source.getFlagForDetection() | measAlg.Flags.STAR)
        psfStars += [source]

    return psfStars, psfCellSet

Variable Documentation

list lsst.meas.algorithms.psfSelectionFromMatchList.args = [None, "MatchList", None]

Definition at line 27 of file psfSelectionFromMatchList.py.