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]

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