我尝试使用interpolate.Rbf外推2d样条的给定范围。
从这里插入str：https://justpaste.it/d7u7z到data_str（不需要''字符串连字符）
我怎么能找到Ti的边界之外的zi的值。我希望能够找到zi时Ti = 0。即计算rbf(0,4500)。
我的代码是：

import io
import numpy as np
import pandas as pd
from scipy.interpolate import griddata, Rbf
from numpy import ma
import matplotlib.pyplot as plt
data_str = # insert the string from the url here
vol = pd.read_csv(io.StringIO(data_str))
vol.set_index('dte',inplace=True)
valid_vol=ma.masked_invalid(vol).T
Ti=np.linspace(float((vol.index).min()),float((vol.index).max()),len(vol.index))
Ki=np.linspace(float((vol.columns).min()),float((vol.columns).max()),len(vol.columns))
Ti,Ki = np.meshgrid(Ti,Ki)
valid_Ti = Ti[~valid_vol.mask]
valid_Ki = Ki[~valid_vol.mask]
valid_vol = valid_vol[~valid_vol.mask]
zi = griddata((valid_Ti, valid_Ki), valid_vol, (Ti, Ki), method='cubic')
points = np.column_stack((Ti.ravel(), Ki.ravel()))
values = zi.ravel()
mask = ~np.isnan(values)
values = values[mask]
points = points[mask]
rbf = Rbf(points[:, 0], points[:, 1], values, function='linear')
interp_value = rbf(0.015, 4545)
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
x = np.linspace(Ti.min(),Ti.max(),100)
y = np.linspace(Ki.min(),Ki.max(),100)
x,y=np.meshgrid(x,y)
z = rbf(x,y)
ax.plot_surface(x, y, z, cmap='viridis')
plt.show()

字符串
我没有在bisplev中使用bisplrep，因为它不能生成样条（我读了文档，没有s给出样条）。