这里有人使用过带旋转经纬度坐标的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 ...
...
当我解析数据集时，我会得到以下错误：KeyError:'perspective_point_height'。由于某种原因，坐标中没有海拔或标高，我想解析函数找不到任何高度坐标？如何将垂直标高添加到坐标？
我还使用了一个数据集，其中高度在坐标中指定，但当我尝试进行交叉分割时，我会得到以下结果：ValueError：未处理的投影：rotated_latitude_t经度。有没有一种方法可以让MetPy使用旋转极坐标？
perspective_point_height的错误是因为(我假设(数据将投影声明为地球静止，但它没有这个属性。。。这似乎很奇怪。
对于Unhandled projection error，这只是因为MetPy还没有将该投影添加为它所理解的投影，并且可以转换为CartoPy。
在这两种情况下，似乎有一些方法可以让MetPy更好地处理这一问题，你能在MetPy的Issue Tracker上打开一个bug吗？根据数据中的具体情况，我们可能会为您提供一些解决方案。