如何在python xarray中存储字典中包含的多维数组？

Question

我有一本字典，其中包含各种大小的numpy服装。 所有数组都有一个共同的轴长（时间），我希望将其存储在其中。

例如：

arr1 = np.random.rand(239, 1)
arr2 = np.random.rand(239, 6)
arr3 = np.random.rand(239, 3, 7)
time = np.random.rand(239, 1)
d = {'A': arr1, 'B': arr2, 'C': arr3, 'time': time}

我需要能够轻松地索引和操作数据，所以我的第一个倾向是使用pandas.Panel存储数据，但是由于尺寸不一致，我一直没有成功。

xarray.Dataset是在此处存储我的数据的正确方法吗？如果可以的话，如何最好地实现呢？

Answer 1

这是使用标准熊猫方法的一种非常简单的方法。 为了简便起见，我将数据缩小并放入数据框，但是概念是相同的。

dr=pd.date_range('1-1-2017', periods=4, freq='d')
df1=pd.DataFrame( np.random.randn(4),   columns=['x'],     index=dr)
df2=pd.DataFrame( np.random.randn(4,2), columns=['y','z'], index=dr)

所以df1和df2看起来像这样：

                   x
2017-01-01 -0.705449
2017-01-02 -0.597631
2017-01-03 -0.844197
2017-01-04 -1.063895
                   y         z
2017-01-01 -0.288822 -0.343934
2017-01-02  1.072678  1.776767
2017-01-03 -0.606593  0.192280
2017-01-04  0.019401  2.007770

像这样重新配置：

df = df1.stack().append(df2.stack()).sort_index()

---

2017-01-01  x   -0.705449
            y   -0.288822
            z   -0.343934
2017-01-02  x   -0.597631
            y    1.072678
            z    1.776767
2017-01-03  x   -0.844197
            y   -0.606593
            z    0.192280
2017-01-04  x   -1.063895
            y    0.019401
            z    2.007770

您甚至可以使用以下命令从此处转换为xarray ：

df.to_xarray()

快速注意事项：

• 不建议使用Panel，而是使用xarrays或multi-indexs。 我采用了上面的多索引方法，但是xarrays是另一个不错的选择
• 有关出色的数据组织理论，请参见Hadley Wickam对“整洁”数据的解释，您可以在此处找到。