IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

Calcul scientifique Python Discussion :

Problème Scipy Griddata : TypeError


Sujet :

Calcul scientifique Python

  1. #1
    Membre à l'essai
    Homme Profil pro
    Inscrit en
    Août 2011
    Messages
    15
    Détails du profil
    Informations personnelles :
    Sexe : Homme

    Informations forums :
    Inscription : Août 2011
    Messages : 15
    Points : 12
    Points
    12
    Par défaut Problème Scipy Griddata : TypeError
    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

  2. #2
    Membre à l'essai
    Homme Profil pro
    Inscrit en
    Août 2011
    Messages
    15
    Détails du profil
    Informations personnelles :
    Sexe : Homme

    Informations forums :
    Inscription : Août 2011
    Messages : 15
    Points : 12
    Points
    12
    Par défaut
    Bon ça y est j'ai résolu mon problème (très bête).
    Dans le début du fichier il faut désactiver la ligne Pylab :
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    from mpl_toolkits.basemap import Basemap
    import matplotlib.pyplot as plt
    from scipy.interpolate import griddata
    #from pylab import *
    import numpy as np
    Et tout baigne.
    Par contre, ne me demandez pas pourquoi ça fait M**der et pourquoi j'ai mis cette foutue ligne...

    Peut-être à bientôt !

+ Répondre à la discussion
Cette discussion est résolue.

Discussions similaires

  1. Réponses: 13
    Dernier message: 15/06/2011, 18h25
  2. Probléme lors d'un Fit sinus avec scipy
    Par astragoth dans le forum Calcul scientifique
    Réponses: 0
    Dernier message: 29/07/2009, 09h49
  3. Problème de hauteur avec un GridData
    Par pathfinder06 dans le forum SWT/JFace
    Réponses: 3
    Dernier message: 18/03/2009, 20h28
  4. Divers problèmes de matrice passant de numpy a scipy
    Par thance dans le forum Calcul scientifique
    Réponses: 25
    Dernier message: 13/06/2007, 00h33
  5. [SWT] Problème de paramètre GridData
    Par yolepro dans le forum SWT/JFace
    Réponses: 4
    Dernier message: 06/12/2002, 11h37

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo