我如何在PythonOCC中使用样条曲线？

Question

我有一个关于如何在pythonOCC中使用样条的双重问题。

首先，我知道可以使用以下方式创建样条：

array = []
array.append(gp_Pnt2d (0,0))
array.append(gp_Pnt2d (1,2))
array.append(gp_Pnt2d (2,3))
array.append(gp_Pnt2d (4,3))
array.append(gp_Pnt2d (5,5))

pt2d_list = point2d_list_to_TColgp_Array1OfPnt2d(array)
SPL1      = Geom2dAPI_PointsToBSpline(pt2d_list).Curve()
display.DisplayShape(make_edge2d(SPL1) , update=True)

我希望您能通过以下方法来计算bspline：

BSPL1      = Geom2dAPI_PointsToBSpline(pt2d_list)

但是我该如何获得：

1. bspline的导数？
2. bspline的结点？
3. 结点是pt2d_list吗？
4. bspline的控制点？
5. 样条函数的系数？

我怎样才能删除或添加bspline的结点？

其次，当在PythonOCC中加载一个CAD图纸 .stp 文件时，可以像这样实现：

from OCC import TopoDS, StlAPI
shape = TopoDS.TopoDS_Shape()
stl_reader = StlAPI.StlAPI_Reader()
stl_reader.Read(shape,str(filename))
display.DisplayShape(shape)

如何从像结、B样条和系数这样的形状中获取数据？

- Steenstrup

---

你在其他地方找到答案了吗？能否请你和我们分享一下呢 :D？ - Jiloc

无望的老旧................. - Dima Tisnek

Answer 1

我已经使用了python-boost来实现这个。请查看这个函数：http://letslearncomputing.blogspot.com/2013/04/c-program-for-cubic-spline-interpolation.html。你可以在这里得到一些你想要的5个值。

你只需要修改代码为C++（而不是C）。

BOOST_PYTHON_MODULE(Spline)
{
    import_array();
    boost::python::numeric::array::set_module_and_type("numpy", "ndarray");
    class_<Spline, Spline*>("Spline", init<>())
        .def("spline", &Spline::spline)
        ;
}

所以在Python中，您可以使用：

from Spline.Spline import *
operation = Spline()
value, error_ = operation.spline(np.array(your_x_array), np.array(your_y_array), 0.01)

Cpp类：

#define NUMBER_OF_SAMPLES 14
class Spline
{
public:
    boost::python::list spline(numeric::array& x_val, numeric::array& y_val, double look_up_val);
};

在boost::python::list Spline::spline(numeric::array& x_val, numeric::array& y_val, double p)函数中，您会获得：

PyArrayObject* x_pyArr = (PyArrayObject*)PyArray_FROM_O(x_val.ptr());

PyArrayObject* y_pyArr = (PyArrayObject*)PyArray_FROM_O(y_val.ptr());
int size = *(x_pyArr->dimensions), i , j;    
double* data_x = (double*)x_pyArr->data;
double* data_y = (double*)y_pyArr->data;
double h[NUMBER_OF_SAMPLES], a, b, c, d, s[NUMBER_OF_SAMPLES] = { 0 }, F[NUMBER_OF_SAMPLES], f[NUMBER_OF_SAMPLES], x[NUMBER_OF_SAMPLES], m[NUMBER_OF_SAMPLES][NUMBER_OF_SAMPLES] = { 0 }, temp;

for (int i = 0; i < size; i++)
{
    x[i] = *(data_x + i);
    f[i] = *(data_y + i);
}

根据链接中的代码，我在Spline::spline函数中返回一个Python列表：

boost::python::list return_val;
// ....
return_val.append(sum);
return_val.append(result);
return return_val;

- Jagoda Gorus

Answer 2

我建议您查看scipy文档，并在其中搜索您要应用的函数。

- ctrl-alt-delete

---

这应该是一条注释，而不是答案。 - Daniel