Discussion :

[hananemto2007]

bonjour

vous avez résolu le problème de spécifier une couleur pour la valeur 0 dans cette discussion : Spécifier une couleur à une valeur

Est ce que vous pouvez m'aider,le problème que j'ai c'est que je veux spécifier pour chaque intervalle une couleur.
Par exemple les valeurs entre 0 et 1 doivent etre en rouge
entre 1 et 5 en bleu

merci d'avance

[Jerome Briot]

Tu pourrais nous dire quels objets tu traces (Line, Surface, Patch) et quelles fonctions tu utilises ?

Montre nous du code

[hananemto2007]

Voila mon code:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
[lon,lat]=meshgrid([ -35.00:0.25:35.00],[ 44.50:-0.25:10.00]);
m_proj('miller','lon',[-35 15],'lat',[10 44]);
temp=load('C:\Documents and Settings\hanane\Bureau\travail_lundi\max_min_12m\mat_max_HS_01.txt');
hold on 
m_pcolor(lon,lat,temp); shading interp;
[c,h]=m_contour(lon,lat,temp,'k');
clabel(c,h,'fontsize',8);
m_coast('patch','k','edgecolor','k','linewi',2);
m_grid('box','fancy','linestyle','none');
ylabel('Latitude (°)','FontSize',11,'FontWeight','demi');
xlabel('Longitude (°)','FontSize',11,'FontWeight','demi');
colorbar

[Jerome Briot]

Est-il possible de joindre un fichier mat_max_HS_01.txt (sous format zip) à ton message ?

[hananemto2007]

Je suis désolé pour le retard.
voila le fichier

merci d'avance

[Jerome Briot]

Comme ceci :

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
42
43
44
[lon,lat]=meshgrid([ -35.00:0.25:35.00],[ 44.50:-0.25:10.00]);
m_proj('miller','lon',[-35 15],'lat',[10 44]);
temp=load('mat_max_HS_01.txt');
 
% Les intervalles
interv=[0 1 5 8];
 
maxtemp=max(temp(:));
% Construction de la nouvelle palette de couleurs (Colormap)
% Choisir N suffisamment grand
N=128;
d=round(diff([interv maxtemp])*N/maxtemp);
% Construction de d développée :
% for n=1:numel(interv)-1
%     d(n)=round((interv(n+1)-interv(n))*N/max(temp(:)));
% end
% d(n+1)=round((max(temp(:))-interv(n+1))*N/max(temp(:)));
 
col=hsv(numel(d));
col=col(end:-1:1,:);
 
newcmap=[];
for n=1:numel(d)    
    newcmap=[newcmap ; repmat(col(n,:),d(n),1)];    
end
 
% Création de l'objet Figure
figure;
% Mise en place de la palette de couleur newcmap
colormap(newcmap)
 
h=m_pcolor(lon,lat,temp); 
shading interp;
 
colorbar
 
hold on 
[c,h]=m_contour(lon,lat,temp,interv,'k');
clabel(c,h,'fontsize',8);
 
m_coast('patch','k','edgecolor','k','linewi',2);
m_grid('box','fancy','linestyle','none');
ylabel('Latitude (°)','FontSize',11,'FontWeight','demi');
xlabel('Longitude (°)','FontSize',11,'FontWeight','demi');

[hananemto2007]

Merci bien pour votre aide.

J'ai testé plusieurs exemples, mais une fois je change interv=[0 1 5 8]
par interv=[0 0.1 0.5 1.25 2.5 4 6 9 14] il me donne des cartes bizar.

le même problème si je change le fichier temp.

merci d'avance

[Jerome Briot]

Le problème vient de la ligne suivante :

Code :

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

m_grid('box','fancy','linestyle','none');

En la mettant en commentaire (avec un % devant), la carte apparait normalement.

Autre chose, chacune des valeurs des intervalles doit être inférieure à la valeur maximale contenue dans temp. Dans l'exemple donné, la valeur 14 est supérieur au max de temp (issu du fichier mat_max_HS_01.txt) qui est 11.

De plus, le code ci-dessus assume que la valeur minimale de temp est 0.

[hananemto2007]

Merci beaucoup pour vorte reponse

ok c'est bon, mais j'ai plusieurs exemples de fichier temp, par exemple pour temp1 les valeurs sont comprisent entre 0 et 4, alors que pour temp2 les valeurs sont comprisent entre 0 et 1.25 et je veux que toute les cartes auront la même colorbar( 0 0.1 0.5 1.25 4 6 9 14) pour faire une comparaison.


Je suis désolé et merci encore une fois.

[Jerome Briot]

Comme ceci :

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
[lon,lat]=meshgrid([ -35.00:0.25:35.00],[ 44.50:-0.25:10.00]);
m_proj('miller','lon',[-35 15],'lat',[10 44]);
temp=load('mat_max_HS_01.txt');

% Les intervalles
interv=[0 0.1 0.5 1.25 2.5 4 6 9];

maxtemp=14;
% Construction de la nouvelle palette de couleurs (Colormap)
% Choisir N suffisamment grand
N=128;
d=round(diff([interv maxtemp])*N/maxtemp);

col=hsv(numel(d));
col=col(end:-1:1,:);

newcmap=[];
for n=1:numel(d)    
    newcmap=[newcmap ; repmat(col(n,:),d(n),1)];    
end

% Création de l'objet Figure
figure;
% Mise en place de la palette de couleur newcmap
colormap(newcmap)

h=m_pcolor(lon,lat,temp); 
shading interp;

hold on 
[c,h]=m_contour(lon,lat,temp,interv,'k');
clabel(c,h,'fontsize',8);
caxis([0 14])
colorbar

m_coast('patch','k','edgecolor','k','linewi',2);
% m_grid('box','fancy','linestyle','none');
ylabel('Latitude (°)','FontSize',11,'FontWeight','demi');
xlabel('Longitude (°)','FontSize',11,'FontWeight','demi');

Il faut aussi chercher un peu dans la documentation pour voir comment son gérer les palettes de couleur (Colormap) et plus généralement la couleur sous MATLAB