OpenMesh面分割

Question

我最近开始在一个项目中使用OpenMesh，在这个项目中，我必须优化网格。我需要使用face split(FaceHandle _fh, Point _p)操作在三角形的质心处插入一个顶点。但是，当我使用这种方法并尝试使用VertexFaceIterator获取新创建的顶点的面时，我总是得到带有顶点索引的无效面，例如(87，87，-1)或(12，12，-1)。似乎OpenMesh在分割后没有更新网格拓扑。

我的代码看起来像这样。faceStartIt为我提供了这些奇怪的索引。

    typedef OpenMesh::TriMesh_ArrayKernelT<> TriMesh;

    TriMesh::FaceIter triangleIt = mesh.faces_begin();

    for( ; triangleIt != mesh.faces_end(); )
    {            
        TriMesh::Point centroid = mesh.calc_face_centroid( *triangleIt );

        if( hasToSplit( centroid ) )
        {               
            TriMesh::VertexHandle centroidHandle = mesh.split( *triangleIt, centroid ); 
            TriMesh::VertexFaceIter faceStartIt = mesh.vf_begin( centroidHandle );
            TriMesh::VertexFaceIter faceEndIt = mesh.vf_end( centroidHandle );
            TriMesh::VertexFaceIter faceIt = faceStartIt; //faceIt++;

            for( ; faceIt != faceEndIt; ++faceIt )
            {                   
                // Do something for each face
            } 
        }
        else
        {
            ++triangleIt;
        }
    }

Answer 1

我找出了我遇到的问题。split方法运行得很好，我做错的地方在// Do something for each face注释中。对于每个面，我都在做边翻转操作，但引用正在丢失，因为翻转会导致后面要评估的面发生变化。解决方案是在一次迭代中，在分割面之后，在std::set中插入我想要翻转的每条边，以确保它们是唯一的。然后在另一次迭代中，在标记的边缘上进行实际翻转。

Answer 2

也许你可以

TriMesh::FaceIter triangleIt = mesh.faces_begin();

使用

TriMesh::FaceIter triangleIt = mesh.faces_sbegin();

garbage_collection()也可以，但是它的coefficient.If很低你使用face_sbegin()，迭代器将跳过无效的faces。