Discussion :

[Jragk]

Bonjour, je souhaiterai tracer un graphique présentant une élévation d'une surface entre deux pas de temps fixés. J'ai les position x, y, z des points pour chaque pas de temps (lus dans un fichier).

Je souhaite tracer un plot "vu de dessus", avec x, y comme axes, et la variation d'altitude en couleur, avec un changement d'échelle de couleur à 0 (certains points montent, d'autres descendent, et je souhaiterai clairement les différencier) Mon programme semble marcher, mais je n'arrive pas à afficher correctement la colorbar avec le changement d'échelle de couleur à 0, et je ne comprend pas pourquoi.
Quelqu'un saurait-il m'aider pour afficher une colorbar correcte ?
Voici mon 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
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
data = load('....surfrun2.dat1');
x1 = data(:,3);
y1 = data(:,4);
z1 = data(:,5);
x1 = x1(z1>10);  %top surface
y1 = y1(z1>10);
z1 = z1(z1>10);
 
 
data = load('....surfrun2.dat101');
x2 = data(:,3);
y2 = data(:,4);
z2 = data(:,5);
x2 = x2(z2>10);  %top surface
y2 = y2(z2>10);
z2 = z2(z2>10);
 
 
dhdt = (z2-z1)/10;
 
 n = 32; % Augmenter n augmente le nombre de teintes
 cmap = [(winter(n));autumn(n)]; 
 colormap(cmap)
 caxis([-3 0.5])
 colorbar;
cm = colormap;
figure; hold on;
m_dhdt = min(dhdt)
M_dhdt = max(dhdt)
for k=1:length(dhdt)
    col(1) = interp1([m_dhdt:(M_dhdt-m_dhdt)/63:M_dhdt],cm(:,1), dhdt(k));
    col(2) = interp1([m_dhdt:(M_dhdt-m_dhdt)/63:M_dhdt],cm(:,2), dhdt(k));
    col(3) = interp1([m_dhdt:(M_dhdt-m_dhdt)/63:M_dhdt],cm(:,3), dhdt(k));
    plot(x1(k)/1000, y1(k)/1000, '.', 'color', col, 'MarkerFaceColor', col)
end
%colorbar;
box on;
        xlabel('x (km)');
        ylabel('y (km)');
        xlim([500 700]);
        ylim([0 50]);

[Jerome Briot]

Est-ce toujours en rapport avec ta première question ?

=> Segmentation de la colorbar avec pcolor

[Jragk]

Correction : mon échelle de couleur ne change pas à 0, mais au milieu de mon interval, je me trompe. Bref, en fait ça ne marche pas du tout...

[Jragk]

C'est lié, mais c'est un autre problème.
Ici, je n'ai pas les résultats sous forme matricielle, donc je ne peux pas utiliser de pcolor...
Cette méthode me semblait la plus simple mais je m'embourbe là, je crois.

[Jragk]

Ok, j'ai de l'amélioration : en modifiant le nombre d'élément que je mets dans mes 'col', j'arrive à ploter le bon graph avec la colorbar à la bonne échelle, et sur le bon graph.

C'est pas très beau comme code, mais ça semble fonctionner. merci quand même pour le passage Dut !

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
dhdt = (z2-z1)/50;
 
 cmap = [(winter(3*32));autumn(0.5*32)]; 
 colormap(cmap)
 
cm = colormap;
figure; hold on;
m_dhdt = min(dhdt)
M_dhdt = max(dhdt)
for k=1:length(dhdt)
    col(1) = interp1([m_dhdt:(M_dhdt-m_dhdt)/111:M_dhdt],cm(:,1), dhdt(k));
    col(2) = interp1([m_dhdt:(M_dhdt-m_dhdt)/111:M_dhdt],cm(:,2), dhdt(k));
    col(3) = interp1([m_dhdt:(M_dhdt-m_dhdt)/111:M_dhdt],cm(:,3), dhdt(k));
    plot(x1(k)/1000, y1(k)/1000, '.', 'color', col, 'MarkerFaceColor', col)
end
 
 
cmap = [(winter(3*32));autumn(0.5*32)]; 
colormap(cmap)
caxis([-3 0.5])
colorbar;
 
box on;
        xlabel('x (km)');
        ylabel('y (km)');
        xlim([500 700]);
        ylim([0 50]);