我使用cgal编写了一个程序如下:
typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef CGAL::Triangulation_vertex_base_2<K> Vb;
typedef CGAL::Triangulation_face_base_with_info_2<FaceInfo2,K> Fbb;
typedef CGAL::Constrained_triangulation_face_base_2<K,Fbb> Fb;
typedef CGAL::Triangulation_data_structure_2<Vb,Fb> TDS;
typedef CGAL::Exact_predicates_tag Itag;
typedef CGAL::Constrained_Delaunay_triangulation_2<K, TDS, Itag> CDT;
typedef CDT::Point Point;
typedef CGAL::Polygon_2<K> Polygon_2;
int main( )
{
//construct two non-intersecting nested polygons
Polygon_2 polygon1;
polygon1.push_back(Point(0,0));
polygon1.push_back(Point(2,0));
polygon1.push_back(Point(2,2));
polygon1.push_back(Point(0,2));
CDT cdt;
insert_polygon(cdt,polygon1);
CDT::Finite_faces_iterator t=cdt.faces_begin();
for (t = cdt.finite_faces_begin(); t != cdt.finite_faces_end(); t++)
{
//how can i achieve that?
}
}
我希望将三角形的顶点从cdt.finite_faces_begin()
变为三角形,t
指向for循环的每次迭代中的三角形。
例如在循环的第一次迭代中,我得到第一个三角形的顶点,在循环的第二次迭代中,我得到第一个+第二个三角形的顶点,在第三个三角形中,我得到顶点第一个+第二个+第三个三角形等等。 我怎么能实现这个目标?
答案 0 :(得分:1)
如果在第i个iteratotion中你只想使用脸部的三个顶点for(inti=0;i<3;++i){t->vertex(i);}
或t->vertex(i)->point()
。
如果您不想获得前一个脸部已经获得的顶点,则必须使用在循环外声明的std::set<Vertex_handle>
。