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

Jaromír BRAMBOR7.3. SKELETON ALGORITHMIQUE DE LA PROPAGATION SIMD EN 4-VOISINAGE(a) Image du masque "*msk" pour lapropagation correspondant à la fonction f deréférence des nivellements(b) Image des valeurs à propager "*val"correspondant à la fonction m du marquer desnivellementsFIG. 7.9 : Les images d’entrée aux fonctions de sur- et sous-nivellements et l’exemple de leur contenu.Une itération de lambda-sous-nivellements est définie par la fonction lolevelSQR4 :lolevelSQR4 :: [Char] → I → I → Ar ( I , I ) (α,α) → Ar ( I , I ) (α,α)lolevelSQR4 axe pvecsz lmb ar = propAlgSQR4 axe pvecsz( λ (vmsk,vval) (emsk,eval) → maxSIMD emsk (minSIMD (vval+lmb) eval) )aroù le terme λ définit le kernel de la propagation directionnelle et le paramètre lmb définit la pente deslambda-sous-nivellements. En spécifiant sa valeur à 0, nous obtenons les sous-nivellements plats. Lasignification des autres paramètres de cette fonction est identique à celle des paramètres pour la fonctiondistance (cf. 7.3.4) : axe définit l’axe de stockage des données dans la mémoire, pvecsz est le nombredes éléments que nous pouvons traiter en même temps par les instructions SIMD de notre architecturemultimédia et ar est l’array d’entrée composé des tuples (masque, valeur).Par analogie, nous définissons une itération des lambda-sur-nivellements par la fonction hilevelSQR4 :hilevelSQR4 :: [Char] → I → I → Ar ( I , I ) (α,α)→ Ar ( I , I ) (α,α)hilevelSQR4 axe pvecsz lmb ar = propAlgSQR4 axe pvecsz( λ (vmsk,vval) (emsk,eval) → minSIMD emsk (maxSIMD (vval−lmb) eval) )aroù le seul changement par rapport à la fonction des lambda-sous-nivellements est dans la définition duterme λ définissant le kernel de la propagation directionnelle.7.3.5.4 Algorithme pour une itération des nivellementsComme nous l’avons déjà expliqué dans la section 7.3.5.2, une itération des nivellements (finaux)sera obtenue par la combinaison des résultats d’une itération des sur- et sous-nivellements (partiels)selon l’équation 7.3.La fonction levelSQR4 définit cette composition :levelSQR4 :: [Char] → I → I → Ar ( I , I ) (α,α)→ Ar ( I , I ) (α,α)levelSQR4 axe pvecsz lmb ar =(mkAr2DFromAr2DPVec axe)◦ listArray (bounds $ lo)map2 ( λ (lomsk,loval) (himsk,hival) → cndmoveSIMD (testSIMD lomsk ( ≥ ) loval ) loval hival )(elems◦ (mkAr2DPVec axe pvecsz)$lo)(elems◦ (mkAr2DPVec axe pvecsz)$hi )wherelo = (lolevelSQR4 axe pvecsz lmb)$arhi = (hilevelSQR4 axe pvecsz lmb)$ar157