这里有人用过旋转纬度经度坐标的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,rlon) float32 ...
snOW_GSP (time,rlon) float32 ...
PRS_GSP (time,rlon) float32 ...
GRAU_GSP (time,rlon) float32 ...
PRG_GSP (time,rlon) float32 ...
PRR_GSP (time,rlon) float32 ...
RELHUM (time,rlon) float32 ...
U (time,srlon) float32 ...
V (time,srlat,rlon) float32 ...
W (time,level1,rlon) float32 ...
...
当我解析然后是数据集时,会出现以下错误: KeyError:“ perspective_point_height”。 由于某种原因,坐标中没有海拔或水平,因此我想解析函数找不到任何高度坐标?如何将垂直标高添加到坐标?
我还使用了在坐标中指定海拔高度的数据集, 但是当我尝试做一个cross_section时,我得到以下信息: ValueError:未处理的投影:rotation_latitude_longitude。 MetPy是否可以使用旋转的极坐标进行工作?