[pyar] griddata with irregular data
Kiko
kikocorreoso en gmail.com
Lun Oct 3 14:13:33 ART 2016
un punto en un plano se define usando coordenada x (lon, en tu caso) e
y (lat, en tu caso). Es por ello que para definir puntos necesitas la
longitud y latitud de cada uno de ellos, de ahí que tus arrays. lat y
lon deban tener la misma longitud.
El 3/10/16, Cristian Vazquez <cristianvazquez4 en gmail.com> escribió:
> buenas, Como va?
>
> Alguien tiene experiencia interpolando datos?
>
> Quiero usar griddata de scipy con los siguientes datos
>
> Run Roms0
> x0: 60 x1: 150 y0 200 y1 350
> len(lon0_z) 90 len(lat0_z) 150
> len(lon1_z) 151 len(lat1_z) 166
> len(z) 2
> len(z[0,:,:]) 150
> z.shape (2, 150, 90)
> lon0_z.shape (90,)
> lat0_z.shape (150,)
> lon1_z.shape (151,)
> lat1_z.shape (166,)
> len(t) 2
>
> k = 1
>
> griddata(lon0_z,lat0_z,z[k,:,:],(lon1_z[None,:],lat1_z[:,None]))
>
> pero no le gusta, primero no puedo agrupar lon0_z con lat0_z y segundo
> parece que las dimenciones de lon0_z y lat0_z tiene que ser identicas.
>
> alguna idea?
>
> --
> Saludos.
> Cristian
>
Más información sobre la lista de distribución pyar