約4000ポイントのデータセット(恒星スペクトル)があります。x
とデータ(同じ長さ)を、 17000ポイント近くあるy
この新しいリストで補間したいと思います。x_new
numpy.interp
区分的線形補間であることを除いて、thasがこの仕事をしていることがわかりました。なめらかにしたいです。
助言がありますか?
約4000ポイントのデータセット(恒星スペクトル)があります。x
とデータ(同じ長さ)を、 17000ポイント近くあるy
この新しいリストで補間したいと思います。x_new
numpy.interp
区分的線形補間であることを除いて、thasがこの仕事をしていることがわかりました。なめらかにしたいです。
助言がありますか?
scipy.interpolate.interp1dには、たとえばkind
に設定できるパラメーターがあります"cubic"
。
補間する x 値をこの関数に直接指定する代わりに、次のようにします。
>>> x = np.arange(0, 10)
>>> y = np.exp(-x / 3.0)
>>> f = interpolate.interp1d(x, y, "cubic")
>>> i = np.arange(0.5, 9.5, 1)
>>> f(i)
array([ 0.85826566, 0.61497421, 0.44064828, 0.31573829, 0.22623637,
0.16210544, 0.11615363, 0.08322771, 0.05963526])
これは役立つかもしれない scipy 補間チュートリアルです。