def tofits(filename, data, hdr=None, clobber=False):
"""simple pyfits wrapper to make saving fits files easier."""
from pyfits import PrimaryHDU, HDUList
hdu = PrimaryHDU(data)
if not hdr is None: hdu.header = hdr
hdulist = HDUList([hdu])
hdulist.writeto(filename, clobber=clobber, output_verify='ignore')
def tofits(filename, data, hdr=None,clobber=False):
"""simple pyfits wrapper to make saving fits files easier."""
from pyfits import PrimaryHDU,HDUList
hdu = PrimaryHDU(data)
if not hdr is None: hdu.header = hdr
hdulist = HDUList([hdu])
hdulist.writeto(filename, clobber=clobber,output_verify='ignore')
ct = new_table([anc,pcc,sxc,syc], tbtype='TableHDU')
ct.name = 'ColdSpot'
ret.append(ct)
for hdu in ret:
hdu.add_checksum()
return ret
if '__main__' == __name__:
import sys, time, datetime
from hashlib import sha1
commit = sha1() # dummy for now
now = datetime.datetime(*time.gmtime()[:6]) # now for now
filename = sys.argv[1]
from pyfits import HDUList
phdu = primary_hdu('GainNoiseColdSpots', now, commit)
pl = gain_noise_payload(linear_fit_gains, median_gains,
overscan_noise, stddev_noise, cold_spots,
analyzed_files, [sha1(f).hexdigest() for f in analyzed_files])
hdul = HDUList([phdu]+pl)
hdul.writeto(filename, checksum=True)
