问题标签 [metpy]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
265 浏览

python - 获取xarray中点到点对角线上的所有网格值

我希望从xarray一个点和另一个点之间的数据集中选择所有数据。这当然很容易,sel例如使用基于纬度/经度网格的数据集从全球数据集到大西洋海岸线的子集我可以做到ds_geo_subset = ds.sel( lat=slice(42, 20), lon=slice(-85,-65)),但我想选择墨西哥湾流前沿的所有数据,基本上是从到点(- 78 lon, 32 lat) 到 (-72 lon, 36 lat) 但我不知道该怎么做。最好通过给出一个完整的点数据集并从每个点获取数据来完成,还是可以将这两个端点提供给 xarray 并返回插值数据?

0 投票
1 回答
113 浏览

metpy - 转换metpy单位GFS precipitation_rate

我按照此处的示例(https://unidata.github.io/python-gallery/examples/Precipitation_Map.html)但是在尝试访问单位时出现维度错误。我假设这与pint它解析单位的方式和降水量是一个速率有关。任何帮助,将不胜感激

0 投票
1 回答
203 浏览

python - 使用简单数据时,metpy interpolate_to_grid 不适合我

我在调整简单数据集以从点生成温度插值图时遇到问题

像使用 colab 笔记本一样安装依赖项

这给出了一个 lat, lon 的数组

这显示数据点正常

在此处输入图像描述

这会产生错误 TypeError: Input z must be at least a (2, 2) 形状的数组,但具有形状 (1, 1)

作为 python 的新手,metpy 一直在努力解决这个问题,感谢您的帮助。

**已编辑:谢谢@dopplerShift,就是这样,这是修复和结果

在此处输入图像描述

0 投票
1 回答
194 浏览

python - 插值后将特征叠加到 matplotlib 图的问题

我不得不承认我仍然无法理解绘图的正确设置和关系以及它与 matplotlib 的部分,仍然混淆 fig 与 plt 与 ax 之间的关系,所以我只是反复试验,文档有时更多让我困惑。:-(

我正在绘制来自 json 的天气值并获得积分。我可以用下面的代码绘制如下图

我可以绘制地图 在此处输入图像描述

然后我使用metpy和这段代码生成插值

我根据需要得到了点的插值,但无法显示特征,该怎么做?如果我取消注释 ax.extent 我看到的只是一个空白的数字。如果我取消注释 ax.features 插值显示为下图而不是地图。

感谢您的帮助和指导。

在此处输入图像描述

0 投票
1 回答
83 浏览

python-xarray - 在 xarray.Dataset 上调用 .metpy.parse_cf() 会给出未定义的“ctables”

在 xarray.Dataset 上调用 .metpy.parse_cf() 会给出未定义的“ctables”

0 投票
3 回答
641 浏览

python - 将网格数据插值到地理位置

我是 MetPy 的忠实粉丝,并查看了它们的插值函数 ( https://unidata.github.io/MetPy/latest/api/generated/metpy.interpolate.html ),但找不到我想要的东西。

我正在寻找一种将网格 2D(经度和纬度)或 3D(经度、纬度和垂直水平)气候数据字段插入特定地理位置(纬度/经度)的功能。

该函数将采用 5 个参数:一个 2D/3D 数据变量和相关的纬度和经度变量,以及两个所需的纬度和经度坐标值。返回的是单个值(对于 2D 字段)或垂直配置文件(对于 3D 字段)。

我基本上是在寻找与旧 Basemap 函数 bm.interp() 等效的函数。Cartopy 没有等价物。CDO(气候数据操作员)操作员 'remapbil,lon=/lat=' 做同样的事情,但直接从命令行处理 netCDF 文件,我正在寻找 Python 解决方案。

我认为这样的功能将是对 MetPy 库的有用补充,因为它允许将网格数据(例如模型或卫星数据)与来自气象站或无线电探空仪剖面(此处仅视为垂直剖面)等点观测数据进行比较。

你能为我指出正确的方向吗?

0 投票
1 回答
126 浏览

python - 使用 Google Colab 时出现 Cartopy (0.14.2) / MetPy (0.12.2) parse(cf) 错误

我有一个学生在使用旧版本的 MetPy 和 Cartopy 时遇到的问题类似于其他人报告的问题。

问题是metpy 下的parse_cf() 命令。它发生在 Google Colaboratory 下,我对此几乎没有经验。在其他平台(anaconda、miniconda、homebrew 等)上一切正常

我之前看到的反馈是指从 0.18.* 回滚 cartopy 库的版本。我也遵循了之前发布的关于 Shapely 的指南,但我们似乎无法在 Google CoLab 中更改 cartopy 或 metpy 版本的版本。

下面是环境和python库的版本信息,后面是代码和错误。

我在这里不知所措,并且尽我所能为我的学生解决了问题。我们俩都将不胜感激。

系统和包信息

要复制的代码

错误信息

0 投票
1 回答
85 浏览

metpy - 示例中的 mpcalc.advection dim_order 和缺失值

在metpy网页上,他们为850hPa温度平流提供的示例如下:

但根据文档,其中指出:

在我看来,示例中的 dim_order 应该是 xy 而不是 yx?当我们给 u 然后是 v,然后是 dx 然后是 dy,而不是反过来呢?

我还发现,即使我的输入温度场在某些网格点有缺失值,metpy 仍然以某种方式返回这些点处的平流值,而不是缺失值。这怎么可能?

0 投票
1 回答
87 浏览

python - 使用 MetPy 生成 Countour 图的问题

我正在尝试使用 MetPy 和 ContourPlot 生成等压线。我使用我自己的 netCDF 文件中的压力数据。背景热图生成正确,但对于轮廓,似乎存在一种插值错误。有谁知道如何解决水平白线的问题?

这是我用于绘图生成的代码:

等高线 MetPy

0 投票
1 回答
201 浏览

python - 在此降水图上显示两个水平颜色条,导致大的空白区域

我正在使用 Metpy 制作一张漂亮的降水图。我认为这更像是一个 matplotlib 问题,但我对如何在这个图中添加两个水平颜色条感到困惑。总体思路是在这张地图上显示降雪率和降雨率,这是我通过在一个图形上使用两个等高线填充来完成的。我想显示两个颜色条,因为我有两个单独的降雪率和降雨率。理想情况下,我希望两个颜色条水平并排。我得到的最接近的方法是将它们堆叠在一起,并在它们之间留出很大的空间。这是我草率的代码。

绘图示例

我已经尝试过 subplots 并且我已经尝试过 gridSpec 但我一直遇到不同的障碍,当我尝试这些不同的方法时会遇到不同的错误。主要是我认为这与我对 Metpy、Cartopy 和 Matplotlib 的入门级知识有关。欣赏此问题的任何见解和可能的解决方案。谢谢!