2013-01-09 12 views
9

私は(x、y)のペアで形成される軌道を持っています。スプラインを使ってこの軌道上の点を補間したいと思います。scipy:補間軌道

どうすればよいですか? scipy.interpolate.UnivariateSplineを使用すると、xでもyも単調でないため動作しません。私はパラメータ化を導入することができます(例えば、軌道に沿って長さd)が、その後、2つの従属変数x(d)y(d)があります。

例:あなたは任意の形状の曲線上で補間することができsplprepを使用

import numpy as np 
import matplotlib.pyplot as plt 
import math 

error = 0.1 
x0 = 1 
y0 = 1 
r0 = 0.5 

alpha = np.linspace(0, 2*math.pi, 40, endpoint=False) 
r = r0 + error * np.random.random(len(alpha)) 
x = x0 + r * np.cos(alpha) 
y = x0 + r * np.sin(alpha) 
plt.scatter(x, y, color='blue', label='given') 

# For this special case, the following code produces the 
# desired results. However, I need something that depends 
# only on x and y: 
from scipy.interpolate import interp1d 
alpha_i = np.linspace(alpha[0], alpha[-1], 100) 
r_i = interp1d(alpha, r, kind=3)(alpha_i) 
x_i = x0 + r_i * np.cos(alpha_i) 
y_i = x0 + r_i * np.sin(alpha_i) 
plt.plot(x_i, y_i, color='green', label='desired') 

plt.legend() 
plt.show() 

example data

+0

@alinsoar:この例の(x、y)データはしませんか? – Nikratio

答えて

11

from scipy import interpolate 
tck,u=interpolate.splprep([x,y],s=0.0) 
x_i,y_i= interpolate.splev(np.linspace(0,1,100),tck) 

これは、与えられたようなプロットを生成しますが、アルファとアルファのパラメタではなく、xとyの点のみを使用します。 Same as yours only using x and y points.

私の元の回答は申し訳ありませんが、私は質問を誤解しました。