我正在读取一个 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 的东西？有什么方法可以让我在 netcdf 文件中正确写入 crs 吗？