-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathfiles.cpp
198 lines (177 loc) · 6.89 KB
/
files.cpp
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
#include "files.h"
#include "cork.h"
#include <iostream>
#include "vec.h"
#include <vcg/complex/trimesh/create/platonic.h>
#include <vector>
using namespace std;
#ifndef uint
typedef unsigned int uint;
#endif
// new version - 2018.11.28 修改 GY (合并原来的 readMeshModel 和 file2corktrimesh)
void MeshModel_to_CorkTriMesh(const MeshModel& mm, CorkTriMesh *out)
{
out->n_vertices = mm.cm.vert.size();
out->n_triangles = mm.cm.face.size();
out->triangles = new uint[out->n_triangles * 3];
out->vertices = new float[out->n_vertices * 3];
for(CMeshO::ConstFaceIterator fi = mm.cm.face.begin(); fi != mm.cm.face.end(); ++fi)
{
int i = fi - mm.cm.face.begin();
out->triangles[3 * i] = (int)(fi->V(0) - &*mm.cm.vert.begin());
out->triangles[3 * i + 1] = (int)(fi->V(1) - &*mm.cm.vert.begin());
out->triangles[3 * i + 2] = (int)(fi->V(2) - &*mm.cm.vert.begin());
}
for(CMeshO::ConstVertexIterator vi = mm.cm.vert.begin(); vi != mm.cm.vert.end(); ++vi)
{
int i = vi - mm.cm.vert.begin();
out->vertices[3 * i] = vi->P().X();
out->vertices[3 * i + 1] = vi->P().Y();
out->vertices[3 * i + 2] = vi->P().Z();
}
}
// new version - 2018.11.28 GY 比上边的不需要变换矩阵的就多了一个变换矩阵和坐标相乘
void MeshModel_to_CorkTriMesh_withMoveMat(const MeshModel& mm, CorkTriMesh *out, const vcg::Matrix44f& move_matrix)
{
out->n_vertices = mm.cm.vert.size();
out->n_triangles = mm.cm.face.size();
out->triangles = new uint[out->n_triangles * 3];
out->vertices = new float[out->n_vertices * 3];
for(CMeshO::ConstFaceIterator fi = mm.cm.face.begin(); fi != mm.cm.face.end(); ++fi)
{
int i = fi - mm.cm.face.begin();
out->triangles[3 * i] = (int)(fi->V(0) - &*mm.cm.vert.begin());
out->triangles[3 * i + 1] = (int)(fi->V(1) - &*mm.cm.vert.begin());
out->triangles[3 * i + 2] = (int)(fi->V(2) - &*mm.cm.vert.begin());
}
for(CMeshO::ConstVertexIterator vi = mm.cm.vert.begin(); vi != mm.cm.vert.end(); ++vi)
{
int i = vi - mm.cm.vert.begin();
Point3f tmpP = move_matrix * vi->P();
out->vertices[3 * i] = tmpP.X();
out->vertices[3 * i + 1] = tmpP.Y();
out->vertices[3 * i + 2] = tmpP.Z();
}
}
// GY 融合两个模型，目前用来融合（所有牙齿）和（附件与倒凹），所以不需要变换矩阵
void CSG_Union(const MeshModel& m1, const MeshModel& m2, MeshDocument* md, const int& which_jaw)
{
CorkTriMesh *ctm1 = new CorkTriMesh(), *ctm2 = new CorkTriMesh(), *ctm_res = new CorkTriMesh();
MeshModel_to_CorkTriMesh(m1, ctm1); // 直接转了，不需要转成 FileMesh 过渡
MeshModel_to_CorkTriMesh(m2, ctm2);
if(isSolid(*ctm1) == false) printf("第一个模型不是 Solid 的!\n");
if(isSolid(*ctm2) == false) printf("第二个模型不是 Solid 的!\n");
computeUnion(*ctm1, *ctm2, ctm_res);
delete ctm1;
delete ctm2;
vector<Point3f> vertices;
vector<Point3i> faces;
for(int i = 0; i < ctm_res->n_vertices; i++)
{
Point3f p;
p.X() = ctm_res->vertices[3 * i];
p.Y() = ctm_res->vertices[3 * i + 1];
p.Z() = ctm_res->vertices[3 * i + 2];
vertices.push_back(p);
}
for(int i = 0; i < ctm_res->n_triangles; i++)
{
Point3i f; // 一个三角形的三个顶点的索引
f.X() = ctm_res->triangles[3 * i];
f.Y() = ctm_res->triangles[3 * i + 1];
f.Z() = ctm_res->triangles[3 * i + 2];
faces.push_back(f);
}
delete ctm_res;
QString qlabel;
if(which_jaw == 0) qlabel = "U_aligner_tmp.obj";
else qlabel = "L_aligner_tmp.obj";
MeshModel* new_mm;
new_mm = md->addNewMesh("", qlabel, false);
vcg::tri::ConcaveCover<CMeshO>(new_mm->cm, vertices, faces);
new_mm->updateDataMask(MeshModel::MM_ALL);
tri::UpdateBounding<CMeshO>::Box(new_mm->cm);
tri::UpdateNormals<CMeshO>::PerVertexNormalizedPerFaceNormalized(new_mm->cm);
// 将CSG结果分解多个连通区域，选择最大的
CMeshO &cm = new_mm->cm;
vector<pair<int, CMeshO::FacePointer> > connectedCompVec;
int numCC = tri::Clean<CMeshO>::ConnectedComponents(cm, connectedCompVec);
printf("Found %d Connected Components\n", numCC);
if(numCC == 1)
{
if(which_jaw == 0)
{
new_mm->setLabel(QString("U_aligner.obj"));
printf("完全连通，U_aligner.obj 有 %d 个点，%d 个面\n", new_mm->cm.vert.size(), new_mm->cm.face.size());
}
else
{
new_mm->setLabel(QString("L_aligner.obj"));
printf("完全连通，L_aligner.obj 有 %d 个点，%d 个面\n", new_mm->cm.vert.size(), new_mm->cm.face.size());
}
return;
}
uint flag_for_seperate = 0;
for(size_t i = 0; i < numCC; ++i)
{
tri::UpdateSelection<CMeshO>::ClearFace(cm);
connectedCompVec[i].second->SetS();
tri::UpdateSelection<CMeshO>::FaceConnectedFF(cm);
tri::UpdateSelection<CMeshO>::ClearVertex(cm);
tri::UpdateSelection<CMeshO>::VertexFromFaceLoose(cm);
// MeshModel *destMesh= md->addNewMesh("",QString("CC %1").arg(i));
MeshModel *destMesh = new MeshModel(md, "", QString("CC %1").arg(i));
// 需要头文件 #include <vcg/complex/trimesh/append.h>
tri::Append<CMeshO, CMeshO>::Mesh(destMesh->cm, cm, true);
if(destMesh->cm.vert.size() < 10000)
{
delete destMesh;
continue;
}
if(which_jaw == 0) destMesh->setLabel(QString("U_aligner.obj"));
else destMesh->setLabel(QString("L_aligner.obj"));
++flag_for_seperate;
md->meshList.push_back(destMesh);
printf("\t删除不连通区域后 %s 剩余 %d 个点，%d 个面\n", (const char*)destMesh->label().toLocal8Bit(), destMesh->cm.vert.size(), destMesh->cm.face.size());
destMesh->updateDataMask(destMesh);
tri::UpdateBounding<CMeshO>::Box(destMesh->cm);
tri::UpdateNormals<CMeshO>::PerVertexNormalizedPerFaceNormalized(destMesh->cm);
}
if(flag_for_seperate != 1)
printf("\n\nError seperating!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n\n");
// 如果选择了最大的连通区域，原来的meshmodel就不要了
md->delMesh(new_mm);
printf("\n");
return;
}
// 融合 单颗牙齿m1 和 对应的附件m2，最后保存在m1中(目前认为一颗牙齿最多只有一个附件)
// 因为牙齿和附件实际坐标是没有变的，只是变换矩阵变，所以需要传入变换矩阵
void CSG_Union_Teeth_and_Att_matrix(const MeshModel& m1, const MeshModel& m2, vcg::Matrix44f move_matrix, MeshModel* res, const int& which_jaw)
{
CorkTriMesh *ctm1 = new CorkTriMesh(), *ctm2 = new CorkTriMesh(), *ctm_res = new CorkTriMesh(); // 不写这个new，生成的是临时变量？
MeshModel_to_CorkTriMesh_withMoveMat(m1, ctm1, move_matrix);
MeshModel_to_CorkTriMesh_withMoveMat(m2, ctm2, move_matrix);
computeUnion(*ctm1, *ctm2, ctm_res);
delete ctm1;
delete ctm2;
vector<Point3f> vertices;
vector<Point3i> faces;
for(int i = 0; i < ctm_res->n_vertices; ++i)
{
Point3f p;
p.X() = ctm_res->vertices[3 * i];
p.Y() = ctm_res->vertices[3 * i + 1];
p.Z() = ctm_res->vertices[3 * i + 2];
vertices.push_back(p);
}
for(int i = 0; i < ctm_res->n_triangles; ++i)
{
Point3i aa;
aa.X() = ctm_res->triangles[3 * i];
aa.Y() = ctm_res->triangles[3 * i + 1];
aa.Z() = ctm_res->triangles[3 * i + 2];
faces.push_back(aa);
}
delete ctm_res;
vcg::tri::ConcaveCover<CMeshO>(res->cm, vertices, faces);
}