Discussion :

[Pelagos]

Bonjour a tous,

J'ai un probleme qui relève plus des maths que de la prog, mais voyons voir tout de même. Je suis relativement débutant en programmation, et je code en F77.

Je cherche a faire un programme d'interpolation lineaire. Comprenez, pour un ensemble de formes geometriques données, (dans mon cas, ce sont des espèces d'ellipsoïde), definies par un ensemble de coordonnées (x,y) je dois calculer les coordonnées de nouveau points avec quelques contraintes (c'est la que ça se complique).

-Le nombre de point doit toujours etre le même quelque soit la forme. C'est a dire, si lance mon programme sur 6 formes, elles doivent toutes avoir a la fin le même nombre K de points interpolés. Le must serait que ce nombre de point soit fixé par l'utilisateur, mais ça, ça change pas grand chose au final.

-Les points doivent tous etre equidistants, y compris le dernier avec le premier. (ce qui veut donc dire que la distance entre chaque point ne sera pas la même d'une forme a l'autre, au final.)

Mon probleme est relativement simple. Prenez un cercle, et essayez de placer 4 points equidistants dessus. vous verrez tres vite que la distance entre ces quatres point n'est absolument pas egale au perimètre du cercle/4. De plus, mes formes geometriques ne sont pas regulières, donc on ne peux pas calculer un rapport entre le perimetre initial et celui interpolé, par exemple.

La seule possibilité que j'ai trouvé jusqu'a present c'est de calculer un premier run avec la distance entre chaque point egale au perimetre divisé par k, et d'affiner la distance par une boucle, cad de faire tourner mon algorythme en reduisant la distance a chaque boucle jusqu'a avoir la bonne.

Mais lorsqu'on a environ 900 formes, avec 300 points a calculer a chaque fois, c'est beaucoup, beaucoup, beaucoup trop long (et souvent hasardeux)

Quelqu'un a une idée ou une suggestion ?
Merci d'avance.

PS: Cela fait deux semaines que je m'arrache les cheveux sur ce p***** d'algorithme, alors même les idées les plus saugrenus sont vraiment, vraiment bienvenues...

[genteur slayer]

tu peux essayer de regarder les travaux de Harlow et Welch (1965) sur la méthode Marker and Cells (Méthode MAC) ou plus adapté à ton problème la méthode des marqueurs de front de Daly et Pracht (1967) ils font de la reconstruction d'interface mais en plaçant des marqueurs (particules sans masse) sur un contour où il faut les re-parsemer au fur et à mesure de l'évolution du contour et idéalement de manière équidistante... il doit y avoir les algo de fourni dans leurs articles: avec un nombre de point suffisant, on interpole n'importe quel contour!

[Pelagos]

Merci beaucoup. Je me mets tout de suite en quete.