所以,我有三个numpy数组,它们在网格上存储纬度,经度和一些属性值 - 也就是说,我有LAT(y,x),LON(y,x),并说温度T(y,x) ),对于x和y的某些限制.网格不一定是规则的 - 事实上,它是三极的.
然后我想将这些属性(温度)值插值到一堆不同的纬度/经度点(存储为lat1(t),lon1(t),大约10,000吨......),这些点不落在实际网格点上.我已经尝试了matplotlib.mlab.griddata,但这需要太长时间(毕竟它并不是为我正在做的事情设计的).我也试过scipy.interpolate.interp2d,但是我得到了一个MemoryError(我的网格大约是400x400).
有没有任何光滑的,最好快速的方式这样做？我不禁想到答案显而易见......谢谢!!
好吧,极端的菜鸟问题在这里.在我的程序中,我生成了一个2D numpy数组,其中一些条目缺失(不是"nan"类型的nonxistant,而是"None"类,或NoneType).我想在这些条目上加上一个掩码,但我似乎遇到了一些麻烦.通常,为了掩盖所有值为2的条目,我会这样做
A = np.ma.masked_where(A [A == 2],A)
在这种情况下,无论我尝试第一个参数,这似乎都不起作用.思考？