如何从 C++ 访问嵌入式 python 中的 numpy 数组？

Question

在 C++ 中访问 2dim numpy 数组的好方法是什么？我已经查看了 numpy/c api 和其他一些帖子，但这并没有让我更进一步。情况如下：

我在一个名为 Testfile.py 的 python 文件中定义了以下 numpy 数组：

import numpy as np
A = np.array([[1, 2, 3], [1, 2, 3], [1, 2, 3]])

现在，我想在 c++ 中访问这个数组，以使用它进行进一步的计算。这是我到目前为止所做的。 注意：为简单起见，我省略了错误处理和引用计数代码 sn-ps。

#define NPY_NO_DEPRECATED_APINPY_1_7_API_VERSION
#define PY_ARRAY_UNIQUE_SYMBOL cool_ARRAY_API
#include <Python.h>
#include <arrayobject.h> // numpy!
#include <iostream>
using namespace std;
#include <stdlib.h>
#include <string>

int main(){

// Name of input-file
char pyfilename[] = "Testfile";

// initilaize python interpreter
Py_Initialize();
import_array(); 

// load input-file
PyObject *pyName = PyUnicode_FromString(pyfilename);
PyObject *pyModule = PyImport_Import(pyName);

// import my numpy array object
char pyarrayname[] = "A";
PyObject *obj = PyObject_GetAttrString(pyModule, pyarrayname);

//------------------------------------------------------------
// The Problem starts here..

// Array Dimensions
npy_intp Dims[] = { PyArray_NDIM(obj) }; // array dimension
Dims[0] = PyArray_DIM(obj, 0); // number of rows
Dims[1] = PyArray_DIM(obj, 1); // number of columns

// PyArray_SimpleNew allocates the memory needed for the array.
PyObject *ArgsArray = PyArray_SimpleNew(2, Dims, NPY_DOUBLE);

// The pointer to the array data is accessed using PyArray_DATA()
double *p = (double *)PyArray_DATA(ArgsArray);

for (int i = 0; i<Dims[0]; i++)
{
    for (int j = 0; j<Dims[1]; j++)
        {
             p[i * Dims[1] + j] = *((int *)PyArray_GETPTR2(obj, i, j));
        }
}
// -----------------------------------------------------------


Py_Finalize();

return 0;
}

我使用 python 3.6 和 MSVC 2015。

编辑：我添加了我使用的标题并稍微改变了问题的表述。

编辑：我添加了 Swift 和 Alan Stokes

提供的建议解决方案策略

Answer 1

通过

访问数组后

double *p = (double*)PyArray_DATA(ArgsArray);
int* pA = (int*)PyArray_DATA(obj);

您可以像使用数组一样使用它。什么是维度？

int height = PyArray_DIM(obj, 0);
int width = PyArray_DIM(obj, 1);

现在您可以使用该指针访问数组中的数据

for ( int i = 0; y<height; y++) 
{ 
   for (int j = 0; x<width; x++) 
   { 
      p[i * width + j] = pA[i * width + j];
   }
}

其实如果只需要复制数组，此时使用memcpy或者std::copy即可。

Tbh，我曾考虑使用 boost.python，但这是我自己的偏好。