Tracé de graphiques 3D

Version imprimable

Bonjour à tous!

j'ai avec moi plusieurs fichiers de données comportant des altitudes 'alt' pour lesquelles on a des valeurs de mesures 'Diff' :
Code:

1 2 3 4 5 6 13.151 .0000001717736 13.250 .0000001581861 13.346 0 13.451 0 13.551 .0000001349428
Le nom des fichiers correspond à la date de mesure, et je veux utiliser ces dates en abscisse d'un graphique, avec les altitudes en ordonnée.
Jusque là j'y arrive bien. Voici mon code :
Code:

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 28 import matplotlib.pyplot as plt import numpy as np import glob tabTime = np.zeros(0) tabAlt = np.zeros(0) tabDiff = np.zeros(0) ext="Diff_*" for files in sorted(glob.glob(ext),reverse=False): filename=(files[5:12]) data = np.loadtxt(files) alt = data[:,0] #données altitude diff = data[:,1] #données difference NumberOfLine = len(alt) #nombre de ligne tabAlt = np.append(tabAlt,alt) tabDiff = np.append(tabDiff, diff) for i in range(len(alt)): tabTime = np.append(tabTime,filename) # remplissage tableau Time avec la date des fichiers plt.plot(tabTime, tabAlt,'+')
Ensuite j'aimerais représenter les valeurs de mesures 'Diff' en utilisant une colorbar (j'ai regardé la fonction ' pcolormesh '). Donc à chaque altitude, la couleur des points varierait en fonction des valeurs de 'Diff'.
Mais j'ai l'impression que pour créer ce type de graphique, il faut absolument que la troisième 'variable' affichée soit calculé en fonction des deux autres... Ce qui n'est pas le cas pour mes valeurs.

Est-ce que vous savez comment utiliser ce type de graphique dans mon cas ?

Merci à vous!

Salut,

Voici un exemple reprenant ton code et les quelques données proposées :
Code:

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 28 data = np.array([ [13.151, 0.0000001717736], [13.250, 0.0000001581861], [13.346, 0.0], [13.451, 0.0], [13.551, 0.0000001349428], ]) tabTime = np.zeros(0) tabAlt = np.zeros(0) tabDiff = np.zeros(0) alt = data[:, 0] #données altitude diff = data[:, 1] #données difference NumberOfLine = len(alt) #nombre de ligne tabAlt = np.append(tabAlt, alt) tabDiff = np.append(tabDiff, diff) for i in range(len(alt)): # Date fictive dummy_date = "23-04-2021" tabTime = np.append(tabTime, dummy_date) # remplissage tableau Time avec la date des fichiers plt.scatter(tabTime, tabAlt, marker='+', c=tabDiff) plt.colorbar() plt.show()
L'astuce est de founir en argument de couleur du graph la liste des différences tabDiff.

Je note quelques points qui devraient être revus à mon sens :

Il y a autant de tracés que de fichiers. C'est peut être voulu. Intuitivement je n'en ferais qu'un, et donc sortirais la création de la figure de la boucle de lecture des fichiers. Dans l'exemple que tu fournis, a chaque boucle un nouveau graph se créé avec de nouvelles données en plus. Seul le dernier, affichant toutes les données, a du sens je pense.
La gestion des dates n'est pas claire. J'en déduis que la date est une chaine de caractères. Visiblement matplotlib s'en fiche un peu, mais sauf erreur de ma part quelque soit l'intervalle de temps entre deux dates, l'espacement des points dans le graph sera identique. Pour matplotlib, ce n'est pas vraiment une date, mais un label. Tu pourrais transformer cette chaine de caractères en objet datetime et ainsi permettre à matplotib de ranger dans l'ordre les points avec un espacement qui a du sens.

Ju

2 pièce(s) jointe(s)

Salut!
Merci beaucoup l'astuce a parfaitement marché!

J'ai sorti la figure de la boucle pour n'afficher que le dernier tracé contenant toutes les données. Voici mon code :

Code:

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
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
 
 
import matplotlib.pyplot as plt
import numpy as np
import glob
from datetime import datetime
 
tabTime = np.zeros(0)
tabAlt = np.zeros(0)
tabDiff= np.zeros(0)
 
 
 
ext="Diff_*"
for files in sorted(glob.glob(ext),reverse=False):
    print(files) 
    filename=(files[5:12])
    print(filename)
 
 
 
    data = np.loadtxt(files)
    alt = data[:,0]  
    diff = data[:,1]
    NumberOfLine = len(alt)
    tabAlt = np.append(tabAlt,alt)
    tabDiff = np.append(tabDiff, diff)
 
 
    a = datetime.strptime(filename, "%y%m%d")
 
    for i in range(len(alt)):
        tabTime = np.append(tabTime,a)
        #tabTime = np.append(tabTime,filename)
 
    #tabTime, tabAlt = np.meshgrid(tabTime, tabAlt)
    #plt.pcolormesh(X)
#tabTime, tabAlt = np.meshgrid(tabTime, tabAlt)
plt.xticks(rotation = 'vertical')   
plt.scatter(tabTime, tabAlt, marker='+', c=tabDiff)    
plt.colorbar()
plt.show()

Concernant la date, en effet il s'agit ici d'une chaîne de caractère extraite à partir du nom de fichier, donc l'espacement entre les tracés est le même. J'ai utilisé datetime, mais j'ai juste un petit soucis au niveau des limites de l'axe temps :

Pièce jointe 596408

Le soucis est le zoom ? Dans ce cas il faut regarder du côté de xlim() :
Code:

1 2 3 4 plt.xlim(( datetime.strptime('22-04-2021', '%d-%m-%Y'), datetime.strptime('26-04-2021', '%d-%m-%Y'), ))

24/04/2021, 10h39
GlmGil

1 pièce(s) jointe(s)
Merci Julien, xlim a bien fonctionné, mais l'affichage des valeurs de 'x' ne semble pas fixe :

Voici les dates que j'ai :

08 08 05

08 08 13

08 08 21

08 08 27

Mais le graphique ne m'affiche pas les dates exactes ( il créé un décalage avec les dates) :

Pièce jointe 596501

J'ai utilisé ce code là :
Code:

1 2 3 4 5 6 plt.xlim(( datetime.strptime("2008-08-04", "%Y-%m-%d"), datetime.strptime("2008-08-28", "%Y-%m-%d"), ))
26/04/2021, 10h22
Julien N
Salut,

Au vu des dates et du graphique, je cois que tout est juste. Je ne vois pas de décalage. La première barre a l'air située au 4 aout, et la dernière au 27. En zoomant à la main dans la fenêtre matplotlib cela devrait être plus visible. Le rendu ne correspond peut-être pas à ce que tu souhaites cela dit, mais le problème semble être autre. Est-ce que ce ne serait pas les petits "traits" de la barre des abscisses (les "ticks") qui ne sont pas alignés avec les points ? Si c'est cela, alors il y a plusieurs moyens de faire ça :
- Forcer à mettre des "ticks" tous les jours
- Placer uniquement les "ticks" pour les dates où l'on a des données
Ju
27/04/2021, 10h44
GlmGil

Salut Julien,

Oui mon problème est bien lié au fait que les ticks ne sont pas alignés avec mes points.
Je ne peux pas mettre les ticks à la main car il y a plusieurs milliers de fichiers au total à traiter, donc il faudrait que j'automatise ça. Je regarde comment mettre les ticks uniquement pour les jours où il y a des données.

Merci!