check bands

checkband.py

import pyfits
from scipy.constants import c
from scipy.integrate import trapz
rimo=pyfits.open('rimo.fits')
for ext in rimo[3:]:
    f = c/ext.data['WAVENUMBER']/1.e6 #GHz
    t = ext.data['TRANSMISSION']
    print "%s = Sum: %.2f, Max: %.2f, Integ: %.2f" % (ext.name, t.sum(), t.max(), trapz(t[1:], f[1:]))