Bonjour,
m'essayant au tracé de graphiques en 2D, j'essaye de faire une interpolation via griddata, en suivant cet exemple :
Le résultat donne ceci :
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()
Donc, je me lance dans un script qui récupère des données : le voici :
Le résultat donne cela :
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()
Donc j'essaye le griddata (je désactive les commentaires) et là, j’obtiens un message d'erreur dans la console :
J'ai beau chercher, je n'arrive à voir la cause de ce message d'erreur.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'
Si une âme charitable passe par là...
Merci
Tintin
Partager