Python baseline correction library

by

I found an answer to my question, just sharing for everyone who stumbles upon this.

There is an algorithm called “Asymmetric Least Squares Smoothing” by P. Eilers and H. Boelens in 2005. The paper is free and you can find it on google.

def baseline_als(y, lam, p, niter=10):
  L = len(y)
  D = sparse.csc_matrix(np.diff(np.eye(L), 2))
  w = np.ones(L)
  for i in xrange(niter):
    W = sparse.spdiags(w, 0, L, L)
    Z = W + lam * D.dot(D.transpose())
    z = spsolve(Z, w*y)
    w = p * (y > z) + (1-p) * (y < z)
  return z

More Related Contents:

  1. How to smooth a curve in the right way?
  2. How to smooth a curve for a dataset
  3. Fitting empirical distribution to theoretical ones with Scipy (Python)?
  4. Histogram Matplotlib
  5. Plotting a fast Fourier transform in Python
  6. How to check BLAS/LAPACK linkage in NumPy and SciPy?
  7. How to make scipy.interpolate give an extrapolated result beyond the input range?
  8. Inverse Distance Weighted (IDW) Interpolation with Python
  9. How to create a density plot in matplotlib?
  10. How to do n-D distance and nearest neighbor calculations on numpy arrays
  11. Relationship between SciPy and NumPy
  12. How do I install SciPy on 64 bit Windows?
  13. Multivariate spline interpolation in python/scipy?
  14. Removing duplicate columns and rows from a NumPy 2D array
  15. scipy curve_fit doesn’t like math module
  16. Numpy Resize/Rescale Image
  17. How to specify upper and lower limits when using numpy.random.normal
  18. How to install numpy and scipy for Ironpython27?
  19. how to annotate heatmap with text in matplotlib
  20. How to filter numpy array by list of indices?
  21. Fitting a closed curve to a set of points
  22. Python/Numpy MemoryError
  23. Numpy: Divide each row by a vector element
  24. new pythonic style for shared axes square subplots in matplotlib?
  25. Compute divergence of vector field using python
  26. Python 4D linear interpolation on a rectangular grid
  27. Consecutive, Overlapping Subsets of Array (NumPy, Python)
  28. How to include constraint to Scipy NNLS function solution so that it sums to 1
  29. Computing cross-correlation function?
  30. Understanding scipy deconvolve

Leave a Comment