27 args = [
None,
"MatchList",
None]
31 """Get a list of suitable stars to construct a PSF."""
39 kernelSize = psfPolicy.get(
"kernelSize")
40 borderWidth = psfPolicy.get(
"borderWidth")
41 sizePsfCellX = psfPolicy.get(
"sizeCellX")
42 sizePsfCellY = psfPolicy.get(
"sizeCellY")
44 mi = exposure.getMaskedImage()
46 if display
and displayExposure:
47 disp = afwDisplay.Display(frame=0)
48 disp.mtv(mi, title=
"PSF candidates")
54 ref, source = val[0:2]
55 if not (ref.getFlagForDetection() & measAlg.Flags.STAR)
or \
56 (source.getFlagForDetection() & measAlg.Flags.BAD):
60 cand = measAlg.makePsfCandidate(source, mi)
64 if cand.getWidth() == 0:
65 cand.setBorderWidth(borderWidth)
66 cand.setWidth(kernelSize + 2*borderWidth)
67 cand.setHeight(kernelSize + 2*borderWidth)
69 im = cand.getMaskedImage().getImage()
71 if not np.isfinite(max):
74 psfCellSet.insertCandidate(cand)
76 if display
and displayExposure:
77 disp.dot(
"+", source.getXAstrom() - mi.getX0(), source.getYAstrom() - mi.getY0(),
78 size=4, ctype=afwDisplay.CYAN)
79 disp.dot(
"o", source.getXAstrom() - mi.getX0(), source.getYAstrom() - mi.getY0(),
80 size=4, ctype=afwDisplay.CYAN)
84 source.setFlagForDetection(source.getFlagForDetection() | measAlg.Flags.STAR)
87 return psfStars, psfCellSet