用SciPy插值不均勻數據

Question

我一直在通過documentation for scipy.interpolate上的例子，我無法弄清楚我將如何採集不均勻間隔的數據並插入它，因爲所有的教程都使用了linspaces - 它們均勻間隔。

例如，我有一些數據傳播，像這樣：

[--1--4-----5-3-22---55-] 

其中每個-代表缺失值。

我將如何去擬合插值函數使用scipy.interpolate？

Answer 1

interpolate.interp1d適用於不均勻間隔的數據。例如，

import re 
import numpy as np 
import scipy.interpolate as interpolate 
import matplotlib.pyplot as plt 

text = '--1--4-----5-3-22---55-' 
parts = [c for c in re.split(r'(-|\d+)', text) if c] 
data = np.array([(x, int(y)) for x, y in enumerate(parts) if y != '-']) 
x, y = data.T 
f = interpolate.interp1d(x, y, kind='cubic') 

newx = np.linspace(x.min(), x.max()) 
newy = f(newx) 
plt.plot(newx, newy) 
plt.scatter(x, y, s=20) 
plt.show() 

產生