Tracé d'une surface 3D

Version imprimable

28/08/2014, 16h23
SCIEBOU

1 pièce(s) jointe(s)

Tracé d'une surface 3D

Bonjour,

J'ai 3 vecteurs X <1x700 double>, Y <1x700 double> et Z <1x700 double> qui m'ont permis d'obtenir le tracé 3D suivant avec 'scatter3(X,Y,Z)':

Pièce jointe 154917

Auriez-vous une idée pour savoir comment je peux tracer une surface passant par tous ces points et faire varier la couleur de cette surface en fonction des valeurs de Z.

Merci beaucoup pour votre aide.
28/08/2014, 16h59
Jerome Briot

Utilise les fonctions griddata ou TriScatteredInterp

Tu trouveras des exemples en effectuant une recherche sur ce forum.

1 pièce(s) jointe(s)

Bonsoir,

Merci pour ton aide. J'ai essayé les deux méthodes griddata et TriScatteredInterp et ai réussi à obtenir le résultat suivant:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 % Tracé 3D: Méthode TriScatteredInterp F1 = TriScatteredInterp(X,Y,Z); [qx,qy] = meshgrid(X,Y); qz = F1(qx,qy); figure(1); mesh(qx,qy,qz); hold on; plot3(X,Y,Z,'o'); %Tracé 3D: Méthode griddata [xg,yg] = meshgrid(X,Y); zg = griddata(X,Y,Z,xg,yg); figure(2); mesh(xg,yg,zg); hold on; plot3(X,Y,Z,'o');
Pièce jointe 154940

Par contre, j'ai encore un souci. C'est qu'à la différence des exemples de tracés 3D sur Matlab, je n'ai pas les vecteurs X et Y répartis de manière régulière. En effet, parmi les 700 points que j'ai, il y a 100 points qui ont les mêmes coordonnées de x. En d'autres termes, pour un x donné, j'ai 100 points définis dans les coordonnées (Y,Z). Du coup, en obtenant la surface passant par tous les points que je souhaite, j'ai eu aussi la surface verticale contenue dans le plan (y,z) (voir figure suivante):

Questions:
1) Pouvez-vous m'aider à définir la fonction griddata (ou une autre) pour ne obtenir que la surface voulue?

2) De plus, avez-vous des idées pour avoir une bonne qualité de visualisation pour cette image?

Merci pour votre aide.

28/08/2014, 20h07
Jerome Briot

Pourrais-tu nous fournir un jeu de données ?
28/08/2014, 21h41
SCIEBOU

1 pièce(s) jointe(s)

Ci-joint mes trois vecteurs X, Y, Z en .mat dans un fichier zip.
J'espère que cela vous ira. Merci pour votre aide.
29/08/2014, 09h21
le fab

Salut

en effet griddata et triscatteredinterp sont un peu bourrin dans leur interpolations
cependant d'après ton graph, il semblerait que tu ne ne leur ai pas demander d'interpoler vraiement : il faudrait que tes qx, qy et xg, yg soient plus définis
là tu les construit à partir des vecteur X et Y de tes données : il faut donner à meshgrid des vecteurs avec plus de points
ensuite griddata et triscatteredinterp feront de l'interpolation

je te conseille aussi de regarder la fonction gridfit dispo sur matlab central
la surface obtenue ne passera pas forcement exactement sur les points contrairement au deux autres fonctions mais le résultat sera peut être plus proche de tes attentes

sinon pour ce type d'affichage mesh est très bien

Fabien

1 pièce(s) jointe(s)

[edit] Grillé par Fabien ;)

Essaie avec la contribution Surface Fitting using gridfit

Code:

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
load X
load Y
load Z
 
% Suppression des NaN
idx = any(isnan([X Y Z]),2);
X(idx) = []; Y(idx) = []; Z(idx) = [];
 
% Suppression des doublons
[~,idx] = unique([X Y Z], 'rows');
X = X(idx); Y = Y(idx); Z = Z(idx);
 
% Approximation de la surface avec gridfit
% http://www.mathworks.com/matlabcentral/fileexchange/8998-surface-fitting-using-gridfit
xi = linspace(min(X), max(X), 50);
yi = linspace(min(Y), max(Y), 50);
zi = gridfit(X, Y, Z, xi, yi);
 
figure
 
plot3(X,Y,Z,'r*')
hold on
surf(xi, yi, zi)
 
axis tight vis3d
grid on

29/08/2014, 12h02
SCIEBOU

1 pièce(s) jointe(s)

Bonjour,

Super! La surface montrée dans la figure est exactement ce que je souhaite obtenir.

Néanmoins, il me semble que mon Matlab (Matlab2012b version 32bits) n'a pas de fonction 'gridfit'. Je me souviens pourtant d'avoir installé toutes les toolbox.

Voici le message que j'ai eu:

Pièce jointe 154967

Comment est ce que je peux installer la fonction 'gridfit'?

Merci beaucoup pour votre aide.
29/08/2014, 12h16
Gooby

Citation:

Envoyé par Dut

Essaie avec la contribution Surface Fitting using gridfit

Par contribution, Dut entend une fonction qui a été proposé par un utilisateur de MATLAB à la communauté. Il faut donc que tu ailles sur le lien proposer par Dut, que tu télécharges ce(tte) fonction/script afin de pouvoir l'utiliser!
29/08/2014, 14h45
le fab

Citation:

Envoyé par SCIEBOU

Comment est ce que je peux installer la fonction 'gridfit'?

y a le lien dans mon msg précédent

Citation:

Envoyé par SCIEBOU

Je me souviens pourtant d'avoir installé toutes les toolbox.

hummm ... soit t'es très riche, soit tu devrais rester plus discret ;)
29/08/2014, 14h48
le fab
Citation:
Envoyé par Dut

Code:

1 2 3 ... xi = linspace(min(X), max(X), 50); yi = linspace(min(Y), max(Y), 50);
ca c'est le truc qui manquait à ton premier code et qui aurait améliorer le résultat via griddata
(mais pas au point de ce que retourne gridfit)
29/08/2014, 15h12
SCIEBOU

Génial!!! Merci beaucoup pour votre aide.
Je suis très content d'avoir une surface si bien dessinée!!!

Bon WE.