Discussion :

[kana83]

Bonjour à tous
J'aimerais avoir vos suggestions sur le problème ci-dessous.

Actuellement sur un projet d'imagerie numérique de fin d'études, nous cherchons à délimiter la zone de surintensité commune à 2 noyaux (image_origine), zone issue du chevauchement de deux noyaux lors de l'acquisition.

Nous arrivons tout à fait à délimiter la ligne médiane entre 2 noyaux de cellules (voir ligne blanche image_resultat) grâce à une méthode basée sur les watershed.
Par contre nous avons du mal à délimiter précisément la zone de surintensité.

Nous avons essayé en seuillant par hystéresis avec l'utilisation de différents seuils bas déterminés par les valeurs des pixels de la ligne de séparation des noyaux (résultat des watershed):
minimum, médian ou moyenne des points de la ligne médiane et les résultats sont moyens.
Le minimum donne de mauvais résultats (sélectionne la totalité des noyaux), le médian est plutôt satisfaisant mais pas optimal et la moyenne ne donne pas de très bon résultats.

Si vous avez d'autres idées, elles sont les bienvenues,soit pour améliorer le seuillage, soit carrément une autre méthode plus efficace.

Merci

[pseudocode]

Hum. C'est pas évident.

Je tenterai une croissance de région en utilisant comme points de départ les pixels qui sont sur la ligne de séparation.

Le critère d'homogénéité va sans doute être difficile à règler. On peut commencer simplement par la valeur de l'intensité (avec un seuil bas).

[kana83]

Le problème est que si nous baissons le seuil, le risque est de trop faire croitre ta région et donc de prendre toute l'image.

[graindecafé]

une croissance contrôlée de région ne permettrait elle pas de limiter le seuillage à une zone proche de celle en surbrillance ?

càd prendre le seuil médian comme seuil de départ et de baisser progressivement ce seuil jusqu'à observer une brusque augmentation du nombre de points pris en compte dans le seuillage.

[souviron34]

moi personnellement je tenterais autre chose :

soit :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
trouver les 2 meilleures ellipses reflétant les 2 coques (la forme obtenue)
calculer mathématiquement l'intersection

soit

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
suivre le contour pour trouver le point d'inflexion d'un coté
trouver le point d'inflexion opposé
calculer l'intensité sur cette ligne
faire une croissance de région à partir de là avec ce paramètre

Une recherche rapide sur la méthode 1 donne :

http://forums.futura-sciences.com/thread12024.html

ou

http://ogigi.polsl.pl/biuletyny/zeszyt_16/z16cz2_1.pdf

ou

http://mathforum.org/library/drmath/view/66877.html

ça me semblerait la méthode la plus adaptée..

[pseudocode]

J'avais aussi envisagée le 1ère méthode donnée par Souviron qui donnera sans doute la solution la plus robuste. Après ça dépend des types de cellules que vous allez rencontrer (courbure/convexité).

Les autres méthodes proposées sont des variations de la croissance de régions.

Pour reprendre l'idée de Souviron, on peut envisager une approche plus statistique:

On se déplace le long de la médiane (step by step) et on analyse les valeurs situées sur la droite perpendiculaire à la médiane. Normalement on devrait avoir des distribution plus ou moins "sympathiques" sur chaque perpendiculaire. A partir de là on devrait pouvoir calculer les points de l'enveloppe de l'intersection. En théorie...