基于numpy python中的条件对二维数组进行子集化

Question

我有一个大小为 3600 * 7200 的 numpy 2D 数组。我有另一个形状相同的数组，我想用作掩码。

问题是当我做这样的事情时：

import numpy as np
N = 10
arr_a = np.random.random((N,N))
arr_b = np.random.random((N,N))
arr_a[arr_b > 0.0]

结果数组不再是二维的，而是一维的。 我如何得到一个二维数组作为回报？

Answer 1

您可以使用 np.where 来保留形状：

np.where(arr_b > 0.0, arr_a, np.nan)

当 arr_b 的值大于 0 时，它将从 arr_a 中获取相应的值，否则将使用 np.nan。

import numpy as np
N = 5
arr_a = np.random.randn(N,N)
arr_b = np.random.randn(N,N)
np.where(arr_b > 0.0, arr_a, np.nan)

Out[107]: 
array([[ 0.5743081 ,         nan, -1.69559034,         nan,  0.4987268 ],
       [ 0.33038264,         nan, -0.27151598,         nan, -0.73145628],
       [        nan,  0.46741932,  0.61225086,         nan,  1.08327459],
       [        nan, -1.20244926,  1.5834266 , -0.04675223, -1.14904974],
       [        nan,  1.20307104, -0.86777899,         nan,         nan]])