Discussion :

[scott88]

Bonjour,

Je dispose d'un profil de faisceau laser (abscisse = distance, de 0 à xx mm - ordonnée = intensité). Ce profil est à révolution de symétrie autour de l'axe des ordonnées.
Je souhaite générer une carte "3D" à partir du profil, en le faisant "tourner" autour de l'axe des ordonnées.
Est-ce possible? Si oui comment ? Peut-on passer par un espace en coordonnées polaire.

Je précise que je ne possède que des points du profil (pas d'expression analytique, lecture des données à partir d'un fichier texte).

En vous remerciant par avance,

[membreComplexe12]

Je pense qu'un petit schema de ce que tu as et de ce que tu veux permettras de nous faire comprendre plus facilement se que tu cherches (perso je n'ai pas vraiment compris)

[scott88]

Ok,

avec la PJ, j'espère que ce sera 1 peu + clair :-).
Les graphes sont donnés à titre d'illustration.

Sur le premier on a le profil du faisceau (ou plutôt demi profil). C'est ce graphe que j'appelle (abusivement peut être) "graphe 2D" (axes X et Z). Ici, le profil est gaussien. => c'est les données que je dispose.
Sur le deuxième graphe, il y le profil "3D", représentant la gaussienne en 3D. Elle est à symétrie de révolution autour de l'axe Z. => c'est ce que je souhaite obtenir.

Mon problème est de savoir s'il est possible de réaliser le 2ème graphe à partir du profil "2D" en faisant une rotation de ce dernier autour de l'axe Z, puisque je considère que mes profils "3D" sont à symétrie de révolution autour de Z.

Merci de votre/vos réponse(s),

[Kangourou]

salut,

un truc du genre :

Code :

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

1
2
3
4
5
6
7
8
9
 
lx = 1:50;
lz = exp(-lx.^2/400);
theta = linspace(0, 2*pi, 65);
x = lx(:)*cos(theta);
y = lx(:) * sin(theta);
z = lz(:) * ones(size(theta));
figure;
surf(x, y, z);

?

[scott88]

Exactement, c'est bien ce que je cherchais. L'astuce était donc bien de passer en coordonnées polaire.

Merci bcp !!!

Question subsidiaire : peut-on revenir, une fois le profil 3D obtenu, dans l'espace cartésien. Cela afin de pouvoir mieux manipuler le profil (là en coord. polaire c'est pas trop pratique). La fonction pol2cart peut-elle servir (là, je n'ai pas encore réussi à l'utiliser correctement).

Merci d'avance !

[scott88]

Au final, j'y suis arrivé via la fonction "griddata". Voici le code (basé sur celui de
Kangourou)

Code :

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

1
2
3
4
5
6
7
8
ti=[-max(lx):1:max(lx)];
[XI,YI] = meshgrid(ti,ti);
 
ZI = griddata(x,y,z,XI,YI);
ZI(isnan(ZI))=0;
 
figure;
mesh(XI,YI,ZI);

[Jerome Briot]

Code :

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

1
2
ti=[-max(lx):1:max(lx)];
[XI,YI] = meshgrid(ti,ti);

s'écrit plus simplement

Code :

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

1
2
ti=-max(lx):max(lx);
[XI,YI] = meshgrid(ti);

Sinon, je pense que INTERP3 conviendrait mieux que GRIDDATA, non ?