查找稀疏csc_matrix中存在非零条目的行的索引

Question

我有一个numpy数组，X：

type(X)
>>> <class 'scipy.sparse.csc.csc_matrix'>

我感兴趣的是找到第0列中存在非零条目的行的索引。我试过了：

getcol =  X.getcol(0)
print getcol

给了我：

(0, 0)  1
(2, 0)  1
(5, 0)  10

这很棒，但我想要的是一个载有0, 2, 5的向量。

我如何获得我正在寻找的指数？

感谢您的帮助。

Answer 1

使用CSC矩阵，您可以执行以下操作：

>>> import scipy.sparse as sps
>>> a = np.array([[1, 0, 0],
...               [0, 1, 0],
...               [1, 0, 1],
...               [0, 0, 1],
...               [0, 1, 0],
...               [1, 0, 1]])
>>> aa = sps.csc_matrix(a)
>>> aa.indices[aa.indptr[0]:aa.indptr[1]]
array([0, 2, 5])
>>> aa.indices[aa.indptr[1]:aa.indptr[2]]
array([1, 4])
>>> aa.indices[aa.indptr[2]:aa.indptr[3]]
array([2, 3, 5])

所以aa.indices[aa.indptr[col]:aa.indptr[col+1]]可以让你得到你想要的东西。