使用Numpy / SciPy进行线条平滑处理

Question

我有一行应该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)

e.g。 (1, 2, 4, 3) -> (1, 1.5, 2, 2.5, 3, 3.5, 4, 3.5, 3)

在Numpy / SciPy中有没有“简单”的方法来实现这个目标？

Answer 1

你可以这样做：

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. ]

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