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

OpenCasCade中曲线分割问题
通过参数进行分割
直接对曲线分割
两种分割方法的效果放在一张图上展示如下
相关代码
OpenCasCade中曲线分割问题
在OpenCasCade中,参数类曲线有两种分割方法(暂不清楚是否还有其他方法),以B样条曲线为例(相关代码附在文字末尾)

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

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

 

 

直接对曲线分割

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

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

 

 

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

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

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

 

 

相关代码

        ...
        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;
            }
        }

 


————————————————
版权声明:本文为CSDN博主「求仁得仁AAA」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/DaGeYong/article/details/81533381

上一篇:Vue中:rules的写法:


下一篇:ggplot2分组箱线图添加均值