1

我一定是犯了某种非常微不足道的错误,但是我正在尝试创建一个带有单个轴名称的结构化数组,例如,我有一个data带有 shape的数组(2, 3, 4),并且我想命名第一个轴以便我可以访问data['a']并且data['b']在这两种情况下都得到(3, 4)成形切片。我试过:

shape = (2, 3, 4)
data = np.arange(np.product(shape)).reshape(shape)

dtype = [(nn, float) for nn in ['a', 'b']]
data = np.array(data, dtype=dtype)

但这似乎将所有数据复制到“a”和“b”中,例如

print(data.shape)
print(data['a'].shape)
> (2, 3, 4)
> (2, 3, 4)

我尝试指定形状(在 dtype 规范中)应该是(3, 4),但重复数据 12 次以上......我尝试将轴顺序更改为(3, 4, 2),但这没有任何作用。任何帮助表示赞赏!

4

1 回答 1

2
In [263]: shape = (2, 3, 4)
     ...: data = np.arange(np.product(shape)).reshape(shape)
     ...: 
     ...: dtype = [(nn, float) for nn in ['a', 'b']]

虽然可以转换data,但更可靠的方法是制作所需的目标数组,并将值复制到其中:

In [264]: res = np.zeros(shape[1:], dtype)
In [265]: res['a'] = data[0]
In [266]: res['b'] = data[1]
In [267]: res
Out[267]: 
array([[( 0., 12.), ( 1., 13.), ( 2., 14.), ( 3., 15.)],
       [( 4., 16.), ( 5., 17.), ( 6., 18.), ( 7., 19.)],
       [( 8., 20.), ( 9., 21.), (10., 22.), (11., 23.)]],
      dtype=[('a', '<f8'), ('b', '<f8')])
In [268]: res['a'].shape
Out[268]: (3, 4)

在这个结构化数组中,一条记录由 2 个浮点数组成,数据缓冲区包含:

In [272]: res.view(float).ravel()
Out[272]: 
array([ 0., 12.,  1., 13.,  2., 14.,  3., 15.,  4., 16.,  5., 17.,  6.,
       18.,  7., 19.,  8., 20.,  9., 21., 10., 22., 11., 23.])

这与data,不同[0,1,2,3,...]。因此,没有任何类型的 reshape 或 view 或 astype 可以将一个转换为另一个。

所以有一个从结构化数组到 (3,4,2) 数组的简单映射,但不是你的源。

In [273]: res.view(float).reshape(3,4,2)
Out[273]: 
array([[[ 0., 12.],
        [ 1., 13.],
        [ 2., 14.],
        [ 3., 15.]],

       [[ 4., 16.],
        [ 5., 17.],
        [ 6., 18.],
        [ 7., 19.]],

       [[ 8., 20.],
        [ 9., 21.],
        [10., 22.],
        [11., 23.]]])
于 2019-01-28T19:50:43.083 回答