Création de Colorbar/Colormap

Version imprimable

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:

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]);

09/06/2011, 11h23
Jerome Briot

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

=> Segmentation de la colorbar avec pcolor
09/06/2011, 11h23
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... :mrgreen:
09/06/2011, 11h24
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.

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:

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]);