用Numpy / SciPy对线进行平滑处理

问题描述

|| 我有一条线应以

scipy.interpolate.splrep

和

scipy.interpolate.splev

平滑。

line = ((x1,y1),(x2,y2),... (xn,yn))
tck = interpolate.splrep(x,y)

我需要为我的x坐标找到更多值，这些值应该均匀排列。

newx = numpy.XXX(x)
newy = interpolate.splev(newx,tck)

例如

(1,2,4,3) -> (1,1.5,2.5,3,3.5,3)

在Numpy / SciPy中是否有一种“简单”的方法来实现这一目标？     

解决方法

        您可以执行以下操作：

import scipy.interpolate as interp
z = arange(0,4)
x = np.array([1,2,4,3])
f = interp.interp1d(z,x)
newx = f(np.linspace(z[0],z[-1],7))

这应该给你

In [40]: print z
[0 1 2 3]

In [41]: print x
[1 2 4 3]

In [42]: print newx
[ 1.   1.5  2.   3.   4.   3.5  3. ]

它将按照在数组中定义的顺序在横坐标点之间线性插值。那是你在想的吗？