numpy.where返回结果不一致求助：元素存在却返回空数组

你遇到的这个问题确实挺诡异的——明明for循环已经定位到目标元素，同一行的np.where能正常返回索引，可下一行的np.where却直接返回空数组，触发索引越界错误。结合你的代码和报错信息，我整理出两个最可能的原因：

1. 变量名拼写错误（直接触发报错的元凶） #

仔细看你的报错行代码：

iind = np.where(progenitors[oind]==-progenitors[out][i])[0][0]

再对比前面能正常工作的print语句：

print "found", p, np.where(progenitors[oind]==-progenitors[ozi][i])[0][0]

这里你把变量ozi写成了out！progenitors[out][i]对应的元素和你要找的-p完全不匹配，np.where自然找不到结果，返回空数组后访问[0][0]就会抛出IndexError。

这应该是导致当前报错的直接原因，先把这个变量名的拼写错误修复了再说。

2. 浮点数精度问题（潜在的隐藏坑） #

如果修复变量名后还是出现类似的奇怪问题，那就要考虑浮点数精度匹配的坑：

• 虽然你的输出里显示的都是整数，但如果progenitors数组的 dtype 是浮点数（比如float64），元素实际存储的是类似76.0的形式。如果数组里的元素因为文件读取或计算产生了微小的精度误差（比如76.0000000001或75.9999999999），直接用==比较就会失败。
• 而for循环里的p == -progenitors[ozi][i]能成功，可能是因为p是从数组里直接取出来的浮点数，和-progenitors[ozi][i]的精度误差完全一致，所以比较相等；但如果后续的比较值是重新计算的（比如整数转浮点数），就可能出现不匹配的情况。

针对这种情况，你可以这么处理：

• 先检查数组类型：用print(progenitors.dtype)查看，如果是浮点数且确定所有元素都是整数，直接转成整数类型：

  progenitors = progenitors.astype(int)
  

• 或者用numpy的isclose函数做近似比较，规避精度问题：

  iind = np.where(np.isclose(progenitors[oind], -progenitors[ozi][i]))[0][0]
  

额外的调试小技巧 #

为了避免这类摸不着头脑的问题，你可以在报错前加一些调试输出，确认比较的两个值是否真的相等：

target_val = -progenitors[ozi][i]
print "Target value:", target_val, "Type:", type(target_val)
print "First element in array:", progenitors[oind][0], "Type:", type(progenitors[oind][0])
# 用for循环找到的索引去验证值是否匹配
match_idx = np.where(progenitors[oind]==p)[0][0]
print "Are target and matched value equal?", target_val == progenitors[oind][match_idx]

这样能直观看到比较的两个值是否存在类型或精度上的差异。

内容的提问来源于stack exchange，提问作者mivkov