Image : ligne de partage des eaux : lpe (win32)

Soyez le premier à donner votre avis sur cette source.

Vue 18 873 fois - Téléchargée 2 134 fois

Description

La LPE est un algorithme en traitement de l'image qui permet de segmenter l'image en région homogène. C'est une technique très puissante simulant l'immersion d'un relief topologique par les eaux.
L'algorithme implémenté est une réponse au principe de S.Beucher: LPE sans biais et avec marqueur.
L'immersion commence par les marqueurs et à chaque fois que deux bassins versant se rencontre, une "digue" est irriguée séparant les deux bassins. L'ensemble des digues correspond à la ligne de partage des eaux.
Les marqueurs permettent de contrôler la segmentation puisqu'il y aura autant de région que de marqueur.

Une application de la LPE: Une image de cellules. Cette technique permet de séparer les cellules qui se touchent et offre une très bonne segmentation.

Ce programme met en jeu des techniques avancées en imagerie :
-Seuillage automatique
-Détection des contours
-Bouchage des trous
-Carte des distances
-Construction des marqueurs
-LPE

2 scriptes exemples ont été implémentés, il suffit de cliquer sur les icones LPE dans la barre d'outil.

Source / Exemple :


//****************************************************************************************
//MORARD Vincent
//vincent.morard@cpe.fr
//pistol.petesampras.free.fr
//****************************************************************************************

#include "CImage.h"

void Message(HWND hList,HWND hDlg,char *Mes,int Buf)
{
	SendMessage(hList,LB_ADDSTRING, 0, (LONG)Mes);
	SendMessage(hDlg,WM_UPDATE,Buf,0);
	Sleep(500);
}
//************************************************************************************
//Thread cells, Ce thread permet de lister toute une série de traitement plus ou moins
//long. A partir de l'image de cellules, ce processus permet d'effectuer une segmentation
//afin d'isoler chaque cellule. La difficulté pour segmenter correctement cette image
//c'est de séparer les cellules qui se touchent. Pour cela, on utilise l'algorithme de 
//la ligne de partage des eaux appliqué sur l'image des distances.
//Cette image des distances permet de créer l'image topologique nécessaire à la LPE et
//de localiser le centre des cellules.
//************************************************************************************
DWORD WINAPI cells(LPVOID lpParam )
{
	HWND hDlg,hOK,hwndList;
	
	hDlg=(HWND) lpParam;
	hwndList=GetDlgItem(hDlg, IDC_LIST1);
	hOK=GetDlgItem(hDlg, IDOK);
	
	bool Noyau[9];
	for(int i=0;i<9;i++) Noyau[i]=1;

	int t = GetTickCount();  //initialisation du temps

	SendMessage(hwndList, LB_ADDSTRING, 0, (LONG)"Début du traitement");

	if(Img[0].LoadBmp("./cell.bmp")==0)return 0;
	Message(hwndList,hDlg,"Ouverture de l'image : Fait",0);
	
	
	if(Img[0].ContrasteAuto(GRAY,&Img[1])==0)return 0;
	Message(hwndList,hDlg,"Contraste automatique : Fait",1);
	

	if(Img[1].Threashold(Img[1].thrFuzzy(),&Img[2])==0)return 0;
	Message(hwndList,hDlg,"Seuillage automatique : Fait",2);

	if(Img[2].NOT(&Img[3])==0)return 0;
	Message(hwndList,hDlg,"Inversion : Fait",3);
	
	
	if(Img[3].morCloseHole(3,&Img[4])==0)return 0;
	Message(hwndList,hDlg,"Bouchage des trous : Fait",4);

	if(Img[4].NOT(&Img[3])==0)return 0;
	if(Img[3].LaplaceI(GRAY,&Img[4])==0)return 0;
	if(Img[4].Threashold(10,&Img[4])==0)return 0;
	Message(hwndList,hDlg,"Détection des contours :Fait",4);
	

	if(Img[4].DistanceMorphologique(3,TRUE,&Img[5])==0)return 0;
	Message(hwndList,hDlg,"Distance morphologique : Fait",5);

	
	if(Img[5].ContrasteAuto(GRAY,&Img[5])==0)return 0;
	if(Img[5].Inversion(&Img[6])==0)return 0;
	Message(hwndList,hDlg,"Contraste auto + inversion : Fait",6);

	if(Img[6].ErosionNG(Noyau,9,0,&Img[7])==0)return 0;
	Message(hwndList,hDlg,"Erosion en niveau de gris: Fait",7);

	if(Img[7].FindMinimum(&Img[8])==0)return 0;
	Message(hwndList,hDlg,"Détermination des marqueur : Fait",8);

	if(Img[6].LPE(9,10,4)==0)return 0;
	Message(hwndList,hDlg,"LPE : Fait",9);

	if(Img[9].lpeFrontiere(&Img[4],&Img[10])==0)return 0;
	if(Img[10].Superposition(&Img[1])==0)return 0;
	Message(hwndList,hDlg,"Superposition : Fait",10);
	
	
	SendMessage(hwndList, LB_ADDSTRING, 0, (LONG)"Fin du traitement, temps (ms) ");
	t = GetTickCount() - t;  
	char Buf[20];
	strcpy(Buf,itoa(t,Buf,10));
	SendMessage(hwndList, LB_ADDSTRING, 0, (LONG)Buf);

	EnableWindow(hOK,TRUE);
	return 1;

}

Conclusion :


Si vous trouvez des bugs, merci de me les signaler.

Codes Sources

A voir également

Ajouter un commentaire Commentaires
j'ai besoin Watersheed en java ou bien en builder C++ !!
Messages postés
1054
Date d'inscription
samedi 2 octobre 2004
Statut
Membre
Dernière intervention
9 juillet 2013
7
Il me semble que oui, c'est la fonction watersheed
A+
Messages postés
1
Date d'inscription
mercredi 4 mai 2011
Statut
Membre
Dernière intervention
4 mai 2011

Ya pas du LPE sur matlab ?!!
Messages postés
1054
Date d'inscription
samedi 2 octobre 2004
Statut
Membre
Dernière intervention
9 juillet 2013
7
Regarde sur mon site Internet, il y a une version plus récente.(htp://ImAnalyse.free.fr) C'est compilé avec VS 2003, 2005 ou 2008.
Les sources sont à cette adresse:

http://ImAnalyse.free.fr/Telechargement/_Sources2.3.zip
A+
Messages postés
5
Date d'inscription
dimanche 10 février 2008
Statut
Membre
Dernière intervention
23 janvier 2010

Bonjour,
Merci pur votre attention
mais ça ne marche pas.
sur quelle version du visual C++ vous travailler?
a+
Afficher les 12 commentaires

Vous n'êtes pas encore membre ?

inscrivez-vous, c'est gratuit et ça prend moins d'une minute !

Les membres obtiennent plus de réponses que les utilisateurs anonymes.

Le fait d'être membre vous permet d'avoir un suivi détaillé de vos demandes et codes sources.

Le fait d'être membre vous permet d'avoir des options supplémentaires.