按非关联坐标剪裁NetCDF文件

Question

查看此NetCDF EURO-CORDEX文件：

<xarray.Dataset>
Dimensions:       (bnds: 2, rlat: 412, rlon: 424, time: 1800)
Coordinates:
    lat           (rlat, rlon) float64 ...
    lon           (rlat, rlon) float64 ...
  * rlat          (rlat) float64 -23.38 -23.26 -23.16 ... 21.61 21.73 21.83
  * rlon          (rlon) float64 -28.38 -28.26 -28.16 ... 17.93 18.05 18.16
  * time          (time) object 2011-01-01 12:00:00 ... 2015-12-30 12:00:00
Dimensions without coordinates: bnds
Data variables:
    pr            (time, rlat, rlon) float32 ...
    rotated_pole  |S1 ...
    time_bnds     (time, bnds) object ...

当我尝试使用xarray（在Python脚本中）或nco（在命令行中）通过lat和lon坐标剪裁它时，它可能不起作用，因为lat和lon不是维度。我读到这些与任何变量都没有关联。我如何关联它们或使它们成为维度？坐标rlat和rlon旁边的星号是否表示我遗漏了什么

注意：我通过xarray.Dataset.where实现剪裁，但是我对操纵尺寸/坐标非常感兴趣，这样通常的剪裁操作就可以工作了。谢谢

Answer 1

您应该能够使用CDO对此进行剪辑：

sellonlatbox,lonmin,lonmax,latmin,latmax infile outfile

CDO将计算出原始数据集中包含指定框的最小可能矩形。这似乎是一个旋转的极轴栅格，因此不会在框外获得一些栅格单元

NCO也可以工作，但需要使用rlon和rlat进行裁剪。因此，您需要选择适当的rlon和rlat，以获得所需的lon和lat