是否有一种方法可以对一个操作进行矢量化,该操作需要几个numpy数组并将它们放入字典列表中?
下面是一个简化的例子。真实的场景可能涉及更多的数组和更多的字典键。
import numpy as np
x = np.arange(10)
y = np.arange(10, 20)
z = np.arange(100, 110)
print [dict(x=x[ii], y=y[ii], z=z[ii]) for ii in xrange(10)]
我可能在xrange
调用中有数千或数十万次迭代。创建x
、y
和z
的所有操作都是矢量化的(我的示例不像上面那么简单)。因此,只剩下1个for循环需要删除,我预计这会带来巨大的速度提升。
我试过使用map
与一个函数来创建字典和各种其他工作。似乎Python的for
循环是慢的部分(像往常一样)。我坚持使用字典是因为已有的API需求。然而,没有字典和记录数组或其他东西的解决方案会很有趣,但最终我认为这不会与现有的API一起工作。
对于您的小示例,我很难获得比列表和字典推导的组合更快的东西
In [105]: timeit [{'x':i, 'y':j, 'z':k} for i,j,k in zip(x,y,z)]
100000 loops, best of 3: 15.5 µs per loop
In [106]: timeit [{'key':{'x':i, 'y':j, 'z':k}} for i,j,k in zip(x,y,z)]
10000 loops, best of 3: 37.3 µs per loop
在分区之前使用数组连接来连接数组的替代方法比较慢。
In [108]: timeit [{'x':x_, 'y':y_, 'z':z_} for x_, y_, z_ in np.column_stack((x,y,z))]
....
10000 loops, best of 3: 58.2 µs per loop
=======================
使用recfunctions
:
In [109]: from numpy.lib import recfunctions
In [112]: M=recfunctions.merge_arrays((x,y,z))
In [113]: M.dtype.names=['x','y','z']
In [114]: M
Out[114]:
array([(0, 10, 100), (1, 11, 101), (2, 12, 102), (3, 13, 103),
(4, 14, 104), (5, 15, 105), (6, 16, 106), (7, 17, 107),
(8, 18, 108), (9, 19, 109)],
dtype=[('x', '<i4'), ('y', '<i4'), ('z', '<i4')])
In [115]: M['x']
Out[115]: array([0, 1, 2, 3, 4, 5, 6, 7, 8, 9])
时间要慢得多,但是如果你想一次访问所有的x
值,这比从所有的字典中获取它们要好得多。
np.rec.fromarrays((x,y,z),names=['x','y','z'])
生成给定名称的重数组。差不多一样的速度
还可以构造一个dtype和形状正确的空数组,并将数组复制到其中。这可能和这个merge
一样快,但描述起来更复杂。
我建议优化数据结构的使用/访问,而不是建设速度。一般只构造一次,然后多次使用。
============
In [125]: dt=np.dtype([('x',x.dtype),('y',y.dtype),('z',z.dtype)])
In [126]: xyz=np.zeros(x.shape,dtype=dt)
In [127]: xyz['x']=x; xyz['y']=y; xyz['z']=z
# or for n,d in zip(xyz.dtype.names, (x,y,z)): xyz[n] = d
In [128]: xyz
Out[128]:
array([(0, 10, 100), (1, 11, 101), (2, 12, 102), (3, 13, 103),
(4, 14, 104), (5, 15, 105), (6, 16, 106), (7, 17, 107),
(8, 18, 108), (9, 19, 109)],
dtype=[('x', '<i4'), ('y', '<i4'), ('z', '<i4')])
这是一个(Num)?神谕的方式:
In [18]: names = np.array(['x', 'y', 'z'])
In [38]: map(dict, np.dstack((np.repeat(names[None, :], 10, axis=0), np.column_stack((x, y, z)))))
Out[38]:
[{'x': '0', 'y': '10', 'z': '100'},
{'x': '1', 'y': '11', 'z': '101'},
{'x': '2', 'y': '12', 'z': '102'},
{'x': '3', 'y': '13', 'z': '103'},
{'x': '4', 'y': '14', 'z': '104'},
{'x': '5', 'y': '15', 'z': '105'},
{'x': '6', 'y': '16', 'z': '106'},
{'x': '7', 'y': '17', 'z': '107'},
{'x': '8', 'y': '18', 'z': '108'},
{'x': '9', 'y': '19', 'z': '109'}]
另外,请注意,如果您不需要一次使用所有的字典,您可以简单地创建一个生成器,并根据需要访问每个条目。
(dict(x=x[ii], y=y[ii], z=z[ii]) for ii in xrange(10))
如果你想要一个嵌套的字典,我建议使用列表推导式:
In [88]: inner = np.dstack((np.repeat(names[None, :], 10, axis=0), np.column_stack((x, y))))
In [89]: [{'connection': d} for d in map(dict, inner)]
Out[89]:
[{'connection': {'x': '0', 'y': '10'}},
{'connection': {'x': '1', 'y': '11'}},
{'connection': {'x': '2', 'y': '12'}},
{'connection': {'x': '3', 'y': '13'}},
{'connection': {'x': '4', 'y': '14'}},
{'connection': {'x': '5', 'y': '15'}},
{'connection': {'x': '6', 'y': '16'}},
{'connection': {'x': '7', 'y': '17'}},
{'connection': {'x': '8', 'y': '18'}},
{'connection': {'x': '9', 'y': '19'}}]
这是使用NumPy
和Pandas
混合的方法-
# Stack into columns & create a pandas dataframe with appropriate col names
a = np.column_stack((x.ravel(),y.ravel(),z.ravel()))
df = pd.DataFrame(a,columns=[['x','y','z']])
# Convert to list of dicts
out = df.T.to_dict().values()
示例运行-
In [52]: x
Out[52]: array([0, 1, 2, 3, 4, 5, 6, 7, 8, 9])
In [53]: y
Out[53]: array([10, 11, 12, 13, 14, 15, 16, 17, 18, 19])
In [54]: z
Out[54]: array([100, 101, 102, 103, 104, 105, 106, 107, 108, 109])
In [55]: out
Out[55]:
[{'x': 0, 'y': 10, 'z': 100},
{'x': 1, 'y': 11, 'z': 101},
{'x': 2, 'y': 12, 'z': 102},
{'x': 3, 'y': 13, 'z': 103},
{'x': 4, 'y': 14, 'z': 104},
{'x': 5, 'y': 15, 'z': 105},
{'x': 6, 'y': 16, 'z': 106},
{'x': 7, 'y': 17, 'z': 107},
{'x': 8, 'y': 18, 'z': 108},
{'x': 9, 'y': 19, 'z': 109}]