CGAL-在3D中进行Delaunay三角剖分后检索到错误的顶点索引

Question

我在VS2012中使用最新的CGAL（4.5）。 在Delaunay三角剖分之后，我将整个四面体网格输出到带有CGAL Delaunay三角剖分的默认流输出的“ testFile1”。 “ testFile1”中的数据正确。（我看到了网格，可以）

然后，我遍历所有四面体并获得顶点的索引，然后将它们输出到“ testFile2”，但是索引与“ testFile1”中的索引完全不同。 我可视化了网格，这完全是一团糟。

我还在自己的迭代中输出了每个四面体的顶点位置，它们与testFile1相同（第一个四面体的第一个顶点的坐标相同，依此类推）。 因此，testFile1和testFile2中的四面体相同，问题在于获取顶点索引。

我的代码：

std::vector<std::pair<Point,int>> V;
V.reserve(pointCount);
//push selected point into V with index info
for(int i=0;i<pointCount;i++)
{
    int id = randomNumbers[i];
    V.push_back(std::make_pair(Point(
        points[id]._p[0], 
        points[id]._p[1], 
        points[id]._p[2]), i));
}

Delaunay T;
T.insert(V.begin(), V.end());

//output vertices position and vertices index and cell neighbor index to testFile1
//data in testFile1 is right
std::ofstream oFileT("testFile1",std::ios::out);
oFileT << T;        



//output indices to testFile2 by my own iteration
std::ofstream oFileT("testFile2",std::ios::out);
int index = 0;
Delaunay::Finite_cells_iterator it;
for(it = T.finite_cells_begin(); it!=T.finite_cells_end(); ++it)
{
    for(int i=0;i<4;i++)
    {
        //here the coord is right 
        float testCoord = T.tetrahedron(it).vertex(i).x();

        //but the index here is totally different with testFile1
        int info = 0;
        if(!T.is_infinite(it->vertex(i)))
            info = it->vertex(i)->info();

        oFileT<<info<<" ";  //output vertices index to testFile2
    }
    oFileT<<std::endl;
    index++;
}

有任何想法吗？

更新：

奇怪的是，例如：

在testFile1中：tet索引：702 25 35 153

在我自己的迭代中，它-> vertex（i）-> point（）。x（）得到702 \\ 25 \\ 35 \\ 153的x坐标，但是它-> vertex（i）-> info（）得到1601 \\ 1352。 ....

Answer 1

沿希尔伯特曲线对顶点进行排序，以加快三角剖分的构造。 如果要使迭代顺序与info（）字段匹配，只需对顶点进行一次迭代并设置info（）