erreur à la recupération de contenu d'un vecteur

**tasnim2010** · 25/09/2012, 18h45

Bonsoir,

Bon, je vais essayer de cous simplifier la description de ma tache. J'ai un objet 3D, chaque point de cet objet est defini par ses coordonnées dans l'espace 3D et un identificateur, chaque 3 points sont reliés entre eux pour former un triangle. j'ai appliqué l'octree pour segmenter mon objet en des regions. Maintenant, je veux recupérer pour chaque région, les coordonnées des points qui y appartient et leurs identificateurs, puis les ecrires dans un fichiers. Pour eviter d'avoir des redondances d'identificateurs dans mon fichier, j'ai crée un vecteur, en parcourant la liste de triangles (nommés face aussi ) de mon objet, si cette face appartient à mon région, je recupère ses trois sommets, je verifie si ces sommets existent déjà dans mon vecteur ou non, le sommet qui n'existe pas dans mon vecteur va etre ajouter à ce vecteur (on ajoute son identificateur) etc.En meme temps, j'ai declaré une variable pour chaque région (idss1 , idss2 ...) dans laquelle , je vais compter le nombre de case du tableau, une autre variable ( f1, f2,...) pour compter le nombre des faces (triangles) dans chaque région. Jusqu'au ici, tout va bien avec ce code:

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
for (int i = 0; i < this->Input->GetNumberOfCells(); i++) {
    vtkIdType v1;
    vtkIdType v2;
    vtkIdType v3;
 
    this->Input->GetFaceVertices(i, v1, v2, v3);
    int dd = PointsNode->GetValue(i);
    switch (dd) {
    case (1):
	{
	    int vv1 = 0;
	    int vv2 = 0;
	    int vv3 = 0;
	    for (int j = 1; j <= ids1.size(); j++) {
		if (ids1[j] == v1)
		    vv1++;
		if (ids1[j] == v2)
		    vv2++;
		if (ids1[j] == v3)
		    vv3++;
	    }
	    if (vv1 == 0) {
		ids1.push_back(v1);
		idss1++;
		cout << "v1 = " << v1 << "n" << endl;
	    }
	    if (vv2 == 0) {
		ids1.push_back(v2);
		idss1++;
		cout << "v2 = " << v2 << "n" << endl;
	    }
	    if (vv3 == 0) {
		ids1.push_back(v3);
		idss1++;
		cout << "v3 = " << v3 << "n" << endl;
	    }
	    cout << "idss1 = " << idss1 << "n" << endl;
 
	    f1++;
	}
    case (2):
	{
	    int vv1 = 0;
	    int vv2 = 0;
	    int vv3 = 0;
	    for (int j = 1; j <= ids1.size(); j++) {
		if (ids1[j] == v1)
		    vv1++;
		if (ids1[j] == v2)
		    vv2++;
		if (ids1[j] == v3)
		    vv3++;
	    }
	    if (vv1 == 0) {
		ids1.push_back(v1);
		idss2++;
	    }
	    if (vv2 == 0) {
		ids1.push_back(v2);
		idss2++;
	    }
	    if (vv3 == 0) {
		ids1.push_back(v3);
		idss2++;
	    }
	    f2++;
	}
    case (3):
	{
	    int vv1 = 0;
	    int vv2 = 0;
	    int vv3 = 0;
	    for (int j = 1; j <= ids3.size(); j++) {
		if (ids3[j] == v1)
		    vv1++;
		if (ids3[j] == v2)
		    vv2++;
		if (ids3[j] == v3)
		    vv3++;
	    }
	    if (vv1 == 0) {
		ids3.push_back(v1);
		idss3++;
	    }
	    if (vv2 == 0) {
		ids3.push_back(v2);
		idss3++;
	    }
	    if (vv3 == 0) {
		ids3.push_back(v3);
		idss3++;
	    }
	    f3++;
	}
    case (4):
	{
	    int vv1 = 0;
	    int vv2 = 0;
	    int vv3 = 0;
	    for (int j = 1; j <= ids4.size(); j++) {
		if (ids4[j] == v1)
		    vv1++;
		if (ids4[j] == v2)
		    vv2++;
		if (ids4[j] == v3)
		    vv3++;
	    }
	    if (vv1 == 0) {
		ids4.push_back(v1);
		idss4++;
	    }
	    if (vv2 == 0) {
		ids4.push_back(v2);
		idss4++;
	    }
	    if (vv3 == 0) {
		ids4.push_back(v3);
		idss4++;
	    }
	    f4++;
	}
    case (5):
	{
	    int vv1 = 0;
	    int vv2 = 0;
	    int vv3 = 0;
	    for (int j = 1; j <= ids3.size(); j++) {
		if (ids5[j] == v1)
		    vv1++;
		if (ids5[j] == v2)
		    vv2++;
		if (ids5[j] == v3)
		    vv3++;
	    }
	    if (vv1 == 0) {
		ids5.push_back(v1);
		idss5++;
	    }
	    if (vv2 == 0) {
		ids5.push_back(v2);
		idss5++;
	    }
	    if (vv3 == 0) {
		ids5.push_back(v3);
		idss5++;
	    }
	    f5++;
	}
    case (6):
	{
	    int vv1 = 0;
	    int vv2 = 0;
	    int vv3 = 0;
	    for (int j = 1; j <= ids6.size(); j++) {
		if (ids6[j] == v1)
		    vv1++;
		if (ids6[j] == v2)
		    vv2++;
		if (ids6[j] == v3)
		    vv3++;
	    }
	    if (vv1 == 0) {
		ids6.push_back(v1);
		idss6++;
	    }
	    if (vv2 == 0) {
		ids6.push_back(v2);
		idss6++;
	    }
	    if (vv3 == 0) {
		ids6.push_back(v3);
		idss6++;
	    }
	    f6++;
	}
    case (7):
	{
	    int vv1 = 0;
	    int vv2 = 0;
	    int vv3 = 0;
	    for (int j = 1; j <= ids7.size(); j++) {
		if (ids7[j] == v1)
		    vv1++;
		if (ids7[j] == v2)
		    vv2++;
		if (ids7[j] == v3)
		    vv3++;
	    }
	    if (vv1 == 0) {
		ids7.push_back(v1);
		idss7++;
	    }
	    if (vv2 == 0) {
		ids7.push_back(v2);
		idss7++;
	    }
	    if (vv3 == 0) {
		ids7.push_back(v3);
		idss7++;
	    }
	    f7++;
	}
    case (8):
	{
	    int vv1 = 0;
	    int vv2 = 0;
	    int vv3 = 0;
	    for (int j = 1; j <= ids8.size(); j++) {
		if (ids8[j] == v1)
		    vv1++;
		if (ids8[j] == v2)
		    vv2++;
		if (ids8[j] == v3)
		    vv3++;
	    }
	    if (vv1 == 0) {
		ids8.push_back(v1);
		idss8++;
	    }
	    if (vv2 == 0) {
		ids8.push_back(v2);
		idss8++;
	    }
	    if (vv3 == 0) {
		ids8.push_back(v3);
		idss8++;
	    }
	    f8++;
	}
 
    }
}

Puis, j'ai crée un fichier pour stocker ces points et la liste des identificateurs de chaque face.Bon, mon problème debute maintenant

. Lorsque j'essaye d'ajouter la liste des sommets de chaque face au fichier, au debut, tout va bien, mais en arrivant à la fin du boucle, je vois que j'ai une liste des sommets qui se repètent telque:

95 115 114
123 124 125
126 127 128
129 130 131
132 133 134
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132
135 133 132

Je sais pas ou est le problème

. ceci est le code:

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
string const faceun("face1.ply");
objetface.open(faceun.c_str(), ios::app);
objetface << "ply" << endl;
objetface << "format ascii 1.0" << endl;
objetface << "element vertex " << idss1 << endl;
objetface << "property float32 x" << endl;
objetface << "property float32 y" << endl;
objetface << "property float32 z" << endl;
objetface << "element face " << f1 << endl;
objetface << "property list uint8 int32 vertex_indices" << endl;
objetface << "end_header" << endl;
for (int j = 1; j <= idss1; j++) {
    double x[3];
    this->Input->GetPointCoordinates(ids1[j], x);
    objetface << x[0] << " " << x[1] << " " << x[2] << endl;
 
}
 
objetface.close();
for (vtkIdType d = 0; d < this->Input->GetNumberOfCells(); d++)	//pour afficher la liste des faces
{
    int dd = PointsNode->GetValue(d);
    vtkIdType k1;
    vtkIdType k2;
    vtkIdType k3;
    int vv1;
    int vv2;
    int vv3;
 
    if (dd == 1) {
	objetface.open(faceun.c_str(), ios::app);
	//recupérer les cellules d'une meme région
	this->Input->GetFaceVertices(d, k1, k2, k3);
	for (int g = 1; g < idss1; g++) {
	    if (ids1[g] == k1)
		vv1 = g;
	    else if (ids1[g] == k2)
		vv2 = g;
	    else if (ids1[g] == k3)
		vv3 = g;
	}
	cout << " ancien sommet du cellule : " << k1 << " " << k2 << " " << k3 << "n" << endl;
	objetface << vv1 << " " << vv2 << " " << vv3 << endl;
	cout << " nouvel sommet du cellule : " << vv1 << " " << vv2 << " " << vv3 << "n" << endl;
	objetface.close();
	cout << "idss1 = " << idss1;
    }
}

Merci de me corriger SVP.
Désolée car mon code est un peu long

**cob59** · 26/09/2012, 10h08

Il manque des break à la fin de chaque bloc case. D'ailleurs ces case m'ont l'air de copiés-collés assez honteux qu'il serait bon de re-factoriser

**tasnim2010** · 26/09/2012, 15h53

Envoyé par cob59

Il manque des break à la fin de chaque bloc case. D'ailleurs ces case m'ont l'air de copiés-collés assez honteux qu'il serait bon de re-factoriser

hh, ah oui, c'est la cause du break. ça fait longtemps avec le switch.. case

. ce copier coller, sert à remplir 8 vecteurs, selon la region de la face. j'ai essayer de diminuer le temps d'execution et la complexité de mon programme par ça, au lieu de faire 8 boucle for. Si vous avez une suggestion, ça sera le bienvenue

**mitkl** · 26/09/2012, 17h11

J'ai pas tout regardé en profondeur mais tu ne serais pas en train de coder toi-même l'enregistrement d'un fichier en PLY ?

Pourquoi ne pas utiliser directement les classes de vtk comme vtkPLYWriter ? Il suffit de former 3 PolyData à partir de tes 3 régions et d'injecter leurs sorties en entrée du Writer et le tour est joué.

**tasnim2010** · 26/09/2012, 18h16

Envoyé par mitkl

J'ai pas tout regardé en profondeur mais tu ne serais pas en train de coder toi-même l'enregistrement d'un fichier en PLY ?

Pourquoi ne pas utiliser directement les classes de vtk comme vtkPLYWriter ? Il suffit de former 3 PolyData à partir de tes 3 régions et d'injecter leurs sorties en entrée du Writer et le tour est joué.

Bien sure j'ai essayé d'utiliser vtkPLYWriter, mais j'ai trouvé un problème, après la génération de mon fichier, lorsque je l'ouvre avec mon application pour le visualiser, il m'affiche, juste un triangle

. c'est pour cela, que j'ai essayé de l'écrire à la main.
Je serais reconnaissante si vous pouvez m'aider à corriger mon code avec vtkPlyWriter:

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
for(int m=1; m<NumberOfSplitNodes; m++)// debut de boucle des régions
{ 
vtkIdType v1;
vtkIdType v2;
vtkIdType v3;
vtkIdList *pts=vtkIdList::New();
vtkPoints* points = vtkPoints::New();
vtkSmartPointer<vtkCellArray> vertices = vtkSmartPointer<vtkCellArray>::New();
vtkSmartPointer<vtkCellArray> triangles = vtkSmartPointer<vtkCellArray>::New();
vtkSmartPointer<vtkCellArray> lines = vtkSmartPointer<vtkCellArray>::New();
vtkSmartPointer<vtkTriangle> triangle = vtkSmartPointer<vtkTriangle>::New();
vtkSmartPointer<vtkLine> line0 = vtkSmartPointer<vtkLine>::New();
vtkSmartPointer<vtkLine> line1 = vtkSmartPointer<vtkLine>::New();
vtkSmartPointer<vtkLine> line2 = vtkSmartPointer<vtkLine>::New();
for ( int k=0; k<this->Input->GetNumberOfPoints();k++)
{
double x[3];
this->Input->GetPointCoordinates(k,x);
vtkIdType pid[1];
pid[0] = points->InsertNextPoint(x[0],x[1],x[2]);
vertices->InsertNextCell(1,pid);
}
for (int k=0; k<this->Input->GetNumberOfCells();k++)
{
vtkIdType pid[1];
this->Input->GetFaceVertices(k,v1,v2,v3);
line0->GetPointIds()->SetId ( 0,v1 );
line0->GetPointIds()->SetId ( 1,v2 );
line1->GetPointIds()->SetId ( 0,v1 );
line1->GetPointIds()->SetId ( 1,v3 );
line2->GetPointIds()->SetId ( 0,v2 );
line2->GetPointIds()->SetId ( 1,v3 );
lines->InsertNextCell(line0);
lines->InsertNextCell( line1 );
lines->InsertNextCell( line2 );					
triangle->GetPointIds()->SetId(v1,v1);
triangle->GetPointIds()->SetId(v2,v2);
triangle->GetPointIds()->SetId(v3,v3);
triangles->InsertNextCell(triangle);
}
 vtkSmartPointer<vtkPolyData> newdata = 
   vtkSmartPointer<vtkPolyData>::New(); 
 
 
newdata->SetPoints(points); 
newdata->SetVerts(vertices);
newdata->SetLines(lines);
newdata->SetPolys(triangles);
 
 
vtkPLYWriter *Writer=vtkPLYWriter::New();
 
Writer->SetInput(newdata);
std::stringstream strfile;
strfile<<"Region"<<m<<".ply";
Writer->SetFileName(strfile.str().c_str());
Writer->Write();
}

Merci d'avance

**Neckara** · 26/09/2012, 18h19

Envoyé par tasnim2010

hh, ah oui, c'est la cause du break. ça fait longtemps avec le switch.. case

. ce copier coller, sert à remplir 8 vecteurs, selon la region de la face. j'ai essayer de diminuer le temps d'execution et la complexité de mon programme par ça, au lieu de faire 8 boucle for. Si vous avez une suggestion, ça sera le bienvenue

Plutôt que d'utiliser 8 vecteurs avec un switch et 8 cases, ne peux-tu pas utiliser un tableau de 8 vecteurs ce qui te permettrais sûrement de diviser ton code par 8?

Exemple :
Remplacer :

Code :

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

1
2
3
4
5
6
7
8
9
switch (n)
{
           case 1 :
                    toto1 = 1;
           break;
           case 2 :
                    toto2 = 1;
           break;
}

Par :

Code :

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

toto[n] = 1;

Je pense que c'est ce que cob59 voulait dire.

**mitkl** · 26/09/2012, 18h44

Envoyé par tasnim2010

Bien sure j'ai essayé d'utiliser vtkPLYWriter, mais j'ai trouvé un problème, après la génération de mon fichier, lorsque je l'ouvre avec mon applicatio pour le visualiser, il m'affiche, juste un triangle

. c'est pour celà, que j'ai essayer de l'ecrire à la main.

même avec le correctif des switch case ?