3 from astropy.io
import fits
10 with fits.open(filename)
as pf:
12 if hdu.name ==
"PRIMARY":
14 elif hdu.name ==
"LDAC_IMHEAD":
17 elif hdu.name ==
"LDAC_OBJECTS":
18 print(
"%d objects" % (len(hdu.data)))
20 ttype = [k
for k, v
in hdu.header.items()
if v ==
"VIGNET"]
22 raise RuntimeError(
"Unable to find a VIGNET column")
23 vignetCol = int(re.search(
r"^TTYPE(\d+)$", ttype[0]).
group(1)) - 1
25 for row
in range(len(hdu.data)):
26 pixelData = hdu.data[row][vignetCol]
27 bad = np.where(pixelData < -1e29)
28 sat = np.where(pixelData > 99e3)
30 mi = afwImage.MaskedImageF(*hdu.data[row][vignetCol].shape)
32 im.getArray()[:] = pixelData
33 msk = mi.getMask().getArray()
34 msk[bad] = afwImage.MaskU.getPlaneBitMask(
"BAD")
35 msk[sat] = afwImage.MaskU.getPlaneBitMask(
"SAT")
36 ds9.mtv(mi, title=row)
41 return psfex.Prefs(filename, md)
def readSExtractor(filename)
def readPrefs(filename, md=None)
table::Key< table::Array< int > > group