使用兩個不同數據集的Python樣條插值

Question

我有一個包含4000個點的數據集（恆星光譜）。我想插入x和y數據（相同長度）與這個新的名單，x_new其中有近17000點。

我找到了numpy.interp thas做這個工作，除了它是分段線性插值。我希望它平穩。

有什麼建議嗎？

Answer 1

該scipy.interpolate.interp1d有一個kind參數，可以設置爲例如。 "cubic"。

提供你希望直接插值到這個函數的x值的Instad，你做這樣的：

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

This是SciPy的插值教程，可能會有所幫助。