CGAL中Polyhedron_3中与半边结构有关的具体使用

CGAL中Polyhedron_3中与半边结构有关的具体使用

 转载自：http://blog.csdn.net/jingwenlai_scut/article/details/2128638

 

 

半边结构可以根据点，半边，面中的任一个来得到另外两个。在我的使用过程中，我觉得半边是联系的纽带。CGAL中有些不方便的是不能用面来直接得到绕面的所有顶点，或者不能根据点来得到一绕该点的一邻域的所有点，这个经常性的操作有时会带来一些不方便。

在使用的过程中，你可以通过使用Source insight 等查看源代码的软件来看其具体包含的函数及使用方法。

<!--[if !supportLists]-->１．<!--[endif]-->据半边来得到顶点。

Halfedge_handle he;

Vertex_handle vh = he->vertex();

2.根据顶点得到其三维坐标　

Point p = v->point();

可查看Point.h看其成员及成员函数，若需转换其成其它形式如CVector3d,可用

Cvector3d point = Cvector3d(v->point().x(),v->point().y(),v->point().z());

<!--[if !supportLists]-->3.     <!--[endif]-->据半边得到面。

Facet_handle f = he->facet();

<!--[if !supportLists]-->4.     <!--[endif]-->据顶点得到半边

Vertex_handle v;

Halfedge_handle he = v->halfedge();

<!--[if !supportLists]-->5.     <!--[endif]-->据顶点得到面

因为没有直接通过顶点得到面的函数，所以先用半边来查找

Halfedge_around_vertex_circulator he = v->vertex_begin();

Facet_handle f = he->facet();

根据上图也可以得到，此面即为与此顶点相关联的面.

<!--[if !supportLists]-->6.     <!--[endif]-->据面得到顶点。

同样地，没有通过面直接得到顶点的函数，则利用半边来实现。

Halfedge_around_facet_circulator he = f->facet_begin();

Vertex_handle v = he->vertex();

需要注意的是，这里并没有说明其遍历的顺序。这里有个不清楚的是，因为一个面相关联的Halfedge有好几个，所以每次遍历的结果有可能不相同。

<!--[if !supportLists]-->7.     <!--[endif]-->据面得到半边

Halfedge_around_facet_circulator he = f->facet_begin();

<!--[if !supportLists]-->8.     <!--[endif]-->据一顶点点得到一邻域的顶点

这也是通过半边来实现的。先得到绕这个点的半边。而这个半边根据上面的图是指向该点的，因此，我们得用该半边的另一半来得到其所指的另一个顶点。

Halfedge_around_vertex_circulator he = v->vertex_begin();

Haledge_handle he_begin = he;

Vertex_handle v;

do{

     v = he->opposite()->vertex();

}while(++he != he_begin);