MetPy: информация о проекции в наборе данных xarray для сечений?
Входные данные для metpy.cross_section() - это набор данных xarray или DataArray, который "должен иметь прикрепленную информацию о проекции".
Чтобы использовать перекрестный раздел MetPy, мне нужно написать свой собственный набор данных xarray, поскольку мои данные NWP с метаданными не могут быть автоматически прочитаны в набор данных xarray. Однако мне неясно, какие метаданные требуются для того, чтобы MetPy мог понять мою проекцию.
Я понимаю, что именно metpy.parse_cf() предоставляет MetPy информацию о проекции, но я не могу найти никакой информации о metpy.parse_cf(), то есть какие переменные являются необходимыми. На данный момент metpy.parse_cf() для меня черный ящик.
1 ответ
"Прикрепленная информация о проекции", необходимая для функции поперечного сечения, имеет crs координировать и axis координаты атрибутов, которые MetPy добавляет в набор данных при использовании ds.metpy.parse_cf() (как описано в руководстве xarray с MetPy).
MetPy создает crs координировать и axis координировать атрибуты с использованием метаданных в наборе данных в соответствии с конвенциями CF. Если ваш считыватель данных GRIB-to-xarray содержит эту информацию для данных на конформной сетке Ламберта, то crs координата должна быть создана правильно (вы можете проверить это, убедившись, что она гласит "Проекция: lambert_conformal_conic"). Однако, если ваш читатель набора данных не включает эти метаданные, возможно, вам придется добавить его самостоятельно.