Discussion :

[matt67]

Bonjour,

J'ai des données de déformations dx pour chaque point de coordonnées [x y]. Cela me donne donc un vecteur de données 3D [x, y, dx] dont je parviens à afficher un nuage de points mais pas de créer une surface passant par chaque point.

J'ai bien essayé meshgrid puis griddata ou TriScatteredInterp mais le résultat à la visualisation est très moche.

Ci-joint mes points de coordonnées.

Merci de toute aide que vous pourrez m'apporter sur ce sujet !

[Jerome Briot]

Essaie 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
15
16
17
18
load Test.mat
 
x = linspace(min(Test(:,1)),max(Test(:,1)),150);
y = linspace(min(Test(:,2)),max(Test(:,2)),150);
 
[xi,yi] = meshgrid(x,y);
 
zi = griddata(Test(:,1),Test(:,2),Test(:,3),xi,yi,'cubic');
 
figure
 
plot3(Test(:,1),Test(:,2),Test(:,3),'r*');
 
hold on
 
surface(xi,yi,zi,'edgecolor','none')
 
axis vis3d

[matt67]

Super merci !

[le fab]

salut

au cas ou, la fonction gridfit est très sympa aussi, et suivant ta répartition de points, peut donner de bien meilleurs résultats que griddata ou TriScatteredInterp

après je trouve TriScatteredInterp très pratique (et élégant) dans sa mise en oeuvre, du coup je le combine souvent à gridfit, meilleur à mon gout en lissage de surface