Fitting de points en 3D

Bonjour,

j'ai cherché dans le forum pendant un moment mais je n'ai rien trouvé, ou alors j'ai mal cherché et je m'en excuse d'avance.
J'ai créé un programme qui va être utilisé pour de la métrologie. Le but est de récupérer des données métrologiques d'une pièce et de les comparer avec les données nominales. J'ai pour commencé utilisé la fonction icp (Iterative Closest Points) qui permettra de trouver la retation et la translation aux points récupéré afin qu'ils soient le plus proches des données nominales. Pour le moment j'aimerais créer des droites et cercles approximés à partir des points. Pour du 2D j'ai trouvé facilement les fonctions à utiliser, mais pour du 3D j'ai un peu de mal.
Connaissez-vous des fonctions qui permettent cela ? Ou dois-je passer par plus d'étape pour le faire ?

Merci d'avance !

salut

quelles fonctions as tu identifier en 2D ?
celles auquelles je pense marchent en 3D tout aussi bien

Fabien

Salut,

Je ne connais rien à rien à cette problématique ==> google "icp matlab" ==> je suis tombé sur cette contribution dans matlab central... cela peut t'être utile, peut-être ?

Cordialement,

Salut

la fonction icp je l'ai déjà implémenté dans Matlab (je vous ai mis mes fichier .m en pièces jointes).

Citation:

---

quelles fonctions as tu identifier en 2D ?
celles auquelles je pense marchent en 3D tout aussi bien

---

En 2D je me souviens que nous avions utilisé en cours une fonction polyfit il me semble afin de retrouver des équations de droites, paraboles ... mais nous avons toujours travailler qu'en 2D. Est celle là à laquelle tu penses ?

Merci

plus personne pour m'aider ? :(

Sinon je continue toujours à chercher. Je sais qu'il existe une multitude de fonction, le problème c'est que je ne sais pas comment m'y prendre pour affiner mes recherches...

mieux que polyfit (et plus générique), l'opérateur \
sinon lsqcurvfit, fmincon ...

Bonjour,

J'ai deux remarques :

1) Je suis désolé, mais tu peux m'expliquer pourquoi tu as besoin de faire du fitting de courbe ? Tu as bien tes deux séries de points, non ?

2) Si tu veux faire du fitting 3D, c'est pareil qu'en 2D. Deux possibilités :
- Soit tu suis ce qui est indiqué au point 4) de la méthode donnée par lefab
- Soit tu appliques ceci : tu choisis l'équation de ta courbe (A) : {Z = f(X,Y)} ou bien (B) : {X = f(t) | Y = g(t) | Z = h(t)}, puis tu prends autant de points significatifs de ta courbe que tu as d'inconnues, puis tu résous les équations pour trouver tes inconnues. (B) est beaucoup plus générique que (A) car tu peux en principe faire des boules avec (B), mais pas avec (A) par exemple (ou alors il faut appliquer 2 fois (A)... bref). A titre informatif, (B) forme l'ensemble des courbes paramétrées 3D. Tu peux d'ailleurs utiliser lsqcurvefit 3 fois pour directement trouver la solution pour la forme (B), de manière numérique (il te suffit juste d'écrire tes équations sous forme d'équations paramétriques, en indiquant tes inconnues dans un vecteur "x").

Cordialement,

Merci pour vos réponses,

en fait je viens de me rendre compte que, en effet, faire du fitting ne me servira à rien. Maintenant que j'ai réussi à recaler les points au plus proches des points nominaux, ce qu'il me reste à faire c'est trouver la distance qu'il y a entre le cercle nominal et chaque point. Je dois fouiller de ce coté là maintenant, j'espère que ça va bien se faire. En tout cas merci pour vos réponses ! Ça fait un moment que je n'ai pas touché à Matlab, et comme nous n'avons jamais reçu de cours là dessus j'ai un peu de mal à trouver ce que je cherche.

Encore merci !

Bonjour,

j'ai un petit peu avancé dans mon projet, j'ai essayé avec un cylindre maintenant, mais maintenant c'est la fonction ICP qui n'a pas l'air de trouver la translation et rotation nécessaire. Est-ce que j'utilise mal cette fonction ? Dans mon code je génère un cylindre nominal, j'en tire des points dits expérimentaux qui subissent une translation/rotation, et normalement après la fonction icp est sensée retrouver la translation/rotation que j'ai appliqués...

Si vous avez un petit moment pour trouver mon erreur, ou m'expliquer en quoi ma démarche n'est pas bonne, je suis ouvert à toutes remarques.

Merci d'avance !