半边数据结构与网格细分算法Loop subdivision（附代码）

网格细分的原理其实并不难理解，它的难点主要在于如何实现。在看过无数有原理无代码的博客后，终于决定写一写我的实现方法，并附上代码供大家参考。c++写的可能比较笨拙，望见谅。

1.半边数据结构

很好理解，就是把网格的每一条边分成两个半边，半边是有方向的同一条边的两个半边方向相反。并且一条边是属于两个面，则半边完全属于一个面。

综合上述，得到半边的数据结构为一条半边的起始顶点origin，这条半边指向的下一条半边next，这条半边对应的另外半边opposite，和这条半边所在的面IncFace。

typedef struct HalfEdge//半边结构

{

        int origin;

        struct HalfEdge* next;

        struct HalfEdge* opposite;

        int IncFace;

}HalfEdge;

 

2.网格细分算法Loop subdivision

我只写了三角形网格的细分，即实现了Loop subdivision算法。

首先，我们应该清楚这个算法是怎样细分的：

​

图中上面两种情况适用于内部点和内部边产生新点的更新，下面两种情况适用于边界点和边界产生新点的更新。

由于我只细分了闭合的三维图形，所以代码只实现了第一种情况，就不讨论第二种了。

细分过程中，

• 每一条边都会产生一个新的顶点，这个顶点由该边所在的两个面中所有顶点决定，即v = 3/8*(v1 + v3) + 1/8*(v0 + v2)
• 每一个原来的顶点都会更新自己的坐标，这个更新由所有与它相邻的原顶点决定，即 ​，其中 ​

所以，每次细分要做的就是处理每一个旧顶点和每一条边，生成新的顶点，再将新的顶点连接成半边和面。

 

算法思路：

1.读入图形文件，分别定义数据结构保存图形的点，线，面。定义半边数据结构，将所有边转换为半边。然后开始细分。

2.更新旧顶点：第一步，遍历所有半边，找到一条以该顶点为起始顶点的半边。第二步，通过这条半边找到它的next半边，next半边的origin为v1。再找next半边的next半边 的opposite半边，进入下一个面。第三步，重复第二步直到找到的origin点vi等于v1，则已经找到所有vi，可计算出旧顶点的更新后坐标。

3.在每条旧边上生成新顶点：遍历每一条半边，通过半边和半边的opposite边的next与origin找到两个面的所有顶点，计算出新的顶点。

4.已经得到了所有新的顶点，连接它们得到新的半边和面。画出新图，细分完成。