15

我有一个大矩阵,我想将其转换为稀疏 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

ETC...

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

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 矩阵,但我想做相反的事情,将这三个向量取出。

我错过了什么?

谢谢你的时间!

4

1 回答 1

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

我相信这些属性是无证的,这可能会使它们发生变化,但我已经在几个版本的 scipy 上使用了它们。

于 2012-04-28T15:30:23.787 回答