delauney 三角化

原创

已于 2022-05-31 21:35:18 修改 · 1.1k 阅读

标签

#算法 #三角网格

收录于

于 2022-05-11 21:58:51 首次发布

欢迎关注更多精彩
关注我，学习常用算法与数据结构，一题多解，降维打击。

给定平面内三角网格进行最优delauney三角化

空圆性

如果一个三角网格中任意一对三角形满足空圆性，则该三角网格为delauney三角化网格。

对其中一个三角形作外接圆，如果圆没有第4个点，则满足空圆。否则不满足。如果一对三角形不满足空圆性，则只要调整中间的分隔边即可。
如上图，右边是不满足空圆性的，调整中间的边，变成左边的图就满足了。

可以证明，如果所有局部满足空圆性，则全局也满足空圆性。

空圆性判断

按照定义可以对三角形作外接圆进行判断，此种方法比较麻烦。

还有一种简便方法。

$可以计算边e的2个对角\alpha_1, \alpha_2, 如果\alpha_1+\alpha_2>\pi, 则不满足空圆性。$

顶点最优化

上述只能保证是一个delauney三角网格，但是网格的分布还不是很均匀。
想要得到均匀的网格，需要对顶点位置进行优化。优化公式如下。

算法实现

代码库：https://github.com/LightningBilly/ACMAlgorithms/tree/master/图形学算法/三角网格算法/MeshWork/src/hw12/

#include "../../PolyMesh/IOManager.h"
#include <string>
#include <cmath>
#include <unistd.h>

using namespace std;

using namespace acamcad;
using namespace polymesh;

PolyMesh* mesh12;

double get_TriFace_Area(MPolyFace* f)
{
   
   
    MHalfedge* he = f->halfEdge();
    MPoint3 v0 = he->fromVertex(