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)