Open3D错误ComputeHalfEdge失败。复制半边

Question

有人知道如何处理三角形网格以允许转换为半边网格吗。我正在尝试使用HalfEdge中的get_边界来检测边，这不是TriangleMesh中包含的函数

最后一行给出了错误

mesh = mesh.remove_duplicated_triangles()
mesh = mesh.remove_degenerate_triangles()
mesh = mesh.remove_duplicated_vertices()
mesh = mesh.remove_non_manifold_edges()
mesh = mesh.remove_unreferenced_vertices()
half_edge_mesh = o3d.geometry.HalfEdgeTriangleMesh.create_from_triangle_mesh(mesh)

Answer 1

在Open3d中，如果要从三角形网格创建半边，请确保网格是流形且没有奇异顶点。如果顶点有多个输出半边，将发生运行时错误。在最新版本的Open3d中，此错误尚未修复。Open3d Link

    for (size_t triangle_index = 0;
     triangle_index < mesh_cpy->triangles_.size(); triangle_index++) {
    const Eigen::Vector3i &triangle = mesh_cpy->triangles_[triangle_index];
    size_t num_half_edges = het_mesh->half_edges_.size();

    size_t he_0_index = num_half_edges;
    size_t he_1_index = num_half_edges + 1;
    size_t he_2_index = num_half_edges + 2;
    HalfEdge he_0(Eigen::Vector2i(triangle(0), triangle(1)),
                  int(triangle_index), int(he_1_index), -1);
    HalfEdge he_1(Eigen::Vector2i(triangle(1), triangle(2)),
                  int(triangle_index), int(he_2_index), -1);
    HalfEdge he_2(Eigen::Vector2i(triangle(2), triangle(0)),
                  int(triangle_index), int(he_0_index), -1);

    if (vertex_indices_to_half_edge_index.find(he_0.vertex_indices_) !=
                vertex_indices_to_half_edge_index.end() ||
        vertex_indices_to_half_edge_index.find(he_1.vertex_indices_) !=
                vertex_indices_to_half_edge_index.end() ||
        vertex_indices_to_half_edge_index.find(he_2.vertex_indices_) !=
                vertex_indices_to_half_edge_index.end()) {
        utility::LogError(
                "ComputeHalfEdges failed. Duplicated half-edges.");//Throw a runtime error.
    }

    het_mesh->half_edges_.push_back(he_0);
    het_mesh->half_edges_.push_back(he_1);
    het_mesh->half_edges_.push_back(he_2);
    vertex_indices_to_half_edge_index[he_0.vertex_indices_] = he_0_index;
    vertex_indices_to_half_edge_index[he_1.vertex_indices_] = he_1_index;
    vertex_indices_to_half_edge_index[he_2.vertex_indices_] = he_2_index;
}

singular vertices 就像红色箭头的顶点一样。在这些顶点中，会发生错误。 如果要查找边界，请尝试Openmesh或Libigl。Openmesh和Libigl都有python版本