Python NumPy - 3D 数组的角度切片
Python NumPy - angled slice of 3D array
在 NumPy 中工作,我了解如何使用 this article 从 3D 数组中切出 2D 数组。
取决于我想要切入的轴:
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]
[24 25 26]]]
切片会给我:
i_slice = array[0]
[[0 1 2]
[3 4 5]
[6 7 8]]
j_slice = array[:, 0]
[[0 1 2]
[9 10 11]
[18 19 20]]
k_slice = array[:, :, 0]
[[0 3 6]
[9 12 15]
[18 21 24]]
但是45度角切片可以吗?如:
j_slice_down = array[slice going down starting from index 0]
[[0 1 2]
[12 13 14]
[24 25 26]]
我能够在所有 3 个轴的上升或下降中实现这一点,甚至一直环绕......在黑名单的日子里有很多 for 循环......但我确定必须有NumPy 中更好的方法。
更新,选择的答案:
我选择了 hpaulj 的答案,用 np.arrange
创建了两个坐标数组。通过一些工作,我能够满足我的需求,即以任何角度、任何轴、3D 阵列的非对称维度和任何位置返回切片,包括一直环绕以使其具有相同的维度作为轴。
为 x
和 y
创建了两个 np.arrange
数组。
np.roll
、递增、np.hstack
和 np.concatenate
等不同方法在 np.arrange
数组 x 轴数组上完成。 y = y[::-1]
为备用角度。
if axis is 'z': #i
slice_notation = np.index_exp[x, y, :]
elif axis is 'y': #k
slice_notation = np.index_exp[x, :, y]
else: #j
slice_notation = np.index_exp[:, x, y]
创建切片表达式,然后我使用 slice_notation
就地执行我需要的操作。
建议的其他方法:np.diagonal
和 np.eye
可能更适合其他人,因为他们的要求可能与我不同。
你可以试试 np.diagonal
:
arr = np.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],
[24 ,25 ,26]]])
np.diagonal(arr).T
array([[ 0, 1, 2],
[12, 13, 14],
[24, 25, 26]])
In [145]: arr[np.arange(3), np.arange(3),:]
Out[145]:
array([[ 0, 1, 2],
[12, 13, 14],
[24, 25, 26]])
在 NumPy 中工作,我了解如何使用 this article 从 3D 数组中切出 2D 数组。
取决于我想要切入的轴:
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]
[24 25 26]]]
切片会给我:
i_slice = array[0]
[[0 1 2]
[3 4 5]
[6 7 8]]
j_slice = array[:, 0]
[[0 1 2]
[9 10 11]
[18 19 20]]
k_slice = array[:, :, 0]
[[0 3 6]
[9 12 15]
[18 21 24]]
但是45度角切片可以吗?如:
j_slice_down = array[slice going down starting from index 0]
[[0 1 2]
[12 13 14]
[24 25 26]]
我能够在所有 3 个轴的上升或下降中实现这一点,甚至一直环绕......在黑名单的日子里有很多 for 循环......但我确定必须有NumPy 中更好的方法。
更新,选择的答案:
我选择了 hpaulj 的答案,用 np.arrange
创建了两个坐标数组。通过一些工作,我能够满足我的需求,即以任何角度、任何轴、3D 阵列的非对称维度和任何位置返回切片,包括一直环绕以使其具有相同的维度作为轴。
为 x
和 y
创建了两个 np.arrange
数组。
np.roll
、递增、np.hstack
和 np.concatenate
等不同方法在 np.arrange
数组 x 轴数组上完成。 y = y[::-1]
为备用角度。
if axis is 'z': #i
slice_notation = np.index_exp[x, y, :]
elif axis is 'y': #k
slice_notation = np.index_exp[x, :, y]
else: #j
slice_notation = np.index_exp[:, x, y]
创建切片表达式,然后我使用 slice_notation
就地执行我需要的操作。
建议的其他方法:np.diagonal
和 np.eye
可能更适合其他人,因为他们的要求可能与我不同。
你可以试试 np.diagonal
:
arr = np.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],
[24 ,25 ,26]]])
np.diagonal(arr).T
array([[ 0, 1, 2],
[12, 13, 14],
[24, 25, 26]])
In [145]: arr[np.arange(3), np.arange(3),:]
Out[145]:
array([[ 0, 1, 2],
[12, 13, 14],
[24, 25, 26]])