我试图找到一个整齐的解决方案,但我正在以相同的方式切片相同形状的几个2D数组。我已经通过定义包含'x,y'中心的列表来尽可能地整理它,例如cpix = [161, 134]
我想要做的是不必像这样写三遍切片:
a1 = array1[cpix[1]-50:cpix[1]+50, cpix[0]-50:cpix[0]+50]
a2 = array2[cpix[1]-50:cpix[1]+50, cpix[0]-50:cpix[0]+50]
a3 = array3[cpix[1]-50:cpix[1]+50, cpix[0]-50:cpix[0]+50]
只是预定义了一些东西(例如面具?),所以我可以做一个
a1 = array1[predefined_2dslice]
a2 = array2[predefined_2dslice]
a3 = array3[predefined_2dslice]
这是numpy支持的吗?
最佳答案
是的,您可以使用 numpy.s_
:
例子:
>>> a = np.arange(10).reshape(2, 5)
>>>
>>> m = np.s_[0:2, 3:4]
>>>
>>> a[m]
array([[3],
[8]])
在这种情况下:
my_slice = np.s_[cpix[1]-50:cpix[1]+50, cpix[0]-50:cpix[0]+50]
a1 = array1[my_slice]
a2 = array2[my_slice]
a3 = array3[my_slice]
您还可以使用
numpy.r_
来将切片对象转换为沿第一个轴的串联。