采用单一维度的numpy视图

时间:2015-11-19 23:54:50

标签: python numpy

设置

我有一个Fortran有序的numpy数组,其中包含单个第一维

In [1]: import numpy as np

In [2]: x = np.array([[1, 2, 3]], order='F')

In [3]: x.strides
Out[3]: (8, 8)

我使用较小尺寸的dtype

来查看此数组
In [4]: y = x.view('i2')

In [5]: y
Out[5]: array([[1, 0, 0, 0, 2, 0, 0, 0, 3, 0, 0, 0]], dtype=int16)

In [6]: y.strides
Out[6]: (8, 2)

问题

首先,我应该关注这一进步的结果吗?下一行有八个以上的字节。

其次,我真的希望这个数组的形状为(4, 3),而不是(12, 3)。通常使用Fortran布局的数组.view将扩展第一维,而不是第二维。

第三,如何将上面的数组(或由类似进程生成的任何数组)转换为(4, 3)数组,这样,在这种情况下,四行中的三行全部为零。我错过了聪明的换位技巧?

版本

In [3]: sys.version_info
Out[3]: sys.version_info(major=3, minor=4, micro=3, releaselevel='final', serial=0)

In [4]: numpy.__version__
Out[4]: '1.10.1'

1 个答案:

答案 0 :(得分:1)

奇怪的是,当我重复你的步骤时,我马上得到(4,3):

In [637]: x=np.array([[1,2,3]],'int64',order='F')
In [638]: x
Out[638]: array([[1, 2, 3]], dtype=int64)
In [639]: y=x.view('i2')
In [640]: y
Out[640]: 
array([[1, 2, 3],
       [0, 0, 0],
       [0, 0, 0],
       [0, 0, 0]], dtype=int16)
In [641]: x.shape
Out[641]: (1, 3)
In [642]: y.shape
Out[642]: (4, 3)
In [643]: x.strides
Out[643]: (8, 8)
In [644]: y.strides
Out[644]: (2, 8)

我在Py3(老式32b机器)上的numpy 1.8.2

订单C

In [655]: x=np.array([[1,2,3]],'int64')
In [656]: y=x.view('i2')
In [657]: y
Out[657]: array([[1, 0, 0, 0, 2, 0, 0, 0, 3, 0, 0, 0]], dtype=int16)
In [658]: y.strides
Out[658]: (24, 2)
In [659]: x.strides
Out[659]: (24, 8)