Discussion :

[nadjib2007]

Bonjour tous le monde; j'ai un petit code en Matlab extrait d'un programme de segmentation d'image par le mod�le de Ron Kimmel)

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
% 'Roughly' correct Phi to be a distance map of its zero set
function u=redistance(Phi);
u=(sign(Phi)+1)*999999; % set to infinity all positive
for i=1:2,
    l2=2;
    if i>1 
        u=(1-sign(Phi))*999999; 
    end % set to infinity all negative
    while l2>1,
          v=Update(u,1);
          l2=sum(sum((u-v).^2));
          u=v;;
    end % while
    if i>1 
       u=u-up; 
   else 
       up=u; 
 
   end % if
end % for
%-----------------------------------------------------------------------------------------------
%Solve |grad u|=F(x,y)  parallel version of the FMM
function res=Update(u,F)
mx=min(u([2:end end],:), u([1 1:end-1],:));
my= min(u(:,[2:end end]), u(:,[1 1:end-1]));
delm=(mx-my);
mask=(abs(delm)<F);
res=min(mask.*(mx+my+sqrt(2.*F.^2-delm.^2))./2+...
   (1-mask).*(min(mx,my)+F),u);

tel que "phi" une matrice carr� quelconque.

j'ai essayer de convertir ce code en C++ Builder (manuellement), mais je ne sais pas si je suis juste ou non, alors s'il y a quelqu'un qui peuvent m'aider je le remercie d'avance.

voici le code que j'ai �crit:

pour la fonction "redistance"

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
78
79
80
81
void __fastcall  TForm1::redistance(float **phi, int H, int W)
 
{
 
float**up=new float*[H];  // la heuteur de l'image
for(int i=0;i<H;i++)
  {
   up[i]=new float[W];   // la largeur de l'image
  }
 
float**v=new float*[H];
for(int i=0;i<H;i++)
  {
   v[i]=new float[W];
  }
 
float l2, l22;
 
for (int i=0;i<H;i++)
    for (int j=0;j<W;j++)
       {
        phi[i][j]=(fctsign(phi[i][j])+1)*999999; // fixé le tous positif
       }
 
        for (int z=1;z<3;z++)
           { 
            l2=2;
            if (z>1)
              {
               for (int i=0;i<H;i++)
               for (int j=0;j<W;j++)
                  {
               phi[i][j]=(1-fctsign(phi[i][j]))*999999; //fixé le tous négative
                  }
              }   // end if
                while (l2>1)
                    {
                    Updt(phi,v,H,W); // fonction update
                     l22=0;
                    for (int i=1;i<H-1;i++)
                      for (int j=1;j<W-1;j++)
                         {
                         l22=l22+pow((phi[i][j]-v[i][j]),2);
                         }
                    l2=l22;
                    for (int i=1;i<H-1;i++)
                    for (int j=1;j<W-1;j++)
                       {
                        phi[i][j]=v[i][j];
                       }
                    } // end while 
 
                if (z>1)
                  {
                   for (int i=0;i<H;i++)
                   for (int j=0;j<W;j++)
                       {
                       phi[i][j]=phi[i][j]-up[i][j];
                       }
                   }
                else
                    {
                     for (int i=0;i<H;i++)
                     for (int j=0;j<W;j++)
                        {
                         up[i][j]=phi[i][j];
                        }
                     }
 
               }// end for z
 
 
for (int i=0;i<H;i++)
   delete[]up[i];
   delete[]up;
 
for (int i=0;i<H;i++)
   delete[]v[i];
   delete[]v;
 
 }

et pour la fonction "Update ":

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
void __fastcall  TForm1::Updt(float **u,float **v, int H, int W)
 
{
 
 float mx,my,diff,minimum;
 
//Résoudre l'équation |grad u|=F(x,y)
 
//  déclaration de la matrice mask
float**mask=new float*[H];
for(int i=0;i<H;i++)
   {
    mask[i]=new float[W];
   }
 
// déclaration de  la matrice produit
float**produit=new float*[H];
for(int i=0;i<H;i++)
   {
    produit[i]=new float[W];
   }
 
for (int i=1;i<H-1;i++)
  for (int j=1;j<W-1;j++)
     {
     my=min(u[i][j+1],u[i][j-1]);
     mx=min(u[i+1][j],u[i-1][j]);
     minimum=min(mx,my);
     diff=mx-my;
 
     if (abs(diff)<1) // abs(mx-my)<1
       {
        mask[i][j]=1;
        produit[i][j]=mask[i][j]*((mx+my+sqrt(2-(diff*diff)))/(2+((1-mask[i][j])*(minimum+1))));
       }
     else
         {
          mask[i][j]=0;
          produit[i][j]=0;
         }
     }
 
for (int i=1;i<H-1;i++)
   for (int j=1;j<W-1;j++)
      {
       v[i][j]=min(produit[i][j],u[i][j]);
      }
 
// Libération de mémoire
 
for(int i=0;i<H;i++)
   delete[]mask[i];
   delete[]mask;
 
for(int i=0;i<H;i++)
   delete[]produit[i];
   delete[]produit;
}

et s'il y a des fautes merci de bien les signal�es.

[Sunchaser]

et s'il y a des fautes merci de bien les signal�es.

Bonsoir,

Ben, pour commencer -> ' signal�es '
Cela s'�crit -> ' signaler '

Voili ... et d'une.



Lol, d�sol�, pour le reste, je ne connais pas.

[nadjib2007]

lol...je crois que j'ai poster mon sujet dans le forum de C++ Builder et non pas dans "lecture et �criture"...

Merci pour la remarque j'ai pas fait attention.

[Crayon]

Salut, j'ai une question qui peut paraitre un peu stupide, est-ce que tu as essayer de mettre des valeurs dans ta fonctions BCB et de comparer les r�sultats avec ta fonction MatLab? Si les r�sultats sont les m�mes il y a des chances que tu as fais un bon travail et sa �vite aux membres du forum de faire la validation

[nadjib2007]

Merci pour la r�ponse j'ai u cette aider mais comment je pourrais faire ca?!

pour mettre des valeurs c'est facile mais comment je peut avoir le r�sultat � la fin, car je doit trouver une matrice !!

Merci encore une fois

[Crayon]

Si c'est une matrice 2D tu peux la mettre dans un tableau, sinon tu peux faire la sommes de tout les �l�ments de la matrice et comparer le r�sultat. Utilises ton imagination, tu devrais trouver une solution.

[nadjib2007]

ok je vais essayer �a manuellement ( stylo et papier) et je vous r�pond ....thnx

[henderson]

Salut !

Tu pourrais utiliser un TStringGrid pour afficher tes valeurs.
A noter que chaque cellule est du type AnsiString (cha�ne de caract�res).
Donc pour afficher une valeur num�rique du type float (ou double) il te faudra faire une conversion (num�rique -> litt�ral) :

Code :

S�lectionner tout - Visualiser dans une fen�tre � part

StringGrid1->Cells[x][y] = FloatToStr(Une_de_tes_valeurs);

A plus !

[nadjib2007]

Merci beaucoup " mod�rateur ", c'est �a qu'on appel => gagner du temps .

[Neilos]

Tu passe en mode debug et tu testes sur une petite matrice.
Tu met un point d'arret a la fin de ton algorithme et tu regardes le contenu de ta matrice...ou alors tu sors tout sur une console, ou ou ou....

Ne connaissant pas bien Matlab je ne m'avancerais pas sur l'exactitude de ta traduction.

[henderson]

Salut !

Merci beaucoup " mod�rateur ", c'est �a qu'on appel => gagner du temps .

Oui parce que tu te sers de BCB pour v�rifier ce qu'il fait (en fait ce que fait ton code) ! D�s que la v�rification est faite (comme le dit Neilos avec une matrice de taille r�duite comportant des valeurs telles que le r�sultat soit conforme � ceux attendus), il ne te reste plus qu'� virer le TStringGrid et passer � la suite !

On gagne aussi beaucoup de temps � �tudier ce que BCB nous offre !

A plus !

[nadjib2007]

Merci beaucoup les amis j'ai r�ussi a trouver l'erreur c'�tait dans les parenth�ses

Merci encore une fois du fond du coeur..

[Neilos]

Un ptit tag alors ?

[nadjib2007]

R�solu....! Ouiii bien sur avec plaisir; et s'il y a un bouton "r�sultat satisfaisante" je le pr�f�re..

[Neilos]

Merci beaucoup et � bient�t sur le forum.

[koala01]

Salut,

Juste une petite remarque au passage...

Sais tu qu'il n'y a rien de plus proche, du point de vue de l'utilisation de la m�moire, du moins, d'un tableau � deux dimensions... qu'un seul tableau contenant (dim1*dim2) �l�ments

La seule chose, c'est qu'au lieu d'acc�der aux �l�ments via

Code :

S�lectionner tout - Visualiser dans une fen�tre � part

tab [i][j]

on y acc�de avec la formule

Code :

S�lectionner tout - Visualiser dans une fen�tre � part

tab[i*jmax+j]

...

Mais cela ouvre un max de possibilit�s...

D'abord, tu peux envisager de travailler avec un pointeur, au lieu d'avec... un pointeur de pointeur

Ensuite, cela te permet de faire l'allocation de m�moire d'une seule traite

Code :

S�lectionner tout - Visualiser dans une fen�tre � part

float* tab=new float[imax*jmax];//ou imax est le nombre de lignes et jmax le nombre de colones

Evidemment, ce qui est vrai pour l'allocation est vrai pour la lib�ration

Code :

S�lectionner tout - Visualiser dans une fen�tre � part

 delete[] tab;

Enfin, cela peut t'ouvrir des horizons, telles que, par exemple, utiliser les conteneurs de la STL (les vector, par exemple)...

Cela n'a l'air de rien, mais cela te lib�re de l'obligation de g�rer la m�moire (allocation/lib�ration), d'autant plus que, tel que tu le fait l�, c'est pas s�curisant (que se passerait-il si le syst�me tombait � court de m�moire pour allouer les W derniers float)

Maintenant, ce que j'en dis, ce n'est que pour �carter un risque de "plantage" qui, d'apr�s la loi de Murphy, appara�tra le jour de la grande d�monstration