我询问了dtype,因为你的例子令人费解。
我可以创建一个包含 3 个元素 (1d) 和 3 个字段的结构化数组:
In [1]: A = np.ones((3,), dtype='i,i,i')
In [2]: A
Out[2]:
array([(1, 1, 1), (1, 1, 1), (1, 1, 1)],
dtype=[('f0', '<i4'), ('f1', '<i4'), ('f2', '<i4')])
我可以通过名称访问一个字段(添加括号不会改变事情)
In [3]: A['f0'].shape
Out[3]: (3,)
但如果我访问 2 个字段,我仍然会得到一个一维数组
In [4]: A[['f0','f1']].shape
Out[4]: (3,)
In [5]: A[['f0','f1']]
Out[5]:
array([(1, 1), (1, 1), (1, 1)],
dtype=[('f0', '<i4'), ('f1', '<i4')])
实际上,如果我查看值,那些额外的括号确实很重要
In [22]: A['f0']
Out[22]: array([1, 1, 1], dtype=int32)
In [23]: A[['f0']]
Out[23]:
array([(1,), (1,), (1,)],
dtype=[('f0', '<i4')])
如果数组是简单的 2d 数组,我仍然没有得到你的形状
In [24]: A=np.ones((3,3),int)
In [25]: A[0].shape
Out[25]: (3,)
In [26]: A[[0]].shape
Out[26]: (1, 3)
In [27]: A[[0,1]].shape
Out[27]: (2, 3)
但是关于确保数组是2d的问题,不管索引返回1d还是2,你的函数基本没问题
def reshape_to_vect(ar):
if len(ar.shape) == 1:
return ar.reshape(ar.shape[0],1)
return ar
您可以测试ar.ndim 而不是len(ar.shape)。但无论哪种方式,它的成本都不高——也就是说,执行时间很短——没有大的数组操作。 reshape 不会复制数据(除非你的步幅很奇怪),所以它只是使用共享数据指针创建一个新数组对象的成本。
查看np.atleast_2d的代码;它测试 0d 和 1d。在第一种情况下,它返回result = ary[newaxis,:]。它首先添加额外的轴,更自然的numpy 添加轴的位置。你在最后添加它。
ar.reshape(ar.shape[0],-1) 是绕过if 测试的巧妙方法。在小时间测试中它更快,但我们谈论的是微秒,即函数调用层的效果。
np.column_stack 是另一个在需要时创建列数组的函数。它使用:
if arr.ndim < 2:
arr = array(arr, copy=False, subok=True, ndmin=2).T