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["Raw Data"] = VALID
OBJECT.show()
OBJECT["Interpolated"] = InterpolatedSpectrum(VALID,label="Interpolated")
OBJECT.show()
OBJECT["Interpolation"] = OBJECT.frame()(wavelengths=WAVELENGTHS_LOWR)
OBJECT.show()
OBJECT["Resampled"] = OBJECT.frame("Interpolated")(wavelengths=WAVELENGTHS_LOWR[1:],resolution=LOWR,method='resample')
OBJECT.show()
plt.legend()
plt.title("Resample vs. Interpolation")
plt.figure(2)
OBJECT["Integrated"] = OBJECT.frame("Interpolated")(wavelengths=WAVELENGTHS[1:],resolution=HIGH_R,method='integrate')
OBJECT.show()
OBJECT["Integrated Quad"] = OBJECT.frame("Interpolated")(wavelengths=WAVELENGTHS[1:],resolution=HIGH_R,method='integrate_quad')
OBJECT.show()
OBJECT["R and Integrated"] = OBJECT.frame("Interpolated")(wavelengths=WAVELENGTHS[1:],resolution=HIGH_R,method='resolve_and_integrate')
OBJECT.show()
plt.legend()
plt.title("High Resolution Integration Methods")
plt.figure(3)
