如何在网格中分割三角形同时维护三角形网格的拓扑关系

                                                    

问题：如何使用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 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 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