如何为多维NumPy数组创建布尔索引数组?

时间:2013-06-22 16:00:25

标签: arrays numpy indexing

我正在看这个问题:Efficient thresholding filter of an array with numpy

我有类似的问题,但有一个二维数组,其中几个坐标和值存储在NumPy数组的行中。我如何能够像上面的问题那样进行类似的过滤?

我的数据如下:

>>> A.shape
(4269862, 5)

>>> A[0]
array([  2.27011719e+02,   0.00000000e+00,   2.88134766e+02,
         2.00000000e+00,   7.69880000e+04], dtype=float32)

这些值对应于X,Y,Z和value1和value2。我想要的是有效地获得例如。所有行的X在300-400,Y在200-250,Z在200-300。

1 个答案:

答案 0 :(得分:2)

您可以创建一个布尔蒙版,在满足所有条件时为真:

idx = ((A[:, 0] > 300) & (A[:, 0] < 400) & 
       (A[:, 1] > 200) & (A[:, 1] < 250) & 
       (A[:, 2] > 200) & (A[:, 2] < 300))

print A[idx]   
# this should give your array rows where idx is True

您可以测试一下:

A = np.random.uniform(150, 500, (200, 5)).astype('i')
idx = ((A[:, 0] > 300) & (A[:, 0] < 400) & 
       (A[:, 1] > 200) & (A[:, 1] < 250) & 
       (A[:, 2] > 200) & (A[:, 2] < 300))
print A[idx]
#[[339 292 231 211 474]
# [371 252 310 281 256]
# [337 263 471 159 397]
# [361 299 383 250 206]
# [360 278 328 194 453]
# [360 258 205 245 427]
# [339 286 331 175 418]]