脾气暴躁的地方

在 您的特定情况下 ，最好的方法 是 将两个条件更改为一个条件：

dists[abs(dists - r - dr/2.) <= dr/2.]

它仅创建一个布尔数组，在我看来是更易于阅读，因为它说，
是

dist

dr

r

r

r = r +dr/2.

问题的答案： 如果您只是想过滤出不符合标准的元素，则
实际上并不需要：

where``dists

dists[(dists >= r) & (dists <= r+dr)]

因为

&

and

或者，如果您

where

 dists[(np.where((dists >= r) & (dists <= r + dr)))]

原因：
不起作用的原因是因为

np.where

and

True

False

a

b

True

aand b

b

[0,1,2] and [2,3,4]

[2,3,4]

In [230]: dists = np.arange(0,10,.5)In [231]: r = 5In [232]: dr = 1In [233]: np.where(dists >= r)Out[233]: (array([10, 11, 12, 13, 14, 15, 16, 17, 18, 19]),)In [234]: np.where(dists <= r+dr)Out[234]: (array([ 0,  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12]),)In [235]: np.where(dists >= r) and np.where(dists <= r+dr)Out[235]: (array([ 0,  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12]),)

您期望比较的只是布尔数组，例如

In [236]: dists >= rOut[236]: array([False, False, False, False, False, False, False, False, False,       False,  True,  True,  True,  True,  True,  True,  True,  True,        True,  True], dtype=bool)In [237]: dists <= r + drOut[237]: array([ True,  True,  True,  True,  True,  True,  True,  True,  True,        True,  True,  True,  True, False, False, False, False, False,       False, False], dtype=bool)In [238]: (dists >= r) & (dists <= r + dr)Out[238]: array([False, False, False, False, False, False, False, False, False,       False,  True,  True,  True, False, False, False, False, False,       False, False], dtype=bool)

现在，您可以调用

np.where

In [239]: np.where((dists >= r) & (dists <= r + dr))Out[239]: (array([10, 11, 12]),)In [240]: dists[np.where((dists >= r) & (dists <= r + dr))]Out[240]: array([ 5. ,  5.5,  6. ])

或者使用花式索引简单地用布尔数组对原始数组进行索引

In [241]: dists[(dists >= r) & (dists <= r + dr)]Out[241]: array([ 5. ,  5.5,  6. ])