def make_hansuchik(self):
     """Make Hansuchik.FITS file"""
     # Hansuchik.py -> HansuchikUVES.FITS
     from Data.Hansuchik import uves_sky
     HansuchikUVES = SpectraStack(filename=self.config["Dirs"]["Data"] + "/UVESSky.fits")
     HansuchikUVES.save(uves_sky.T,"UVESSky")
     HansuchikUVES.write(clobber=True)
 def make_palsky(self):
     """Make palsky.FITS file"""
     # palsky_100318.dat -> Palsky.FITS
     palskydata = np.genfromtxt(self.config["Dirs"]["Data"] + "/palsky_100318.dat").T
     palSKY = SpectraStack(filename=self.config["Dirs"]["Data"] + "/PalSky.fits")
     palSKY.save(palskydata,"PalSky")
     palSKY.write(clobber=True)
 def make_atmosphere(self):
     """Make Atmosphere.FITS file"""
     # atmosphere.py -> atmosphere.FITS
     from Data.atmosphere import palextinct
     atmosphereEXT = SpectraStack(filename=self.config["Dirs"]["Data"] + "/atmosphere.fits")
     atmosphereEXT.save(palextinct.T,"Atmosph")
     atmosphereEXT.write(clobber=True)
 def make_quimby(self):
     """Make Quimby.FITS file"""
     # Quimby.py -> QuimbySky.FITS
     from Data.Quimby import quimby_sky
     QuimbySKY = SpectraStack(filename=self.config["Dirs"]["Data"] + "/QuimbySky.fits")
     QuimbySKY.save(quimby_sky.T,"QuimbySky")
     QuimbySKY.write(clobber=True)
 def make_turnrose(self):
     """Make Turnrose.FITS file"""
     # Turnrose.py -> Turnrose.FITS
     from Data.Turnrose import skyspec
     TurnroseSKY = SpectraStack(filename=self.config["Dirs"]["Data"] + "/TurnroseSKY.fits")
     TurnroseSKY.save(skyspec.T,"TurnroseSKY")
     TurnroseSKY.write(clobber=True)
 def make_massey(self):
     """Make the Massey.FITS file"""
     # Massey.py -> MasseySky.FITS
     from Data.massey import skyab
     # Convert some sky spec from AB to fLambda
     # This sky spectrum is imported from the massey module, and uses the mag conversion functions above.
     skyflam = np.array([skyab[:,0], abmag_to_flambda(skyab[:,1], skyab[:,0])])
     MasseySky = SpectraStack(filename=self.config["Dirs"]["Data"] + "/MasseySky.fits")
     MasseySky.save(skyab.T,"SkyAB")
     MasseySky.save(skyflam,"SkyFL")
     MasseySky.write(clobber=True)
Exemple #7
0
 def make_atmosphere(self):
     """Make Atmosphere.FITS file"""
     # atmosphere.py -> atmosphere.FITS
     from Data.atmosphere import palextinct
     atmosphereEXT = SpectraStack(filename=self.config["Dirs"]["Data"] +
                                  "/atmosphere.fits")
     atmosphereEXT.save(palextinct.T, "Atmosph")
     atmosphereEXT.write(clobber=True)
Exemple #8
0
 def make_quimby(self):
     """Make Quimby.FITS file"""
     # Quimby.py -> QuimbySky.FITS
     from Data.Quimby import quimby_sky
     QuimbySKY = SpectraStack(filename=self.config["Dirs"]["Data"] +
                              "/QuimbySky.fits")
     QuimbySKY.save(quimby_sky.T, "QuimbySky")
     QuimbySKY.write(clobber=True)
Exemple #9
0
 def make_turnrose(self):
     """Make Turnrose.FITS file"""
     # Turnrose.py -> Turnrose.FITS
     from Data.Turnrose import skyspec
     TurnroseSKY = SpectraStack(filename=self.config["Dirs"]["Data"] +
                                "/TurnroseSKY.fits")
     TurnroseSKY.save(skyspec.T, "TurnroseSKY")
     TurnroseSKY.write(clobber=True)
Exemple #10
0
 def make_hansuchik(self):
     """Make Hansuchik.FITS file"""
     # Hansuchik.py -> HansuchikUVES.FITS
     from Data.Hansuchik import uves_sky
     HansuchikUVES = SpectraStack(filename=self.config["Dirs"]["Data"] +
                                  "/UVESSky.fits")
     HansuchikUVES.save(uves_sky.T, "UVESSky")
     HansuchikUVES.write(clobber=True)
Exemple #11
0
 def make_palsky(self):
     """Make palsky.FITS file"""
     # palsky_100318.dat -> Palsky.FITS
     palskydata = np.genfromtxt(self.config["Dirs"]["Data"] +
                                "/palsky_100318.dat").T
     palSKY = SpectraStack(filename=self.config["Dirs"]["Data"] +
                           "/PalSky.fits")
     palSKY.save(palskydata, "PalSky")
     palSKY.write(clobber=True)
Exemple #12
0
 def make_massey(self):
     """Make the Massey.FITS file"""
     # Massey.py -> MasseySky.FITS
     from Data.massey import skyab
     # Convert some sky spec from AB to fLambda
     # This sky spectrum is imported from the massey module, and uses the mag conversion functions above.
     skyflam = np.array(
         [skyab[:, 0],
          abmag_to_flambda(skyab[:, 1], skyab[:, 0])])
     MasseySky = SpectraStack(filename=self.config["Dirs"]["Data"] +
                              "/MasseySky.fits")
     MasseySky.save(skyab.T, "SkyAB")
     MasseySky.save(skyflam, "SkyFL")
     MasseySky.write(clobber=True)
Exemple #13
0
from AstroObject.AstroObjectLogging import *
from AstroObject.AnalyticSpectra import InterpolatedSpectrum,GaussianSpectrum,FlatSpectrum,BlackBodySpectrum,UnitarySpectrum,Resolver
from AstroObject.AstroSpectra import SpectraStack

LOG = logging.getLogger('AstroObject')
LOG.configure()
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["Post Interpolation"] = OBJECT.frame()(wavelengths=WAVELENGTHS)
OBJECT.show()
OBJECT["Resampled"] = OBJECT.frame("Interpolated")(wavelengths=WAVELENGTHS_LOWR[1:],resolution=LOWR,method='resample')
OBJECT.show()
plt.legend()
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()