Discussion :

[loranne]

Bonsoir,
J'ai besoin pour un projet de dessiner une rose des vents à partir de données météorologiques de vent
J'utilise :

Code :

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

1
2
figure
rose(force,direction)

Et j'obtiens :
Pièce jointe 206178
Mon problème est que je ne sais pas utiliser les légendes et la configuration des axes,
J'aimerais que le 0 (représentant le Nord) soit en haut puis pouvoir indiquer les points cardinaux, et ensuite colorier le centre des triangles et changer la couleur en fonction de l'intensité avec une échelle
J'ai cherché sur l'aide de matlab mais je ne trouve pas ce que je voudrais, et je ne peux pas utiliser les fonctions déjà faites comme 'Wind_Rose' car je ne comprend pas tout...
Merci de votre aide !

[Jerome Briot]

J'aimerais que le 0 (représentant le Nord) soit en haut puis pouvoir indiquer les points cardinaux

Tu peux t'inspirer des solutions données dans cette discussion : Graphe en polaire décalé

Pour le reste, peux-tu nous donner un fichier .mat (mis dans une archive zip ou rar) contenant les variables force et direction ?

[loranne]

Merci beaucoup j'ai essayé avec la 2ème façon du lien que vous m'avez donné (un peu plus simple ^^) et ça fonctionne très bien, mon 0 se trouve en haut de ma figure ! Par contre, je n'ai pas su trouver la ligne qui me permettrait d'indiquer 'Nord' devant ou à la place du 0° ou 'Sud' à la place de 180°

Je vous met donc en pièce jointes les données que j'utilise, c'est une moyenne pour chaque demi-heure de la force et de la direction du vent en Décembre afin de caractériser le comportement du vent sur une journée
Je vous remercie de votre aide !

[Jerome Briot]

je n'ai pas su trouver la ligne qui me permettrait d'indiquer 'Nord' devant ou à la place du 0° ou 'Sud' à la place de 180°

Par 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
load forcedirection
 
figure
 
rose(forcedirection(:,1), forcedirection(:,2))
 
set(gca, 'CameraPosition', get(gca, 'CameraPosition').*[1 1 -1])
rotate(findall(gca),[0 0 1],-90)
 
t = findall(0,'type', 'text');
 
set(t(1),'string', 'Nord')
set(t(2),'string', 'Sud')
set(t(7),'string', 'Ouest')
set(t(8),'string', 'Est')
 
delete(t([3:6 9:12]))

[Jerome Briot]

colorier le centre des triangles et changer la couleur en fonction de l'intensité avec une échelle

Sur ce genre de graphique, la couleur ne représente pas l'intensité mais le nombre d'éléments comptés (comme avec un histogramme).

Bref, voici comment ajouter de la couleur :

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
28
29
30
31
32
33
34
35
36
load forcedirection
 
figure
forcedirection(:,2) = forcedirection(:,2)*pi/180;
 
[tout,rout] = rose(forcedirection(:,2), 100);
h = polar(tout,rout);
 
x = get(h, 'XData');
x = reshape(x, 4, []);
y = get(h, 'YData');
y = reshape(y, 4, []);
rout = reshape(rout, 4, []);
delete(h)
 
hold on
colormap(parula(max(rout(:))))
k = 0;
for n = 1:size(x,2)
    if any(x(:,n))
        k = k+1;
        h(k) = patch('Vertices', [x(:,n) y(:,n)], 'Faces', [1 2 3 4], 'FaceVertexCData', rout(2,n), 'FaceColor', 'flat');
    end
end
 
set(gca, 'CameraPosition', get(gca, 'CameraPosition').*[1 1 -1])
rotate(findall(gca),[0 0 1],-90)
 
t = findall(0,'Type', 'text');
 
set(t(1),'string', 'Nord')
set(t(2),'string', 'Sud')
set(t(7),'string', 'Ouest')
set(t(8),'string', 'Est')
 
delete(t([3:6 9:12]))

Le code "transforme" les objets Line en objet Patch

[loranne]

Merci beaucoup !!
Je l'ai adapté !
En fait je m'étais trompée au début quand je traçais

Code :

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

1
2
figure
rose(force,direction)

en fait avec la fonction rose, on ne peut pas indiquer la force du vent sur la graphique c'est cela ?
Donc mon premier graphique ne représente pas les bons résultats de la direction du vent ? mais la force je ne sais pas trop par quel moyen c'est ça ?
La direction moyenne vraie du vent est celle du dernier graphique que vous m'avez envoyé ?
Je n'ai pas la fonction 'parula', si je remplace par 'spring' par exemple c'est bon ? C'est juste un changement de couleur c'est ça ?
Je vous remercie !!

[Jerome Briot]

Je n'ai pas la fonction 'parula', si je remplace par 'spring' par exemple c'est bon ? C'est juste un changement de couleur c'est ça ?

Oui, tu peux aussi utiliser la palette de couleur jet :

Code :

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

colormap(jet(max(rout(:))))

en fait avec la fonction rose, on ne peut pas indiquer la force du vent sur la graphique c'est cela ?
Donc mon premier graphique ne représente pas les bons résultats de la direction du vent ? mais la force je ne sais pas trop par quel moyen c'est ça ?
La direction moyenne vraie du vent est celle du dernier graphique que vous m'avez envoyé ?

La fonction rose trace un histogramme en coordonnée polaire.
Sur la précédente image que j'ai postée, on peut simplement dire qu'un maximum de vents avaient pour direction celle du secteur jaune.
Le graphique permet de visualiser correctement la direction des vents mais il n'y a aucune infirmation sur leurs intensités.