convertir code C en MATLAB

**zied113** · 21/05/2014, 22h19

Bonsoir,

SVP je veux exécute le code ci -joint en C sous Matlab . Merci bcp

**Jerome Briot** · 21/05/2014, 23h17

Pour commencer, il est plus respectueux de citer ses sources : Perception and Neurodynamics Laboratory (PNL) et plus précisément : http://www.cse.ohio-state.edu/~dwang...are/wang-nc97/

Le code est suffisamment simple pour être facilement transformé en MEX :

Code C :

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
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
       /*****************************************************************/
       /* This program is for 2D segmentation using the LEGION network  */
       /* The equations were modeled as a PDP algorithm for efficiency  */
       /* by DeLiang Wang in 7/94. For a detailed description of the    */
       /* algorithm see D.L. Wang and D. Terman: "Image segmentation    */
       /* based on oscillatory correlation", Neural Computation, Vol. 9 */
       /* pp. 805-836, 1997 (For errata see Neural Computation 9,       */
       /* 1623-1626, 1997), on page 822.                                */
       /*****************************************************************/
 
#include "mex.h"
#include "segment.h"
 
#define MAXVAL 32768.0
               /* 2147483648.0 or 32768.0,machine (HP or SUN) dependable*/
#define MAX 32767 /* 2147483647 or 32767, machine dependable */
#define GRAY 255  /* Maximum value of a pixel */
 
int bi();
int jump[2], z, Num = 1;             /* so we can change step size at will */
                               /* Num: number of the segments */
FILE *outfile, *fopen();
float x[N1][N2][2], lateral[N1][N2];
                      /* lateral: the sum of input from lateral connections */
float W[N1][N2][8];   /* strength of the interactions for 4 directions */
int stim[N1][N2], matrix[N1][N2];
int y[N1][N2], cycle[N1][N2], trigger[N1][N2];
     /* y: integer of x; cycle: test for 1st cycle; trigger: jumping signal */
 
 
void inputs(void)
{
   int i, j;
   int width, height, c;
   char magic[25];
   float stimulus = 0.2;
   FILE *IN;
 
   for (i = 0; i < N1; i++) 
     for (j = 0; j < N2; j++)
       matrix[i][j] = 0;
 
   IN = fopen("lake1.pgm", "r");
   fscanf(IN,"%s\n%d %d\n%d\n",magic,&width,&height,&c);
   fprintf(stderr,"[Type %s %dx%d colours %d]\n",magic,width,height,c);
 
   for(i=0;i<160;i++)  
      for(j=0;j<160;j++)
         stim[i][j] = getc(IN);
 
   fclose(IN);
}
 
void connections(void)
{
   int i, j, k;
   float baseline = 0.0, Wtotal = 6.0, sum;
 
   for (i = 0; i < N1; i++) /* for eight nearest neighbors */
     for (j = 0; j < N2; j++) {
       W[i][j][0] = W[i][j][1] = W[i][j][2] = W[i][j][3] = 0.0;
       W[i][j][4] = W[i][j][5] = W[i][j][6] = W[i][j][7] = 0.0;
     }
 
   for (i = 0; i < N1; i++) 
     for (j = 0; j < N2; j++) {
	 sum = 0.0;
 
         if (j-1 >= 0)
	    W[i][j][0] = (float) GRAY/(1+abs(stim[i][j]-stim[i][j-1]));
         if (i-1 >= 0 && j-1 >= 0)
	    W[i][j][1] = (float) GRAY/(1+abs(stim[i][j]-stim[i-1][j-1]));
         if (i-1 >= 0)
	    W[i][j][2] = (float) GRAY/(1+abs(stim[i][j]-stim[i-1][j]));
         if (i-1 >= 0 && j+1 < N2)
	    W[i][j][3] = (float) GRAY/(1+abs(stim[i][j]-stim[i-1][j+1]));
         if (j+1 < N2)
	    W[i][j][4] = (float) GRAY/(1+abs(stim[i][j]-stim[i][j+1]));
         if (i+1 < N1 && j+1 < N2) 
	    W[i][j][5] = (float) GRAY/(1+abs(stim[i][j]-stim[i+1][j+1]));
         if (i+1 < N1)
	    W[i][j][6] = (float) GRAY/(1+abs(stim[i][j]-stim[i+1][j]));
         if (i+1 < N1 && j-1 >= 0)
	    W[i][j][7] = (float) GRAY/(1+abs(stim[i][j]-stim[i+1][j-1]));
 
         for (k = 0; k < 8; k++) lateral[i][j] += W[i][j][k];
    }
}
 
void initiate(void)
{
   int i,j,t, temp;  /* labels for setting up display intervals */
 
   for (i = 0; i < N1; i++)   /* randomize the initial phases */
     for (j = 0; j < N2; j++) {
        y[i][j] = rand();
        x[i][j][0] = y[i][j]/(MAXVAL + 1.0);
	y[i][j] -= 2;
	if (y[i][j] < 0)
	  y[i][j] = 0;
	trigger[i][j] = 0;
	cycle[i][j] = 0;
     }
   z = 0; 
 
 /* initial display values for 2D */
   for (i = 0; i < N1; i++)
      for (j = 0; j < N2; j++) {
        temp = x[i][j][0]*RP;       
        fprintf(outfile, "%d ", temp);
      }
}
 
void select(void)
{
   int i, j, max, imax=0, jmax=0;
   float baseline = 1200.0, temp;
/* baseline 1200 -> 17 regions; 1000 -> 23 regions. Both reasonable */
 
   max = -1;
   for (i = 0; i < N1; i++)
     for (j = 0; j < N2; j++) {
       trigger[i][j] = 0;
       if (lateral[i][j] > baseline && y[i][j] != MAX && y[i][j] > max) {
	 max = y[i][j];
	 imax = i;
	 jmax = j;
       }
     }
 
   if (cycle[imax][jmax]) {     /* find out jumpers */
     for (i = 0; i < N1; i++)
       for (j = 0; j < N2; j++)
	 if (y[i][j] == max && cycle[i][j])
	   trigger[i][j] = 1;
   } else
     trigger[imax][jmax] = 1;
 
   for (i = 0; i < N1; i++)   /* advance oscillators on the left branch */
     for (j = 0; j < N2; j++) 
       if (y[i][j] != MAX && (cycle[i][j] || lateral[i][j] > baseline)) {
	 y[i][j] = y[i][j] + (MAX - max);
         x[i][j][0] = y[i][j]/(MAXVAL + 1.0);
       }
}
 
void evaluate(void)
{
   float temp, thetae = 2.0, temp2, temp3;
   int i,j,t, cat, count;  /* labels for setting up display intervals */
 
   jump[0] = 0;
 
   select();
   for ( t = 0; t < WIDTH-1; t++) {
      count = 0;
      for (i = 0; i < N1; i++)
	for (j = 0; j < N2; j++) { /* compute net input */
	  if (x[i][j][0] > thetae && jump[0] == 0) {  /* jumping down */
             x[i][j][1] = LP;
	     y[i][j] = -1;
	     cycle[i][j] = 1;
	     count++;
	  } else if (x[i][j][0] < thetae) {  /* summing up overall input */
             temp =  -20.0 * bi((float) z,0.5);  /* used to be 50 */
 
/* the following is based on traversing through strongest link */
             temp2 = 0.0;
             if (j-1 >= 0) {
                temp3 = W[i][j][0] * bi(x[i][j-1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i-1 >= 0 && j-1 >= 0) {
                temp3 = W[i][j][1] * bi(x[i-1][j-1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i-1 >= 0) {
                temp3 = W[i][j][2] * bi(x[i-1][j][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i-1 >= 0 && j+1 < N2) {
                temp3 = W[i][j][3] * bi(x[i-1][j+1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (j+1 < N2) {
                temp3 = W[i][j][4] * bi(x[i][j+1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i+1 < N1 && j+1 < N2) {
                temp3 = W[i][j][5] * bi(x[i+1][j+1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i+1 < N1) {
                temp3 = W[i][j][6] * bi(x[i+1][j][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i+1 < N1 && j-1 >= 0) {
                temp3 = W[i][j][7] * bi(x[i+1][j-1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             temp += temp2;
 
	     if (trigger[i][j] || temp > 0.25) {        /* jumping up */
	        x[i][j][1] = RP; 
		y[i][j] = MAX;
	        jump[1] = 1;
	        z++;
	     } else
	        x[i][j][1] = x[i][j][0];
	  } else 
	     x[i][j][1] = x[i][j][0];
	}
      z -= count;
 
      if (jump[1] == 0 && z == 0) 
        select();
 
    /* store for display in 2d case*/
      if (jump[1] == 0 && jump[0] == 1) {
mexPrintf("we are here \n");
	Num++;
        for (i = 0; i < N1; i++)
          for (j = 0; j < N2; j++) {
             if (z > 0) {
	          if (x[i][j][1] < thetae)
                     cat = 0;       
		  else cat = x[i][j][1];       
	     } else  
		  cat = x[i][j][1];       
             fprintf(outfile, "%d ", cat);
 
             if (Num < 25 && cat > 0)
                matrix[N1-1-i][j] = 26-Num;  /* 17 regions, or 23 */
	  }
      }
 
 
      for (i = 0; i < N1; i++)
         for (j = 0; j < N2; j++)  /* replace OLD by NEW */
           x[i][j][0] = x[i][j][1];
      jump[0] = jump[1];
      jump[1] = 0;
   }
}
 
void store(void)
{
   int i, j;
   FILE *out, *out1, *fopen();
 
   mexPrintf("The number of the segments/cycles is %d\n", Num);
 
   out = fopen("NUMBER", "w");
   out1 = fopen("Result", "w");
   fprintf(out, "%d ", Num);
   for (i = 0; i < N1; i++)
     for (j = 0; j < N2; j++)
       fprintf(out1, "%d ", matrix[i][j]);
   fclose(out);
   fclose(out1);
}	
 
void mexFunction(int nlhs, mxArray *plhs[],
        int nrhs, const mxArray *prhs[])
{
    int *pr_seed;
 
    pr_seed = (int *)mxGetData(prhs[0]);
 
    srand(*pr_seed);
 
    outfile = fopen("Out", "w");
 
    inputs();
    connections();
    initiate();
    evaluate();
    store();
    fclose(outfile);
}
 
int bi(x, threshold)
   float x, threshold;
{
    if (x > threshold) return (1);
    else return (0);
}

A compiler comme ceci :

Code :

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

mex segment.c

Et à utiliser comme ceci :

Code :

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

1
2
seed = int32(5);
segment(seed);

**zied113** · 22/05/2014, 13h43

Merci Mr.Jérôme Briot pour votre aide , mais je vous demande plus de détailles car je suis débutant en programmation ,

est ce que je dois convertie le code .txt en code .c puis en MEX ??? Merci d’avance

**Jerome Briot** · 22/05/2014, 16h39

Comme mentionné dans la FAQ, il te faut tout d'abord un compilateur pour compiler le fichier MEX.

Quelle version de MATLAB utilises-tu ? Avec quel système d'exploitation ? 32 ou 64 bits ?

Voir la FAQ Quels sont les différents compilateurs supportés pour compiler un fichier MEX ? pour choisir un compilateur.

**zied113** · 22/05/2014, 16h49

MATLAB 7.10.0 (R2010a) , windows 7 ,32 bits

**Jerome Briot** · 22/05/2014, 17h24

OK, donc tu suis les instructions données ici Comment choisir le compilateur à utiliser pour créer un fichier MEX ? et tu choisis le compilateur Lcc dans la liste proposée.

**zied113** · 22/05/2014, 21h22

Code :

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

1
2
3
4
Please verify your choices: 
 
Compiler: Lcc-win32 C 2.4.1 
Location: C:\PROGRA~1\MATLAB\R2010a\sys\lcc

J'ai essai , matlab m'affiche ce message :

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
Are these correct [y]/n? y 
 
Trying to update options file: C:\Users\zied\AppData\Roaming\MathWorks\MATLAB\R2010a\mexopts.bat 
From template:              C:\PROGRA~1\MATLAB\R2010a\bin\win32\mexopts\lccopts.bat 
 
Done . . . 
 
************************************************************************** 
  Warning: The MATLAB C and Fortran API has changed to support MATLAB 
           variables with more than 2^32-1 elements.  In the near future 
           you will be required to update your code to utilize the new 
           API. You can find more information about this at: 
           http://www.mathworks.com/support/sol...ution=1-5C27B9 
           Building with the -largeArrayDims option enables the new API. 
**************************************************************************

**Jerome Briot** · 22/05/2014, 22h05

Donc maintenant tu enregistres le code C que j'ai fourni précédemment dans un fichier segment.c

Tu copies le fichier segment.h dans le même dossier que segment.c puis tu fais :

Code :

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

mex segment.c

**zied113** · 22/05/2014, 23h07

j'ai copier les deux fichier segment .c + segment.h a "C:\PROGRA~1\MATLAB\R2010A\BIN\"
j'ai executé la commande mex segment.c

Matlab affiche ce message

C:\PROGRA~1\MATLAB\R2010A\BIN\MEX.PL: Error: 'segment.c' not found.

??? Error using ==> mex at 222
Unable to complete successfully.

:''(

**Jerome Briot** · 23/05/2014, 08h13

Envoyé par zied113

j'ai copier les deux fichier segment .c + segment.h a "C:\PROGRA~1\MATLAB\R2010A\BIN\"

Premièrement, il ne faut pas utiliser le fichier segment.c que tu as téléchargé sur le site du "Perception and Neurodynamics Laboratory" mais celui que j'ai modifié dans ce message . As-tu au moins remarqué les quelques différences entre les deux codes… j'en doute.

Deuxièmement, il ne faut pas copier les fichiers dans le dossier d'installation de MATLAB.
Place les, par exemple, dans le dossier dans lequel MATLAB démarre.

Pour finir, il faut absolument que le dossier où se trouvent ces deux fichiers soit le dossier courant de MATLAB avant de les compiler avec mex

**zied113** · 23/05/2014, 12h02

Bonjour , normalement c'est bon , matlab m’affiche le message suivant

Code :

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

1
2
3
4
5
6
7
>> mex segment.c
seed = int32(5);
segment(seed);
we are here 
.......
we are here 
The number of the segments/cycles is 348

Mais la question qui ce pose comment modifier l'Input de segment.c avec mon Input ??????? !!!!!

en attendant monsieur votre réponse veuillez agréer mes salutations les meilleures.

**Jerome Briot** · 23/05/2014, 12h16

Envoyé par zied113

Mais la question qui ce pose comment modifier l'Input de segment.c avec mon Input ??????? !!!!!

C'est à dire ?

Tu veux appliquer le code sur une autre image que celle mentionnée ici ?

Code C :

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

IN = fopen("lake1.pgm", "r");

**zied113** · 23/05/2014, 12h47

oui exactement c'est ca

**Jerome Briot** · 23/05/2014, 13h48

Voici une version :

Code C :

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
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
       /*****************************************************************/
       /* This program is for 2D segmentation using the LEGION network  */
       /* The equations were modeled as a PDP algorithm for efficiency  */
       /* by DeLiang Wang in 7/94. For a detailed description of the    */
       /* algorithm see D.L. Wang and D. Terman: "Image segmentation    */
       /* based on oscillatory correlation", Neural Computation, Vol. 9 */
       /* pp. 805-836, 1997 (For errata see Neural Computation 9,       */
       /* 1623-1626, 1997), on page 822.                                */
       /*****************************************************************/
 
#include "mex.h"
#include "segment.h"
 
#define MAXVAL 32768.0
               /* 2147483648.0 or 32768.0,machine (HP or SUN) dependable*/
#define MAX 32767 /* 2147483647 or 32767, machine dependable */
#define GRAY 255  /* Maximum value of a pixel */
 
int bi();
int jump[2], z, Num = 1;             /* so we can change step size at will */
                               /* Num: number of the segments */
FILE *outfile, *fopen();
float x[N1][N2][2], lateral[N1][N2];
                      /* lateral: the sum of input from lateral connections */
float W[N1][N2][8];   /* strength of the interactions for 4 directions */
int stim[N1][N2], matrix[N1][N2];
int y[N1][N2], cycle[N1][N2], trigger[N1][N2];
     /* y: integer of x; cycle: test for 1st cycle; trigger: jumping signal */
 
 
int inputs(char *filename)
{
   int i, j;
   int width, height, c;
   char magic[25];
   float stimulus = 0.2;
   FILE *IN;
 
   for (i = 0; i < N1; i++) 
     for (j = 0; j < N2; j++)
       matrix[i][j] = 0;
 
   IN = fopen(filename, "r");
   if(IN==NULL)
       return -1;
 
   fscanf(IN,"%s\n%d %d\n%d\n",magic,&width,&height,&c);
   fprintf(stderr,"[Type %s %dx%d colours %d]\n",magic,width,height,c);
 
   for(i=0;i<160;i++)  
      for(j=0;j<160;j++)
         stim[i][j] = getc(IN);
 
   fclose(IN);
 
   return 0;
}
 
void connections(void)
{
   int i, j, k;
   float baseline = 0.0, Wtotal = 6.0, sum;
 
   for (i = 0; i < N1; i++) /* for eight nearest neighbors */
     for (j = 0; j < N2; j++) {
       W[i][j][0] = W[i][j][1] = W[i][j][2] = W[i][j][3] = 0.0;
       W[i][j][4] = W[i][j][5] = W[i][j][6] = W[i][j][7] = 0.0;
     }
 
   for (i = 0; i < N1; i++) 
     for (j = 0; j < N2; j++) {
	 sum = 0.0;
 
         if (j-1 >= 0)
	    W[i][j][0] = (float) GRAY/(1+abs(stim[i][j]-stim[i][j-1]));
         if (i-1 >= 0 && j-1 >= 0)
	    W[i][j][1] = (float) GRAY/(1+abs(stim[i][j]-stim[i-1][j-1]));
         if (i-1 >= 0)
	    W[i][j][2] = (float) GRAY/(1+abs(stim[i][j]-stim[i-1][j]));
         if (i-1 >= 0 && j+1 < N2)
	    W[i][j][3] = (float) GRAY/(1+abs(stim[i][j]-stim[i-1][j+1]));
         if (j+1 < N2)
	    W[i][j][4] = (float) GRAY/(1+abs(stim[i][j]-stim[i][j+1]));
         if (i+1 < N1 && j+1 < N2) 
	    W[i][j][5] = (float) GRAY/(1+abs(stim[i][j]-stim[i+1][j+1]));
         if (i+1 < N1)
	    W[i][j][6] = (float) GRAY/(1+abs(stim[i][j]-stim[i+1][j]));
         if (i+1 < N1 && j-1 >= 0)
	    W[i][j][7] = (float) GRAY/(1+abs(stim[i][j]-stim[i+1][j-1]));
 
         for (k = 0; k < 8; k++) lateral[i][j] += W[i][j][k];
    }
}
 
void initiate(void)
{
   int i,j,t, temp;  /* labels for setting up display intervals */
 
   for (i = 0; i < N1; i++)   /* randomize the initial phases */
     for (j = 0; j < N2; j++) {
        y[i][j] = rand();
        x[i][j][0] = y[i][j]/(MAXVAL + 1.0);
	y[i][j] -= 2;
	if (y[i][j] < 0)
	  y[i][j] = 0;
	trigger[i][j] = 0;
	cycle[i][j] = 0;
     }
   z = 0; 
 
 /* initial display values for 2D */
   for (i = 0; i < N1; i++)
      for (j = 0; j < N2; j++) {
        temp = x[i][j][0]*RP;       
        fprintf(outfile, "%d ", temp);
      }
}
 
void select(void)
{
   int i, j, max, imax=0, jmax=0;
   float baseline = 1200.0, temp;
/* baseline 1200 -> 17 regions; 1000 -> 23 regions. Both reasonable */
 
   max = -1;
   for (i = 0; i < N1; i++)
     for (j = 0; j < N2; j++) {
       trigger[i][j] = 0;
       if (lateral[i][j] > baseline && y[i][j] != MAX && y[i][j] > max) {
	 max = y[i][j];
	 imax = i;
	 jmax = j;
       }
     }
 
   if (cycle[imax][jmax]) {     /* find out jumpers */
     for (i = 0; i < N1; i++)
       for (j = 0; j < N2; j++)
	 if (y[i][j] == max && cycle[i][j])
	   trigger[i][j] = 1;
   } else
     trigger[imax][jmax] = 1;
 
   for (i = 0; i < N1; i++)   /* advance oscillators on the left branch */
     for (j = 0; j < N2; j++) 
       if (y[i][j] != MAX && (cycle[i][j] || lateral[i][j] > baseline)) {
	 y[i][j] = y[i][j] + (MAX - max);
         x[i][j][0] = y[i][j]/(MAXVAL + 1.0);
       }
}
 
void evaluate(void)
{
   float temp, thetae = 2.0, temp2, temp3;
   int i,j,t, cat, count;  /* labels for setting up display intervals */
 
   jump[0] = 0;
 
   select();
   for ( t = 0; t < WIDTH-1; t++) {
      count = 0;
      for (i = 0; i < N1; i++)
	for (j = 0; j < N2; j++) { /* compute net input */
	  if (x[i][j][0] > thetae && jump[0] == 0) {  /* jumping down */
             x[i][j][1] = LP;
	     y[i][j] = -1;
	     cycle[i][j] = 1;
	     count++;
	  } else if (x[i][j][0] < thetae) {  /* summing up overall input */
             temp =  -20.0 * bi((float) z,0.5);  /* used to be 50 */
 
/* the following is based on traversing through strongest link */
             temp2 = 0.0;
             if (j-1 >= 0) {
                temp3 = W[i][j][0] * bi(x[i][j-1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i-1 >= 0 && j-1 >= 0) {
                temp3 = W[i][j][1] * bi(x[i-1][j-1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i-1 >= 0) {
                temp3 = W[i][j][2] * bi(x[i-1][j][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i-1 >= 0 && j+1 < N2) {
                temp3 = W[i][j][3] * bi(x[i-1][j+1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (j+1 < N2) {
                temp3 = W[i][j][4] * bi(x[i][j+1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i+1 < N1 && j+1 < N2) {
                temp3 = W[i][j][5] * bi(x[i+1][j+1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i+1 < N1) {
                temp3 = W[i][j][6] * bi(x[i+1][j][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             if (i+1 < N1 && j-1 >= 0) {
                temp3 = W[i][j][7] * bi(x[i+1][j-1][0], thetae);
                if (temp2 < temp3) temp2 = temp3; }
             temp += temp2;
 
	     if (trigger[i][j] || temp > 0.25) {        /* jumping up */
	        x[i][j][1] = RP; 
		y[i][j] = MAX;
	        jump[1] = 1;
	        z++;
	     } else
	        x[i][j][1] = x[i][j][0];
	  } else 
	     x[i][j][1] = x[i][j][0];
	}
      z -= count;
 
      if (jump[1] == 0 && z == 0) 
        select();
 
    /* store for display in 2d case*/
      if (jump[1] == 0 && jump[0] == 1) {
mexPrintf("we are here \n");
	Num++;
        for (i = 0; i < N1; i++)
          for (j = 0; j < N2; j++) {
             if (z > 0) {
	          if (x[i][j][1] < thetae)
                     cat = 0;       
		  else cat = x[i][j][1];       
	     } else  
		  cat = x[i][j][1];       
             fprintf(outfile, "%d ", cat);
 
             if (Num < 25 && cat > 0)
                matrix[N1-1-i][j] = 26-Num;  /* 17 regions, or 23 */
	  }
      }
 
 
      for (i = 0; i < N1; i++)
         for (j = 0; j < N2; j++)  /* replace OLD by NEW */
           x[i][j][0] = x[i][j][1];
      jump[0] = jump[1];
      jump[1] = 0;
   }
}
 
void store(void)
{
   int i, j;
   FILE *out, *out1, *fopen();
 
   mexPrintf("The number of the segments/cycles is %d\n", Num);
 
   out = fopen("NUMBER", "w");
   out1 = fopen("Result", "w");
   fprintf(out, "%d ", Num);
   for (i = 0; i < N1; i++)
     for (j = 0; j < N2; j++)
       fprintf(out1, "%d ", matrix[i][j]);
   fclose(out);
   fclose(out1);
}	
 
void mexFunction(int nlhs, mxArray *plhs[],
        int nrhs, const mxArray *prhs[])
{
    int *pr_seed, status;
    mwSize buflen;
    char *buf;
 
    pr_seed = (int *)mxGetData(prhs[0]);
 
    buflen = mxGetM(prhs[1])*mxGetN(prhs[1]);
    buf = mxCalloc(buflen, sizeof(char));
 
    mxGetString(prhs[1], buf, buflen+1);
 
    srand(*pr_seed);    
 
    status = inputs(buf);
 
    if(status) {        
        mxFree(buf);
        mexErrMsgTxt("Pb with file");
    }        
 
    outfile = fopen("Out", "w");
    connections();
    initiate();
    evaluate();
    store();
    fclose(outfile);
 
    mxFree(buf);
}
 
int bi(x, threshold)
   float x, threshold;
{
    if (x > threshold) return (1);
    else return (0);
}

A utiliser comme ceci :

Code :

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

1
2
seed = int32(5);
segment(seed,'lake1.pgm');

**zied113** · 23/05/2014, 16h22

mon entrée dans la fonction segmentation est une image de "Canal de corrélation croisée" d'un signal mixer, le Canal corrélation croisée est définie dans matlab par cette fonction :

Code :

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

cross = crossChanCorr(c)

jai tapé la code suivant

Code :

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

1
2
3
mex segment.c
seed = int32(5);
segment(seed,cross);

j'ai obtenu le message suivant :

Code :

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

1
2
we are here 
The number of the segments/cycles is 695

je sais pas est ce que juste ou nn !!!!!!

vous trouvez ci joint limage de "Canal de corrélation croisée" avec le système.

**Jerome Briot** · 23/05/2014, 18h28

La ligne suivante n'a pas besoin d'être répétée à chaque fois :

Code :

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

mex segment.c

Sauf si tu modifies le contenu du fichier segment.c

Pour le reste, tout ce que je peux dire c'est que le code génère trois fichiers : Out, NUMBER et Result

**zied113** · 25/05/2014, 13h18

Bonjour M. Jérôme Briot , j'ai téléchargé le code de groupe.c et groupe.h par le lien : http://www.cse.ohio-state.edu/pnl/sh...tnn99/dynamic/

J'ai essayé de le exécuté sur matlab par la code :

Code :

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

mex group.c

matlab affiche l'erreur suivante

Code :

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

1
2
3
4
5
6
7
8
9
10
mex group.c
c:\users\zied\appdata\local\temp\mex_kf~1\group.obj: multiple definition of _main 
first definition is in file c:\progra~1\matlab\r2010a\sys\lcc\lib\libcrt0.obj 
Specified export _mexFunction is not defined 
Missing exports. Aborting 
 
  C:\PROGRA~1\MATLAB\R2010A\BIN\MEX.PL: Error: Link of 'group.mexw32' failed. 
 
??? Error using ==> mex at 222
Unable to complete successfully.

Invité · 25/05/2014, 13h41

Envoyé par Dut

Premièrement, il ne faut pas utiliser le fichier segment.c que tu as téléchargé sur le site du "Perception and Neurodynamics Laboratory" mais celui que j'ai modifié dans ce message . As-tu au moins remarqué les quelques différences entre les deux codes… j'en doute.

Différences sûrement toujours pas remarquées...
Il serait temps de lire minutieusement le contenu de ce lien afin de bien comprendre la logique des fichiers MEX qui ne passent pas par une fonction d'entrée main() mais mexFunction() !

**zied113** · 26/05/2014, 17h18

bonjour , j'ai pas trouvé la bonne explication dans ce lien http://matlab.developpez.com/faq/?page=mex

**Jerome Briot** · 26/05/2014, 21h38

Ouvre les deux fichiers segment.c (l'original et celui que j'ai modifié) côte à côte dans un éditeur de texte.

Compare les lignes par lignes pour bien identifier les différences, dont la plus notable est, comme l'a fait remarqué Winjerome, l'absence de main dans le MEX

convertir code C en MATLAB

MATLAB

Discussions similaires

Partager

Partager