Numpy.cumsum相反

时间:2015-02-19 21:38:11

标签: python numpy

这是前进方向的cumsum:

> import numpy as np
> np.arange(10)
array([0, 1, 2, 3, 4, 5, 6, 7, 8, 9])
> np.cumsum(np.arange(10))
array([ 0,  1,  3,  6, 10, 15, 21, 28, 36, 45])

我想反向执行cumsum,这会给我

array([45, 45, 44, 42, 39, 35, 30, 24, 17, 9])

最简单,最有效的方法是什么?

2 个答案:

答案 0 :(得分:5)

最简单的我能想到并产生你的结果

import numpy as np
x = np.arange(10)
x[::-1].cumsum()[::-1]

给出了

array([45, 45, 44, 42, 39, 35, 30, 24, 17,  9])

编辑:正如 dg99 指出的,还有a post about the efficiency of reversing an array。因此,[::-1]似乎是最好的。因此,x[::-1].cumsum()[::-1]似乎也是最有效的反向积累方式。

第二次编辑: 为了完整起见,如果您有一个多维数组,您可以通过以下方式获得最内层维度的反向cumsum:

x[...,::-1].cumsum(axis=-1)[...,::-1]

例如,

x = np.array(((range(10), range(10)), (range(10), range(10))))
print(x)

打印

array([[[0, 1, 2, 3, 4, 5, 6, 7, 8, 9],
    [0, 1, 2, 3, 4, 5, 6, 7, 8, 9]],

   [[0, 1, 2, 3, 4, 5, 6, 7, 8, 9],
    [0, 1, 2, 3, 4, 5, 6, 7, 8, 9]]])

 x[...,::-1].cumsum(axis=-1)[...,::-1]

返回

array([[[45, 45, 44, 42, 39, 35, 30, 24, 17,  9],
    [45, 45, 44, 42, 39, 35, 30, 24, 17,  9]],

   [[45, 45, 44, 42, 39, 35, 30, 24, 17,  9],
    [45, 45, 44, 42, 39, 35, 30, 24, 17,  9]]])

答案 1 :(得分:1)

怎么样

np.cumsum(np.arange(10)[::-1])

参见讨论逆转numpy数组的this question