lsst.meas.extensions.psfex.utils Namespace Reference

Functions
def	readSExtractor (filename)
def	readPrefs (filename, md=None)

Function Documentation

readPrefs()

def lsst.meas.extensions.psfex.utils.readPrefs	(		filename,
			md = None
	)

Definition at line 42 of file utils.py.

def readPrefs(filename, md=None):
    return psfex.Prefs(filename, md)

readSExtractor()

def lsst.meas.extensions.psfex.utils.readSExtractor

(

filename

)

Definition at line 11 of file utils.py.

def readSExtractor(filename):
    with fits.open(filename) as pf:
        for hdu in pf:
            if hdu.name == "PRIMARY":
                pass
            elif hdu.name == "LDAC_IMHEAD":
                hdr = hdu.data[0][0]    # the fits header from the original fits image
                print(hdr[3])
            elif hdu.name == "LDAC_OBJECTS":
                print("%d objects" % (len(hdu.data)))
                # Find the VIGNET column
                ttype = [k for k, v in hdu.header.items() if v == "VIGNET"]
                if not ttype:
                    raise RuntimeError("Unable to find a VIGNET column")
                vignetCol = int(re.search(r"^TTYPE(\d+)$", ttype[0]).group(1)) - 1
 
                for row in range(len(hdu.data)):
                    pixelData = hdu.data[row][vignetCol]
                    bad = np.where(pixelData < -1e29)
                    sat = np.where(pixelData > 99e3)
                    pixelData[bad] = 0.0
                    mi = afwImage.MaskedImageF(*hdu.data[row][vignetCol].shape)
                    im = mi.getImage()
                    im.getArray()[:] = pixelData
                    msk = mi.getMask().getArray()
                    msk[bad] = afwImage.MaskU.getPlaneBitMask("BAD")
                    msk[sat] = afwImage.MaskU.getPlaneBitMask("SAT")
                    afwDisplay.Display().mtv(mi, title=row)
                    input("Next ")