Discussion :

[csdty]

Bonjour à tous,
Voila mon problème : J'ai des points (de coordonnées x,y,z ) repartis de manière éparses (nuage de points) dans [0,1]^3. Je peux les afficher avec plot3d; mais j'aimerais tracer une surface fermée "au plus près" autour de l'ensemble qu'il forme.
J'ai essayé avec convhull mais cela me crée un volume convexe passant par un maximum de points . Le problème c'est que cette surface est très "globale". Il faudra au moins que le volume ainsi créé soit connexe par arcs.
Je ne suis pas très bon en MATLAB, mais je pense que ce problème est assez compliqué.
Avez-vous une idée ?
Merci d'avance

P.S : n’hésitez pas à demander si vous voulez que je vous envoie les coordonnées des points

[le fab]

hello
matlab propose griddata et TriScatteredInterp pour ça
moi j'utilise plutôt gridfit qui donne des meilleurs résultats à mes yeux, surtout quand les points sont pas si reprochés que ça

fab

[csdty]

Merci pour ta reponse !
C'est tres impressionant cette fonction ( elle marche mieux que griddata) ;
mais elle ne fait pas ce que je veux : je me suis mal exprimee...je vais essayer de mieux m'exprimer en expliquant mon probleme en 2D :
mes vecteur de coordonnes (x,y) sont repartis comme on peut le voir sur l'image . J'aimerais creer un contour ( surface en 3d ) qui fasse le tour du nuage de points . Comme le fait convhull (trait rouge sur l'image); mais je voudrai que ca soit plus "pres" des donnees donc pas de facon convexe , mais plutot convexe par arcs ou autre chose
Vous voyez ce que je veux dire ?


D'avance merci
P.S : j'ai mis les donnees pour que vous puissez les reprensenter en 3D

[le fab]

ah ca je te confirme qu'elle fonctionne mieux que griddata !!

sinon pour ton probleme, je le comprend, mais je sais pas faire comme ca
j'y réfléchi et si j'ai une idée ...

[Jerome Briot]

Tu peux essayer l'algorithme Alpha Shapes en téléchargeant par exemple cette contribution : Alpha shapes by Jonas Lundgren

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.mat
 
[V,S] = alphavol(X1, 150);
 
figure('doublebuffer','on')
 
plot3(X1(:,1),X1(:,2),X1(:,3),'r*')
hold on
trisurf(S.bnd,X1(:,1),X1(:,2),X1(:,3), ...
    'FaceColor','b', ...
    'FaceAlpha',1)
 
axis vis3d
grid on

[csdty]

C'est super impressionnant comme fonction !!
Ca a l'air de bien marcher !!

Merci beaucoup (a tout les deux) !