制作较小的3D数组：

In [155]: a = np.arange(24).reshape(2,3,4)
In [158]: a
Out[158]: 
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]]])

选择两个“行”（在中间尺寸上）：

In [159]: a[:,[0,2],:]
Out[159]: 
array([[[ 0,  1,  2,  3],
        [ 8,  9, 10, 11]],

       [[12, 13, 14, 15],
        [20, 21, 22, 23]]])

如果我们使用相同形状的2个列表（或数组），我们最终会从[159]中选择2个“行”：

In [160]: a[[0,1],[0,2],:]
Out[160]: 
array([[ 0,  1,  2,  3],
       [20, 21, 22, 23]])

如果第一个列表/数组是“列向量”，我们选择（2，2）“块”：

In [161]: a[[[0],[1]],[0,2],:]
Out[161]: 
array([[[ 0,  1,  2,  3],
        [ 8,  9, 10, 11]],

       [[12, 13, 14, 15],
        [20, 21, 22, 23]]])

ix_可用于创建相同的2个数组：

In [162]: np.ix_([0,1],[0,2])
Out[162]: 
(array([[0],
        [1]]),
 array([[0, 2]]))

因此，使用ix_数组：

In [163]: I,J = np.ix_([0,1],[0,2])
In [164]: a[I,J,:]
Out[164]: 
array([[[ 0,  1,  2,  3],
        [ 8,  9, 10, 11]],

       [[12, 13, 14, 15],
        [20, 21, 22, 23]]])

当我说它们相互广播时，我的意思与在加法或乘法期间广播的意思相同：

In [165]: I*10 + J
Out[165]: 
array([[ 0,  2],
       [10, 12]])

参考：https://numpy.org/doc/stable/user/basics.indexing.html#advanced-indexing

编辑

In [166]: np.ix_([7,10,12,13,14], [1,4])
Out[166]: 
(array([[ 7],
        [10],
        [12],
        [13],
        [14]]),
 array([[1, 4]]))

关于您的错误：

In [167]: np.ix_([7,10,12,13,14], [1,4],:,:)
  Input In [167]
    np.ix_([7,10,12,13,14], [1,4],:,:)
                                  ^
SyntaxError: invalid syntax

ix_是一个函数。在函数调用中不允许使用“：”。它只在索引中有效，在索引中它被转换为slice。这就是为什么会出现语法错误。