Discussion :

[SCIEBOU]

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)':



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.

[Jerome Briot]

Utilise les fonctions griddata ou TriScatteredInterp

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

[SCIEBOU]

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 :

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

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');

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.

[Jerome Briot]

Pourrais-tu nous fournir un jeu de données ?

[SCIEBOU]

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.

[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

[Jerome Briot]

[edit] Grillé par Fabien

Essaie avec la contribution Surface Fitting using gridfit

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
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

[SCIEBOU]

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:



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

Merci beaucoup pour votre aide.

[Gooby]

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!

[le fab]

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

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

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

hummm ... soit t'es très riche, soit tu devrais rester plus discret

[le fab]

Code :

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

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)

[SCIEBOU]

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

Bon WE.