用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) 

例如， (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. ] 

这将只是他们在阵列中定义的顺序横坐标点之间的线性插值。那是你在想什么？