61
61
# write out the fits images, the data numbers are still JUST the RGB
62
62
# scalings; don't use for science
63
red = pyfits.PrimaryHDU()
64
red.header.update('LATOBS', "32:11:56") # add spurious header info
65
red.header.update('LONGOBS', "110:56")
63
red = pyfits.PrimaryHDU(data=npr)
64
red.header['LATOBS'] = "32:11:56" # add spurious header info
65
red.header['LONGOBS'] = "110:56"
67
66
red.writeto('red.fits')
68
green = pyfits.PrimaryHDU()
69
green.header.update('LATOBS', "32:11:56")
70
green.header.update('LONGOBS', "110:56")
68
green = pyfits.PrimaryHDU(data=npg)
69
green.header['LATOBS'] = "32:11:56"
70
green.header['LONGOBS'] = "110:56"
72
71
green.writeto('green.fits')
73
blue = pyfits.PrimaryHDU()
74
blue.header.update('LATOBS', "32:11:56")
75
blue.header.update('LONGOBS', "110:56")
73
blue = pyfits.PrimaryHDU(data=npb)
74
blue.header['LATOBS'] = "32:11:56"
75
blue.header['LONGOBS'] = "110:56"
77
76
blue.writeto('blue.fits')