Bonjour,
m'essayant au tracé de graphiques en 2D, j'essaye de faire une interpolation via griddata, en suivant cet exemple :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
import numpy as np
from scipy.interpolate import griddata
import matplotlib.pyplot as plt
import numpy.ma as ma
from numpy.random import uniform, seed
# make up some randomly distributed data
seed(1234)
npts = 4
x = uniform(-2,2,npts)
y = uniform(-2,2,npts)
z = x*np.exp(-x**2-y**2)
# define grid.
xi = np.linspace(-2.1,2.1,100)
yi = np.linspace(-2.1,2.1,100)
# grid the data.
zi = griddata((x, y), z, (xi[None,:], yi[:,None]), method='cubic')
 # contour the gridded data, plotting dots at the randomly spaced data points.
CS = plt.contour(xi,yi,zi,15,linewidths=0.5,colors='k')
CS = plt.contourf(xi,yi,zi,15,cmap=plt.cm.jet)
plt.colorbar() # draw colorbar
# plot data points.
plt.scatter(x,y,marker='o',c='b',s=5)
plt.xlim(-2,2)
plt.ylim(-2,2)
plt.title('griddata test (%d points)' % npts)
 
plt.show()
Le résultat donne ceci :


Donc, je me lance dans un script qui récupère des données : le voici :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from scipy.interpolate import griddata
from pylab import *
import numpy as np
 
#lecture donnees
z = np.loadtxt('donnee/donnee.dat')
x = np.loadtxt('donnee/lon.dat')
y = np.loadtxt('donnee/lat.dat')
 
# define grid.
xi = np.linspace(-3.9,6.1,10)
yi = np.linspace(41,51,10)
# grid the data.
#zi = griddata((x, y),z,(xi[None,:], yi[:,None]),method='cubic')
# contour the gridded data, plotting dots at the randomly spaced data points.
#CS = plt.contour(xi,yi,zi,15,linewidths=0.5,colors='k')
#CS = plt.contourf(xi,yi,zi,15,cmap=plt.cm.jet)
#plt.colorbar() # draw colorbar
plt.scatter(x,y,marker='o',c='b',s=5)
plt.xlim(-3.9,6.1)
plt.ylim(41,51)
 
plt.show()
Le résultat donne cela :


Donc j'essaye le griddata (je désactive les commentaires) et là, j’obtiens un message d'erreur dans la console :
Traceback (most recent call last):
File "aleatoire.py", line 16, in <module>
zi = griddata((x, y),z,(xi[None,:], yi[:,None]),method='cubic')
TypeError: griddata() got an unexpected keyword argument 'method'
J'ai beau chercher, je n'arrive à voir la cause de ce message d'erreur.
Si une âme charitable passe par là...
Merci
Tintin