从2D NumPy数组中查找一维数组?

时间:2016-08-11 20:14:34

标签: python python-2.7 numpy

假设我有

[[array([x1, y1]), z1]
 [array([x2, y1]), z2]
 ......
 [array([xn, yn]), zn]
]

我希望找到array([x5, y5])的索引。如何才能有效地使用NumPy?

1 个答案:

答案 0 :(得分:0)

首先,由于混合数据格式,我认为你不能以矢量化的方式提取数组。因此,您可以使用loop comprehension从每个列表元素中提取与数组对应的第一个元素作为2D数组。所以,假设A是输入列表,我们会有 -

arr = np.vstack([a[0] for a in A])

然后,只需使用NumPy's broadcasting feature以矢量化方式进行比较,因为它将沿所有行广播该比较,并使用np.all(axis=1)查看所有匹配的行。最后,使用np.flatnonzero获取最终索引。因此,谜题的最终和平将是 -

idx = np.flatnonzero((arr == search1D).all(1))

您可以阅读this post的答案,了解在1D数组问题中搜索此类2D数组中的索引的其他替代方法。

示例运行 -

In [140]: A
Out[140]: 
[[array([3, 4]), 11],
 [array([2, 1]), 12],
 [array([4, 2]), 16],
 [array([2, 1]), 21]]

In [141]: search1D = [2,1]

In [142]: arr = np.vstack([a[0] for a in A]) # Extract 2D array

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

In [144]: np.flatnonzero((arr == search1D).all(1)) # Finally get indices
Out[144]: array([1, 3])