Discussion :

[svagrim]

Bonjour,

je reprend depuis peu le C++ après de longues années où je n'en avais pas fait et j'ai un peu du mal.

depuis peu j'ai décidé d'utiliser des classes pour faire ce que j'avais commencé.
il m'a semblé que ça me faciliterais grandement la vie pour le développement futur du programme.
mais pour le moment ça me le facilite pas du tout ...

alors mon problème:
j'ai une classe que j'appelle Images.
je veux faire un programme membre dans cette classe qui "lit" les images pour les stocker dans un tableau dynamique de 5 dimensions (x, y, z, temps, énergie).

du coup dans mon main.cpp j'ai mis ça pour le moment :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
#include <cstdlib>
#include <iostream>
 
#include "read_vff.h"
#include "test_class.h"
 
using namespace std;
 
int main(int argc, char *argv[])
{      
Images.reading();
std::cin.get();
}

dans mon test_class.h j'ai ça :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
#ifndef TEST_CLASS_H
#define TEST_CLASS_H
 
/*reading the vff files, header and images, and sort them in the variable ... */
 
class Images
{   
public:
   void reading(vector<std::string> &files);        
   vector<std::string> &ttype;
   //déclatration tableau dynamique 5 dimensions
 
};

et dans mon test_class.cpp j'ai ça :

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
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
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <iostream>
#include <fstream>
#include <string>
#include <sstream>
#include <sys/types.h>
#include <dirent.h>
#include <errno.h>
#include <vector>
#include <conio.h>
#include <algorithm>
#include <functional>
#include <iterator>
 
#include "test_class.h"
using namespace std;
 
/*function getdir2*/
int getdir2 (string dir, vector<string> &files)
{
  DIR *dp;
  struct dirent *dirp;
  if((dp = opendir(dir.c_str())) == NULL) {
  cout << "Error(" << errno << ") opening " << dir << endl;
  return errno;
}
 
while ((dirp = readdir(dp)) != NULL) 
{
  files.push_back(string(dirp->d_name));
}
 
closedir(dp);
return 0;
}
// end function getdir2
 
//function main
void Images::reading()
{
vector<string> files = vector<string>();
string input = ""; 
 
// How to get a string/sentence with spaces
cout << "Enter the path of the folder where there is only the projections :\n>";
getline(cin, input);
 
getdir2(input,files);
 
//int const tailleTableau(20);
//cout << tailleTableau << endl;
for (unsigned int i = 2;i < files.size();i++) 
{
  cout << files[i] << endl; 
}
 
char *ttype[files.size()];
char *tsize[files.size()];
char *tbits[files.size()];
char *trank[files.size()];
char *tbands[files.size()];
char *ty_bin[files.size()];
char *tz_bin[files.size()];
char *trfan_y[files.size()];
char *trfan_z[files.size()];
char *tmin[files.size()];
char *tmax[files.size()];
char *twater[files.size()];
char *tair[files.size()];
char *tbonehu[files.size()];
char *tangle_increment[files.size()];
char *tcenter_of_rotation[files.size()];
char *tcentral_slice[files.size()];
char *tspacing[files.size()];
char *torigin[files.size()];
char *tformat[files.size()];
char *telementsize[files.size()];
char *tmodality[files.size()];
char *tdate[files.size()];
char *ttitle[files.size()];
char *tsubject[files.size()];
char *treverse_order[files.size()];
char *trecon_sysid[files.size()];
char *tmacaddr[files.size()];
char *tbinning[files.size()];
char *trawsize[files.size()];
char *tcurrent[files.size()];
char *tcurrentUnit[files.size()];
char *tdataSourceJID[files.size()];
char *texam[files.size()];
char *texposure[files.size()];
char *tframesAverage[files.size()];
char *tgantryPosition[files.size()];
char *tincrementAngle[files.size()];
char *toffset[files.size()];
char *tproject[files.size()];
char *tprotocol[files.size()];
char *tscanTechnique[files.size()];
char *tscanType[files.size()];
char *tsequence[files.size()];
char *tserialNumber[files.size()];
char *tseries[files.size()];
char *tseriesVersion[files.size()];
char *ttablePosition[files.size()];
char *tvoltage[files.size()];
char *tvoltageUnit[files.size()];
 
for (unsigned int i = 2;i < files.size();i++) 
{
 
  /* Read projection file*/
 
  string sep = "\\";
  string file_to_open = input+sep+files[i];
 
  ifstream input_file;
 
  input_file.open (file_to_open.c_str(), ifstream::in);
  input_file.seekg(0,ios::beg);
 
  if(input_file == NULL)
  {
    std::cerr<<"Cannot open input file \n";
  }
 
  /* Read the header */
  std::string line;
 
  cout << endl;
  cout << file_to_open << endl;
  cout << "line :" << std::getline(input_file, line )<< endl;
 
  string test ("=");
  size_t found;
 
 
  while ( std::getline(input_file, line ) )
  {
    found = line.find(test);
    cout << "found : " << int(found) << endl;
 
    string name_info;
    string info;
    name_info = line.substr(0,found);
    info = line.substr(found+1,(line.length()-found-2));
    cout << (line.length()-found-2) << endl;
    cout << name_info << endl;
    cout << "info :" << info << endl;
 
 
 
    if(!strcmp(name_info.c_str(), "type" ))
    {
      cout << "1"; //à la place ça enregistrait l'info dans le tableau dynamique déclaré dans la classe.
    }
    else if(!strcmp(name_info.c_str(), "size" ))
    {
      cout << "2";
    }
    else if(!strcmp(name_info.c_str(), "bits" ))
    {
      cout << "3";
    }
    else if(!strcmp(name_info.c_str(), "rank" ))
    {
      cout << "4";
    }
    else if(!strcmp(name_info.c_str(), "bands" ))
    {
      cout << "5";
    }
    else if(!strcmp(name_info.c_str(), "y_bin" ))
    {
      cout << "6";
    }
    else if(!strcmp(name_info.c_str(), "z_bin" ))
    {
      cout << "7";
    }
    else if(!strcmp(name_info.c_str(), "rfan_y" ))
    {
      cout << "8";
    }
    else if(!strcmp(name_info.c_str(), "rfan_z" ))
    {
      cout << "9";
    }
    else if(!strcmp(name_info.c_str(), "min" ))
    {
      cout << "10";
    }
    else if(!strcmp(name_info.c_str(), "max" ))
    {
      cout << "11";
    }
    else if(!strcmp(name_info.c_str(), "water" ))
    {
      cout << "12";
    }
    else if(!strcmp(name_info.c_str(), "air" ))
    {
      cout << "13";
    }
    else if(!strcmp(name_info.c_str(), "bonehu" ))
    {
      cout << "14";
    }
    else if(!strcmp(name_info.c_str(), "angle_increment" ))
    {
      cout << "15";
    }
    else if(!strcmp(name_info.c_str(), "center_of_rotation" ))
    {
      cout << "16";
    }
    else if(!strcmp(name_info.c_str(), "central_slice" ))
    {
      cout << "17";
    }
    else if(!strcmp(name_info.c_str(), "spacing" ))
    {
      cout << "18";
    }
    else if(!strcmp(name_info.c_str(), "origin" ))
    {
      cout << "19";
    }
    else if(!strcmp(name_info.c_str(), "format" ))
    {
      cout << "20";
    }
    else if(!strcmp(name_info.c_str(), "elementsize" ))
    {
      cout << "21";
    }
    else if(!strcmp(name_info.c_str(), "modality" ))
    {
      cout << "22";
    }
    else if(!strcmp(name_info.c_str(), "date" ))
    {
      cout << "23";
    }
    else if(!strcmp(name_info.c_str(), "title" ))
    {
      cout << "24";
    }
    else if(!strcmp(name_info.c_str(), "subject" ))
    {
      cout << "25";
    }
    else if(!strcmp(name_info.c_str(), "reverse_order" ))
    {
      cout << "26";
    }
    else if(!strcmp(name_info.c_str(), "recon_sysid" ))
    {
      cout << "27";
    }
    else if(!strcmp(name_info.c_str(), "macaddr" ))
    {
      cout << "28";
    }
    else if(!strcmp(name_info.c_str(), "binning" ))
    {
      cout << "29";
    }
    else if(!strcmp(name_info.c_str(), "rawsize" ))
    {
      cout << "30";
    }
    else if(!strcmp(name_info.c_str(), "current" ))
    {
      cout << "31";
    }
    else if(!strcmp(name_info.c_str(), "currentUnit" ))
    {
      cout << "32";
    }
    else if(!strcmp(name_info.c_str(), "dataSourceJID" ))
    {
      cout << "33";
    }
    else if(!strcmp(name_info.c_str(), "exam" ))
    {
      cout << "34";
    }
    else if(!strcmp(name_info.c_str(), "exposure" ))
    {
      cout << "35";
    }
    else if(!strcmp(name_info.c_str(), "framesAverage" ))
    {
      cout << "36";
    }
    else if(!strcmp(name_info.c_str(), "gain" ))
    {
      cout << "37";
    }
    else if(!strcmp(name_info.c_str(), "gantryPosition" ))
    {
      cout << "38";
    }
    else if(!strcmp(name_info.c_str(), "ncrementAngle" ))
    {
      cout << "39";
    }
    else if(!strcmp(name_info.c_str(), "offset" ))
    {
      cout << "40";
    }
    else if(!strcmp(name_info.c_str(), "project" ))
    {
      cout << "41";
    }
    else if(!strcmp(name_info.c_str(), "protocol" ))
    {
      cout << "42";
    }
    else if(!strcmp(name_info.c_str(), "scanTechnique" ))
    {
      cout << "43";
    }
    else if(!strcmp(name_info.c_str(), "scanType" ))
    {
      cout << "44";
    }
    else if(!strcmp(name_info.c_str(), "sequence" ))
    {
      cout << "45";
    }
    else if(!strcmp(name_info.c_str(), "serialNumber" ))
    {
      cout << "46";
    }
    else if(!strcmp(name_info.c_str(), "series" ))
    {
      cout << "47";
    }
    else if(!strcmp(name_info.c_str(), "seriesVersion" ))
    {
      cout << "48";
    }
    else if(!strcmp(name_info.c_str(), "tablePosition" ))
    {
      cout << "49";
    }
    else if(!strcmp(name_info.c_str(), "voltage" ))
    {
      cout << "50";
    }
    else if(!strcmp(name_info.c_str(), "voltageUnit" ))
    {
      cout << "51";
    }
 
    else
    { 
      cout << "qwerty";
    }
 
    // afficher la ligne à l'écran
    std::cout << line << std::endl;   
  }
 
  cout << "blurp" << endl;
  std::cin.get();
}
}

et j'ai un problème de compilation.
il ne veut pas de mon void reading(vector<std::string> &files);
(il veut pas du vector<std::string> &files) plus exactement.

quelqu'un verrait le problème ?
(et si quelqu'un a une idée de comment faire un tableau dynamique 5 dimensions au passage ...)

Merci d'avance.

[Bousk]

Bonsoir,

vector, tout comme string, se trouve dans le namespace std.

[svagrim]

tu parles pour le point h c'est ça ?
j'ai vu sur les tutoriels C++ qu'il fallait pas mettre de std dans le fichier de la classe.

Dans les .h, il est recommandé de ne jamais mettre la directive using namespace std; car cela pourrait avoir des effets néfastes, par la suite, lorsque vous utiliserez la classe.
Par conséquent, il faut rajouter le préfixe std:: devant chaque string du .h. Sinon, le compilateur vous sortira une erreur du type string does not name a type.

du coup tu crois que c'est vraiment ça le problème ?
j'ai un peu peur de potentiel futur effet nefastes. c'est quoi ces effets néfastes dont ils parlent ?

[Médinoc]

La réponse n'est pas mettre un using namespace, mais écrire std::vector<std::string>

[svagrim]

ah oki j'avais pas compris.
merci j'essaye ça.

[svagrim]

nan, toujours la même erreur :

In file included from main.cpp:5:
test_class.h:9: error: variable or field `reading' declared void
test_class.h:9: error: expected `;' before '(' token

main.cpp: In function `int main(int, char**)':

main.cpp:12: error: expected primary-expression before '.' token

[Neckara]

Bonjour,

Il me semble que personne n'a encore répondu à cette question que tu te posais :

j'ai un peu peur de potentiel futur effet nefastes. c'est quoi ces effets néfastes dont ils parlent ?

Si tu met un using namespace X dans un header, tous les fichiers incluant cet header auront le using namespace X.
Or il arrive que dans les bibliothèques ou dans le code qu'on écrit, de se retrouver avec des noms de classes identiques (ex : X::String et Y::String respectivement contenu dans la bibliothèque X et Y).
Et si tu fais un using namespace X et un using namespace Y, comment savoir si String provient de la bibliothèque X ou de la bibliothèque Y ?

Ainsi si tu met un using namespace X dans un header, tu ne pourras plus faire de using namespace Y (ou de header le faisant) dans les fichiers qui inclueront ton header.

[svagrim]

merci Neckara pour les infos.

et merci JolyLoic pour l'aide.

concernant le reading.
Je suis d'accort avec toi sur l'idée mais j'ai fait que suivre la façon dont dans mon labo ils nomment leur fonction.
ils sont canadiens anglophones (je suis sur vancouver) donc j'imagine que c'est pas un souci "reading" vu qu'ils les appellent comme ça.
en fait je pense que c'est plus pour être sur que c'est une fonction programmer à la main et pas un appel d'une fonction de la bibliothèque par defaut de matlab/c++ ou autre langage qu'ils utilisent.

pour mon problème oui j'avais vu que j'avais zappé de définir le constructeur.

sinon effectivement je reprend pas le c comme ça.
j'ai fait une bonne partie des tutorials des sites du zero.

mais je me melange encore les pinceaux.
faire ce programme reprèsente environ 10% de mon taf.
du coup je commence, je peux pas m'y remettre pendant un bon moment.
j'y reviens et j'ai oublié des trucs que j'ai appris parce que pas utilisé de suite après.
bref pas l'idéal pour un langage que j'ai pas utilisé depuis longtemps.
mais ça commence à revenir petit à petit.

merci pour l'avertissement pour les **

sinon oui tu as bien compris pour mon tableau.
ils ont la même taille. et l'indice correspond à l'indice d'une image en fait. les tableaux sont les tableaux avec les informations sur l'acquisition et la reconstruction de l'image.

pour le tableau à 5 dimensions.
les indices sont tous des entiers (sous matlab j'avais des indices entiers et un vecteur pour une correspondance indice entier/valeur réelle si c'est pas un entier) c'était le plus simple.

par contre je vois pas ce que tu veux dire par creux ou plein.
Il n'y a aucune partie vide si c'est ce que tu veux savoir.
il y a des valeurs dans tous les éléments de la matrice.

[JolyLoic]

par contre je vois pas ce que tu veux dire par creux ou plein.
Il n'y a aucune partie vide si c'est ce que tu veux savoir.
il y a des valeurs dans tous les éléments de la matrice.

C'était bien ça. Si je demande, c'est parce que même s'il n'y a que 100 indices par dimension, le nombre total de cases augmente très vite (100^5 == 10^10 = environ 75Go si tu travaille avec des doubles = il va falloir faire un peu attention à la manière dont la mémoire est gérée).

Et donc, si la grande majorité des cases étaient avec une même valeur, il pourrait être intéressant de stocker les informations dans un std::map<tuple<int, int, int, int, int>, double>, en ne codant que les cases n'ayant pas cette valeur.

Sinon, si tu est certain d'avoir assez de mémoire contigüe pour tout stocker, un très bourrin vector<double> avec l'accès à une case du type v[x + y*nbX + z *nbX*nbY + t*nbX*nbY*nbZ + e*nbX*nbY*nbZ*nbT]Après, tu peux fragmenter un peu, par exemple un vector<vector<double>>, le premier niveau représentant e et t, le second x,y,z.

Tous ça dépend des calculs que tu vas vouloir faire. Je te conseille en particulier de te débrouiller pour que tes boucles internes soient de préférence sur des données stockées physiquement proches les unes des autres dans ton tableau. C'est important pour que le cache joue bien son rôle, et d'autant plus si tu dois paralléliser ton algorithme un jour.

Et j'ai failli oublié : Comme tu vois, il y a plein de variantes possible, selon les niveaux de découpage, les collections choisies à chaque niveau, l'ordre des indices..., et tu va probablement devoir jouer avec tout ça. Donc n'utilise pas directement les structures que je t'ai conseillées, mais wrap les dans une classe avec une interface unique, mais qui te permettra de tester diverses combinaisons.