如何找到网格边与TopoDS_Edge的对应关系？

Question

我对TopoDS_Face进行了三角剖分形成了3D网格，我现在可以知道网格中那些边是TopoDS_Face的边缘，但是我不知道网格中的边和TopoDS_Face的边的对应关系。
我说的TopoDS_Face的边指的是TopoDS_Edge。
所以我的问题就是如何把三角化以后的边与TopoDS_Edge对应起来？
Open CASCADE:7.5.0
C++:17
编译器：VS2022社区版
操作系统：Windows11

Answer 1

“Devil组”引证GPT后的撰写：

用BRep_Tool类的函数来获得与TopoDS_Edge相对应的Geom_Curve对象。然后，你可以使用Geom_Curve对象的函数来获得其曲线参数，并与网格边上的点进行比较，以确定哪些网格边与TopoDS_Edge相对应。

以下是一个示例代码片段，其中使用了BRep_Tool类的函数和Geom_Curve对象的函数来找到TopoDS_Face的边缘和网格中的边之间的对应关系：

#include <BRep_Tool.hxx>
#include <Geom_Curve.hxx>
#include <TopoDS.hxx>
#include <TopoDS_Edge.hxx>
#include <TopoDS_Face.hxx>

// Assume 'triangles' is a vector containing the triangles of the mesh
// Assume 'face' is the TopoDS_Face object corresponding to the mesh
// Assume 'tolerance' is a tolerance value for comparing points

for (const auto& triangle : triangles) {
    for (int i = 0; i < 3; i++) {
        const auto& p1 = triangle.vertices[i];
        const auto& p2 = triangle.vertices[(i + 1) % 3];
        
        // Check if the edge between p1 and p2 is on the edge of the face
        TopoDS_Edge edge;
        if (!BRep_Tool::FindEdge(face, gp_Pnt(p1.x, p1.y, p1.z), gp_Pnt(p2.x, p2.y, p2.z), edge)) {
            // The edge is not on the edge of the face
            continue;
        }
        
        // Get the curve of the edge
        Handle(Geom_Curve) curve = BRep_Tool::Curve(edge);
        
        // Check if the points on the mesh edge match the curve of the edge
        if (curve->Value(curve->FirstParameter()).Distance(gp_Pnt(p1.x, p1.y, p1.z)) > tolerance
            || curve->Value(curve->LastParameter()).Distance(gp_Pnt(p2.x, p2.y, p2.z)) > tolerance) {
            // The mesh edge does not match the curve of the edge
            continue;
        }
        
        // The mesh edge matches the curve of the edge
        // Do something with the edge and its corresponding mesh edge
    }
}

假设'triangles'是一个包含网格三角形的向量，'face'是一个TopoDS_Face对象，'tolerance'是用于比较点的公差值。对于每个网格三角形的每条边，该代码检查该边是否在TopoDS_Face的边缘上，如果是，则找到与TopoDS_Edge对应的网格边并检查其是否与TopoDS_Edge的曲线匹配。如果是，则可以将它们彼此对应起来。