问题:如何使用OpenMesh完成三角形单元的分割,同时保护原始的拓扑信息。
解决方法一:
采用OpenMesh中封装的接口,mesh.split()可直接完成面片的分割
解决方法二:
将原始的三角形删除,然后再添加分割后的三角形。
特别注意!!!
使用方法一,之前的三角形并没有被删除,只是在此基础上添加了新的三角形。若只有一个三角形,将其分割成两个
三角形,最终会得到三个三角形,并且拓扑关系维护有问题,这留作进一步探讨!!!建立大家使用第二种方法。
方法一主要代码
关于OpenMesh常用的数据定义,请看:OpenMesh番外篇:头文件
//方法一
//例子,将两个相邻的三角形分割,分割后得到四个三角形
MyMesh mesh;
//首先添加4个顶点
VertexHandle v[4];
v[0] = mesh.add_vertex(Point(0,0,0));
v[1] = mesh.add_vertex(Point(100,0,0));
v[2] = mesh.add_vertex(Point(100,100,0));
v[3] = mesh.add_vertex(Point(0,100,0));
//添加面片
FaceHandle face[2];
vector<VertexHandle> vechandles;
vechandles.push_back(v[0]);
vechandles.push_back(v[1]);
vechandles.push_back(v[2]);
face[0] = mesh.add_face(vechandles);
vechandles.clear();
vechandles.push_back(v[0]);
vechandles.push_back(v[2]);
vechandles.push_back(v[3]);
face[1] = mesh.add_face(vechandles);
vechandles.clear();
//写出测试
OpenMesh::IO::write_mesh(mesh,"before.off");
//取共用边的中点作为分割点
VertexHandle vertex_split = mesh.add_vertex(Point(50,50,0));
//对两个边进行分割,split分割函数默认将多边长分割成三角形的形式
mesh.split(face[0],vertex_split);
mesh.split(face[1],vertex_split);
//写出测试
OpenMesh::IO::write_mesh(mesh,"after.off");
return;
方法二主要代码:
将原始的三角形删除,然后添加分割后的三角形,分割后的三角形并不难求得
//方法二
MyMesh mesh;
mesh.request_face_status();
mesh.request_edge_status();
mesh.request_vertex_status();
//首先添加4个顶点
VertexHandle v[4];
v[0] = mesh.add_vertex(Point(0,0,0));
v[1] = mesh.add_vertex(Point(100,0,0));
v[2] = mesh.add_vertex(Point(100,100,0));
v[3] = mesh.add_vertex(Point(0,100,0));
//添加面片
FaceHandle face[2];
vector<VertexHandle> vechandles;
vechandles.push_back(v[0]);
vechandles.push_back(v[1]);
vechandles.push_back(v[2]);
face[0] = mesh.add_face(vechandles);
vechandles.clear();
vechandles.push_back(v[0]);
vechandles.push_back(v[2]);
vechandles.push_back(v[3]);
face[1] = mesh.add_face(vechandles);
vechandles.clear();
//写出测试
OpenMesh::IO::write_mesh(mesh,"before.off");
//取共用边的中点作为分割点
VertexHandle vertex_split = mesh.add_vertex(Point(50,50,0));
//删除第一个三角形,然后添加两个三角形
mesh.delete_face(face[0]);
vechandles.push_back(v[0]);
vechandles.push_back(v[1]);
vechandles.push_back(vertex_split);
mesh.add_face(vechandles);
vechandles.clear();
vechandles.push_back(v[2]);
vechandles.push_back(v[1]);
vechandles.push_back(vertex_split);
mesh.add_face(vechandles);
vechandles.clear();
mesh.delete_face(face[1]);
vechandles.push_back(v[0]);
vechandles.push_back(v[3]);
vechandles.push_back(vertex_split);
mesh.add_face(vechandles);
vechandles.clear();
vechandles.push_back(v[2]);
vechandles.push_back(v[3]);
vechandles.push_back(vertex_split);
mesh.add_face(vechandles);
vechandles.clear();
//写出测试
OpenMesh::IO::write_mesh(mesh,"after.off");
示例
before
after
译者说明:因科研需要,避免重复造轮子,近来借用OpenMesh,因其结构封装严谨,灵活可靠,且十分方便,故抽出时间对官网的说明做点翻译。英文水平不高,如翻译有误请多多包涵,也请您多多指正。
英文原文地址:http://www.openmesh.org/Documentation/OpenMesh-Doc-Latest/index.html