Discussion :

[cindy63]

Bonjour,
Je dois segmenter une image 3D contenant 3 régions. L'histogramme du volume peut être modélisé par un mélange de gaussiennes. Je cherche, donc, à estimer l'amplitude, la moyenne et l'écart-type de chaque gaussienne pour classifier les voxels.
L'algorithme EM et les moindres carrés non linéaires (méthode de Levenberg Marquardt) appliqués sur l'histogramme ne donnent pas de bons résultats même avec une initialisation proche de la solution optimale ou ne convergent pas.

Auriez-vous une idée pour déterminer les paramètres des gaussiennes?

Je travaille en C# et je ne peux utiliser de librairie extérieure.
J'ai mis en fichier joint un exemple d'histogramme. Les données sont bruitées même si cela n'est pas visible sur l'histogramme. Un lissage ne suffit pas à débruiter suffisamment l'histogramme pour que l'algorithme EM converge vers la solution optimale.

Merci de votre aide.