matplotlib.mlab.griddata muy lento y devuelve una matriz de nan cuando se ingresan datos válidos
Estoy tratando de asignar un conjunto de datos de cuadrícula irregular (datos de satélite sin procesar) con latitudes y longitudes asociadas a un conjunto de latitudes y longitudes cuadriculadas regularmente dado porbasemap.makegrid()
. Estoy usandomatplotlib.mlab.griddata
conmpl_toolkits.natgrid
instalado. A continuación se muestra una lista de las variables que @ @ utiliza como saliwhos
en ipython y algunas estadísticas sobre las variables:
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
ydatalon
son las coordenadas de datos originales
gridlat
ygridlon
son las coordenadas para interpolar a
var
contiene los datos reales
Utilizando estas variables, cuando llamo agriddata(datalon, datalat, var, gridlon, gridlat)
ha tardado hasta 20 minutos en completarse y devuelve una matriz denan
. Al observar los datos, las latitudes y longitudes parecen ser correctas con las coordenadas originales que se superponen a una parte de la nueva área y algunos puntos de datos que se encuentran fuera de la nueva área. ¿Alguien tiene alguna sugerencia? Los valores nan sugieren que estoy haciendo algo estúpido ...