三个 2d 轮廓到 3d 阵列 [重复]

Question

numpy 中是否有可用的 staking 方法将三个 2d 配置文件转换为 3d 数组？其中 x 是直视视图，y 是水平视图，z 是鸟瞰图。因此，如果在 3d 空间中，x、y 和 z 中的每个值都等于 1，则保持 1 的值，否则给出 0 的值。所以在这个小例子中，我相信应该有 108 个唯一值（9*12） ...

import numpy as np

x = np.array([[1, 1, 1],
              [0, 1, 1],
              [0, 1, 1]])

y = np.array([[1, 1, 1],
             [ 1, 1, 1],
             [ 0, 1, 1],
             [ 0, 0, 1]])

z = np.array([[0, 1, 1],
             [ 0, 1, 1],
             [ 0, 0, 1],
             [ 0, 0, 0]])

即如果我们从x 数组在[0, 0] 开始，值是1，在位置y[0,0] 值也是1 但是在z[0, 0] 值是0 所以@987654330 的3d 数组@ 应该是 0 [0,0,0]。使用x 和y 的相同位置，但将z 更改为z[0, 1] 也是1 的值，因此xyz 位置[0, 0, 1] 的3d 数组应为1。所以它应该是一些排列，本质上是创建一个矩形 3D 网格。

Answer 1

您可以结合广播和按位/逻辑运算来执行此操作：

>>> x[np.newaxis, :, :] & y[:, np.newaxis, :] & z[:, :, np.newaxis]
array([[[0, 0, 0],
        [0, 1, 1],
        [0, 1, 1]],

       [[0, 0, 0],
        [0, 1, 1],
        [0, 1, 1]],

       [[0, 0, 0],
        [0, 0, 0],
        [0, 1, 1]],

       [[0, 0, 0],
        [0, 0, 0],
        [0, 0, 0]]])
>>>