DCEL C++实现

已于 2022-12-25 22:43:47 修改
阅读量340 收藏 1
点赞数
分类专栏： 几何计算 文章标签： c++ 开发语言
于 2022-12-24 23:01:25 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/yhaida/article/details/128431601
版权

DECL

DECL结构是比较通用的来表示几何的数据结构
主要包括三种数据结构：
1.点 vertex
2.半边 edge
3.面 face

点

存储点的位置以及任意一条由该点为起点向量的半边

半边

其中几何边由两条半边组成
在这里插入图片描述
这也是一种表示方式
在这里插入图片描述

面

存储一个面的外边界和内边界，如果外边界或者内边界不存在那么直接等于nullptr
在这里插入图片描述
像下面这种情况内部空洞需要用一个vector来存储了
在这里插入图片描述

通用的表达方式

在这里插入图片描述

数据结构简单实现（不包括具体操作）

template<typename type = float, size_t dim = DIM3>
	struct VertexDCEL
	{
		Vector<float, dim> point;// 1.Coordinates
		EdgeDCEL<type, dim>* incident_edge = nullptr;// 2.以point任意为起点的半边

		VertexDCEL(Vector<type, dim>& _point) :point(_point) {}

		//......
	};

	template<typename type = float, size_t dim = DIM3>
	struct EdgeDCEL
	{
		//1.起始点
		VertexDCEL<type, dim>* origin = nullptr;

		//2.与该半边相关联的三条半边
		EdgeDCEL<type, dim>* twin = nullptr;
		EdgeDCEL<type, dim>* next = nullptr;
		EdgeDCEL<type, dim>* prev = nullptr;

		//3.该半边对应的面
		FaceDCEL<type, dim>* incident_face = nullptr;

		int id = -1;
		EdgeDCEL() :id(-1) {}
		EdgeDCEL(VertexDCEL<type, dim>* _origin) :origin(_origin)
		{
			id++;
		}

		VertexDCEL<type, dim>* destination()
		{
			return twin->origin;
		}

		//......
	};

	template<typename type = float, size_t dim = DIM3>
	struct FaceDCEL
	{
		EdgeDCEL<type, dim>* outer = nullptr;
		std::vector< EdgeDCEL<type, dim>*>inner;//可能会有空洞
		
		void print()//只输出外边界
		{
			if (outer)
			{
				auto edge_ptr = outer;
				auto next_ptr = outer->next;

				edge_ptr->origin->print();
				while (next_ptr != edge_ptr) {
					next_ptr->origin->print();
					next_ptr = next_ptr->next;
				}
			}
		}
		//......
	};

	template<typename type, size_t dim = DIM3>
	class PolygonDCEL
	{
		typedef Vector<type, dim> VetorNf;
		std::vector< VertexDCEL<type, dim>* > vertex_list;
		std::vector< EdgeDCEL<type, dim>* > edge_list;
		std::vector< FaceDCEL<type, dim>* > face_list;

		//保存空的边，没有组成面->边集合
		EdgeDCEL<type, dim>* empty_edge = new EdgeDCEL<type, dim>();
	public:
		explicit PolygonDCEL(std::vector< VetorNf>&);
		
		bool split(VertexDCEL<type, dim>* _v1, VertexDCEL<type, dim>* _v2);

		bool join(VertexDCEL<type, dim>* _v1, VertexDCEL<type, dim>* _v2);

		std::vector<VertexDCEL<type, dim>*> getVertexList();

		std::vector<FaceDCEL<type, dim>*> getFaceList();

		std::vector<EdgeDCEL<type, dim>*> getEdgeList();

		VertexDCEL<type, dim>* getVertex(VectorNf&);

		void getEdgesWithSamefaceAndGivenOrigins(VertexDCEL<type, dim>* _v1, VertexDCEL<type, dim>* _v2,
			EdgeDCEL<type, dim>** edge_leaving_v1, EdgeDCEL<type, dim>** edge_leaving_v2);
	};

构造简单多边形的函数

template<typename type, size_t dim>
	inline PolygonDCEL<type, dim>::PolygonDCEL(std::vector<VetorNf>& _points)
	{
		int size = _points.size();
		if (size < 3)
			return;
		//1.根据所给点，创建vertex
		for (size_t y = 0; i < _points.size(); ++i)
		{
			vertex_list.push_back(new VertexDCEL<type, dim>(_points[i]));
		}
		//2.根据创建vertex，创建半边
		for (size_t i = 0; i <= vertex_list.size() - 2; ++i)
		{
			auto hfedge = new EdgeDCEL<type, dim>(vertex_list[i]);
			vertex_list[i]->incident_dege = hfedge;

			auto edge_twin = new EdgeDCEL<type, dim>(vertex_list[i + 1]);

			hfedge->twin = edge_twin;
			edge_twin->twin = hfedge;

			edge_list.push_back(hfedge);
			edge_list.push_back(edge_twin);
		}
		//最后一条边
		auto hfedge = new EdgeDCEL<type, dim>(vertex_list.back());
		vertex_list[vertex_list.size() - 1]->incident_dege = hfedge;

		auto edge_twin = new EdgeDCEL<type, dim>(vertex_list.front());

		hfedge->twin = edge_twin;
		edge_twin->twin = hfedge;

		edge_list.push_back(hfedge);
		edge_list.push_back(edge_twin);

		//设置next，prev指针
		//edge_list里面是根据半边的顺序存储,所以注意是每隔2个处理一次
		for (size_t i = 2; i <= edge_list.size() - 3; ++i)
		{
			if (i % 2 == 0)//设置半边
			{
				edge_list[i]->prev = edge_list[i - 2];
				edge_list[i]->next = edge_list[i + 2];
			}
			else//设置twin半边
			{
				edge_list[i]->prev = edge_list[i + 2];
				edge_list[i]->next = edge_list[i - 2];
			}
		}
		edge_list[0]->prev = edge_list[edge_list.size() - 2];
		edge_list[0]->next = edge_list[2];
		edge_list[1]->prev = edge_list[3];
		edge_list[1]->next = edge_list[edge_list.size() - 1];

		edge_list[edge_list.size() - 2]->prev = edge_list[edge_list.size() - 4];
		edge_list[edge_list.size() - 2]->next = edge_list[0];
		edge_list[edge_list.size() - 1]->prev = edge_list[1];
		edge_list[edge_list.size() - 1]->next = edge_list[edge_list.size() - 3];
		
		//3.新建face
		FaceDCEL<type, dim>* f1 = new FaceDCEL<type, dim>();//简单多边形的内半边
		FaceDCEL<type, dim>* f2 = new FaceDCEL<type, dim>();//简单多边形的外半边

		f1->outer = edge_list[0];
		f2->inner.push_back(edge_list[1]);

		//设置整个内部face f1外边界的半边对应incident_face
		f1->outer->incident_face = f1;
		EdgeDCEL<type, dim>* edge = f1->outer->next;
		while (edge!= f1->outer)
		{
			edge->incident_face = f1;
			edge = edge->next;
		}
		 
		//设置整个外部face f2内边界的半边对应incident_face
		f2->inner[0]->incident_face = f2;
		edge = f2->inner[0]->next;
		while (edge != f2->inner[0])
		{
			edge->incident_face = f2;
			edge = edge->next;
		}
	}

重要操作一 找到一个vertex的所有出度边

在这里插入图片描述

template<typename type, size_t dim>
	inline void PolygonDCEL<type, dim>::getEdgesWithSamefaceAndGivenOrigins(VertexDCEL<type, dim>* _v1, VertexDCEL<type, dim>* _v2, EdgeDCEL<type, dim>** edge_leaving_v1, EdgeDCEL<type, dim>** edge_leaving_v2)
	{
		std::vector<EdgeDCEL<type, dim>*> edges_with_v1_ori, edges_with_v2_ori;
		auto v1_inci_edge = _v1->incident_edge;
		edges_with_v1_ori.push_back(v1_inci_edge);

		auto next_edge = v1_inci_edge->twin->next;
		while (next_edge != v1_inci_edge)
		{
			edges_with_v1_ori.push_back(next_edge);
			next_edge = next_edge->twin->next;
		}

		auto v2_inci_edge = _v2->incident_edge;
		edges_with_v2_ori.push_back(v2_inci_edge);

		auto next_edge = v2_inci_edge->twin->next;
		while (next_edge != v2_inci_edge)
		{
			edges_with_v2_ori.push_back(next_edge);
			next_edge = next_edge->twin->next;
		}

		//edges_with_v1_ori与edges_with_v2_ori里面存的半边对应的incident_face是不一样的
		for (auto& ev1 : edges_with_v1_ori)
		{
			for (auto& ev2 : edges_with_v2_ori)
			{
				if (ev1->incident_face->outer != nullptr) {//找内部的面
					if (ev1->incident_face == ev2->incident_face) {//找相同的内部面
						//这两条半边必定是相反的方向！！！！可以根据其遍历两个部分
						**edge_leaving_v1 = ev1;
						*edge_leaving_v2 = ev2;
						return;
					}
				}
			}
		}
	}

重要操作二 split

在这里插入图片描述

template<typename type, size_t dim>
	inline bool PolygonDCEL<type, dim>::split(VertexDCEL<type, dim>* _v1, VertexDCEL<type, dim>* _v2)
	{
		EdgeDCEL<type, dim>* edge_oriV1;
		EdgeDCEL<type, dim>* edge_oriV2;
		getEdgesWithSamefaceAndGivenOrigins(_v1, _v2, &edge_oriV1, &edge_oriV2);
		if (edge_oriV1->id == -1 || edge_oriV2->id == -1)
			return false;
		//判断是否是邻边
		if (edge_oriV1->next->origin == _v2 || edge_oriV1->prev->origin == _v2)
			return false;

		//创建新对边
		auto half_edge1 = new EdgeDCEL<type, dim>(_v1);
		auto half_edge2 = new EdgeDCEL<type, dim>(_v2);
		//设置新半边
		half_edge1->twin = half_edge2;
		half_edge2->twin = half_edge1;

		half_edge1->next = edge_oriV2;
		half_edge2->next = edge_oriV1;

		half_edge1->next->prev = half_edge1;
		half_edge2->next->prev = half_edge2;

		half_edge1->prev = edge_oriV1->prev;
		half_edge2->prev = edge_oriV2->prev;

		half_edge1->prev->next = half_edge1;
		half_edge2->prev->next = half_edge2;

		//新建两个面，并删除旧的面，设置面
		FaceDCEL<type, dim>* new_face1 = new FaceDCEL<type, dim>();
		new_face1->outer = half_edge1;
		half_edge1->incident_face = new_face1;
		auto temp_edge = half_edge1->next;
		while (temp_edge != half_edge1) {
			temp_edge->incident_face = new_face1;
			temp_edge = temp_edge->next;
		}

		FaceDCEL<type, dim>* new_face2 = new FaceDCEL<type, dim>();
		new_face2->outer = half_edge2;
		half_edge2->incident_face = new_face2;
		temp_edge = half_edge2->next;
		while (temp_edge != half_edge2) {
			temp_edge->incident_face = new_face2;
			temp_edge = temp_edge->next;
		}

		face_list.push_back(new_face1);
		face_list.push_back(new_face2);

		//删除原先的面
		FaceDCEL<type, dim>* previous_face = edge_oriV1->incident_face;
		auto itr = std::find(face_list.begin(), face_list.end(), previous_face);
		if (itr != face_list.end()) {
			face_list.erase(itr);
			delete previous_face;
		}

		return true;
	}

单调多边形

在这里插入图片描述
正例
在这里插入图片描述
反例
在这里插入图片描述

单调多边形拆分

在这里插入图片描述
在这里插入图片描述

vertex种类

1.开始顶点
2.结束顶点
3.普通顶点
4.分裂顶点
5.合并顶点
在这里插入图片描述
在这里插入图片描述

开始顶点与结束顶点是一对，面向多边形内部的角度是锐角
在这里插入图片描述
普通顶点，两个端点是上下一侧
在这里插入图片描述
分裂顶点和合并顶点是需要实际操作的，单调多边形拆分主要是拆这两个顶点
在这里插入图片描述
分裂顶点是与多边形上边内部点相连内对角线，进行拆分。
合并顶点是与多边形下 边内部点相连内对角线，进行拆分。

Plane sweep algorithm

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

伪代码

在这里插入图片描述
在这里插入图片描述

辅助顶点

在这里插入图片描述
如何去决定内对角线应该和谁相连呢？
在这里插入图片描述

对于不同vertex不同的操作

start vertex

在这里插入图片描述

end顶点

在这里插入图片描述
翻译：
如果边的辅助点的类别在这个顶点结束是一个归并顶点（helper(ei)->merge vertex)，那么我们可以在这个辅助点和当前顶点之间添加一条对角线
从T中删除顶点的端点

分裂顶点

在这里插入图片描述

合并顶点

在这里插入图片描述

普通顶点

在这里插入图片描述

地理信息系统中叠置分析数据结构：topology模型DCEL
XiaoCai
02-17 884
叠置分析(2):数据结构：topology模型/DCEL 转载地址： http://www.whudj.cn/?p=824 1. 简介 叠置分析算法的基础-数据结构，即拓扑(topology)数据模型，在不同的领域/软件中，对拓扑数据模型的实现有不一样的叫法，比如geos中称之为planargraph(平面图)或者geomgraph,计算几何(Computational Geometry)的书籍以及CGAL中将其称为DCEL(Double Connected Edge List)，ArcGIS早期的数据模型
__stdcall、Q_DECL_IMPORT、Q_DECL_EXPORT
weixin_30342827的博客
06-13 527
https://www.cnblogs.com/songfeixiang/p/3733661.html 在C语言中，假设我们有这样的一个函数：intfunction(inta,intb)调用时只要用result=function(1,2)这样的方式就可以使用这个函数。但是，当高级语言被编译成计算机可以识别的机器码时，有一个问题就凸现出来：在CPU中，计算机没有办法知道一个函数调用需要多...
DCEL数据结构
baidu_34931359的博客
04-11 707
虽然https://blog.csdn.net/baidu_34931359/article/details/129962825?已经实现了平面点云的三角化，但是其在效率上还存在较大的优化空间，在非法边的检索过程中会变已生成三角形进行遍历，若可提前直到非法边所对应的三角形，可大幅提高检索效率，将该过程的时间复杂度由n降至常数。
计算几何课堂：DCEL（双向边链表）
主要研究图形学相关领域
07-30 2054
计算几何课堂：DCEL（双向边链表）
DCEL(DoubleConnectEdgeList 双向链连接表)
ZZLCSDN2017的博客
11-27 909
DCEL(DoubleConnectEdgeList 双向链连接表) DCEL用于表示平面划分: 中心线、德劳内三角划分、三角划分、维诺图等大多数平面划分结构 它连接这三种元素的关系：顶点（点Point3d）、边（有序的两个点Tuple<Point3d, Point3d>）、面（List<Tuple<Point3d, Point3d>>（可以转化为Polyline）） 对于每一个面，内部的vector顺序都是逆时针的，最外面的外包框是顺时针的。 对于一个生成的DCEL可
利用随机增量法与DCEl数据结构完成Delaunay算法（C#）
04-22
利用随机增量法与DCEl数据结构完成Delaunay算法（C#），程序内附大量注释，能够帮助完成对数据结构及Delaunay算法的了解。
C++ 新特性 | C++ 11 | decltype 关键字
个人技术博客
08-23 1594
decltype 是 C++11 新增的一个用来推导表达式类型的关键字。和 auto 的功能一样，用来在编译时期进行自动类型推导。引入 decltype 是因为 auto 并不适用于所有的自动类型推导场景，在某些特殊情况下 auto 用起来很不方便，甚至压根无法使用。
C/C++中decltype关键字用法总结
最新发布
红军不怕远征难，万水千山只等闲
04-06 3578
decltype可以作用于变量、表达式及函数名。①作用于变量直接得到变量的类型；②作用于表达式，结果是左值的表达式得到类型的引用，结果是右值的表达式得到类型；③作用于函数名会得到函数类型，不会自动转换成指针。 总的来说，decltype在C++中是一个非常强大的工具，它允许开发者以更简洁和灵活的方式处理复杂的类型问题。
c++处理类型（auto、decltyppe）
dzy521的博客
01-28 320
一、类型别名 1.使用关键字typedef定义类型别名 2.c++11中可以使用别名声明来定义类型别名 using ll = long long; 3.错误理解：把使用类型别名的部分用原来的样子来理解。本例中pstring本来是指针类型，但是替换后数据类型变成了char，*成了声明符的一部分 typedef char *pstring; const pstring pp=0; //指向char的常量指针 const char *pp = 0; //指向const char的指针 二、au
c++11中的declval和decltype
fpcc的专栏
12-08 2707
对declval和decltype进行分析说明，并举例在实际中应用
C++11中std::declval的实现机制浅析
01-20
本文主要给大家介绍了关于C++11中std::declval实现机制的相关内容，分享出来供大家参考学习，下面来一起看看详细的介绍： 在vs2013中，declval定义如下 template <_Ty> typenamea dd_rvalue_reference<_Ty>::type declval() _noexcept; 其中，add_rvalue_reference为一个traits，定义为 template <_Ty> struct add_rvalue_reference { typedef _Ty&& type; } 可见，declval被定义为一个函数，并且只有申明，没有
delaunay三角生成算法VC++实现版
02-24
delaunay三角生成算法是一种在...总之，Delaunay三角生成算法的VC++实现版提供了一个独立的工具，便于理解和应用这种重要的几何算法，对于学习和研究图形算法、进行相关项目开发的人员来说，这是一个非常有价值的资源。
VC++独立实现Delaunay三角算法详解
本文将介绍如何使用VC++（即Visual C++，是微软公司的一个集成开发环境）来实现Delaunay三角生成算法，具体是通过MFC（Microsoft Foundation Classes）框架来创建一个可以随机生成三角网的应用程序，并包含了DCEL...
【计算几何软件工具与库】CGAL库功能：用于解决几何问题的强大C++库
[【计算几何软件工具与库】CGAL库功能：用于解决几何问题的强大C++库](https://user-images.githubusercontent.com/147084302/277174476-0caab2dc-6d40-41d9-908b-1299299e248b.jpg) # 1. 计算几何与CGAL简介 计算...
轻量级3D曲面参数化编辑器：深入代码与功能解析
1. vcproj\DCEL：该项目包含了DCEL数据结构的实现，包括双连通边列表的构建和操作。 2. vcproj\geodesic：该项目负责计算表面上的测地线，为曲面参数化提供数学基础。 3. vcproj\OBJModel：负责处理OBJ模型的基本...
一个对象中如何保存拓扑关系和引用关系
Ausma的博客
05-06 390
比如，一个实体对象可以保存一个由点、线、面组成的拓扑结构，这个拓扑结构可以通过在实体对象的 TDF_Label 中保存一个 TopoDS_Shape 对象来表示；另外，一个实体对象可能会引用一组点、线、面等，这些被引用的对象可以通过在实体对象的 TDF_Label 中保存一组 TDF_Reference 对象来表示。一般来说，拓扑关系指的是对象之间的拓扑结构，比如点、线、面等之间的包含、相交等关系；而引用关系则是指一个对象引用另一个对象的方式，比如一个实体可能会引用一组点、线、面等构成自己的形状。
DCEL：如何连接和添加边
主要研究图形学相关领域
07-16 903
DCEL：如何连接和添加边
c++11 std::declval 详解
linuxheik的专栏
12-29 1696
c++11 std::declval 详解 函数模板 std::declval (c++11 only) template 　　typename add_rvalue_reference::type declval() noexcept; 功能描述:  返回一个类型的右值引用，不管是否有没有默认构造函数或该类型不可以创建对象。（可以用于抽象基类）; 参数： 无
c++11之for( decl : coll)循环赋值
lee2813的博客
07-26 436
1.语法结构： for( decl : coll){ statement; } 2.示例代码： for(int i : {1,2,3,4,5,6,7,8,9}{ std::cout<<i<<" "; } 其输出结果为： 1 2 3 4 5 6 7 8 9 3.作用 满足特殊情况下的输出
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

yhaida

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值