python - MetPy 和 COSMO 的旋转纬度经度网格？

Question

这里有人使用带有旋转纬度经度坐标的 MetPy 吗？是否可以？我想要完成的是制作数据的垂直横截面。我用 xr.open_dataset 打开了 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 ...
...

当我解析数据集时，出现以下错误:关键错误:“透视点高度”。由于某种原因，坐标中没有海拔或水平，我想解析函数无法找到任何高度坐标？如何将垂直级别添加到坐标？

我还使用了一个数据集，其中高度在坐标中指定，但是当我尝试做 cross_section 时，我得到以下结果:ValueError:未处理的投影:rotated_latitude_longitude。MetPy 有没有办法可以使用旋转的极坐标？

Answer 1

perspective_point_height 的错误是因为(我假设)数据将投影声明为对地静止，但它没有该属性......这看起来很奇怪。

出现未处理的投影错误，这只是因为 MetPy 尚未将该投影添加为它可以理解并可以转换为 CartoPy 的投影。

在这两种情况下，MetPy 似乎都有办法更好地处理这个问题，你可以在 MetPy's Issue Tracker 上打开一个错误吗？ ？根据数据中具体发生的情况，我们可能会向您指出一些解决方法。