Sizin durumunuzda en iyi yol , iki kriterinizi bir kriter olarak değiştirmektir:
dists[abs(dists - r - dr/2.) <= dr/2.]
Sadece bir boolean dizi oluşturur ve bence, diyor çünkü okumak daha kolaydır olan dist
bir mesafede dr
ya r
? ( r
Başlangıçtan ziyade ilgi alanınızın merkezi olarak yeniden tanımlasam da r = r + dr/2.
) Ancak bu, sorunuzun yanıtını vermiyor.
Sorunuzun cevabı: Sadece kriterlerinize uymayan unsurları filtrelemeye çalışıyorsanız,
aslında ihtiyacınız yok:where
dists
dists[(dists >= r) & (dists <= r+dr)]
Çünkü &
size bir elementwise verecektir and
(parantezler gereklidir).
Veya where
herhangi bir nedenle kullanmak istiyorsanız , şunları yapabilirsiniz:
dists[(np.where((dists >= r) & (dists <= r + dr)))]
Nedeni:
Çalışmamasının nedeni np.where
, bir boole dizisi değil, bir dizin listesi döndürmesidir. and
Elbette beklediğiniz True
/ False
değerlerine sahip olmayan iki sayı listesi arasına girmeye çalışıyorsunuz . Eğer a
ve b
ikisi de True
değerler, daha sonra a and b
döner b
. Yani böyle bir şey söylemek [0,1,2] and [2,3,4]
size verecektir [2,3,4]
. İşte eylemde:
In [230]: dists = np.arange(0,10,.5)
In [231]: r = 5
In [232]: dr = 1
In [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]),)
Karşılaştırmayı beklediğiniz şey basitçe boole dizisiydi, örneğin
In [236]: dists >= r
Out[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 + dr
Out[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)
Şimdi np.where
birleşik boole dizisini çağırabilirsiniz :
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. ])
Veya basit bir şekilde orijinal diziyi süslü indekslemeyi kullanarak boole dizisiyle indeksleyin
In [241]: dists[(dists >= r) & (dists <= r + dr)]
Out[241]: array([ 5. , 5.5, 6. ])
()
etrafta dolaşmak zorundasın(ar>3)
ve(ar>6)
?