[SciPy-User] Splines in scipy.signal vs scipy.interpolation
Thu Jan 28 08:20:20 CST 2010
On Jan 20, 11:56 pm, Tony S Yu <ton...@MIT.EDU> wrote:
> I'm having trouble making splines from scipy.signal work with those in scipy.interpolation.
> Both packages have functions for creating (`signal.cspline1d`/`interpolate.splrep`) and evaluating (`signal.cspline1d_eval`/`interpolate.splev`) splines. There are, of course, huge differences between these functions, which is why I'm trying to get them to talk to each other.
> In particular, I'd like to create a smoothing spline using `cspline1d` (which allows easier smoothing) and evaluate using `splev` (which allows me to get derivatives of the spline).
bouncing between two murky packages doesn't sound as though it'll
interpolate though has both smoothing and derivs --
interpolator = interpolate.UnivariateSpline( x, y, k=3, s=s )
# s=0 interpolates
yy = interpolator( xx )
y1 = interpolator( xx, 1 ) # deriv
Just curious, are your real knots uniform, how many ?
See also http://projects.scipy.org/scipy/ticket/864
"The documentation for class scipy.interpolate.UnivariateSpline? is
misleading, and maybe completely wrong.
UnivariateSpline? behaves in ways that are unpredictable ...
(Fitpack is just a big dense package => big dense doc.)
More information about the SciPy-User