Problème d'initialisation du signe

**nadjib2007** · 22/07/2007, 22h50

Bonjour tous le monde;

Je fais une étude sur les méthodes de segmentation d'images; et pour cela il faut initialiser un contour (moi j'utilise un cercle comme une courbe initiale); et il faut initialiser aussi les différents signes:

si on considère la matrice contour mc[i][j]

à l'intérieur du cercle : mc[i][j]=-1

à l'extérieur du cercle : mc[i][j]=+1

sur le cercle ( c-a-d : les pixels appartenant aux contour): mc[i][j]=0

comment je peut écrire un programme sous Builder C++ qui me permet de faire cette opération

Merci d'avance.

**hiko-seijuro** · 22/07/2007, 22h55

euh j'ai pas bien compris tu veux simplement initialiser la matrice ?

**nadjib2007** · 22/07/2007, 23h21

oui c'est ca en mettant la condition si le pixel est à l'intérieur, éxterieur ou sur le cercle.

Merci

**befalimpertinent** · 23/07/2007, 09h39

Il y a surement plusieurs façon de procéder. En voici deux. Je suppose que ton cercle est représenté par son centre et son rayon en pixel.

Soit tu parcours tous tes pixels et tu calcules à chaque fois sa position par rapport au cercle.
Soit tu commences par trouver les pixels sur le contour du cercle, tu les initialises dans ta matrices à m[i][j]=0. Puis tu parcours tout tes pixels dans un l'ordre (de gauche à droite et de haut en bas par exemple) et tu regarde si le pixel courant est déjà initialisé à 0 dans ce cas tu incrémente un compteur de 1. Tous les pixels situés sur la même ligne seront des m[i][j]==-1 jusqu'à ce que ton compteur passe à 2.
Je te passes les détails de ré-initialisation du compteur et des cas au bord de l'image.

Inutile de dire que niveau perfs la deuxième méthode est largement meilleur que la première.

Bon courage

**nadjib2007** · 23/07/2007, 19h28

Envoyé par befalimpertinent

Il y a surement plusieurs façon de procéder. En voici deux. Je suppose que ton cercle est représenté par son centre et son rayon en pixel.
[LIST=1][*]Soit tu parcours tous tes pixels et tu calcules à chaque fois sa position par rapport au cercle.[*]Soit tu commences par trouver les pixels sur le contour du cercle, tu les initialises dans ta matrices à m[i][j]=0. Puis tu parcours tout tes pixels dans un l'ordre (de gauche à droite et de haut en bas par exemple) et tu regarde si le pixel courant est déjà initialisé à 0 dans ce cas tu incrémente un compteur de 1. Tous les pixels situés sur la même ligne seront des m[i][j]==-1 jusqu'à ce que ton compteur passe à 2.
Je te passes les détails de ré-initialisation du compteur et des cas au bord de l'image.

Inutile de dire que niveau perfs la deuxième méthode est largement meilleur que la première.

Bon courage

Merci beaucoup pour la réponse et vous avez raison la 2ème est beaucoup mieux que la première, mais toujours problème d'implémentation, je vais essayer ici de poster mon programme et a vous de m'aider a corriger les erreurs:

Code :

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

1
2
3
4
5
6
7
 
par exemple pour le cercle je peut le dessiner en rouge et je met la condition:
 
for (int i=0;i<H;i++)
    for (int j=0;j<W;j++)
         if (Image1->Canvas->Pixels[j][i]==clRed)
           mc[i][j]=0; //comme ca on a trouver les pixels du contour

mais après ça je ne sais pas quoi faire car il faut que à l'intérieur soit -1 et à l'extérieur soit +1, et si on arrive a trouver l'un des deux on peut facilement trouver l'autre car l'image contient : pixels internes+pixels externes+contour; donc si on connait deux on peut tirer le troisième.

j'attends vos réponse et Merciiiiiiii

**befalimpertinent** · 24/07/2007, 09h45

Une fois que tu as initialiser les pixels sur contour du cercle avec ta boucle :

Envoyé par nadjib2007

Code :

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

1
2
3
4
5
6
7
 
par exemple pour le cercle je peut le dessiner en rouge et je met la condition:
 
for (int i=0;i<H;i++)
    for (int j=0;j<W;j++)
         if (Image1->Canvas->Pixels[j][i]==clRed)
           mc[i][j]=0; //comme ca on a trouver les pixels du contour

Moi je verrais bien un truc du style :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
 
bool inside = false;
//tu re-parcours les pixels de ton image
for (int i=0;i<H;i++)
   {
    for (int j=0;j<W;j++)
       {
         if (mc[i][j]==0)
           inside=!inside;
         else
           {
             if(inside) //pixels dans le cercle
               mc[i][j]=-1;
             else //pixels a l exterieur
               mc[i][j]=1;
           }        
       }
   }

C'est très simplifié car il y a surement des petites subtilités si le cercle touche le bord ou dépasse le cadre de l'image. ça dépend aussi de la façon dont tu initialise la matrice de contour à sa création ( 1 ?) pour que la condition mc[i][j]==0 corresponde réellement aux pixels du cercle.

**yan** · 24/07/2007, 11h26

Je pense que le plus rapide est d'initialiser ton image a -1, et de traiter ton cercle
un truc du genre

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
InitImage(mc,-1);
int r_r = rayon*rayon;
for (int i=yc-r;i<yc+r;i++)
   {
   //si l'indice de ligne est hors de l'image on passe a la ligne suivante
   if (i<0 || i>=H) continue;
 
   //calcule des deux indices de colonne sur le cercle pour la ligne i
   int dy =(i-yc);
   int x_cercle = sqrt(r_r-dy*dy);
   int xcmin =xc - x_cercle;
   int xcmax =xc + x_cercle;
 
   for (int j=Xc-r;j<Xc+r;j++)
       {
       //si l'indice de colone est hor de l'image on passe au pixel  suivant
       if (j<0 || j>=W) continue;
 
       if (j==xcmin ||j==xcmax)
		{
		//pixel sur le cercle
		mc[i][j]==0)
		}	
	else if (j>xcmin &&j<xcmax)
		{
		//pixels dans le cercle
                mc[i][j]=1;
		}
         else
               {
               //pixels hors du cercle
               mc[i][j]=-1;
               }        
       }
   }

tu doit pouvoir diviser par 4 le traitement du cercle, en traitant 1/4 du cercle et en répercutant le résultat sur les autres partie par "miroire"( je trouve plus le mot

)

**befalimpertinent** · 24/07/2007, 11h43

Envoyé par Mongaulois

tu doit pouvoir diviser par 4 le traitement du cercle, en traitant 1/4 du cercle et en répercutant le résultat sur les autres partie par "miroire"( je trouve plus le mot

)

Symétrie je crois

**yan** · 24/07/2007, 11h46

Envoyé par befalimpertinent

Symétrie je crois

ouai merci

**nadjib2007** · 25/07/2007, 00h14

Envoyé par Mongaulois

Je pense que le plus rapide est d'initialiser ton image a -1, et de traiter ton cercle
un truc du genre

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
InitImage(mc,-1);
int r_r = rayon*rayon;
for (int i=yc-r;i<yc+r;i++)
   {
   //si l'indice de ligne est hors de l'image on passe a la ligne suivante
   if (i<0 || i>=H) continue;
 
   //calcule des deux indices de colonne sur le cercle pour la ligne i
   int dy =(i-yc);
   int x_cercle = sqrt(r_r-dy*dy);
   int xcmin =xc - x_cercle;
   int xcmax =xc + x_cercle;
 
   for (int j=Xc-r;j<Xc+r;j++)
       {
       //si l'indice de colone est hor de l'image on passe au pixel  suivant
       if (j<0 || j>=W) continue;
 
       if (j==xcmin ||j==xcmax)
		{
		//pixel sur le cercle
		mc[i][j]==0)
		}	
	else if (j>xcmin &&j<xcmax)
		{
		//pixels dans le cercle
                mc[i][j]=1;
		}
         else
               {
               //pixels hors du cercle
               mc[i][j]=-1;
               }        
       }
   }

tu doit pouvoir diviser par 4 le traitement du cercle, en traitant 1/4 du cercle et en répercutant le résultat sur les autres partie par "miroire"( je trouve plus le mot

)

Merci beaucoup pour la réponse je vais tester ce programme et si j'aurai des résultats je les poster ....

et merci beaucoup a vous " befalimpertinent " pour tes réponses ça ma fait plaisir

**nicroman** · 25/07/2007, 00h53

A la seule lecture du code , il me semble que ca ne marchera pas sur les parties "horizontales" du cercle...
(edit) Et le premier test me semble louche: si le cercle est plus grand d'entrée de jeu, on fait rien et on laisse à -1 comme si c'était hors du cercle

Choix 1: Refaire la boucle dans le sens vertical (sans mise à +1 à l'interieur cette fois).
Choix 2: Utiliser à mort les symétries infinies du cercle, notamment autour des diagonales (pour éviter le problême évoqué).

**yan** · 25/07/2007, 09h23

C'est juste une ébauche. Y as sûrement quelque erreur

Envoyé par nicroman

A la seule lecture du code , il me semble que ca ne marchera pas sur les parties "horizontales" du cercle...

comment ca?

Envoyé par nicroman

(edit) Et le premier test me semble louche: si le cercle est plus grand d'entrée de jeu, on fait rien et on laisse à -1 comme si c'était hors du cercle

euh. ben si une parti du cercle est hors de l'image, on peut pas faire grand chose, donc ou passe a une autre partie

Envoyé par nicroman

Choix 1: Refaire la boucle dans le sens vertical (sans mise à +1 à l'intérieur cette fois).

pas terrible pour de l'image. Les parcours horizontaux sont plus rapide que les parcours verticaux. Bien sur, uniquement si l'image est implémenté suivant le sens horizontale

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
if (j==xcmin ||j==xcmax)
		{
		//pixel sur le cercle
		mc[i][j]==0)
		}	
	else if (j>xcmin &&j<xcmax)
		{
		//pixels dans le cercle
                mc[i][j]=1;
		}
         else
               {
               //pixels hors du cercle
               mc[i][j]=-1;
               }
}

peut etre remplacer par

Code :

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

1
2
3
4
5
6
7
8
9
10
11
if (j==xcmin ||j==xcmax)
		{
		//pixel sur le cercle
		mc[i][j]==0)
		}	
	else if (j>xcmin &&j<xcmax)
		{
		//pixels dans le cercle
                mc[i][j]=1;
		}
}

puise que l'image est déjà initialisé a -1

mais bon ca ne se vera pas

**yan** · 25/07/2007, 09h25

Aprés si le cercle fait une bonne partie de l'image, autant traiter le cercle en même temp que l'initialisation, mais y aura un plus de test

**JolyLoic** · 25/07/2007, 22h29

Je pense que la lecture de http://www.cs.unc.edu/~mcmillan/comp...e7/circle.html pourra être assez instructive, sur la façon de tracer la frontière du cercle, ensuite, en parcourant le tableau et en comptant le nombre de fois où l'on traverse la frontière, ça ne devrait pas être trop difficile.

**nicroman** · 26/07/2007, 00h25

Excellent lien...
Ce qui est marrant ... les deux premiers essais font exactement la même erreur que Mongaulois.
Mais la dernière, utiliser la symétrie diagonale... par contre....

**yan** · 26/07/2007, 09h38

Envoyé par nicroman

Excellent lien...
Ce qui est marrant ... les deux premiers essais font exactement la même erreur que Mongaulois.
Mais la dernière, utiliser la symétrie diagonale... par contre....

ah ouai je vien de comprendre le problème.
En utilisant les xcmax xcmin de la ligne précédente,
y as sûrement moyen de palier au problème
genre

Code :

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

1
2
3
4
5
6
if (j==xcmin || (j>xcmin_old &&j<=xcmin) 
     ||j==xcmax || (j>=xcmax &&j<xcmax_old)  )
		{
		//pixel sur le cercle
		mc[i][j]==0)
		}

**nadjib2007** · 27/07/2007, 15h24

ohh ...vraiment c'est extra ordinaire vos réponses et vos aidés merci beaucoup à vous tousssssssssssssssss

Problème d'initialisation du signe

C++

Discussions similaires

Partager

Partager