Discussion :

[nadjib2007]

Bonjour tous le monde;

Je suis entrain de programmer un logiciel pour segmenter une image à l'aide de la méthode des ensembles de niveaux (Level Sets); et pour cela j'utilise comme fonction distance "phi" un cercle (contour initial).

j'ai ecrit ce petit programme:

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
// Initialisation de la matrice "Phi"
 
// choix du centre
float x0=76;
float y0=91;
// choix du rayon
float R=10;
 
 
// dessiner le cercle
 
for(int i=0;i<H;i++)
  for(int j=0;j<W;j++)
    {
    phi[i][j]=sqrt(((x0-j)*(x0-j))+((y0-i)*(y0-i)))-R;
    }

Avec le Matlab tous marche bien; j'ai eu un cercle bien dessiner, tel que avec l'utilisation de l'équation du cercle j'aurai des distance négatifs à l'intérieur du cercle et positifs à l'extérieur et le "zéro" sera pour les pixels qui appartient aux cercle => c'est mon but.

mais le problème, comme "i" et "j" sont des entier et la matrice phi,x0,y0 sont des "float" on aura pas un 0 dans tous les points du cercle, peut être on trouve 0.000001 mais pas un 0.

alors comment on résout ce problème es qu'on change "float " ou int ou ...je ne sais pas...

Merci de bien vouloir répondre

[rtg57]

Bonjour,

C'est toujours le même problème avec les nombres à virgule flottante: la précision.
Même si l'on utilise des float de grande précision, de part leur principe de codage, il y a toujours une marge d'erreur.

Personnellement, j'ai résolu le problème en n'utilisant plus ces nombres
J'utilise des type int ou long, voir double et pourquoi pas long double
Je fais tous mes calculs avec cela, puis je divise finalement le résultat par cent ou mille ou plus avec forcage du nombre de décimales.
Exemple: si je fais un programme gérant de l'argent, un fonctionnement au centime près est satisfaisant.
Donc je fais tous mes calculs en centimes, avec des variables de type int. Sachant qu'elles sont codées sur 32 bits, cela m'autorise des montants de -2 147 483 648 centimes jusqu'à 2 147 483 647 centimes. Ce qui fait -21 474 836,48€ à -21 474 836,47€
C'est largement suffisant pour calculer mon salaire
L'affichage des résultats en euros se fait par une division par 100, avec arrondissement à la 2ème décimale.
L'avantage de cette solution est que le 0 centime existe, et donnera 0 euros
Peut être que cette astuce vous permettra de résoudre votre problème ...

[henderson]

Salut !

Il y a de fortes chances pour que seuls quatre points satisfassent directement à la relation : d - r = 0 !

Donc on peut faire ainsi :

Code :

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

1
2
 
float d;

Code :

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

1
2
3
4
5
 
d = sqrt(...) - R;
//A toi de voir pour l'intervalle...
if((d > -0.5) && (d < 0.5)) d = 0.0;
Tab[..][..] = d;

A plus !

[Sunchaser]

La solution antique mais fiable pour les cercles :