三角网格下的Half-Edge数据结构实现方法

Half-Edge结构即半边结构是图形学中一种十分重要的数据结构，许多算法都是基于这种结构实现的，从数学角度看这种结构十分容易，很好理解，而实际编码的时候似乎不是很轻松，之前写过两次都感觉写的不是很好，而网上也没有太多讲解实际编码实现的文章，所以记录一下简易的实现的步骤。

本文基于三角网格实现简易的Half-Edge，并能够读取.obj的模型数据

 

一、结构特点

最质朴的顶点+索引方式存储网格数据不能实现高效地完成一些功能的需求，例如：访问某个点的所有相邻点、访问某条边的所有邻边、访问包含某个点的所有面，而这些需求往往是一些图形算法的基础，所以很有必要寻求一种更加优秀的数据结构来实现功能。所谓的半边结构就是将每一条边拆成两条有方向的半边，如图所示。如果对其结构不是很熟悉可以参看其他的文章。

可知对于一个半边网格，需要实现三种结构：顶点、边、面。

struct HE_Vertex {
	HE_Edge * v_edge; //以该顶点为首顶点的出边(任意一条)
	HE_Face * v_face; //包含该顶点的面(任意一张)
	float x, y, z; 
};

struct HE_Edge {
	HE_Edge * e_pair;		//对偶边
	HE_Edge * e_succ;		//后继边
	HE_Vertex * e_vert;	        //首顶点
	HE_Face * e_face;		//右侧面
};

struct HE_Face {
	std::vector<HE_Vertex*> f_verts;
	HE_Edge * f_edge;
};

这里的实现简化了一些内容，譬如顶点的纹理坐标、法线等，如果有需求可以自行添加。

二、构建过程

这里的构建完全基于.obj的模型数据。

void HE_Mesh::LoadFromObj( const char * file )
{
	ifstream fs;
	fs.open(file);
	while (!fs.eof())
	{
		char line[100];
		fs.getline(line, 100);
		if (line[0] == 'v' && line[1] == ' ')
		{
			char * str = strtok(line, " ");
			str = strtok(NULL," ");
			float x = atof(str);
			str = strtok(NULL , " ");
			float y = atof(str);