Algorithmes de la morphologie mathématique pour - Pastel - HAL

Algorithmes de la morphologie mathématique pour les architectures orientées fluxJaromír BRAMBORnous décrit la définition du skeleton algorithmique ngbGAlgoGen destiné pour le travail géodésique surle voisinage.Algorithme 5.9 : ngbGAlgoGen Skeleton algorithmique généralisé de travail géodésique sur levoisinage1 ngbGAlgoGen :: [Streamize α] → ( [ExtrNgb α], [SampFnc α])→ [NgbGOp α]2 → ( Ar ( I , I ) α, Ar ( I , I ) α) → Ar ( I , I ) α3 ngbGAlgoGen strms (extrs N ,extrs M ) ops (ar , msk) = array4 (bounds ar)5 (compute strms (extrs N ,extrs M ) ops ar msk)6 where7 compute [ ] [ ] [ ] ar msk = [ ]8 compute (s:ss) ( (e N :es N ),(e M ,es M )) (o:os) ar msk =9 ( (zip (s ar))10 $ (map2 o11 ( (map (e N ar ) ) $ (s ar))12 ( (map (e M msk))$ (s msk))13 )14 )15 +(compute ss (es N ,es M ) os ar msk)La fig. 5.13 nous illustre sur le diagramme du flux la structure et la fonction de ce skeleton dans le casspécial où nous ne fractionnons pas l’image aux différentes zones et où nous avons, par conséquent, sonfonctionnement décrit par une seule fonction de parcours de l’image, une seule fonction d’extraction desvoisins et une seule fonction de kernel opérant sur le voisinage et sur le masque (strms, extsN, extsP,ops sont les listes contenant un seul élément).ar(Arrayd’entrée)boundsstrm 0 extrN 0Msk(Arrayde masque)computestrm 0extrM 0NgbOp 0zip ++[]Création de l’arrayde sortiearrayArrayde sortieFIG. 5.13 : Fonctionnement du skeleton algorithmique ngbGAlgoGen de traitement du voisinage avec lemasque dans le cas spécial où nous ne fractionnons pas l’image aux différentes zones5.3.2 Itérations, fin de propagationCe skeleton est élémentaire et il ne décrit qu’une itération de la propagation géodésique. Pour l’utilisercorrectement, l’élément structurant employé doit obligatoirement être de la taille unité. La propagationgéodésique elle-même sera obtenue par l’incorporation de ce skeleton dans un skeleton ou unalgorithme gérant la propagation. Il existe, en effet, deux types de fonctions géodésiques :• fonctions géodésiques avec le nombre d’itérations donné,• fonctions géodésiques itérant jusqu’à l’idempotence.Pour les premières, le nombre d’itérations est connu d’avance est la construction de l’algorithme peuts’appuyer sur ce fait. Pour les deuxièmes, nous devons tester l’arrêt de la propagation (test d’idempotencede la fonction géodésique sur une image donnée). Ce test compare les valeurs de l’image d’entrée (avant116