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

這將只是他們在陣列中定義的順序橫座標點之間的線性插值。那是你在想什麼？