Discussion :

[Ensitof]

Bonjour,

Je ne suis pas habitué à poster des questions sur des forums, mais là, je bloque sur un problème qui, j’en suis certain, tient à un détail qui m’échappe.

J’aimerai représenter les résultat d’un levé (x,y,z) en colorant mes points en fonction de leur hauteur (valeur de z). La fonction scatter3D de matplotlib est faite pour ça visiblement. Tout marche du tonnerre pour la colorbar (jet) qui est bien limitée aux valeurs extrêmes de z mais... rien de chez rien pour la coloration de mes données (couleur uniforme alors que je demande de colorer en fonction de z via mon argument c=z).

Voici le bout de code relatif:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
from mpl_toolkits.mplot3d.axes3d import *
import pylab as plt
x=y=z=range(10)
fig = plt.figure()
ax = Axes3D(fig)
col=ax.scatter3D(x,y,z,c=z,cmap=plt.cm.jet)
cbar = fig.colorbar(col,shrink=0.9,extend='both')
cbar.ax.set_ylabel('axis 3 data values')
plt.show()

Si je ne représente pas la colorbar sur mon graphique de sortie (code un peu plus simple), j'ai logiquement (et malheureusement) toujours le même résultat, à savoir une couleur uniforme pour mes points:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
 
from mpl_toolkits.mplot3d.axes3d import *
import pylab as plt
x=y=z=range(10)
fig = plt.figure()
ax = Axes3D(fig)
ax.scatter3D(x,y,z,c=z,cmap=plt.cm.jet)
plt.show()

Voyez-vous ce qui cloche dans mes arguments?

Merci beaucoup d’avance !

EDIT....: Il s'agit d'un bug, qui a été pris en compte et sera (ou a déjà été) corrigé.

[Aspleet]

Bonjour,

J'utilise python et matplotlib depuis bientôt 2 mois dans le cadre de traitement de données de calcul. Je n'ai donc pas une grande expérience en la matière. Comprendre : merci de me ménager et de parler lentement...

J'aurais besoin de réaliser exactement la fonction décrite dans ce post. A vrai dire, je préfèrerais même une carte colorée à plat si possible. J'ai vérifié avec la dernière version de matplotlib à jour, et ce bug ne semble pas avoir été résolu pour le moment. Par ailleurs, je travaille sur environ 35000 points répartis de manière irrégulière et non une jolie grille de données. A ma connaissance, aucune fonction de matplotlib ne peut plotter cela à part scatter.

Quelqu'un pourrait-il m'aiguiller sur une solution à ce problème ?

Il reste la possibilité bien sûr d'interpoler une surface entre mes 35000 points, puis d'en retirer une grille de points qu'un simple pcolor m'affiche superbement. Ceci étant, mon ordinateur est tout à fait incapable de générer une interpolation aussi lourde, tout du moins avec un simple scipy.interpolate.bisplrep ou scipy.interpolate.BivariateSpline. Je devrais alors coder une fonction qui découpe mon espace en petites zones interpolables séparément, et tenter de recoller ensuite toutes les interpolations ... Et cette perspective ne m'enchante guère.

[Matthieu Brucher]

Surprenant, cela fonctionnait dans une version très ancienne de Matplotlib, OK, c'était avant la réécriture de la partie 3D.
Tu utilises quelle version ?

[Aspleet]

Je suis sous 0.99.1.1

[Matthieu Brucher]

Passe en 1.0.0, j'ai testé il y a quelques jours, chez moi c'est bon.