Tracer un disque - conversion de code C

Version imprimable

Bonjour !

Voila j'ai un petit probleme : je cree des images en C, ce sont des images avec un fond de 0 et un disque d'intensite fixe.

J'essaie ensuite de recreer la meme image avec Matlab. Mais impossible !

Voici ma fonction C :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 void Circ(float **data, int N, int M, int centerx, int centery, int radius,float Isig) { int n,m; float temp; for(n=0;n<=N-1;n++) { for(m=0;m<=M-1;m++) { temp=sqrt(pow((double)((n-N/2)-centerx),2)+pow((double)((m-M/2)-centery),2)); if(temp<=radius) { data[n][m]=Isig; } else data[n][m]=0; } } }
Et la meme avec matlab :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 function s = circulaire(N,M,centerx,centery,radius,intensity) s = zeros(N,M); for n=1:N for m=1:M temp = sqrt(((n-N/2)-centerx-1)^2+((m-M/2)*1.0-centery-1)^2); if (temp<=radius) s(n,m) = intensity; else s(n,m) = 0; end end end s = reshape(s,N*M,1);
Le seul point pour lequel ca marche c'est quand le centre du disque est a [0,0].

Merci de votre aide pour trouver pourquoi ca ne marche pas, la vraiment je seche :aie:

cecile

04/03/2011, 23h42
Invité

bonsoir,

La ligne

Code:

s = reshape(s,N*M,1);

est de trop, elle fait passer la variable s d'une matrice à un vecteur, tu ne peux plus faire ensuite

Code:

imshow(s)

De plus,

Code:

temp = sqrt((n-centerx-1)^2+(m-centery-1)^2);

me semblerait plus juste.
05/03/2011, 00h03
Hew
Merci pour ta reponse. En fait la ligne reshape est bien la pour stocker la matrice s de la meme facon que je stocke mes matrices en c. Pour faire l'affichage, je refais un reshape.

La ligne
Code:

1 2 temp = sqrt((n-N/2-centerx-1)^2+(m-M/2-centery-1)^2);
Est la parce que centerx et centery peuvent etre negatifs. Le point [0,0] est donc au centre de l'image.

Je cherche encore !
05/03/2011, 00h13
Hew

1 pièce(s) jointe(s)

Voila, pour vous aider je rajoute une image du probleme :
le disque cree par ma fonction C, celui cree par ma fonction matlab et enfin la difference des deux.

Citation:

Envoyé par Hew

je cree des images en C, ce sont des images avec un fond de 0 et un disque d'intensite fixe.

J'essaie ensuite de recreer la meme image avec Matlab. Mais impossible !

Pourquoi ne pas passer par un fichier C-MEX Circ.c :

Code:

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
#include "mex.h"
#include "math.h"
 
void Circ(float *, int, int, int, int, int, float);
 
void mexFunction(
        int nlhs, mxArray *plhs[],
        int nrhs, const mxArray *prhs[]) {
 
    int centerx, centery, radius;
    float *data, Isig;
    mwSize M, N;
 
    M = (mwSize)mxGetScalar(prhs[0]);
    N = (mwSize)mxGetScalar(prhs[1]);
    centerx = (int)mxGetScalar(prhs[2]);
    centery = (int)mxGetScalar(prhs[3]);
    radius = (int)mxGetScalar(prhs[4]);
    Isig = (float)mxGetScalar(prhs[5]);
 
    plhs[0] = mxCreateNumericMatrix(M, N, mxSINGLE_CLASS, mxREAL);
    data = mxGetData(plhs[0]);
 
    Circ(data, N, M, centerx, centery, radius, Isig);
 
}
 
void Circ(float *data, int N, int M, int centerx, int centery, int radius, float Isig) {
    int n, m;
    float temp;
 
    for(n=0;n<N;n++) {
        for(m=0;m<M;m++) {
            temp = sqrt(pow(((n-N/2)-centerx), 2)+pow(((m-M/2)-centery), 2));
            if(temp<=(float)radius)
                data[m+n*M] = Isig;
        }
    }
}

A compiler comme ceci :

Code:

mex Circ.c

Puis sous MATLAB :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
N = 20;
M = 20;
centerx = 2;
centery = 4;
radius = 4;
Isig = pi;
 
s = Circ(N, M, centerx, centery, radius, Isig);
 
figure(1);
image(s);
axis image

Il faudra bien entendu ajouter quelques tests dans le code pour améliorer sa robustesse...

Merci Dut !
C'est exactement ce qu'il me fallait. Je ne sais pas trop m'y prendre avec les mex, c'est pour ca que j'avais repousse l'utilisation jusqu'a maintenant mais c'est vrai que c'est l'deal.

Voici mon fichier c-mex :

Code:

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
 
#include "mex.h"
#include "math.h"
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <limits.h>
 
void CircMex(float **, int, int, int, int, int, float);
 
void mexFunction(
        int nlhs, mxArray *plhs[],
        int nrhs, const mxArray *prhs[]) {
 
    int centerx, centery, radius;
    float **data, Isig;
    mwSize M, N;
 
    M = (mwSize)mxGetScalar(prhs[0]);
    N = (mwSize)mxGetScalar(prhs[1]);
    centerx = (int)mxGetScalar(prhs[2]);
    centery = (int)mxGetScalar(prhs[3]);
    radius = (int)mxGetScalar(prhs[4]);
    Isig = (float)mxGetScalar(prhs[5]);
 
    plhs[0] = mxCreateNumericMatrix(M, N, mxSINGLE_CLASS, mxREAL);
    data = mxGetData(plhs[0]);
 
    CircMex(data, N, M, centerx, centery, radius, Isig);
 
}
 
void CircMex(float **data, int N, int M, int centerx, int centery, int radius,float Isig)
{
  int n,m;
  float temp;
  for(n=0;n<=N-1;n++)
  {
	  for(m=0;m<=M-1;m++)
	  {
		  temp=sqrt(pow((double)((n-N/2)-centerx),2)+pow((double)((m-M/2)-centery),2));
		  if(temp<=radius)
		  {
			  data[n][m]=Isig;
		  }
		  else
			  data[n][m]=0;
      }
  }
 
}

Ecris comme ca, ca compile sans erreur et me donne un fichier circmex.mexglx.
J'appelle ensuite ma fonction sous matlab comme ceci :

Code:

1
2
3
 
circ3 = zeros(32,32);
circ3 = CircMex(circ3,32,32,posx(498),posy(498),rad(498),isig(498));

Je recois l'erreur suivante :

Code:

1
2
3
 
??? Undefined function or method 'CircMex' for input arguments of type
'double'.

En effet, ma fonction est definie pour un float, je reessaie donc comme ceci :

Code:

1
2
3
 
circ3 = zeros(32,32);
circ3 = CircMex(float(circ3),32,32,posx(498),posy(498),rad(498),isig(498));

et je recois l'erreur suivante :

Code:

1
2
3
4
5
6
7
8
9
 
??? Operands to the || and && operators must be convertible to logical
scalar values.
 
Error in ==> float at 56
  elseif (TotalBits == 64) && (ExpBits == 11)
 
Error in ==> testmex at 11
circ3 = CircMex(float(circ3),32,32,posx(498),posy(498),rad(498),isig(498));

Bon la je ne sais pas trop ce que ca veut dire :aie:

Je vais reessayer en creant ma fonction pour un double et pas un float et voir ce que ca donne.
Merci pour l'aide et d'autres suggestions sont les bienvenues !

cecile

07/03/2011, 19h35
Hew
J'ai fait un petit changement dans mon test :
Code:

1 2 3 4 circ = zeros(32,32,'double'); circ = reshape(circ,32*32,1); circ3 = circmex(circ,32,32,posx(498),posy(498),rad(498),isig(498));
Cette fois le programme s'execute sans planter mais la matrice circ3 est une matrice vide :(

Je seche maintenant :aie:

Voila comment je ferais :

Code:

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
#include "mex.h"
#include "math.h"
 
void CircMex(float *, int, int, int, int, int, float);
 
void mexFunction(
        int nlhs, mxArray *plhs[],
        int nrhs, const mxArray *prhs[]) {
 
    int centerx, centery, radius;
    float *data, Isig;
    mwSize M, N;
 
    data = mxGetData(prhs[0]);
 
    M = (mwSize)mxGetScalar(prhs[1]);
    N = (mwSize)mxGetScalar(prhs[2]);
    centerx = (int)mxGetScalar(prhs[3]);
    centery = (int)mxGetScalar(prhs[4]);
    radius = (int)mxGetScalar(prhs[5]);
    Isig = (float)mxGetScalar(prhs[6]);
 
    CircMex(data, N, M, centerx, centery, radius, Isig);
 
}
 
void CircMex(float *data, int N, int M, int centerx, int centery, int radius, float Isig) {
    int n, m;
    double temp;
    for(n=0;n<N;n++) {
        for(m=0;m<M;m++) {
            temp=sqrt(pow((double)((n-N/2)-centerx), 2)+pow((double)((m-M/2)-centery), 2));
 
            if(temp<=(double)radius)
                data[n+m*N]=Isig;
 
        }
    }
 
}

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
clc
N = 20;
M = 10;
centerx = -1;
centery = 1;
radius = 3;
Isig = pi;
 
circ3 = zeros(M,N,'single');
Circ(circ3,N, M, centerx, centery, radius, Isig);
 
figure(1);
image(circ3);
axis image

Utilise l'indexage linéaire dans les fichiers MEX, c'est comme ça que sont stockées les valeurs sous MATLAB (voir mxArray)... ou bien des directives define comme ici : how to access multi-dimensional array input in C mex file

Par contre, je ne suis pas certain que la ligne suivante fasse exactement ce que tu souhaites :

Code:

temp=sqrt(pow((double)((n-N/2)-centerx), 2)+pow((double)((m-M/2)-centery), 2));

08/03/2011, 00h27
Hew
Merci Dut, je comprends mieux le coup de l'indexage lineaire. Ca a l'air de marcher : ma matrice n'est plus vide maintenant. Par contre toujours pas de disque.
Est-ce que tu penses que cette fonction :
Code:

1 2 temp=sqrt(pow((double)((n-N/2)-centerx), 2)+pow((double)((m-M/2)-centery), 2));
ne fait pas ce que je veux en c-mex ? Ca marche plutot bien en c...

Je vais investiguer de ce cote la.
08/03/2011, 03h45
Hew
Hello !

Voila, j'ai resolu mon probleme : en appelant ma fonction comme ceci :
Code:

1 2 circ3 = circmex(32,32,posx(498),posy(498),rad(498),isig(498));
Merci beaucoup de ton aide Dut, au moins j'en sais un peu plus sur les mex files maintenant ! :ccool: