Question

我正在尝试对点集进行delaunay三角剖分，找到最接近输入点的点，然后得到它的事件顶点，但不知何故下面的代码不起作用。

#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Delaunay_triangulation_2.h>
#include <fstream>
typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef CGAL::Delaunay_triangulation_2<K> Triangulation;
typedef Triangulation::Edge_iterator Edge_iterator;
typedef Triangulation::Point Point;
typedef Triangulation::Vertex_handle Vertex_handle;
typedef Triangulation::Vertex_circulator Vertex_circulator;
int main( )
{
  std::ifstream in("data.txt");
  assert(in);
  std::istream_iterator<Point> begin(in);
  std::istream_iterator<Point> end;
  Triangulation T;
  T.insert(begin, end);
  std::cout << T.number_of_vertices() <<std::endl;
  Vertex_handle handle = T.nearest_vertex(Point(1, 1));
  std::cout << handle->point() << std::endl;
  std::cout<<"incidents: \n" << std::endl;
  Vertex_circulator circulator = T.incident_vertices(handle), done(circulator);
  do
    {
      std::cout << circulator->point() << std::endl;
    } while(++circulator != done);
  return 0;
}

例如，如果data.txt是

2 3
1 1
1 0

输出

3
1 1
incidents:

1 0
2 3
2.02461e-307 6.94896e-308

为什么我有最后一行？

Answer 1

CGAL的2D三角剖分具有连接到凸包的所有顶点的无限顶点（有关详细信息，请参阅user manual）。

您可以使用is_infinite函数来测试单形是否是无限的。

如何与cgal循环器一起工作？

1 个答案: