我一直在尝试执行一个简单的操作,但我似乎找不到使用Numpy函数的简单方法,而不会创建不必要的数组副本。
假设我们有以下三维数组:
In [171]: x = np.arange(24).reshape((4, 3, 2))
In [172]: x
Out[172]:
array([[[ 0, 1],
[ 2, 3],
[ 4, 5]],
[[ 6, 7],
[ 8, 9],
[10, 11]],
[[12, 13],
[14, 15],
[16, 17]],
[[18, 19],
[20, 21],
[22, 23]]])
以下数组:
In [173]: y = np.array([0, 1, 1, 0])
我想在x
中为每一行选择最后一个维度的值,其索引是y
中的对应元素。换句话说,我想:
array([[ 0, 2, 4],
[ 7, 9, 11],
[13, 15, 17],
[18, 20, 22]])
我现在唯一的解决方案是在x
和y
的第一维上使用for循环,如下所示:
z = np.zeros((4, 3), dtype=int)
for i, row in enumerate(x):
z[i, :] = row[:, y[i]]
有没有办法在这里避免使用for循环,使用numpy函数或花哨的索引?
谢谢!
答案 0 :(得分:5)
棘手的方面是你不希望每个切片的所有第0维,你希望切片对应到第0维中的每个元素。所以你可以这样做:
>>> x[np.arange(x.shape[0]), :, y]
array([[ 0, 2, 4],
[ 7, 9, 11],
[13, 15, 17],
[18, 20, 22]])
答案 1 :(得分:1)
花式索引:
__future__
给出:
x[np.arange(y.size),:,y]