python中多维数组的最佳实践

Question

假设我有一个我希望保存在 python 中的对象集合，比如说，一个数字列表：[0.12, 0.85, 0.11, 0.12], [0.23, 0.52, 0.10, 0.19] 等。进一步假设这些对象由 3 个属性索引，例如“origin”、“destination”和“month”。我希望将这些对象存储在一个可以轻松切片的类数组对象中，最好使用数字索引或名称。

所以，即，

obj[2,1,7] # might return: [0.23, 0.52, 0.10, 0.19]

或者，

obj['chicago','new york','jan'] # might return: [0.12, 0.85, 0.11, 0.12]

还有，

obj[:,'new york','jan'] # would return data with first index = any.

我正在寻找在 python 中实现这一目标的最佳实践。我确实找到了this 帖子，这似乎很合适，但它似乎需要一些开销，并且几乎没有讨论替代方案。我还发现了一个名为xarray 的包，尽管它似乎不那么受欢迎。我正在从 R 转换，我将在其中执行 array() 函数，它为任何类似向量的结构添加多维索引。

Answer 1

经过一番摸索，看来 xarray 适合我的需要。不幸的是，由于我缺乏经验，我无法谈论与其他软件包的兼容性或性能。

import numpy as np
import xarray as xr
cityOrig = ['chicago','new york', 'boston']
cityDest = ['chicago','new york', 'boston']
month = ['jan','feb','mar','apr']
data = np.random.rand(4,3,3,4)

myArray = xr.DataArray(data,
                       dims=['dat','orig','dest','month'],
                       coords = {'orig':cityOrig,'dest':cityDest,'month':month})

print(myArray[:,1,2,1].data)
[0.64  0.605 0.445 0.059]
print(myArray.loc[:,'chicago','new york','jan'].data)
[0.64  0.605 0.445 0.059]