Ejemplo n.º 1
0
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
Ejemplo n.º 2
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