Python 将多项式拟合到大型CMIP6 NetCDF xarray的块？

我的目标是计算各种CMIP6模型的piControl运行的“zos”（海平面）数据的时间维度上的二次拟合，并从其他运行中减去。在此特定示例中，数据如下所示：

Out[108]: 
<xarray.Dataset>
Dimensions:    (bnds: 2, ncells: 830305, time: 3, vertices: 16)
Coordinates:
  * time       (time) object 2401-01-16 12:00:00 ... 2401-03-16 12:00:00
    lon        (ncells) float64 dask.array<chunksize=(5000,), meta=np.ndarray>
    lat        (ncells) float64 dask.array<chunksize=(5000,), meta=np.ndarray>
Dimensions without coordinates: bnds, ncells, vertices
Data variables:
    time_bnds  (time, bnds) object dask.array<chunksize=(3, 2), meta=np.ndarray>
    lon_bnds   (ncells, vertices) float64 dask.array<chunksize=(5000, 16), meta=np.ndarray>
    lat_bnds   (ncells, vertices) float64 dask.array<chunksize=(5000, 16), meta=np.ndarray>
    zos        (time, ncells) float32 dask.array<chunksize=(3, 5000), meta=np.ndarray>

---

仅供参考，

polyfit

最近才被合并，尚未发布到xarray版本中（截至2020年5月2日），因此我猜测您为什么会得到

AttributeError

。但是，一旦xarray发布了新版本，

polyfit

将对数据集对象进行操作，并将使用

cftime.datetime

坐标！谢谢与此同时，当我将dask=“parallelized”更改为dask=“allowed”时，我发现有些东西对我来说是可行的。仅供参考，polyfit最近才合并，尚未发布到xarray版本中（截至2020年5月2日），所以我猜这就是为什么会出现

AttributeError

。但是，一旦xarray发布了新版本，

polyfit

将对数据集对象进行操作，并将使用

cftime.datetime

坐标！谢谢同时，当我将dask=“parallelized”更改为dask=“allowed”时，我发现事情对我来说是可行的。