网格简化算法（二）

网格简化算法（二）边折叠

效果展示：

可见，折叠边(u, v)，三角形1和2消失

算法思路：

具体解释

1、对每个顶点定义一个4x4的对称误差矩阵。
定义顶点到共享自己的面的距离的平方和为顶点到这些面的误差。

由于初始情况下，顶点都是在那些面上的，所以初始误差均为0。Q矩阵就是点所有相邻面的二次误差矩阵的和。

2、收缩一条边（v1, v2），则需要计算新的顶点位置。收缩一条边的Q可以确定：Qbar = Q1+ Q2。但是位置的选择需要计算，一种选择是从v1，v2，(v1 + v2)/2这三个位置上找，分别计算他们的Δ(v)。最小的那个就是我们要寻找的办法。另一种选择就是通过数值计算Δ(vbar)最小来寻找位置，由于Δ的表达式是一个二次项形式，因此令一阶导数为0

等价于求解：

其中q_ij为矩阵Qbar中对应的元素。如果系数矩阵可逆，那么通过求解上述方程就可以得到新顶点vbar的位置，如果系数矩阵不可逆，就通过第一种简单策略来得到新顶点vbar的位置。

3、由此得到所有边的cost。维护一个最小堆，每次取出堆顶元素，删除边表中该边，将新的顶点加入顶点表，当然，如果新的顶点在顶点表中，则不需要加入了。同时更新一下所有顶点的Q矩阵。

整理一下伪代码

1、建立边表和顶点表、面表
2、使用点到面的距离公式，建立Q矩阵
3、对所有的边计算cost，建立（维护）最小堆
4、pop出堆顶，删除边，更新顶点表、面表、Q矩阵
5、重复3 - 4直到面数小于阈值

Reference

[1] 拉风小宇博主
[2] https://www.cnblogs.com/shushen/p/5311828.html
[3] 百度百科