生成三角网算法java_C++ 基于凸包的Delaunay三角网生成算法

Delaunay三角网，写了用半天，调试BUG用了2天……醉了。

基本思路比较简单，但效率并不是很快。

1. 先生成一个凸包；

2. 只考虑凸包上的点，将凸包环切，生成一个三角网，暂时不考虑Delaunay三角网各种规则。将生成的三角形放进三角形集合 Triangles 中；

3.将其它非凸包上的点全都插入。每插入一个点 ptA，都要判断被插入的点存在于 Triangles 集合的哪个三角形(trianA)之内，并将 ptA 与此三角形的三个点进行连接，删除 trianA，并将新生成的三角形加入到集合 Triangles 中。初始三角网生成结束；

3.1. 若 ptA 处在三角形DEF的DE边上，那么只连接点F与ptA ；

4.遍历三角形集合 Triangles(曾考虑用邻接矩阵，但是使用矩阵的复杂度反而会更高)，每遍历到一个三角形DEF，都要遍历DEF的三条边DE,EF, FD，并分别寻找另外一个三角形，此三角形与DEF存在一个公共边；

5.使用LOP(Local Optimization Procedure： 局部优化处理)处理步骤4返回的两个三角形 与DEG。

5.1. 做三角形DEF的外接圆cirDEF，如果点G在圆cirDEF之内，则从Triangles 集合中删除三角形DEF与DEG，并加入两个新三角形：FGD与FGE。实际上，是将两个三角形组成的四边形的对角线对调了。

6.使用新生成的两个三角形FGD与FGE执行步骤5，递归。递归出口为，对三角形进行LOP处理时，没有出现5.1的情况。

Note:

1.暂时存在的问题：

a. 当点数过多时，生成过程中会出现“栈溢出”情况。测试时，2000个点之内是成功的。此情况是由于三角形过多导致的递归过深，需要重新组织算法结构才能避免。

b. 生成时间太长。可采取的办法1是完全采用不同的生成算法，比如分块等等。其次是由于寻找三角形的复杂度过高，达到了O(n)，利用空间换时间的办法能接近O(1)，只不过多一个成员来存储三角形的拓扑结构，这个比较简单但优化效果有限。

2. 一旦遇到了一些算法上的问题，那么搬出数学知识往往非常有效(尽管我的数学水平一望见底)。使用矩阵的方式能够让算法的实现更加清晰，但复杂度有可能会升高。

File: Delaunay.h

1 #pragma once

2

3 #include "ConvexHull.h"

4 #include "Triangle.h"

5 #include "Circle.cpp"

6

7 classTrianIndex final8 {9 public:10 TrianIndex()11 {12 _isVaild = false;13 }14

15 ~TrianIndex()16 {17 _isVaild = false;18 }19

20 TrianIndex(unsigned int iA, unsigned int iB, unsigned intiC)21 {22 Init(iA, iB, iC);23 }24

25 TrianIndex(std::arraypts) :26 TrianIndex(pts[0], pts[1], pts[2])27 {28 }29

30 TrianIndex(const TrianIndex&other)31 {32 this->_isVaild =other._isVaild;33 this->_ptIdx =other._ptIdx;34 }35

36 unsigned int& Get(inti)37 {38 if (i < 0 || i > 2 || !_isVaild)39 {40 ErrorThrow("Error Triangle Point Index[0, 2] Or Invaild Triangle:" +std::to_string(i));41 return _ptIdx[0];42 }43

44 return_ptIdx[i];45 }46

47 unsigned int& operator[](inti)48 {49 returnGet(i);50 }51

52 const unsigned int Get(int i) const

53 {54 TrianIndex* pThis = const_cast(this);55 return pThis->Get(i);56 }57

58 const unsigned int operator[](int i) const

59 {60 TrianIndex* pThis = const_cast(this);61 return (*pThis)[i];62 }63

64 TrianIndex& operator=(const TrianIndex&other)65 {66 this->_isVaild =other._isVaild;67 this->_ptIdx =other._ptIdx;68 return *this;69 }70

71 bool IsVaild() const

72 {73 return_isVaild;74 }75

76 void SetVaild(boolisVaild)77 {78 _isVaild =isVaild;79 }80

81 void Init(unsigned int iA, unsigned int iB, unsigned intiC)82 {83 _ptIdx[0] =iA;84 _ptIdx[1] =iB;85 _ptIdx[2] =iC;86 _isVaild = true;87 }88

89 private:90 bool_isVaild;91 std::array_ptIdx;92 };93

94 //CHECK IT: It Has Been Abandoned.95 //96 classTrianMatrix final97 {98 public:99 enum IsConnected : bool

100 {101 NotConnected = false,102 Connected = true

103 };104

105 TrianMatrix():106 TrianMatrix(0)107 {108 }109

110 TrianMatrix(size_t ptCount)111 {112 _sign.clear();113 _sign.resize(ptCount);114 for (unsigned int x = 0; x < ptCount; ++x)115 {116 _sign[x].resize(ptCount);117 for (unsigned int y = 0; y < ptCount; ++y)1