Python稀疏矩阵删除重复索引除了一个?

时间:2017-03-31 18:01:17

标签: python matrix scipy sparse-matrix

我正在计算向量矩阵之间的余弦相似度,我得到的结果是稀疏矩阵,如下所示:

  
      
  • (0,26)0.359171459261
  •   
  • (0,25)0.121145761751
  •   
  • (0,24)0.316922015914
  •   
  • (0,23)0.157622038039
  •   
  • (0,22)0.636466644041
  •   
  • (0,21)0.136216495731
  •   
  • (0,20)0.243164535496
  •   
  • (0,19)0.348272617805
  •   
  • (0,18)0.636466644041
  •   
  • (0,17)1.0
  •   

但是有重复的例子:

  

(0,24)0.316922015914和(24,0)0.316922015914

我想做的是通过indice删除它们(如果我有(0,24)那么我不需要(24,0)因为它是相同的)只留下其中一个并删除第二,对于矩阵中的所有向量。 目前我有以下代码来创建矩阵:

vectorized_words = sparse.csr_matrix(vectorize_words(nostopwords,glove_dict))
cos_similiarity = cosine_similarity(vectorized_words,dense_output=False)

总而言之,我不想删除所有重复项,我希望只使用其中一个使用pythonic方式。

提前谢谢!

1 个答案:

答案 0 :(得分:2)

我认为最简单的方法是获得coo格式矩阵的上三角形:

首先制作一个小的对称矩阵:

In [876]: A = sparse.random(5,5,.3,'csr')
In [877]: A = A+A.T
In [878]: A
Out[878]: 
<5x5 sparse matrix of type '<class 'numpy.float64'>'
    with 11 stored elements in Compressed Sparse Row format>
In [879]: A.A
Out[879]: 
array([[ 0.        ,  0.        ,  0.81388978,  0.        ,  0.        ],
       [ 0.        ,  0.        ,  0.73944395,  0.20736975,  0.98968617],
       [ 0.81388978,  0.73944395,  0.        ,  0.        ,  0.        ],
       [ 0.        ,  0.20736975,  0.        ,  0.05581152,  0.04448881],
       [ 0.        ,  0.98968617,  0.        ,  0.04448881,  0.        ]])

转换为coo,并将较低三角形数据值设置为0

In [880]: Ao = A.tocoo()
In [881]: mask = (Ao.row>Ao.col)
In [882]: mask
Out[882]: 
array([False, False, False, False,  True,  True,  True, False, False,
        True,  True], dtype=bool)
In [883]: Ao.data[mask]=0

转换回0,并使用eliminate_zeros修剪矩阵。

In [890]: A1 = Ao.tocsr()
In [891]: A1
Out[891]: 
<5x5 sparse matrix of type '<class 'numpy.float64'>'
    with 11 stored elements in Compressed Sparse Row format>
In [892]: A1.eliminate_zeros()
In [893]: A1
Out[893]: 
<5x5 sparse matrix of type '<class 'numpy.float64'>'
    with 6 stored elements in Compressed Sparse Row format>
In [894]: A1.A
Out[894]: 
array([[ 0.        ,  0.        ,  0.81388978,  0.        ,  0.        ],
       [ 0.        ,  0.        ,  0.73944395,  0.20736975,  0.98968617],
       [ 0.        ,  0.        ,  0.        ,  0.        ,  0.        ],
       [ 0.        ,  0.        ,  0.        ,  0.05581152,  0.04448881],
       [ 0.        ,  0.        ,  0.        ,  0.        ,  0.        ]])

coocsr格式都采用就地eliminate_zeros方法。

def eliminate_zeros(self):
    """Remove zero entries from the matrix

    This is an *in place* operation
    """
    mask = self.data != 0
    self.data = self.data[mask]
    self.row = self.row[mask]
    self.col = self.col[mask]

您可以将此代码作为仅消除lower_triangle值的模型,而不是使用Ao.data[mask]=0