LOG = logging.getLogger('AstroObject')
LOG.configure_from_file('Examples/config.yml')
LOG.start()
WAVELENGTHS = ((np.arange(98)+1)/2.0 + 1.0) * 1e-7
HIGH_R =WAVELENGTHS[1:]/np.diff(WAVELENGTHS)
WAVELENGTHS_LOWR = ((np.arange(23)+0.25)*2.0 + 1.0) * 1e-7
LOWR = WAVELENGTHS_LOWR[1:]/np.diff(WAVELENGTHS_LOWR)/2
VALID = np.array([(np.arange(50) + 1.0) * 1e-7, np.sin(np.arange(50)/2.0)+2.0 + np.arange(50)/25.0])
OBJECT = SpectraStack()
OBJECT.read("Examples/SNIa.R1000.dat")
OBJECT["Interpolateable"] = InterpolatedSpectrum(OBJECT.d,"Interpolateable")
wl, rs = get_resolution_spectrum(np.min(OBJECT.f.wavelengths),np.max(OBJECT.f.wavelengths),200)
OBJECT["Raw Data"] = OBJECT.f(wavelengths = wl, resolution = rs, method = 'resample')
OBJECT.show()
wl, rs = get_resolution_spectrum(np.min(OBJECT.f.wavelengths),np.max(OBJECT.f.wavelengths),50)
OBJECT["Low Res Data"] = OBJECT["Interpolateable"](wavelengths = wl, resolution = rs, method = 'resample')
OBJECT.show()
for line in OBJECT.info():
print line
print "Valid:",OBJECT.valid()
OBJECT["Raw Data"].logarize()
OBJECT["Logarized"] = OBJECT["Raw Data"]
for line in OBJECT.info():
print line
print "Valid:",OBJECT.valid()
OBJECT.show()
OBJECT["Raw Data"].linearize()
OBJECT["Linearized"] = OBJECT["Raw Data"]
