过滤numpy数组以仅保留给定值的一行

Question

我有一个较大的nx 2 numpy数组，其格式设置为（x，y）坐标。我想过滤此数组，以便：

1. 标识具有重复x值的坐标对。
2. 仅保留y值最高的重复项的坐标对。

例如，在以下数组中：

arr = [[1, 4]
       [1, 8]
       [2, 3]
       [4, 6]
       [4, 2]
       [5, 1]
       [5, 2]
       [5, 6]]

我希望结果是：

arr = [[1, 8]
       [2, 3]
       [4, 6]
       [5, 6]]

Ive探索了np.unique和np.where，但无法弄清楚如何利用它们来解决此问题。非常感谢！

Answer 1

这是一种基于np.maximum.reduceat-

def grouby_maxY(a):
    b = a[a[:,0].argsort()] # if first col is already sorted, skip this
    grp_idx = np.flatnonzero(np.r_[True,(b[:-1,0] != b[1:,0])])
    grp_maxY = np.maximum.reduceat(b[:,1], grp_idx)
    return np.c_[b[grp_idx,0], grp_maxY]

另外，如果你要带np.unique，我们可以用它来寻找grp_idx有np.unique(b[:,0], return_index=1)[1]。

样品运行-

In [453]: np.random.seed(0)

In [454]: arr = np.random.randint(0,5,(10,2))

In [455]: arr
Out[455]: 
array([[4, 0],
       [3, 3],
       [3, 1],
       [3, 2],
       [4, 0],
       [0, 4],
       [2, 1],
       [0, 1],
       [1, 0],
       [1, 4]])

In [456]: grouby_maxY(arr)
Out[456]: 
array([[0, 4],
       [1, 4],
       [2, 1],
       [3, 3],
       [4, 0]])