Discussion :

[conanran]

Bonjour,

J’ai codé un petit code qui permet de tracer des cercles à différents diamètres.
Mon objectif consiste à dégrader la couleur de la surface selon le diamètre du disque et d’afficher un color bar à coté avec la dégradation de couleur et la valeur de diamètre.

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
VThetaDeg = 0:1:360;
VTheta = VThetaDeg *pi / 180;
X_centre=1:10:150;
Y_centre=20*ones(1,length(X_centre));
Rmax=10;
Rmin=3;
Rayon=zeros(1,length(X_centre));
x_depart=X_centre(1);
X_fin=X_centre(end);
for n=1:length(X_centre)
X=X_centre(n);
Rayon(n)=(((X-x_depart)*(Rmax-Rmin))/(X_fin-x_depart))+Rmin;
end
 
for k=1:length(X_centre)
XCercle = X_centre(k) + Rayon(k) * cos(VTheta);
YCercle = Y_centre(k) + Rayon(k) * sin(VTheta);
fill(XCercle,YCercle,'r')
hold on
end

Merci

[Jerome Briot]

Tu veux plutôt parler de "dégradé", non ?

Une première idée à partir de ces données simplifiées (je présume) :

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
VThetaDeg = 0:1:360;
VTheta = VThetaDeg *pi / 180;
X_centre = 1:10:150;
Y_centre = 20*ones(1,length(X_centre));
Rmax = 10;
Rmin = 3;
Rayon = zeros(1,length(X_centre));
x_depart = X_centre(1);
X_fin = X_centre(end);
for n = 1:length(X_centre)
    X=X_centre(n);
    Rayon(n) = (((X-x_depart)*(Rmax-Rmin))/(X_fin-x_depart))+Rmin;
end
 
[a,b,~] = unique(Rayon);
 
cmap = hsv(max(b));
set(gcf, 'colormap', cmap)
 
hold on
 
for k = 1:length(X_centre)
    XCercle = X_centre(k) + Rayon(k) * cos(VTheta);
    YCercle = Y_centre(k) + Rayon(k) * sin(VTheta);
    h = fill(XCercle, YCercle, 'r');
    set(h,'facecolor', cmap(b(k),:))
end
axis equal
 
c = colorbar;
yl = get(c, 'ylim');
yt = linspace(yl(1), yl(2), 2*numel(a)+1);
set(c, 'ytick', yt(2:2:end-1), 'yticklabel', num2str(a(:)))

Donne nous des données réelles pour pouvoir te fournir une solution plus précise.

[conanran]

Merci bien, c'est génial

[conanran]

Bonjour,
En effet, quant j'ajoute une autre cercle avec un diamètre plus petit ou deux cercles de deux diamètres différents, elle ne donne pas le résultat souhaité.
Merci

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
clear all, clc
%%
VThetaDeg = 0:1:360;
VTheta = VThetaDeg *pi / 180;
X_centre = 1:10:151;
Y_centre = 20*ones(1,length(X_centre));
Rmax = 10;
Rmin = 3;
Rayon = zeros(1,length(X_centre));
x_depart = X_centre(1);
X_fin = X_centre(end-1);
for n = 1:length(X_centre)-1
    X=X_centre(n);
    Rayon(n) = (((X-x_depart)*(Rmax-Rmin))/(X_fin-x_depart))+Rmin;
end
 Rayon(length(X_centre))=2; %Rayon ajouté 
 
 
[a,b,~] = unique(Rayon);
 
cmap = parula(max(b));
set(gcf, 'colormap', cmap)
 
hold on
 
for k = 1:length(X_centre)
    XCercle = X_centre(k) + Rayon(k) * cos(VTheta);
    YCercle = Y_centre(k) + Rayon(k) * sin(VTheta);
    h = fill(XCercle, YCercle, 'r');
    set(h,'facecolor', cmap(b(k),:))
end
axis equal
 c=colorbar
% c=colorbar; colormap(jet)
yl = get(c, 'ylim');
yt = linspace(yl(1), yl(2), 2*numel(a)+1);
set(c, 'ytick', yt(2:2:end-1), 'yticklabel', num2str(a(:)))

Je pense à cause de la fonction unique, par exemple quant je fais Rmax=Rmin=10, normalement, je trouve des cercle avec la même couleur.
Merci pour votre aide?

[Jerome Briot]

Il faut utiliser le 3ème argument de sortie de unique et non pas le 2ème :

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
clear all, clc
%%
VThetaDeg = 0:1:360;
VTheta = VThetaDeg *pi / 180;
X_centre = 1:10:151;
Y_centre = 20*ones(1,length(X_centre));
Rmax = 10;
Rmin = 3;
Rayon = zeros(1,length(X_centre));
x_depart = X_centre(1);
X_fin = X_centre(end-1);
for n = 1:length(X_centre)-1
    X=X_centre(n);
    Rayon(n) = (((X-x_depart)*(Rmax-Rmin))/(X_fin-x_depart))+Rmin;
end
 Rayon(length(X_centre))=2; %Rayon ajouté 
 
 
[a,~,c] = unique(Rayon);
 
cmap = parula(max(c));
set(gcf, 'colormap', cmap)
 
hold on
 
for k = 1:length(X_centre)
    XCercle = X_centre(k) + Rayon(k) * cos(VTheta);
    YCercle = Y_centre(k) + Rayon(k) * sin(VTheta);
    h = fill(XCercle, YCercle, 'r');
    set(h,'facecolor', cmap(c(k),:))
end
axis equal
 c=colorbar;
% c=colorbar; colormap(jet)
yl = get(c, 'ylim');
yt = linspace(yl(1), yl(2), 2*numel(a)+1);
set(c, 'ytick', yt(2:2:end-1), 'yticklabel', num2str(a(:)))

[conanran]

Merci pour votre réponse,
Dans le cas ou j'ai des centaines des cercles, je ne veux pas afficher toutes les valeurs de diamètres sur le colorbar. Je veux mettre le min et le max et quelques valeurs tout au long du colorbar selon un pas déterminé. Comment je peux le faire ?
Merci

[Jerome Briot]

Modifie la fin du code comme ceci :

Code :

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

1
2
3
4
5
6
yt = linspace(yl(1), yl(2), 2*numel(a)+1);
yt = yt(2:2:end-1);
 
idx = [1 numel(yt)];
 
set(c, 'ytick', yt(idx), 'yticklabel', num2str(a(idx).'))

À toi d'adapter les valeurs dans idx selon ton besoin

[conanran]

Merci Jerome Briot,
par contre il n'affiche pas des valeurs au milieu. je veux ajouter des valeurs au niveau des points rouge par exemple.

[Jerome Briot]

C'est bien ce que j'ai écris dans la dernière phrase de mon précédent message.

Si tu veux un point sur deux :

Code :

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

idx = 1:2:numel(yt);