MetPy и повернутая сетка широты и долготы от COSMO? - PullRequest
Новая
       29

MetPy и повернутая сетка широты и долготы от COSMO?

0 голосов
/ 07 ноября 2019

Кто-нибудь здесь использовал MetPy с координатами повернутой широты и долготы? Является ли это возможным? Я пытаюсь сделать вертикальные сечения моих данных. Я открыл свой файл netCDF с помощью xr.open_dataset: 

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:'спектива_point_height '. По какой-то причине в координатах нет высоты или уровня, и я представляю, что функция синтаксического анализа не может найти координату высоты? Как я могу добавить вертикальные уровни к координатам?

Я также использовал набор данных, в котором высота указана в координатах, но затем, когда я пытаюсь выполнить перекрестное сечение, я получаю следующее: ValueError: Unhandledпроекция: rotated_latitude_longitude. Есть ли способ, чтобы MetPy мог работать с повернутыми полюсными координатами?

1 Ответ

0 голосов
/ 08 ноября 2019

Ошибка с perspective_point_height заключается в том, что (я предполагаю) данные объявляют проекцию как геостационарную, но у нее нет этого атрибута ... что кажется странным.

С Unhandled projection error это потому, что MetPy не добавил эту проекцию, поскольку он понимает ее и может конвертировать в CartoPy.

В обоих случаях кажется, что есть способы, которыми MetPy может лучше справиться с этим, вы можете открытьошибка на Отслеживание проблем MetPy ? В зависимости от того, что именно происходит в данных, вероятно, будут некоторые обходные пути, на которые мы можем вам указать.

...