Discussion :

[MariettaUARK]

Bonjour,

J'ai une figure dont les coordonnées des abscisses vont de 0 à 360 et je voudrais la centrer sur 0 et non sur 180.
Le souci c'est que je n'arrive pas à décaler la matrice de points en même temps que les valeurs de l'axe des abscisses et donc les valeurs de la carte ne sont plus en adéquation avec ses coordonnées.
Ce n'est pas forcément clair, voilà le code :

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
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
 
data = 'goossens_datamap.csv' 
 
# importation de la matrice de données
contour_data = pd.read_csv(data, skiprows=0) # importation du fichier de données au format 3 colonnes x,y,z
contour_data.astype(float) 
Mat_d = contour_data.pivot_table(index='x', columns='y', values='z').T.values # conversion colonne => matrice 'numpy.ndarray'
 
# Création des valeurs de longitude et latitude
x = np.linspace(0, 360, np.shape(Mat_d)[1])
y = np.linspace(-90, 90, np.shape(Mat_d)[0])
 
#Affichage
implot = plt.imshow(Mat_d, extent=[x[0],x[-1],y[0],y[-1]], origin={'lower', 'lower'}, cmap='jet', interpolation='blackman')
plt.xlabel('Longitude [°]')
plt.ylabel('Latitude [°]')
 
plt.show()

Est-il possible de décaler toute la carte de 180x ? (cela revient concrètement à inverser les deux moitiés droite et gauche de la carte, ce qui peut être un début de solution : bidouiller les données de départ, pas très élégant mais incontestablement efficace)

Vous noterez que je crée un peu 'artificiellement' les x et y alors qu'il sont aussi contenus dans les 1ère et 2ème colonnes du fichier d'origine.

Merci.
- M.

[robinechuca]

Pouvez-vous joindre le fichier csv, de façon à ce que l'on puisse tester votre code pour mieux comprendre ce qu'il fait.
Les gens forts n'en n'ont pas forcément besoin mais personnellement, ça me serai bien utile, ne serait-ce que pour comprendre votre question.

[MariettaUARK]

Voilà le lien du fichier (je ne peux pas le déposer directement ici, le site n'accepte pas les csv, et même en .txt ou .dat de toute façon le fichier est trop lourd)

https://owen.link/05e6977431avc64dk9...063b6#bZpJWD7H

[danielhagnoul]

Et comme 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
data = 'Tests/en_cours/goossens_datamap.csv'
 
# importation de la matrice de données
# importation du fichier de données au format 3 colonnes x,y,z
contour_data = pd.read_csv(data, skiprows=0)
contour_data.astype(float)
 
# conversion colonne => matrice 'numpy.ndarray'
Mat_d = contour_data.pivot_table(index='x', columns='y', values='z').T.values
 
# Création des valeurs de longitude et latitude
x = np.linspace(-180, 180, np.shape(Mat_d)[1])
y = np.linspace(-90, 90, np.shape(Mat_d)[0])
 
# Affichage
plt.imshow(Mat_d, extent=[x[0], x[-1], y[0], y[-1]],
           origin={'lower', 'lower'}, cmap='jet', interpolation='blackman')
 
plt.xlabel('Longitude [°]')
plt.ylabel('Latitude [°]')
 
plt.xlim(-50, 100)
plt.ylim(-40, 60)
 
plt.show()

[robinechuca]

Je ne maîtrise pas bien les matrices numpy, mais en passant par les listes je m'en sort.
Je suis d'accord avec la solution de danielhagnoul mais si tu veux garder un axe allant de 0 a 360, en inversant comme tu dis les 2 moitiers.
Il y a un moyen qui est celui là:

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
 
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
 
data = 'goossens_datamap.csv' 
 
# importation de la matrice de données
contour_data = pd.read_csv(data, skiprows=0) # importation du fichier de données au format 3 colonnes x,y,z
contour_data.astype(float) 
Mat_d = contour_data.pivot_table(index='x', columns='y', values='z').T.values # conversion colonne => matrice 'numpy.ndarray'
 
# Création des valeurs de longitude et latitude
x = np.linspace(0, 360, np.shape(Mat_d)[1])
y = np.linspace(-90, 90, np.shape(Mat_d)[0])
 
#Affichage
Mat_d = [list(ligne) for ligne in Mat_d]  # je transforme l'array en liste
Mat_d = [ligne[len(ligne)//2:] + ligne[:len(ligne)//2] for ligne in Mat_d] # ensuite je bidouille les liste pour inverser les moitiers
implot = plt.imshow(Mat_d, extent=[x[0], x[-1], y[0], y[-1]], origin={'lower', 'lower'}, cmap='jet', interpolation='blackman')
plt.xlabel('Longitude [°]')
plt.ylabel('Latitude [°]')
 
plt.show()

[Beginner.]

Salut,

Waw, comment vous avez fait pour télécharger le fichier ?
En ce qui me concerne le site se fiche de moi...

J'ai cliqué puis il me dit :

Félicitations !

Vos fichiers ont bien été téléchargés, vous les retrouverez dans votre dossier de téléchargement habituel.

Mais j'ai rien !!!

[robinechuca]

Waw, comment vous avez fait pour télécharger le fichier ?

Je n'ai personnellement aucun souci (firefox linux).
Voila un lien weetransfère si ça peu aider: https://we.tl/t-0Kax7UM7YK

[Beginner.]

Salut,
Merci cette fois c'est bon...

[MariettaUARK]

Merci à tout le monde pour votre aide.

@danielhagnoul, ton code fait un zoom sur une région de la carte (celle avec le spot) : en fait quand j'ai dit que je voulais centrer ma carte, ce n'était pas de la centrer sur le spot mais sur le méridien zéro.

@robinechuca, oui ton code permet bien d'inverser les deux moitiés de ma carte donc je vais faire comme ça en attendant merci. Ce qui m'embête c'est qu'ici j'ai juste de la chance d'avoir besoin d'un décalage de 180° ce qui correspond à l'inversion des deux moitiés de la carte mais pour des décalages de degrés différents j'aurai à nouveau le souci... Mais à chaque jour suffit sa peine =) Je prends ta solution.

@Begginer, désolée pour ta difficulté à télécharger les data, je ne sais pas d'où ça venait. Merci robinechuca d'avoir remis un lien.

[danielhagnoul]

Si le zoom ne vous intéresse pas, il suffit de le supprimer en enlevant ces deux lignes :

Code :

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

1
2
plt.xlim(-50, 100)
plt.ylim(-40, 60)

[MariettaUARK]

Si le zoom ne vous intéresse pas, il suffit de le supprimer en enlevant ces deux lignes :

Code :

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

1
2
plt.xlim(-50, 100)
plt.ylim(-40, 60)

Oui mais hormis ça, ici seules les valeurs de l'axe des longitudes ont été changées, et la figure reste à la même place, c'était bien mon problème.