如何使用CGAL检测网格中的孤立面片（或面片组）？

要解决这个问题，核心是识别网格的连通分量——你的闭合实体是一个大的连通分量，孤立面片组是独立的小连通分量。以下是具体实现步骤：

1. 选择适配的CGAL网格类型 #

优先使用CGAL::Surface_mesh，它能很好配合CGAL的连通性分析工具，也兼容STL/OBJ文件的读取。如果习惯用Polyhedron_3也可以，后续函数调用略有不同。

2. 读取网格文件 #

你已经掌握了STL读取，这里给出完整的初始化代码（以Surface_mesh为例）：

#include <CGAL/Surface_mesh.h>
#include <CGAL/IO/read_STL.h>
#include <CGAL/boost/graph/connected_components.h>
#include <iostream>
#include <algorithm>

typedef CGAL::Simple_cartesian<double> Kernel;
typedef Kernel::Point_3 Point;
typedef CGAL::Surface_mesh<Point> Mesh;

int main() {
    Mesh mesh;
    if (!CGAL::IO::read_STL("input_mesh.stl", mesh)) {
        std::cerr << "Failed to read STL file." << std::endl;
        return EXIT_FAILURE;
    }
    // 连通分量分析代码放在这里
}

OBJ文件只需替换read_STL为CGAL::IO::read_obj即可。

3. 计算所有连通分量 #

使用CGAL::connected_components函数，给每个面分配所属的连通分量ID，然后统计每个分量的面片数量：

// 为每个面存储连通分量ID
std::vector<int> component_ids(num_faces(mesh));
auto face_index_map = get(CGAL::face_index, mesh);

// 计算连通分量总数
int total_components = CGAL::connected_components(
    mesh,
    boost::make_iterator_property_map(component_ids.begin(), face_index_map)
);

// 统计每个分量的面片数
std::vector<int> component_sizes(total_components, 0);
for (int id : component_ids) {
    component_sizes[id]++;
}

4. 识别孤立面片组 #

假设主闭合实体是面片数最多的分量，剩下的分量就是孤立面片组：

// 找到主分量的ID
int main_component_id = std::max_element(component_sizes.begin(), component_sizes.end()) - component_sizes.begin();

// 遍历所有分量，输出孤立组信息
for (int i = 0; i < total_components; ++i) {
    if (i != main_component_id && component_sizes[i] > 0) {
        std::cout << "Found isolated patch with " << component_sizes[i] << " faces." << std::endl;
        // 若需要操作这些面片，可以通过遍历面并匹配component_ids[i]来获取
    }
}

关键细节 #

• 连通性判断：默认基于面共享边的规则，完全符合你的需求——孤立面片和主实体无共享边，会被分到不同分量。
• 若使用Polyhedron_3，需改用CGAL::Polygon_mesh_processing::connected_components（包含头文件CGAL/Polygon_mesh_processing/connected_components.h），迭代器和属性映射的用法类似。
• 可选验证：如果需要确认主实体是闭合的，可调用CGAL::Polygon_mesh_processing::is_closed函数，传入主分量的ID即可。

内容的提问来源于stack exchange，提问作者egev