Discussion :

[pizzsto]

Bonjour,

J'ai un nuage de points qui représentent un os. Ce nuage de points est sous forme de matrice dont la première colonne reprend les coordonnées X la deuxième les Y et la troisième les Z (d'un axe orthonormé).

Je voudrais générer une surface passant par tous les points de la matrice.

J'ai déjà essayé les fonctions meshgrid puis surf mais sans résultat (je crois que c'est dû au faite que mon vecteur des Z est un vecteur colonne).

Avez-vous une solution à mon problème?

merci.

[Jerome Briot]

Est-ce que ce nuage de points représente une surface ouverte (une partie de l'os) ou bien une surface fermée (l'os en entier) ?

[pizzsto]

le nuage reprend l'entièreté de l'os (surface fermée).
Maintenant il est possible de ne prendre qu'une demi matrice de points pour n'avoir que la moitié de l'os.

[Jerome Briot]

A ma connaissance, il n'y a pas de code MATLAB tout fait.

Après, tout dépend du type de surface que tu souhaites. Maillage ? Surface spline ? Alpha shape ? ...

Par simple curiosité, de quel os s'agit-il ?

[pizzsto]

il s'agit du métacarpe du pouce.

Ce que je désirais c'est, d'avoir l'équation de la surface qui interpolerait au mieux tout les points du nuage de points.
Je m'étais dit qu'avec la fonction surf(x,y,z) ca pouvait se faire mais pour eux le z utilisé dans la fonction est une fonction de x et y alors que dans mon cas j'ai deja un vecteur pour z (qui n'est pas une matrice de dimension demandée).

[paradize3]

Comme dut, je ne suis pas vraiment sur de ce que tu veux... (maillage, ... ?)

Si ton but est de de ploter la surface tu peux essayer ceci:

Code :

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

1
2
3
4
5
6
7
8
9
 
% Plot the test surface
hf = figure('renderer','zbuffer','Tag','UTGR');
surface(X,Y,Z, RGB, ...
    'FaceColor', 'texturemap', 'EdgeColor', 'none', 'CDataMapping', 'direct', ...
    'LineStyle', 'none','tag','smoothmap_surfaceplot')
title('mon image')
axis image
axis vis3d

Ceci dit, je n'ai jamais essaye de le faire avec des surfaces fermees.

Bonne chance,

Greg

[pizzsto]

En fait, pour recapituler, j'ai un nuage de points qui represente un os et j'aurais voulu faire passer une surface par tous les points dans le but de créer un volume qui represente en réalité l'os considéré.

Mais matlab peut pas le faire, enfin je crois. La raison est la suivante: il y a moyen d'interpoler une fonction par matlab qui fait passer la fonction par tous les points du nuage et puis matlab crée la surface a partir de la fonction. mais ici c'etait pas possible car pour des valeur de x,y données il y a au moins deux coordonnées z et donc pas possible de créer une fonction.

J'ai trouvé une solution, je coupe mon os en deux pour avoir en tous points une seule coordonnée z pour chaque doublet (x,y). puis j'utilise ceci:
[XI,YI]=meshgrid(xv,yv);
ZI=griddata(x,y,z,XI,YI,'cubic');
figure
surf(XI,YI,ZI)

(pour la signification des variables tels que XI,...allez voir l'help de matlab c'est super bien fait )


mais bon dans cette solution vous avez l'os coupé en deux ne tiens qu'a vous de reprendre des points sur la surface interpolée et de les mettre dans une nouvelle matrice et de regénérer un volume avec une matrice concaténée de matrices réalisées sur les deux parties d'os.

Mais comment reprend ton des points sur la surface interpolée??????? ca je vous le demande ????

[FR119492]

Salut !

Avant de chercher une fonction toute faite dans MatLab ou ailleurs, tu devrais commencer par décider quel algorithme tu vas utiliser. Comme Dut, je pense d'ailleurs que tu ne trouveras pas de fonction toute faite. Si j'étais confronté à ton problème, j'essaierais avec un spline "tricubique" et je procéderais comme suit:

1. Je commencerais par chercher une paire de sous-programmes pour les splines bicubiques (j'ai ça en fortran, mais c'est dans une bibliothèque propriétaire, donc je ne peux pas te l'envoyer, mais le problème est suffisamment classique pour que tu puisses trouver quelque-chose ailleurs).
2. Je décortiquerais le code pour bien comprendre ce qu'il fait.
3. Je le réécrirais pour l'étendre de 2 à 3 dimensions.

Enfin, ce que je ne ferais pas, mais que tu peux faire si tu es obligé d'utiliser MatLab ou si tu y tiens absolument, c'est d'écrire toi-même le toolbox MatLab en traduisant simplement les sous-programmes que tu as développés en Fortran.

Jean-Marc Blanc