Discussion :

[Seli-vince]

Bonsoir,

J'ai à disposition :
- Un vecteur r de 1x50 représentant les rayons
- Un vecteur alpha de 1x50 représentant les angles
- Un vecteur z de 1x50 représentant la hauteur pour chaque couple (r,alpha). (zi est la valeur de (ri, alpha i))

J'aurais souhaité représenter en 3D une surface à partir de mes 50 points (r,alpha,z) dans l'espace.

Mon problème étant que surf et équivalent fonctionnent en coordonnées cartésiennes !

En me renseignant sur internet, j'ai vu deux possibilitées, toutes 2 vaines :
- La première me faisait transformer tout mes points du plan (r, theta) en points d'un plan (x,y) via pol2cart. Mais ensuite en passant à surf, il aurait fallu une matrice z avec une valeur pour chaque couple (x,y) possible incluant (xi,yj) alors que mon vecteur z ne propose qu'une valeur pour les couples (xi,yi).
- La deuxième méthode commencait par un meshgrid, mais nécéssitait ensuite une fonction de r et theta pour trouver les Z (permettant ainsi d'avoir les z pour tout les couples comme expliqué ci-dessus).

Je ne voit donc pas comment faire, si l'un de vous aurait une idée je suis preneur !

Merci d'avance, et n'hésitez pas à me le faire savoir si je ne suis pas clair !

Vincent.

[Jerome Briot]

Pourrais-tu nous fournir les données dans un fichier mat (compressé en zip ou rar) ?

[Seli-vince]

Hum un peu dur le fichier original comme y a des choses importées depuis mon ordinateur.

J'ai recréé ci-dessous un modèle qui reproduit environ ce que j'ai à disposition à la fin.

Cylindrique.zip

Ce que je cherche à représenter est une demie-sphère (globe coupé à l'équateur avec dans le plan de l'equateur r la distance depuis le centre et alpha l'angle par rapport à une direction arbitraire, et z la hauteur).

N'hésitez pas à me poser d'autres questions si ce n'est pas encore clair !

[Jerome Briot]

Hum un peu dur le fichier original comme y a des choses importées depuis mon ordinateur.

Ne peux-tu pas faire simplement ceci ?

Code :

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

save('data.mat', 'r', 'alpha', 'z')

Et compresser le fichier mat en zip ou rar ?

[Seli-vince]

Je n'avais jamais fais cela et ne connaissait pas désolé !

Voilà donc mes 3 colonnes, la seule différence avec l'énoncé de mon dernier message est que le globe n'est pas coupé à l'équateur mais à un parrallèle autre, et que je n'attend pas une sphère parfaite mais biscornue !

Les 2 exemples étant vraiment similaires, une réponse sur l'un ou l'autre me débloquera de toute manière alors n'hésitez pas à prendre celui simplifié si cela peut vous aider !

data.zip

[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
load data
 
r = reshape(r, 73, []);
alpha = reshape(alpha, 73, []);
z = reshape(z, 73, []);
 
figure
plot3(r(:),alpha(:),z(:), 'r*');
hold on
surf(r,alpha,z)
 
axis equal

[Seli-vince]

Merci pour votre réponse.

Votre façon de faire utilise (r,alpha,z) en tant que coordonées cartésiennes (x,y,z), alors que mes données sont en coordonnées polaire (r,angle,z).

Comme expliqué ci-dessus, mon graphe final doit avoir la forme d'un globe coupé à un parallèle (ou dans l'exemple simplifié coupé à l'équateur : une demi-sphère).

Donc les r indiques la distance par rapport au centre du repère, alpha l'angle fait par rapport à l'axe des x dans le plan (x,y) et z la hauteur du point.

[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
load data
 
r = reshape(r, 73, []);
alpha = reshape(alpha, 73, []);
z = reshape(z, 73, []);
 
x = r.*cosd(alpha);
y = r.*sind(alpha);
 
figure
surf(x,y,z)
 
axis equal vis3d
grid on