def May14(): clf() s = RamanSpectrum("/home/chris/Documents/DataWeiss/150514/150514_12.txt") a = RamanSpectrum("/home/chris/Documents/DataWeiss/150514/150514_13.txt") b = RamanSpectrum("/home/chris/Documents/DataWeiss/150514/150514_14.txt") c = add_RamanSpectra(a, b) # c = add_RamanSpectra(s,c) d = RamanSpectrum("/home/chris/Documents/DataWeiss/150514/150514_15.txt") v = RamanSpectrum("/home/chris/Documents/DataWeiss/150514/150514_16.txt") z = 0.15 e = subtract_RamanSpectra(c, d * z) # e.smooth() l = subtract_RamanSpectra(c, v * z) # l.smooth() e = RamanSpectrum(e.append(l)) e = autobaseline(e, (180, 277), order=4) e = autobaseline(e, (277, 1700), order=4, join="start") # j = subtract_RamanSpectra(c,RamanSpectrum(pandas.Series(f[2],f[1]))) e.plot() # d.plot() # v.plot() # CdMeOTPRef.index = array(CdMeOTPRef.index)-5 (CdMeOTPRef / 120).plot() (MeOTPRef / 240).plot() return 0
def Feb10(): figure() a = RamanSpectrum('/home/chris/Documents/DataWeiss/150210/43_ long scan.SPE') adding = pandas.Series([NaN]*len(arange(300,1500,0.5)),arange(300,1500,0.5)) d=a.append(adding) d = d.interpolate(method='index') d = d[arange(300,1500,0.5)] e = FourierFilter(d,width = 1100) e.plot() b = RamanSpectrum('/home/chris/Documents/DataWeiss/150210/44.SPE') c=b+a a.autobaseline((300,1600),order = 4) #a = smooth(a,window_len=9) b.autobaseline((300,1600),order = 4) #b=smooth(b,window_len=9) #a.plot() #b.plot() c= autobaseline(c,(300,1600),order = 4) c = smooth(c, window_len=9) #c.plot() legend(['a','d','e']) return d