使用scipy sqlrep进行定期插值

Question

我无法让scipy.interpolate.splrep使用定期数据点... 我想我没有正确处理'per'参数。到目前为止，我使用此代码：

from scipy.interpolate import splrep
import numpy as np
x = np.concatenate((np.linspace(-180, 180, 720),np.linspace(180, -180, 720)))
y = np.concatenate((np.linspace(-180, 180, 720),np.linspace(-180, 180, 720)))
spl = splrep(x, y, per=True)

但是它不起作用，请引发“输入数据错误”

Answer 1

如the documentation所述，数据点x和y需要定义一条曲线y = f（x）。这意味着：对于x中的每个点，在y中必须恰好一个点。

例如sin(x)的以下数据：

x中的每个点在y中都只有一个点，并且1d样条拟合效果很好。

现在比较一下您建议的数据：

x中的每个点出现两次，甚至更糟的是，在y中出现两个不同的值。因此，此数据无效。