Delaunay 三角测量 3d

如何解决Delaunay 三角测量 3d

我正在尝试编写一个程序来进行 delaunay 三角剖分。这是我目前编写的代码：

using namespace std;
#include <iostream>
#include <fstream>
#include <vector>
#include <stdlib.h>
#include <Cgal/Exact_predicates_inexact_constructions_kernel.h>
#include <Cgal/delaunay_triangulation_3.h>
#include <Cgal/delaunay_triangulation_cell_base_3.h>
#include <Cgal/point_generators_3.h>
#include <Cgal/Triangulation_vertex_base_3.h>
#include <Cgal/Random.h>
#include <Cgal/point_generators_3.h>
#include <Cgal/draw_triangulation_3.h>

int main()
{
#ifdef Cgal_LINKED_WITH_TBB
  typedef Cgal::Exact_predicates_inexact_constructions_kernel K;
  // delaunay T3
  typedef Cgal::Triangulation_data_structure_3<
            Cgal::Triangulation_vertex_base_3<K>,Cgal::delaunay_triangulation_cell_base_3<K>,Cgal::Parallel_tag>                               Tds;
  typedef Cgal::delaunay_triangulation_3<K,Tds>              Triangulation;
  typedef Triangulation::Point                                Point;
  const int NUM_INSERTED_POINTS = 5000;
  typedef Cgal::Creator_uniform_3<double,K::Point_3> Creator;

  ifstream leer;
  ofstream escribir;
  int i;
  double x,y,z;
  vector<Point> V;
  V.reserve(NUM_INSERTED_POINTS);

cout << "Hola";

  // Construct the locking data-structure,using the bounding-Box of the points
  Triangulation::Lock_data_structure locking_ds(
      Cgal::BBox_3(-1.,-1.,1.,1.),50);


  // abrir el fichero

  leer.open("datos.txt");
 
  // Comprobar que el fichero de enTrada se ha abierto correctamente

  if (!leer)
    cout << "Error abriendo el fichero" << endl;

  else
  {

    // leer los datos del fichero mientras no se haya acabado el fichero

    while (!leer.eof())
    {
        leer >> x;
        leer >> y;
        leer >> z;

        for (int i = 0; i != NUM_INSERTED_POINTS; ++i)
          V.push_back(Point(x,z));
    }

    // cerrar el fichero

    leer.close();
  }

    // Construct the triangulation in parallel
  Triangulation T(V.begin(),V.end(),&locking_ds);
  assert(T.is_valid());

  //Me queda pasar a un fichero de salida las posiciones de las esferas o los datos que me de la triangulacion. No sé cuál es la salida que me da la
  //triangulación

  escribir.open("triangulacion.txt");
  // writing file output (simpletriangulation.cpp);
  escribir << T;

  //Dibujar la triangulación
  Cgal::draw(T);

#endif //Cgal_LINKED_WITH_TBB



  return 0;
}

我的家伙是我不太确定这段代码的“输出”是什么。我的意思是，我在“triangulacion.txt”文件中获得的是球体的位置或它们的半径？如果没有，我如何获得这些？

谢谢。