=======
Copyright (2015-2020) C. Le Losq.
lelosq@ipgp.fr
Rampy is a Python library that aims at helping processing spectroscopic data, such as Raman, Infrared or XAS spectra. It offers, for instance, functions to subtract baselines as well as to stack, resample or smooth spectra. It aims at facilitating the use of Python in processing spectroscopic data. It integrates within a workflow that uses Numpy/Scipy as well as optimisation libraries such as lmfit or emcee, for instance.
The /examples/ folder contain various examples.
Rampy is tested on Python 3.6 (see Travis badge; no garantee that it works on other Python versions)
The following libraries are required and indicated in setup.py:
Optional dependencies:
Installation of gcvspline is necessary for use of the rampy.rameau() class.
Installation of cvxpy is necessary for use of the rampy.mixing() class.
Additional libraries for model fitting may be wanted:
Install with pip:
pip install rampy
If you want to use gcvspline and cvxpy, also install it:
pip install gcvspline
pip install cvxpy
Given a signal [x y] containing a peak, and recorded in a text file myspectrum.txt.
You can import it, remove a automatic background, plot the result, and print the centroid of the peak as:
import matplotlib.pyplot as plt
import numpy as np
import rampy as rp
spectrum = np.genfromtxt("myspectrum.txt")
bir = np.array([[0,100., 200., 1000]]) # the frequency regions devoid of signal, used by rp.baseline()
y_corrected, background = rp.baseline(spectrum[:,0],spectrum[:,1],bir,"arPLS",lam=10**10)
plt.figure()
plt.plot(spectrum[:,0],spectrum[:,1],"k",label="raw data")
plt.plot(spectrum[:,0],background,"k",label="background")
plt.plot(spectrum[:,0],y_corrected,"k",label="corrected signal")
plt.show()
print("Signal centroid is %.2f" % rp.centroid(spectrum[:,0],y_corrected))See the /example folder for further examples.
rampy can be used also to analyse the output of the RADIS package.
See for instance https://github.com/charlesll/rampy/issues/13
Updated January 2020