从x数组数据集创建NetCDF时纬度经度错误

Question

我想用xarray创建一个NetCDF文件，并尝试理解此处有关“创建数据集”的文档。

下面是示例中的代码（将ds保存到NetCDF）：

import xarray as xr
import numpy as np
import pandas as pd

temp = 15 + 8 * np.random.randn(2, 2, 3)
precip = 10 * np.random.rand(2, 2, 3)
lon = [[-99.83, -99.32], [-99.79, -99.23]]
lat = [[42.25, 42.21], [42.63, 42.59]]

ds = xr.Dataset({'temperature': (['x', 'y', 'time'],  temp),
                 'precipitation': (['x', 'y', 'time'], precip)},
                coords={'lon': (['x', 'y'], lon),
                        'lat': (['x', 'y'], lat), 
                        'time': pd.date_range('2014-09-06', periods=3),
                        'reference_time': pd.Timestamp('2014-09-05')})

ds.to_netcdf('C:\\Temp\\test.nc')

从上面的例子中，我希望得到一个具有三个维度（x，y，时间）的两个变量（温度，降水）的NetCDF。我希望尺寸在 x 方向上为 2，在 y 方向上为 2，在时间方向上为 3。根据@Bart的测试（评论），NetCDF就是这种情况。因此，当在QGIS 3.4（EPSG 4326）中打开NetCDF“temp”或“precip”变量时，我希望NetCDF数据为：

< li >由2x2网格组成(' x-y '方向) < li >具有3个时间步长(“时间”或“z”方向) < li >位于用于创建数据集的经纬度位置

相反，数据是:

< li >由2x3网格制成 < li >使用两个波段(而不是三个时间步长) < li >位于纬度/经度0/0，每个单元格大1度:

请在这里找到NetCDF的可视化。红色方块用温度变量的蓝色阴影标记六个单元格，左上角点的纬度/经度(0/0)和两个“带”

因此，xarray NetCDF数据格式似乎与QGIS对该格式的解释不同。有人可以编辑该示例，使其在QGIS中的“正确”位置生成一个2x2网格，或者提供/向我指出一个简单的示例，以从xarray数据集创建一个正确的地理参考NetCDF吗？

Answer 1

我不得不交换尺寸的顺序，我必须摆脱嵌套的纬度列表才能正常工作。在此处将 xarray 示例修改为以下代码可以工作：

import xarray as xr
import numpy as np
import pandas as pd

temp = 15 + 8 * np.random.randn(3, 3, 2)  # time, y, x
precip = 10 * np.random.rand(3, 3, 2)  # time, y, x
lat = [100, 101, 102]
lon = [10, 11]

ds = xr.Dataset({'temperature': (['time', 'y', 'x'],  temp),
                 'precipitation': (['time', 'y', 'x'], precip)},
                coords={'x': (['x'], lon),
                        'y': (['y'], lat),
                        'time': pd.date_range('2014-09-06', periods=3)})

ds.to_netcdf('C:\\Temp\\test_xarray.nc')

请注意，我按照问题中的要求将“2x2 网格”更改为“2x3 网格”，以确保创建 NetCDF 时 x-y 的顺序很重要。