2012-08-22 40 views
5

我有一个矩阵,它应该在对角线上有一个矩阵,但列是混合起来的。根据其对角线排列一个numpy矩阵

Messed up matrix

但我不知道怎么回事,没有明显的循环,有效地互换行以获得在对角线上的统一。我甚至不知道我会通过什么样的关键排序。

有什么建议吗?

+2

因为他们'重新漂浮,你能保证他们的列正好是1.0和唯一的吗? – wim

+0

@wim:存在浮点错误的问题。但是,沿着对角线的每个条目都保证在该行中具有最高值。 – mac389

回答

6

您可以使用numpy的的argmax确定了球门立柱订购和使用argmax结果列索引重新排列矩阵:

>>> z = numpy.array([[ 0.1 , 0.1 , 1. ], 
...     [ 1. , 0.1 , 0.09], 
...     [ 0.1 , 1. , 0.2 ]]) 

numpy.argmax(z, axis=1) 

>>> array([2, 0, 1]) #Goal column indices 

z[:,numpy.argmax(z, axis=1)] 

>>> array([[ 1. , 0.1 , 0.1 ], 
...  [ 0.09, 1. , 0.1 ], 
...  [ 0.2 , 0.1 , 1. ]]) 
+0

非常好。对于1000x1000阵列,您的运行时间为0.1s,我的运行时间为10s。 – Snowball

3
>>> import numpy as np 
>>> a = np.array([[ 1. , 0.5, 0.5, 0. ], 
...    [ 0.5, 0.5, 1. , 0. ], 
...    [ 0. , 1. , 0. , 0.5], 
...    [ 0. , 0.5, 0.5, 1. ]]) 
>>> np.array(sorted(a, cmp=lambda x, y: list(x).index(1) - list(y).index(1))) 
array([[ 1. , 0.5, 0.5, 0. ], 
     [ 0. , 1. , 0. , 0.5], 
     [ 0.5, 0.5, 1. , 0. ], 
     [ 0. , 0.5, 0.5, 1. ]]) 

它实际上按行排序,而不是列(但结果是相同的)。它的工作原理是由列的1是指数排序

+0

标准库再次获胜。任何方式来做到这一点,而不必使用numpy列表?我可以想象这对于大型阵列来说相当缓慢。 – jozzas

+0

@jozzas:它实际上并不依赖于numpy。如果你拿出'np.array'部件,它会工作得很好。 – Snowball

+0

对不起,我的意思是Python列表。是否有一个只有numpy的解决方案? – jozzas