Open CASCADE学习|通过指定点的曲线

在OpenCASCADE中，如果你想通过一系列指定的点来创建一条曲线，你可以使用Geom2dAPI_Interpolate类来实现二维曲线的插值，或者使用GeomAPI_Interpolate类来实现三维曲线的插值。这些类允许你定义一条B样条曲线，这条曲线将精确地通过你提供的点集。

Geom2dAPI_Interpolate类用于通过一系列点插值生成B样条曲线。这个类可以处理C2连续性（如果没有切线要求）或者C1连续性（如果有点上的切线要求）。如果需要周期性，曲线将闭合，并且连续性将为C1。使用这个类可以定义通过点的B样条曲线，也可以通过提供每个点的参数值和切向量来进一步约束曲线。

#include <Geom2dAPI_Interpolate.hxx>
#include <BRepBuilderAPI_MakeEdge.hxx>
#include <GeomAPI.hxx>
#include <Geom2d_BSplineCurve.hxx>
#include <TopoDS_Edge.hxx>
#include <math.h>
#include"Viewer.h"


int main(int argc, char* argv[])
{
    const int N = 100;
    Handle(TColgp_HArray1OfPnt2d) aPoints = new TColgp_HArray1OfPnt2d(1, N);
    for (int i = 0; i < N; i++)
    {
        double x = 3 * cos(M_PI / 50 * i);
        double y = 3 * sin(M_PI / 50 * i);
        aPoints->SetValue(i+1, gp_Pnt2d(x, y));
    }

    Geom2dAPI_Interpolate aInterpolater(aPoints, Standard_False, Precision::Approximation());
    aInterpolater.Perform();
    Handle(Geom2d_BSplineCurve) aBSplineCurve;
    aBSplineCurve = aInterpolater.Curve();
    //std::cout << "ok";
    gp_Pln plane = gp_Pln(gp::Origin(), gp::DZ());
    TopoDS_Edge s = BRepBuilderAPI_MakeEdge(GeomAPI::To3d(aBSplineCurve, plane));


    Viewer vout(50, 50, 500, 500);
    vout << s;
    vout.StartMessageLoop();
    return 0;
}