Python MetPy和COSMO旋转的经纬度网格?

Python MetPy和COSMO旋转的经纬度网格?,python,metpy,Python,Metpy,这里有人用过旋转经纬度坐标的MetPy吗?可能吗?我试图完成的是制作数据的垂直横截面。我使用xr.open\u数据集打开了我的netCDF文件: Dimensions: (bnds: 2, level: 80, level1: 81, rlat: 127, rlon: 162, srlat: 127, srlon: 161, time: 48) Coordinates: * time (time) datetime64[ns] 2019-03-07T08:00:

这里有人用过旋转经纬度坐标的MetPy吗?可能吗?我试图完成的是制作数据的垂直横截面。我使用xr.open\u数据集打开了我的netCDF文件:

Dimensions:       (bnds: 2, level: 80, level1: 81, rlat: 127, rlon: 162, srlat: 127, srlon: 161, time: 48)
Coordinates:
  * time          (time) datetime64[ns] 2019-03-07T08:00:00 ... 2019-03-07T19:45:00
    lon           (rlat, rlon) float32 ...
    lat           (rlat, rlon) float32 ...
    slonu         (rlat, srlon) float32 ...
    slatu         (rlat, srlon) float32 ...
    slonv         (srlat, rlon) float32 ...
    slatv         (srlat, rlon) float32 ...
  * rlon          (rlon) float64 -1.385 -1.375 -1.365 ... 0.205 0.215 0.225
  * rlat          (rlat) float64 -0.995 -0.985 -0.975 ... 0.245 0.255 0.265
  * srlon         (srlon) float64 -1.38 -1.37 -1.36 -1.35 ... 0.19 0.2 0.21 0.22
  * srlat         (srlat) float64 -1.0 -0.99 -0.98 -0.97 ... 0.23 0.24 0.25 0.26
Dimensions without coordinates: bnds, level, level1
Data variables:
    time_bnds     (time, bnds) datetime64[ns] ...
    rotated_pole  int32 ...
    vcoord        (level1) float32 ...
    TOT_PREC      (time, rlat, rlon) float32 ...
    TOT_PR        (time, rlat, rlon) float32 ...
    SNOW_GSP      (time, rlat, rlon) float32 ...
    PRS_GSP       (time, rlat, rlon) float32 ...
    GRAU_GSP      (time, rlat, rlon) float32 ...
    PRG_GSP       (time, rlat, rlon) float32 ...
    PRR_GSP       (time, rlat, rlon) float32 ...
    RELHUM        (time, level, rlat, rlon) float32 ...
    U             (time, level, rlat, srlon) float32 ...
    V             (time, level, srlat, rlon) float32 ...
    W             (time, level1, rlat, rlon) float32 ...
...
当我分析then dataset时,我得到以下错误: KeyError:“透视图\点\高度”。 由于某种原因,坐标中没有高度或标高,我认为解析函数无法找到任何高度坐标?如何将垂直标高添加到坐标

我还使用了一个数据集,其中高度在坐标中指定, 但是当我尝试做一个横截面时,我得到了以下结果: ValueError:未处理的投影:旋转的经纬度。
有没有一种方法可以让MetPy处理旋转的极坐标?

透视点高度的错误是因为(我假设)数据声明投影是地球静止的,但它没有这个属性…这似乎很奇怪

对于
未处理的投影错误
,这只是因为MetPy没有将该投影添加为它能够理解的投影,并且可以转换为CartoPy


在这两种情况下,MetPy似乎都有办法更好地处理这个问题,你能在上面打开一个bug吗?根据数据中到底发生了什么,我们可能会为您指出一些解决方法。

透视点高度的错误是因为(我假设)数据声明投影是地球静止的,但它没有该属性…这似乎很奇怪

对于
未处理的投影错误
,这只是因为MetPy没有将该投影添加为它能够理解的投影,并且可以转换为CartoPy


在这两种情况下,MetPy似乎都有办法更好地处理这个问题,你能在上面打开一个bug吗?取决于数据中到底发生了什么,我们可能会向您指出一些解决方法。

我认为我们需要更多关于您试图通过MetPy实现的目标的信息,以了解问题所在。我认为我们需要更多关于您试图通过MetPy实现的目标的信息,以了解问题所在。