hw8_sol - write results using 32-bit floats...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
#! /usr/bin/env python # Joseph Harrington # AST 5765/4762 # Homework 8 # 21 October 2007 # Updated 27 Oct 2008 by Kevin Stevenson # Updated 1 Oct 2009 by Joseph Harrington # Problem 1 from hw7_sol import * # Problem 2 # see routine normmedcomb() in medcombine.py import medcombine as mc # Problem 3 # get away from the edges, particularly the junk in the upper right normregion = ((225, 225), (-225, -225)) # ((y1, x1), (y2, x2)) # do the normalized median combination normsky, normfact = mc.normmedcomb(objdata, normregion)
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: # write results, using 32-bit floats pf.writeto('sky_13s_mednorm'+fext, np.float32(normsky), objhead, clobber=True) # Problem 4 # see routine skycormednorm() in medcombine.py # Problem 5 for k in np.arange(nobj): objdata[k] = mc.skycormednorm(objdata[k], normsky, normregion) # save # make the data float32, not float64, to save space pf.writeto(objfile[-1]+'_nosky'+fext, np.float32(objdata[-1]), objhead, clobber=True)...
View Full Document

This note was uploaded on 11/09/2009 for the course AST 4762 taught by Professor Harrington during the Fall '09 term at University of Central Florida.

Ask a homework question - tutors are online