我有一个相当大的一维numpy数组,我想在其中的一个切片上应用一些排序,并检索置换向量以进行其他处理。
但是,ndarray.sort()方法(这是一个就地操作)不返回此向量,我可以使用ndarray.argsort()方法获取置换向量并使用它对切片进行置换。但是,我不知道如何在原地做。
Vslice = V[istart:istop] # This is a view of the slice
iperm = Vslice.argsort()
V[istart:istop] = Vslice[iperm] # Not an inplace operation...
附属问题:为什么在我们研究v的视图时,下面的代码不修改v?
Vslice = Vslice[iperm]
最美好的祝福!
弗朗索瓦
最佳答案
要回答“为什么要查看的作业不修改原始作业”的问题,请执行以下操作:
您需要将Vslice = Vslice[iperm]
更改为Vslice[:] = Vslice[iperm]
,否则将为Vslice
分配一个新值,而不是更改Vslice
中的值:
>>> a = np.arange(10, 0, -1)
>>> a
array([10, 9, 8, 7, 6, 5, 4, 3, 2, 1])
>>> b = a[2:-2]
>>> b
array([8, 7, 6, 5, 4, 3])
>>> i = b.argsort()
>>> b[:] = b[i] # change the values inside the view
>>> a # note `a` has been sorted in [2:-2] slice
array([10, 9, 3, 4, 5, 6, 7, 8, 2, 1])