Discussion :

[Winounet]

Bonjour à tous,

Je suis encore novice sur Matlab mais étant limité pour de la représentation graphique sous Excel je n'ai pas trouvé de meilleur moyen.
Je cherche à obtenir un graphique 3D avec une base en 2D type circulaire et en Z des valeurs positives et négatives. (Idéalement c'est un peu le logo Matlab que je cherche à obtenir dans la forme).
J'ai eu des valeurs en x et en y qui me représentent un "quasi-cercle" (enfin là n'est pas le problème) et des valeurs associés aux pixels à l'intérieur de ce cercle.
L'objectif est de voir les variations de contraste en positif ou négatif.
J'ai pensé à revenir à l'équation paramétrique et tracer le cercle mais j'ai peur de ne pas pouvoir tracer ensuite les Z.
Donc pour l'instant j'ai beau essayer toutes sortes de possibilités, les seules représentations que j'arrive à obtenir donnent des résultats horribles.

Si quelqu'un a une piste ou a déjà eu ce genre de programme à faire..
Merci,

[Jerome Briot]

J'ai eu des valeurs en x et en y qui me représentent un "quasi-cercle" (enfin là n'est pas le problème) et des valeurs associés aux pixels à l'intérieur de ce cercle.
L'objectif est de voir les variations de contraste en positif ou négatif.

Donne nous un jeu de valeurs, stocké dans un fichier mat (mis dans une archive zip ou rar)

[Winounet]

Merci de la réponse rapide,

Voilà en pièces jointes mes 3 variables.
Pour le moment je suis en train de tracer un cercle en plot, et j'essaye de pouvoir cumuler avec un autre graphique représentant les Z mais ça ne me donne pas de surface :/

[Jerome Briot]

Une première idée :

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
minx = min(VarName1);
maxx = max(VarName1);
 
miny = min(VarName2);
maxy = max(VarName2);
 
cx = minx + (maxx-minx)/2;
cy = miny + (maxy-miny)/2;
 
[th,r] = cart2pol(VarName1-cx, VarName2-cy);
 
thi = linspace(0, 2*pi, 40);
ri = linspace(0, max(r(:)), 25);
 
[thi, ri] = meshgrid(thi, ri);
 
[xi,yi] = pol2cart(thi,ri);
 
xi = xi + cx;
yi = yi + cy;
 
 
F = TriScatteredInterp(VarName1, VarName2, VarName3);
 
zi = F(xi, yi);
 
 
figure(1)
clf
plot3(VarName1, VarName2, VarName3, 'r+');
hold on
% plot3(xi, yi, zi, 'b*')
surf(xi, yi, zi, 'edgecolor', 'none')
 
axis vis3d square tight
grid

[Winounet]

Pour une première idée je pouvais pas espérer mieux x) c'est exactement ce genre de graphique que j'essayais d'obtenir ! Je connaissais pas le TriScatteredInterp et je n'avais pas trouvé quelque chose de semblable !

merci merci !

[Winounet]

Enfin je viens de voir que l'axe des Y est à l'envers enfin c'est pas grand chose car c'est pas le principal mais je ne vois pas comment inverser les valeurs ><

[Jerome Briot]

Ajoute ceci à la fin du code :

Code :

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

set(gca, 'YDir', 'reverse')

[Winounet]

Merci je venais de trouver