Python xarray 将grib文件保存为Xarray中的NetCDF，crs信息丢失

我正在读取grib文件，计算一个变量的平均值，并将平均值写入netcdf文件。这是我正在使用的代码：

import xarray as xr
import pandas as pd
    
ds = xr.open_dataset('gribfile.grib',engine='cfgrib',chunks={"time": 10})
da = ds['ws'].mean(dim='time')

da.to_netcdf('mean.nc')

代码段设法计算平均值并将其写入文件。但是，写入的地理参考数据不足。Qgis无法从文件中读取正确的crs

到目前为止，我所做的是将原始grib文件的crs添加到netcdf文件中，但两个文件的crs似乎并不相同

---

xarray和/或cfgrib是否做了任何改变文件crs的事情？我是否有办法将crs正确写入netcdf文件？

请共享grib文件，以便社区能够重现您的问题。到目前为止，我想您将向我们展示数据集的属性就足够了。请记住，您存储的是数据数组，而不是数据集。也许您必须将变量强制转换为Dataset

xarray.Dataset（da）

，然后再次从源文件添加CRS。如果这不起作用，我将尝试使用

proj

从源文件读取crs，并在计算过程后将其添加到da中。