这似乎是一个已知问题,并且在该问题的线程和最近合并的 PR中进行了相当多的讨论。
.encoding
基本上,数据集在属性中携带原始分块。因此,当您调用第二次写入操作时,ds[var].encoding['chunks']
(如果存在)中定义的块将用于写入var
zarr。
根据 GH issue 中的对话,目前最好的解决方案是手动删除有问题的变量的块编码:
for var in ds:
del ds[var].encoding['chunks']
但是,应该注意的是,这似乎是一个不断变化的情况,最好检查进度以适应最终解决方案。
这是一个展示问题和解决方案的小示例:
import xarray as xr
# load data and write to initial chunking
x = xr.tutorial.load_dataset("air_temperature")
x.chunk({"time":500, "lat":-1, "lon":-1}).to_zarr("zarr1.zarr")
# display initial chunking
xr.open_zarr("zarr1.zarr/").air
# rechunk
y = xr.open_zarr("zarr1.zarr/").chunk({"time": -1})
# display
y.air
#write w/o modifying .encoding
y.to_zarr("zarr2.zarr")
# display
xr.open_zarr("zarr2.zarr/").air
# delete encoding and store
del y.air.encoding['chunks']
y.to_zarr("zarr3.zarr")
# display
xr.open_zarr("zarr3.zarr/").air