OpenCasCade —— 分割一条曲线并返回分割点的坐标及参数

OpenCasCade中曲线分割问题

在OpenCasCade中，参数类曲线有两种分割方法（暂不清楚是否还有其他方法），以B样条曲线为例（相关代码附在文字末尾）

通过参数进行分割

分别获得曲线的 FirstParameter 和 LastParameter ，然后对参数进行分割，获得n个ui，并对每个ui调用D0（获得这个点的坐标值）或D1（获得这个点的坐标值和切向量）。这个方法的优点是，简单易行，好操作，缺点均分参数轴获得的曲线的分割并不是均匀的，图例如下：

直接对曲线分割

通过类 GCPnts_UniformAbscissa 实现对一个Geom_Curve对象的平均分割，获得分割点的坐标及对应的参数，将坐标和参数存入一个列表中，供其他操作使用，分割曲线的图例如下：

两种分割方法的效果放在一张图上展示如下

这样会更明显一点，黄色点位采用参数分割获得的分割点，白色点为直接对曲线进行分割获得的分割点：

相关代码

        ...
        Handle(Geom_BSplineCurve) curve = new Geom_BSplineCurve(Poles, knots, mults, degree);
        ...
        Standard_Integer splitN = 10;
        Standard_Real firstParam = curve->FirstParameter();
        Standard_Real lastParam = curve->LastParameter();
        for(Standard_Integer i=0; i<splitN+1; ++i)
        {
            Standard_Real ui = i * (lastParam - firstParam) / splitN;
            gp_Pnt pi;
            gp_Vec veci;
            curve->D1(ui, pi, veci);
            TopoDS_Vertex verti = BRepBuilderAPI_MakeVertex(pi);
            Handle(AIS_Shape) vertishape = new AIS_Shape(verti);
            myAISContext->SetColor(vertishape,Quantity_NOC_YELLOW,Standard_False); 
            myAISContext->Display(vertishape,Standard_False);
        }

直接分割曲线对应的代码段：

        ...
        Handle(Geom_BSplineCurve) curve = new Geom_BSplineCurve(Poles, knots, mults, degree);
        ...
        Standard_Integer splitN = 10;
        GeomAdaptor_Curve GAC(curve);
        GCPnts_UniformAbscissa UA(GAC, splitN+1);
        gp_Pnt* pnts = new gp_Pnt[splitN+1];
        Standard_Real* params = new Standard_Real[splitN+1];
        if(UA.IsDone())
        {
            Standard_Real n = UA.NbPoints();
            Standard_Integer index = 0;
            for(; index<n+1; ++index)
            {
                Standard_Real parami = UA.Parameter(index+1);
                params[index] = parami;
                gp_Pnt tpnt;
                curve->D0(parami, tpnt);
                pnts[index] = tpnt;
            }
        }