Mthode des Ensembles de Niveaux par Elments Finis
Méthode des Ensembles de Niveaux par Eléments Finis P 1 Jérôme Piovano Stage de DEA Sous l’encadrement de Théodore Papadopoulo INRIA – Sophia Antipolis Projet Odyssée
1. Introduction • Segmentation d’image • Trouver des regions d’images selon certaines caractéristiques • Ensembles de niveaux ou « Levels Sets » • • Modélise l’évolution d’une hypersurface à travers une fonction continue Application à la segmentation d’image • Éléments finis • Méthode d’approximation discrète de fonctions continues Implémenter la méthode des ensembles de niveaux à l’aide des éléments finis
Plan Définitions n n Ensemble de niveaux Éléments finis Modélisation n Équations d’évolution Discrétisation temporelle Discrétisation spatiale Algorithmique n n n Notion de bande Évolution de la bande Simplification des équations Résultats n n Évolution par courbure moyenne Évolution pour contour géodésiques Conclusion
Définitions Ensembles de niveaux ou Levels Sets Interface représentée par le niveau 0 d’une « fonction distance » Évolution de l’interface par l’intermédiaire de la fonction distance t + |r | = 0; n n n = 1 = Détection de contours géodésiques. Schémas par « bande » à bases de différences finies instabilité
Définitions Méthode des éléments finis Approximation discrète d’une fonction continue 1. 2. 3. Partitionnement de l’espace en éléments formant un maillage Calcul des valeurs de aux sommets du maillage Représentation de par interpolation linéaire de ses valeurs aux sommet
Plan Définitions n n Ensemble de niveaux Éléments finis Modélisation n Équations d’évolution Discrétisation temporelle Discrétisation spatiale Algorithmique n n n Notion de bande Évolution de la bande Simplification des équations Résultats n n Évolution par courbure moyenne Évolution pour contour géodésiques Conclusion
Modélisation Calcul de la fonction distance grâce aux éléments finis Soit u l’approximation de la fonction distance par élément finis u définie par 2 facteurs : 1. Espacement constant entre ses différents niveaux (rxu)2 - 1 = 0 2. Vitesse d’évolution ut - = 0 1. Calculs des valeurs de u sur les sommets du maillage en minimisant une énergie associée a ces 2 termes
Modélisation Discretisation temporelle Exprimer l’évolution de u sous forme discrète dans le temps. v = “Pas d’évolution” u(x, t + t) = u(x, t) + v(x, t) = u(x, t + t) - u(x, t) v(x, t) = t ut(x, t) (rxu)2 - 1 = 0 ut - = 0 ( rxu + rxv)2 - 1 = 0 v - t = 0
Modélisation Discretisation spatiale Discretisation de Galerkin : n Utilise des fonctions de « bases » comme des fonctions tests mesurant la déviation au voisinage du sommet auquel elles sont attachées ( rxu + rxv)2 - 1 = 0 v - t = 0 s ((ru + r v)2 - 1) i = 0 8 i 2 1 … n s (v - t ) i = 0 8 i 2 1 … n Résolution d’un système de 2 n équations à n inconnues qui est donc surdéterminé Résolution par moindres carrés
Modélisation Reformulation des équations On peut exprimer les équations précédente en fonction des valeurs aux sommets du maillage de u et v s ((ru + r v)2 - 1) i s (v - t ) i (u + v)TQi(u + v) - s i Pi v - ts i Les vitesses nécessitent le calcul d’une dérivée seconde théoreme de Green
Plan Définitions n n Ensemble de niveaux Éléments finis Modélisation n Équations d’évolution Discrétisation temporelle Discrétisation spatiale Algorithmique n n n Notion de bande Évolution de la bande Simplification des équations Résultats n n Évolution par courbure moyenne Évolution pour contour géodésiques Conclusion
Algorithmique Adaptation du problème à un maillage 2 D régulier
Algorithmique La fonction distance n’est pas calculé sur la totalité de l’espace, mais au voisinage du niveau 0 n n Ajout des éléments proches du niveau 0 Suppression des éléments éloignés du niveau 0
Algorithmique Dynamique d’évolution
Algorithmique Les equations d’evolutions peuvent se simplifier sur des maillages 2 D réguliers à des calculs par difference finies. n Sur un maillage de type: (u + v)TQ 0(u + v) 1/6((S 2 - S 3)2 + (S 3 - S 4)2 + (S 5 - S 6)2 + (S 6 - S 1)2) + 1/3((S 1 - S 0)2 + (S 2 - S 0)2 + (S 4 - S 0)2 + (S 5 - S 0)2 )
Algorithmique Avantage de la methode Les equations d’evolutions peuvent se simplifier sur des maillages 2 D réguliers à des calculs par difference finies. n Sur un maillage de type:
Fin
- Slides: 17