python - 了解 xarray groupby

Question

我正在尝试计算每个组中的成员数量，类似于pandas.DataFrame.groupby.count。但是，它似乎不起作用。这是一个例子:

In [1]: xr_test = xr.DataArray(np.random.rand(6), coords=[[10,10,11,12,12,12]], dims=['dim0'])
        xr_test
Out[1]: <xarray.DataArray (dim0: 6)>
        array([ 0.92908804,  0.15495709,  0.85304435,  0.24039265,  0.3755476 ,
                0.29261274])
        Coordinates:
          * dim0     (dim0) int32 10 10 11 12 12 12

In [2]: xr_test.groupby('dim0').count()
Out[2]: <xarray.DataArray (dim0: 6)>
        array([1, 1, 1, 1, 1, 1])
        Coordinates:
          * dim0     (dim0) int32 10 10 11 12 12 12

但是，我期望这样的输出:

Out[2]: <xarray.DataArray (dim0: 3)>
        array([2, 1, 3])
        Coordinates:
          * dim0     (dim0) int32 10 11 12

发生什么事了？

换句话说:

In [3]: xr_test.to_series().groupby(level=0).count()
Out[3]: dim0
        10    2
        11    1
        12    3
        dtype: int64

Answer 1

这是一个错误! Xarray 当前做出(在本例中是错误的)假设，即与维度对应的坐标具有所有唯一值。这通常是个好主意，但不应该是必需的。如果您创建另一个坐标，这应该可以正常工作，例如， xr_test = xr.DataArray(np.random.rand(6), coords={'aux': ('x', [10,10,11,12,12,12])}, dims=['x'] ) xr_test.groupby('aux').count()