Discussion :

[ERCO503]

Bonjour,

J'affiche des données avec "pcolor" où à chaque (theta,rho) est associé une valeur, ce qui donne donc un disque coloré.

Cependant, les axes sont cartésiens, ce qui n'est pas très utile. J'aimerais pouvoir avoir des axes polaires afin de pouvoir estimer rapidement et visuellement le (theta,rho) ayant une valeur élevée.

Merci de votre aide,

Éric

[Jerome Briot]

Un exemple :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
th = linspace(0,2*pi,30);
r = linspace(0,10,20);
c = rand(20,30);
 
[th,r] = meshgrid(th,r);
 
[x,y] = pol2cart(th,r);
 
figure
 
pcolor(x,y,c)
axis equal tight

[ERCO503]

Merci pour le réponse, mais malheureusement c'est toujours sur cela que je tombe alors que ce n'est pas là mon problème.

Comme j'ai dit plus haut, j'affiche des données polaires (au même titre que l'exemple que vous avez donné), mais les axes sont cartésiens (comme dans votre exemple).

Si je vous demande d'approximer visuellement les coordonnées polaires (theta,rho) du rectangle plutôt orange entre les rectangles cyans près de (x,y) = (5,-7), on en a aucune idée puisque les axes sont cartésiens et non polaires. Ce pourquoi j'aimerais que les axes soient polaires .

Merci,

Éric

[ERCO503]

Je suis en train d'y parvenir en tentant de combiner "polar" et "pcolor", mais si quelqu'un l'a déjà fait proprement, ce me serait grandement utile.

Voici une image jointe qui illustre mieux ce que je veux (il y a encore quelques failles dans ma combinaison, comme le chiffre "10" qui est le seul rayon affiché [et le colorbar trop près est seulement dû à "axis tight", sinon c'est bon]).



Merci,

Éric

[Jerome Briot]

Une idée à améliorer :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
th = linspace(0,2*pi,30);
r = linspace(0,10,20);
c = rand(20,30);
 
[th,r] = meshgrid(th,r);
 
[x,y] = pol2cart(th,r/max(r(:)));
 
figure
 
polar(nan,nan)
hold on
pcolor(x,y,c)
axis equal tight

C'est une base qui peut être améliorée

Tu peux aussi t'inspirer de cette contribution Graphique de type Radar (toile d'araignée) pour générer la grille

[ERCO503]

Ah, merci, simple et efficace.

Savez-vous pourquoi "r" doit être ramené à 1 ? Parce que si j'enlève la division par max(r(:)), ça ne fonctionne plus.

EDIT : Ça fonctionne si j'utilise polar(nan,max(r(:)));, mais il doit y avoir un problème d'affichage puisque je perds les données des rayons sauf le dernier...

Merci,

Éric

[Jerome Briot]

Savez-vous pourquoi "r" doit être ramené à 1 ? Parce que si j'enlève la division par max(r(:)), ça ne fonctionne plus.

Parce que les limites des axes de l'objet Axes sont [-1 1] comme on peut le voir avec ce code minimal :

Code :

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

1
2
3
polar(nan,nan)
 
axis on

[Invité]

Bonsoir,

Pour contourner cela, on peut tracer les données r et th et immédiatement supprimer le tracé du polar :

Code :

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

delete( polar(th,r) )

[ERCO503]

Merci, il ne reste que mon problème de rayons non affichés sauf le dernier (comme sur l'image que j'avais jointe).

EDIT : Bon, voilà :

Code :

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

1
2
3
4
5
6
7
8
9
10
[th,r] = meshgrid(hc{2}/180*pi,hc{1});
[x,y] = pol2cart(th,r);
figure;
polar(nan,max(r(:)));
hold all;
h = pcolor(x,y,hn);colorbar;
set(h,'EdgeColor','none');
set(h,'Layer','top');
axis equal;
hold off;

Merci, je crois que c'est résolu.

EDIT : En fait, non, ça m'envoie une erreur, mais ça fonctionne tout de même... Étrange ?!

Il faut que j'enlève "axis equal" :

Code :

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

1
2
3
4
5
6
7
8
[th,r] = meshgrid(hc{2}/180*pi,hc{1});
[x,y] = pol2cart(th,r);
figure;
polar(nan,max(r(:)));
hold all;
h = pcolor(x,y,hn);colorbar;
set(h,'EdgeColor','none');
hold off;

[Invité]

La fonction pcolor trace chaque couleur entre les points (x,y) désignés. Il y a donc numel(x)-1 éléments.
Pour une taille m x n on peut faire :

Code :

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

1
2
3
4
5
6
7
th = linspace(0,2*pi,n+1);
r = linspace(0,10,m+1);
c = rand(m,n);

% On double les dernières ligne et colonne
c = [c ; c(end,:)];
c = [c c(:,end)];