print 'pixel scale value not found, setting to one'
scale = 1.
try:
zpt = float(elements['zeropoint'][0][1])
except:
print 'zeropoint not found, setting to 25.'
zpt = 25.
try:
for t in elements['array']:
if t[0]['name'] == 'ept':
pts = []
for z in t[2]['axis']:
pts.append(map(float, z[1].split('\n')))
tmp = numarray.array(pts)
pts = numarray.swapaxes(tmp, 1, 0)
break
else:
print 'no ept data in .xml file - aborting'
sys.exit()
except:
raise
print 'problem with data in .xml file - aborting'
sys.exit()
data = []
# storing log r (arcsecs), app mag, erase flag, number of pixels (for
# error), prf intensity, exp intensity
top_lum = pts[-1][1] + zpt
imin = 0
hdr = fitsobj[0].header
pix = fitsobj[0].data
fitsobj.close()
doc = minidom.parse(sys.argv[-1 + last].split('.')[0] + '.xml')
rootNode = doc.documentElement
elements = xml_read(rootNode).walk(rootNode)
xsky = float(elements['sky'][0][1])
for t in elements['array']:
if t[0]['name'] == 'marks':
marks = []
for z in t[2]['axis']:
marks.append(map(int, z[1].split('\n')))
tmp = numarray.array(marks)
marks = numarray.swapaxes(tmp, 1, 0)
break
for t in elements['array']:
if t[0]['name'] == 'prf': # ellipse data
prf = []
head = []
pts = []
for z in t[2]['axis']:
head.append(z[0]['name'])
pts.append(map(float, z[1].split('\n')))
for z in range(len(pts[0])):
tmp = []
for w in head:
tmp.append(pts[head.index(w)][z])
prf.append(tmp)