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

Télécharger le projet
Commenter

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.
Télécharger le projet

Codes Sources

A voir également

S'inscrire maintenant

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.

Du même auteur (Pistol_Pete)