Question

I have this data https://drive.google.com/file/d/15ugQJq_85n--eSEs6TnrFltPpzIugaJq/view?usp=sharing

我正在尝试使用<代码>scipy.interpolate上的网格数据。缺乏数据,而且――不准确的数值需要忽略。我的表面上有200个不同的T s和Ks。这是我尝试的:

import pandas as pd
import numpy as np
from scipy.interpolate import griddata

imp_vol=pd.read_csv("data.csv",index_col= dte )

T=imp_vol.index
K=imp_vol.columns

Ti=np.linspace(float(T.min()),float(T.max()),200)
Ki=np.linspace(float(K.min()),float(K.max()),200)

Ti,Ki = np.meshgrid(Ti,Ki)
T,K = np.meshgrid(T,K)

iv_interpol=griddata((T[~np.isnan(imp_vol)],K[~np.isnan(imp_vol)]), imp_vol[~np.isnan(imp_vol)].T, (Ti, Ki),method= linear )

fig = plt.figure()
ax = fig.add_subplot(111, projection= 3d )
ax.plot_surface(Ti, Ki, iv_interpol)

但是,我仍然有不同的方面错误。即使我把阵列变成了1-D,它还是赢得了打工。

我发现的错误是:

“IndexError: boolean指数与面值0的指数阵列不相匹配;尺寸为10,但对应的生物量为15”

这是我不实的数据框架的形象。

https://i.stack.imgur.com/EGWGm.png” rel=“nofollow noreferer”>array

Answer 1

我不得不在发现问题所在之处以及相关阵列的规模方面显示:

印有<条码>。由于这些阵列的面积是20个,因此,等等可能不会成为问题。

But T/K 涉及从档案中填满的数据,而到<> 网格数据的论据有适当的 b索引,例如。

T[~np.isnan(imp_vol)]

T is 2d from meshgrid. What s the shape? dtype?

imp_vol是一个框架,形状不明(我不看相关图像,除非......)。 <代码>np.isnan,大概有2个阵列。看一看这种阵列与<代码”之间的大小有误。 T (or K or ...)。

我赢得了更多的愤怒。你们应该做些事情,不要 random然地tries开方案规划的墙壁, something着一些东西。

edit

从图像的快速可见度来看,当羊毛(15,10)时,我看到15行,10列。 T/K from meshgrid will also (15,10) or (10,15)。从错误中我说(10,15)。您是否有意忽视<代码>indexing=ij ?

Three important debugging steps:

读到整个错误信息,并试图使其感觉到
印刷一些测试价值,特别是形状
重新印发有关职能文件。

我在像这样的情况下因缺乏细致的细节而被迫结束投票。这将拯救我们两个时间。

edit

友情链接