Discussion :

[kayser33]

salut a tous : j'ai essayé de programmé la méthode de binarisation by OTSU , j'implémente le programme en C et je l'envoie ensuite sur une cible (carte DSP) mais l'image binarisée que j'obtiens ne correspond pas du tout à l'image initiale.

quelqu'un pourrait il m'aider à trouver le bug SVP

voici le code

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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
 
void Bin_Otsu3( void *_imgsrc, void *_imgdst, uint32 _largeur, uint32 _hauteur )
{    
 
	uint16		histogram[256];
	uint16		omega[256]; 
	uint32		myu[256];   
 
	double	max_sigma;
	double	sigma;          //inter-class variance 
 
	uint8	i, threshold; 	// threshold for binarization 
	uint16	n;
	double	k; 				
 
 
	uint8	*psrc, *pdst, *pend;
 
	pend    = (uint8*)_imgsrc + _largeur * _hauteur;
	n		= _largeur * _hauteur; 	// n = 88 *72 = 6336 pixels 
 
 
	// Histogram generation 
 
	for (i = 0; i < 256; i++ )	
		histogram[ i ] = 0;
 
	for (psrc = (uint8*)_imgsrc; psrc < pend; psrc++)
		histogram[ *psrc ]++;
 
 
 
	// On calcule les deux premiers moments 
	// myu=moy et omega= ecart-type			
 
	omega[0] = histogram[0];		
	myu[0]   = 0;					
 
	for (i = 1; i < 256; i++) 
	{
		omega[i] = omega[i-1] + histogram[i]; 
 
		myu[i]   = myu[i-1] + i * histogram[i];
	}
 
 
	// sigma maximization
 
	threshold = 0;
	max_sigma = 0.0;
 
	for (i = 0; i < 256; i++) 
	{
		if ((omega[i] == 0) || (omega[i] == n))
			sigma = 0.0;
		else
		{
			k = myu[255] * omega[i] - n * myu[i];	
 
 
			sigma =  (k * k)/ (n * n * omega[i] * (n - omega[i]));	
		}
 
		if (sigma > max_sigma) 
		{
			max_sigma = sigma;
			threshold = i;
		}
	}
 
	// binarization with calculated threshold 
 
	for (psrc = (uint8*)_imgsrc, pdst = (uint8*)_imgdst; psrc < pend; psrc++, pdst++)
	{
		*pdst = (*psrc < threshold) ? 0 : 255;
	}
}

MERCI

[Invité(e)]

Bonjour,

Pour commencer, les boucle for me paraissent louches :

Code :

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

1
2
3
4
5
 
uint8 i;
for (i = 1; i < 256; i++) {
  ...
}

et gcc est d'accord avec moi :

Code :

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

warning: comparison is always true due to limited range of data type

puisque i varie entre 0 et 255, il ne sera jamais supérieur ou égal à 256, on ne sertira donc jamais de le boucle.

De plus, juste pour simplifier le code, plutôt qu'écrire

Code :

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

1
2
for (i = 0; i < 256; i++ )	
  histogram[ i ] = 0;

pour quoi ne pas initialiser histogram lors de sa déclaration ?

Code :

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

uint16 histogram[256] = { 0 };

[kayser33]

jte remercie je vais modifier l'initialisation et voir le résultat.

[Invité(e)]

Bonjour,

jte remercie je vais modifier l'initialisation et voir le résultat.

Ma remarque sur l'initialisation était secondaire, le problème vient de la déclaration de i. Il faut que i puisse prendre une valeur supérieure à 255. Il faut donc que i ne soit pas de type uint8, mais au moins uint16.

[kayser33]

ah ok! en fait comme i ne prenait que 255 valeurs je me disais qu'un uint8 était largement suffisant mais l'histoire de la boucle infinie je n'y avais pas fait attention.

j'espère que ça va enfin régler mon problème

merci à toi et à tous ceux qui prennent le temps de répondre !

[kayser33]

en fait tu l'a compilé avec quel logiciel?

[kayser33]

je t'envoie un exemple d'image sur laquelle je cherche a isoler la personne présente . La binarisation est correcte.

Quand je transfère le programme sur cible l'image binarisée que j'obtiens n'a pas la meme forme.

[kayser33]

désolé j'ai oublié de joindre le fichier ;-)

[stephl]

Je ne vois toujours pas le fichier. L'as-tu joint?

[kayser33]

desolé au moment ou je téléchargé les pièces la page internet s'est fermée ! les voici

[stephl]

Pas de problème.

Merci. Cependant, ce n'est pas vraiment ce que j'attendais (je me vois mal reconstuire le tableau d'octets correspondant à l'image rien qu'en la regardant).
En fait, il me faudrait un fichier contenant les octets de l'image source tels que ta fonction Bin_Otsu3() les reçoit, ainsi qu'un autre fichier contenant les octets de l'image de destination après binarisation (les valeurs correctes calculées par matlab). Bien entendu, il faut qu'il s'agisse d'une image source avec laquelle tu as constaté des résultats différents sous matlab et avec ton programme c tournant sur le dsp.

[kayser33]

autant pour moi! je t'envoie les données hexadecimales le plus vite possible!merci