ndarray的稀疏邻接矩阵

Question

假设我们有一个尺寸为(Nx, Ny, Nz)的3维网格，对应于一个ndarray：

A = np.zeros((Nx,Ny,Nz))

如果对A进行矢量化处理，我们将得到一个尺寸为(1, Nx * Ny * Nz)的相应数组：

vec_A = A.reshape((1,Nx x Ny x Nz))

我很好奇一个人如何在scipy中为网格的矢量化元素构造一个稀疏的邻接矩阵。也就是说，如果vec_A的元素m具有坐标(i, j, k)，则{-1,0,1}x{-1,0,1}{-1,0,1}中向量的所有可能组合将具有关联的矢量化坐标n，然后对{{1} }包含在稀疏矩阵中。

在我的问题(n,m)，Nx和Ny中，大约有200个，给出了长度为8000000的矢量化版本。