使用 xarray to_netcdf 写入磁盘会丢失坐标变量

Question

我正在尝试使用 xarray to_netcdf 保存数据集，但是当我将数据集写入磁盘然后读回时，坐标已更改。这只是一个偶然的问题，我无法重现它。

编辑：我已经发现问题与附加到变量的编码有关，现在可以用一个虚拟示例重现该问题：

lon = [[-99.83, -99.32], [-99.79, -99.23]]
lat = [[42.25, 42.21], [42.63, 42.59]]
time = pd.date_range("2014-09-06", periods=3)

ds = xr.Dataset(
    data_vars=dict(
        temperature=(["x", "y", "time"], temperature),),
    coords=dict(
        lon=(["x", "y"], lon),
        lat=(["x", "y"], lat),
        time=time,),
    attrs=dict(description="Weather related data."),)


ds.temperature.encoding={'zlib': False,
 'shuffle': False,
 'complevel': 0,
 'fletcher32': False,
 'contiguous': True,
 'chunksizes': None,
 'source': 'fileloc/file.nc',
 'original_shape': (3, 2, 2),
 'coordinates': 'time, lat, lon'}

print(ds) 给出了我期望的坐标，但是 ds.to_netcdf('dummy.nc')、ds2=xr.open_dataset('dummy.nc') 丢失 'lon' 作为坐标并设置它作为数据变量。

有谁知道这种编码中的什么导致“lon”坐标发生变化？我看不到任何与“lon”具体相关的内容。

我已经更新到最新版本的 xarray 并尝试了decode_coords='all'，但仍然看到这个问题。