Question

numpy中是否有一种很好的方法可以获得array1中每个元素位于array2的元素指数？

一个例子：

array1 = np.array([1, 3, 4])
array2 = np.arange(-2, 5, 1, dtype=np.int)

np.where(array1[0] == array2)
# (array([3]),)
np.where(array1[1] == array2)
# (array([5]),)
np.where(array1[2] == array2)
# (array([6]),)

我想做

np.where(array1 == array2)
# (array([3 5 6]),)

这样的事情可能吗？我们保证array1中的所有条目都可以在array2中找到。

Answer 1

方法＃1：在那里使用np.in1d获取匹配发生位置的掩码，然后np.where获取这些索引位置 -

np.where(np.in1d(array2, array1))

方法＃2：使用np.searchsorted -

np.searchsorted(array2, array1)

请注意，如果array2未排序，我们需要使用附加的可选参数sorter。

示例运行 -

In [14]: array1
Out[14]: array([1, 3, 4])

In [15]: array2
Out[15]: array([-2, -1,  0,  1,  2,  3,  4])

In [16]: np.where(np.in1d(array2, array1))
Out[16]: (array([3, 5, 6]),)

In [17]: np.searchsorted(array2, array1)
Out[17]: array([3, 5, 6])

运行时测试 -

In [62]: array1 = np.random.choice(10000,1000,replace=0)

In [63]: array2 = np.sort(np.random.choice(100000,10000,replace=0))

In [64]: %timeit np.where(np.in1d(array2, array1))
1000 loops, best of 3: 483 µs per loop

In [65]: %timeit np.searchsorted(array2, array1)
10000 loops, best of 3: 40 µs per loop

查找两个numpy数组之间的索引映射

1 个答案: