你能在python中增加一个切片吗?

时间:2012-04-16 18:54:54

标签: python numpy

我正在通过移植我在matlab / octave中编写的一些代码来尝试numpy。在matlab中,我可以定义python切片的等价物,然后根据需要增加它。例如,在我的matlab代码中我有

HXx_range = 1:NHXx;
HXy_range = 1:NHXy;

blah blah blah

Hx(HXx_range, HXy_range) = Da(Hx_media(HXx_range, HXy_range)).*Hx(HXx_range, HXy_range) + Db(Hx_media(HXx_range, HXy_range)).*(Ez(HXx_range,HXy_range) -   Ez(HXx_range,**HXy_range+1**));
Hy(HYx_range, HYy_range) = Da(Hy_media(HYx_range, HYy_range)).*Hy(HYx_range, HYy_range) + Db(Hy_media(HYx_range, HYy_range)).*(Ez(**HYx_range+1**,HYy_range) - Ez(HYx_range,HYy_range));
Ez(EZx_range, EZy_range) = Ca(Ez_media(EZx_range, EZy_range)).*Ez(EZx_range, EZy_range) + Cb(Ez_media(EZx_range, EZy_range)).*(Hy(EZx_range,EZy_range) - Hy(**EZx_range-1**,EZy_range) + Hx(EZx_range,**EZy_range-1**) - Hx(EZx_range,EZy_range));

'**'中的术语(如'HXy_range + 1')是关键部分; HXy_range + 1等于2:(NHXy + 1)。在python中,我可以用类似的方式定义切片:

HXx_range = slice(0, NHXx)

然而,HXx_range + 1给了我一个错误。当然,我可以为此做一个新的切片,但它不是那么干净。有办法解决这个问题吗?

感谢。

2 个答案:

答案 0 :(得分:3)

如果将HXy_range定义为numpy数组,则可以根据需要增加它。如果您希望将其用作切片,可以形成slice(*HXy_range)

In [26]: HXy_range = np.array([1,10])

In [27]: HXy_range+1
Out[27]: array([ 2, 11])

In [28]: slice(*(HXy_range+1))
Out[30]: slice(2, 11, None)

答案 1 :(得分:1)

不,Python slice个实例是immutable。要使用标准slice实例,您必须每次都创建一个新实例。不幸的是,you can't subclass slice either