来自scipy CSR稀疏矩阵的访问值,列索引和row_ptr数据

时间:2012-04-28 03:02:12

标签: python scipy sparse-matrix csr

我有一个大矩阵,我想将其转换为稀疏的CSR格式。

当我这样做时:

import scipy as sp
Ks = sp.sparse.csr_matrix(A)

print Ks

A密集,我得

 (0, 0) -2116689024.0
 (0, 1) 394620032.0
 (0, 2) -588142656.0
 (0, 12)    1567432448.0
 (0, 14)    -36273164.0
 (0, 24)    233332608.0
 (0, 25)    23677192.0
 (0, 26)    -315783392.0
 (0, 45)    157961968.0
 (0, 46)    173632816.0

等...

我可以使用以下方法获取行索引,列索引和值的向量:

Knz = Ks.nonzero()
sparserows = Knz[0]
sparsecols = Knz[1]

#The Non-Zero Value of K at each (Row,Col) 
vals = np.empty(sparserows.shape).astype(np.float)
for i in range(len(sparserows)):

    vals[i] = K[sparserows[i],sparsecols[i]]

但是有可能提取稀疏CSR格式中包含的向量(值,列索引,行指针)吗?

SciPy的文档解释说,可以从这三个向量生成CSR矩阵,但我想反过来,将这三个向量输出。

我错过了什么?

感谢您的时间!

1 个答案:

答案 0 :(得分:16)

value = Ks.data
column_index = Ks.indices
row_pointers = Ks.indptr

我认为这些属性没有记录,可能会使它们发生变化,但我已经在几个版本的scipy中使用它们。