我正在尝试将具有相关纬度和经度的不规则网格数据集(原始卫星数据)映射到由basemap.makegrid()
. 我正在使用matplotlib.mlab.griddata
安装mpl_toolkits.natgrid
。whos
以下是ipython中用作输出的变量列表以及变量的一些统计信息:
Variable Type Data/Info
-------------------------------
datalat ndarray 666x1081: 719946 elems, type `float32`, 2879784 bytes (2 Mb)
datalon ndarray 666x1081: 719946 elems, type `float32`, 2879784 bytes (2 Mb)
gridlat ndarray 1200x1000: 1200000 elems, type `float64`, 9600000 bytes (9 Mb)
gridlon ndarray 1200x1000: 1200000 elems, type `float64`, 9600000 bytes (9 Mb)
var ndarray 666x1081: 719946 elems, type `float32`, 2879784 bytes (2 Mb)
In [11]: var.min()
Out[11]: -30.0
In [12]: var.max()
Out[12]: 30.0
In [13]: datalat.min()
Out[13]: 27.339874
In [14]: datalat.max()
Out[14]: 47.05302
In [15]: datalon.min()
Out[15]: -137.55658
In [16]: datalon.max()
Out[16]: -108.41629
In [17]: gridlat.min()
Out[17]: 30.394031556984299
In [18]: gridlat.max()
Out[18]: 44.237140350357713
In [19]: gridlon.min()
Out[19]: -136.17646180595321
In [20]: gridlon.max()
Out[20]: -113.82353819404671
datalat
是原始datalon
数据坐标
gridlat
并且gridlon
是要插值的坐标
var
包含实际数据
使用这些变量,当我调用griddata(datalon, datalat, var, gridlon, gridlat)
它时,它需要长达 20 分钟才能完成并返回一个nan
. 通过查看数据,纬度和经度似乎是正确的,原始坐标与新区域的一部分和位于新区域之外的一些数据点重叠。有没有人有什么建议?nan 值表明我在做一些愚蠢的事情......