Question

我想使用CGAL执行以下操作：相交表面网格;切网;与另一个表面网格;切割网;将切割的网格划分为2个不同的网格。

为此，我使用了CGAL Mesh Slicer，它为我提供了一组折线。不幸的是，切片器不提供每条折线所属的切割网格的哪个面的信息。此外，它有效地实现了切片机动作。如何检索此信息？

我需要这些信息来沿着折线执行切割网格的细分，然后在2个单独的网格中进行划分。

在我的上下文中，切割网格是这样的表面网格：

这是切割网格+切割网格：

在我的代码中，我从切割网格的每个面生成一个平面，并使用切割网格执行切片操作。

这是我的代码：

    //Import off files and instantiate meshes
    const char* filename1 = "test.off";
    const char* filename2 = "cutMesh2.off";
    std::ifstream input(filename1);
    Mesh cuttedMesh, cutMesh;
    if (!input || !(input >> cuttedMesh))
    {
        std::cerr << "First mesh is not a valid off file." << std::endl;
        return 0;
    }
    input.close();
    input.open(filename2);
    if (!input || !(input >> cutMesh))
    {
        std::cerr << "Second mesh is not a valid off file." << std::endl;
        return 0;
    }
    input.close();
    // AABB Slicer constructor from the cutted mesh
    AABB_tree tree(edges(cuttedMesh).first, edges(cuttedMesh).second, cuttedMesh);
    tree.accelerate_distance_queries();
    CGAL::Polygon_mesh_slicer<Mesh, K> slicer_aabb(cuttedMesh, tree);

    std::cout << cutMesh.num_vertices()<< std::endl;
    // For each face of the cut mesh
    BOOST_FOREACH(face_descriptor f, faces(cutMesh))
    {
        std::cout << "Face " << f << std::endl;
        Point points [3];
        int i = 0;
        //for each point of the current face
        BOOST_FOREACH(vertex_descriptor v, CGAL::vertices_around_face(cutMesh.halfedge(f), cutMesh))
        {
            points[i]= cutMesh.point(v);
            ++i;
        }
        Polylines polylines;
        // Perform the slice between the current face of the cut mesh and the cutted mesh
        slicer_aabb(K::Plane_3(points[0],points[1],points[2]), std::back_inserter(polylines));
        std::cout << "the slicer intersects " << polylines.size() << " polylines" << std::endl;

        //for each polyline computed by this face of the cutmesh
        BOOST_FOREACH(Polyline_type polyline,polylines)
        {
            std::cout << "Polyline : " << polyline.size() << " points"<< std::endl;

            BOOST_FOREACH(Point point, polyline)
            {
                std::cout << "Point : " << point << std::endl;
            }
        }
        std::cout << std::endl;
        polylines.clear();
    }

感谢您的帮助。

Answer 1

该帖子的答案已由 sloriot 提出：

尝试使用位于Polygon_mesh_processing / include / CGAL / Polygon_mesh_processing /internal/clip.h中的未记录的clip（）函数

CGAL切片机：如何进行网格切割？

1 个答案: