对给定列的numpy矩阵排序

2024-10-01 09:32:13 发布

您现在位置:Python中文网/ 问答频道 /正文

我尝试使用Ned Batchelder代码按人的顺序对NumPy矩阵进行排序,如下文所述:

Sort numpy string array with negative numbers?

代码在一维数组上运行,命令是:

print (sorted(a, key=natural_keys))

现在,我的问题是我的数据是一个10列的矩阵,我想根据一列对它进行排序(比如MyColumn)。我找不到一种方法来修改代码以打印按此列排序的整个矩阵。我能想到的就是:

^{pr2}$

但是,当然,只有MyColumn出现在输出中,尽管它是正确排序的。。。在

有办法打印整个矩阵吗?

下面是我用来加载数组的命令(我将原来的imputfile简化为一个3列数组):

data = np.loadtxt(inputfile, dtype={'names': ('ID', 'MyColumn', 'length'),
'formats': ('int32', 'S40', 'int32')},skiprows=1, delimiter='\t')

ID  MyColumn    length
164967  BFT_job13_q1_type2  426
197388  BFT_job8_q0_type2   244
164967  BFT_job13_q0_type1  944
72406   BFT_job1_q0_type3   696

以下是理想情况下的输出:

ID  MyColumn    length
72406   BFT_job1_q0_type3   696
197388  BFT_job8_q0_type2   244
164967  BFT_job13_q0_type1  944
164967  BFT_job13_q1_type2  426

Tags: 代码命令id排序矩阵数组lengthq1
1条回答
网友
1楼 · 发布于 2024-10-01 09:32:13

如果您有一个np.matrix,名为m

col = 1
m[np.array(m[:,col].argsort(axis=0).tolist()).ravel()]

如果您有一个np.ndarray,称为a

^{pr2}$

如果您有一个具有命名列的结构化数组:

def mysort(data, col_name, key=None):
    d = data.copy()
    cols = [i[0] for i in eval(str(d.dtype))]
    if key:
        argsort = np.array([key(i) for i in d[col_name]]).argsort()
    else:
        argsort = d[col_name].argsort()
    for col in cols:
        d[col] = d[col][argsort]
    return d

对于您的具体情况,您需要the following ^{} function

def key(x):
    x = ''.join([i for i in x if i.isdigit() or i=='_'])
    return '{1:{f}{a}10}_{2:{f}{a}10}_{3:{f}{a}10}'.format(*x.split('_'), f='0', a='>')

d = mysort(data, 'MyColumn', key)

相关问题 更多 >