试图将自定义C ++矩阵传递给numpy数组

时间:2016-11-14 10:11:39

标签: python c++ arrays numpy wrapper

我正在尝试使用Python包装来使用自定义C ++的东西。我们使用的主要类型是2D灰度图像类型,其数据在1D缓冲区中分配。我尝试以这种方式包装它(遵循示例in an ubuntu forum):

PyObject* to_python_object(const custom2DImage& img) {
  int type_num = (int)NPY_UBYTE;
  long int dims[2] = {img.nr(), img.nc()};

  uchar** tmp_img = new uchar*[img.nr()];
  tmp_img[0] = new uchar[img.nr() * img.nc()];
  for (int i = 1; i < img.nr(); ++i)
    tmp_img[i] = tmp_img[0] + img.nc();
  for (int i = 0; i < img.nc(); ++i)
    memcpy(tmp_mat[i], &img(i, 0), img.nc() * sizeof(uchar));
  PyObject* py_img = PyArray_SimpleNewFromData(2, dims, type_num, tmp_img[0]);

  Py_INCREF(py_img);
  delete[] tmp_img[0];
  delete[] tmp_img;

  PyObject *repr = PyObject_Repr(py_img);
  const char *s = PyString_AsString(repr);
  cout << s << endl;
  Py_XDECREF(repr);

  return py_img;
}

我在C ++函数中的Python对象表示是可以的,但是一旦我尝试在我的python main中打印它,就会出现段错误(即使形状很好)。 Python代码如下:

img = cst.read_image(filename);
if img is None:
    print("Can not load image from", filename)
    sys.exit(-1)
print(img)

你知道为什么我会遇到这个问题吗?

1 个答案:

答案 0 :(得分:1)

好的,我认为PyArray_SimpleNewFromData功能正在复制数据,但似乎并非如此。我的错误是释放tmp_img 2D数组。