问题标签 [metpy]

当我打算从数组数据集计算温度平流时，出现了一个错误：

TypeError：无效的索引器数组，没有整数 dtype：array(None, dtype=object)

示例：我读取 WRF 输出，xarray 数据数组是：

xarray.DataArray'ua' (时间: 365south_north: 156west_east: 108)

坐标：

XLONG (south_north, west_east) float32 ...

XLAT (south_north, west_east) float32 ...

interp_level () float64 850.0

时间（Time）对象 1997-01-01 00:00:00 ...

我跑

错误信息是：

谢谢你

遵循地理彩色卫星图像教程时，我在 MetPy 中遇到错误。具体来说，标题为“使用 Cartopy 地球同步投影绘图”的部分。这种损坏发生在大约两周前，功能尚未恢复。考虑以下代码：

此块在功能上类似于 MetPy geocolor 卫星教程中提供的块。直到最近它工作得很好。现在出现以下错误：

这里似乎有什么问题？有可用的解决方法吗？

我正在使用来自 NCAR RDA 的 ERA5 再分析数据。我正在使用metpy dewpoint_from_relative_humidity() 来计算露点，然后将其转换为DataArray 对象并分配与我拥有ERA5 数据的Dataset 对象匹配的坐标和尺寸。我想在metpy 中使用precipitable_water() 来计算可降水量，但是我正在努力找出执行此计算的最有效方法，因为它旨在处理数组（探测）而不是数据网格。我已经查看了 xr.apply_ufunc() 但我不确定这是否是该 xarray 功能的适当应用程序，或者坦率地说我将如何执行它。我可以编写嵌套的 for 循环/while 循环，但我认为这不是最好的方法。有没有人有什么建议？我尝试了下面的行，但它'

我正在尝试使用一个nc文件创建降水图，类似于我在此处找到的 NWS 示例。

不过，就我而言，我的降水数据已经在mm. 如何保持相同的单位？我确实阅读了以下内容，

创建一个numpy.ma.MaskedArray附加单位。这是一个薄包装器numpy.ma.masked_array()，可确保单元正确附加到结果（否则单元会默默丢失）。单位取自 data_units 参数，或者如果是None，则使用数据上的单位。

我为我的文件遵循了给定的参数（masked_array(data, data_units=None, **kwargs)），但是

1. kwargs 没有定义，并且
2. 当我不包括 kwargs 时，我得到

“AttributeError：‘MaskedArray’对象没有属性‘units’”。

我是初学者，所以请温柔。感谢您的指导！这是我的代码...

使用“contourf”方法时 ArtistAnimation 是否有效？我正在尝试构建模型预测的 MP4 动画。

我正在使用一种类似于 MetPy Monday 视频系列中显示的方法，在该视频系列中，您将绘图附加到一个空的艺术家数组（用于雷达图）。雷达/卫星图似乎成功地构建了动画，因为它们使用了“pcolormesh”，但对于我使用“contourf”的模型图来说，情况并非如此。

在 JupyterLab 中执行此代码时，

我在堆栈跟踪的末尾看到了这一点：

AttributeError：“QuadContourSet”对象没有属性“set_visible”

stackoverflow 上的一个旧线程尝试解决此问题： https ://github.com/matplotlib/matplotlib/issues/6139

实施这些修复时，set_visible 属性错误不再出现，而是出现：

AttributeError：“QuadContourSet”对象没有属性“get_figure”

不知道从这里去哪里。如果有好消息，我会在地图上一次看到每个时间步的所有数据（因此数据检索是可以的），但是 ArtistAnimation 以及它如何与 contourf 方法一起工作。任何指导或有用的提示都会很棒！

我正在尝试使用metpy.calc 函数，但是每次我尝试使用metpy.units 分配单位时，它都不会运行。如果我允许它会坐在那里继续运行几个小时。过去，metpy 函数在使用相同数据集时对我来说工作得非常好，但由于某种原因将不再工作。我已经尝试更新metpy和spyder，并且还使用了多台笔记本电脑以确保它不仅仅发生在一台笔记本电脑上。我还尝试了 metapy.units 两种不同的方式。它使用单位（'kg/kg'）运行，但实际上不会将单位分配给数组。关于如何让它发挥作用的任何想法？

这是我的代码：

长期程序员，但对 python 和 metpy 不熟悉。我正在尝试通过遵循 unidata 网站和 youtube (metpy mondays) 上的示例来学习

我正在尝试构建一个简单的 skew-t 课程并遇到问题。当我尝试像这样使用导入时：

我收到此错误：

我之前对这个导入没有任何问题，所以想知道这里发生了什么？我的环境中的包之间是否存在一些不一致，或者我缺少什么？我已经尝试更新所有内容以确保我没有使用任何旧的东西。

作为参考，我将 conda 版本 4.8.3 与 python 3.7.7 和metpy 0.12.1 一起使用。如果我需要提供任何其他可能有用的信息，请告诉我。

在此先感谢您的帮助/见解！

新手用户在这里。

我在尝试使用 parcel_profile 时遇到了一个问题，它不断弹出一个错误，指出该变量没有与之关联的单位：

回溯（最近一次通话最后）：文件“Advanced_Sounding_3Dnetcdf2.py”，第 165 行，在 prof = mpcalc.parcel_profile(p, T[0], Td[0]).to('degC')

有问题的代码片段：

当我打印“p”时，它如下所示：

数据读入如下：

我这里有另一个编码错误还是有其他问题？

谢谢。

我刚刚注意到，当我尝试在 NCAR 的 Cheyenne 超级计算机上导入metpy.plots 时，它在使用metpy 0.10.0（使用Cartopy 0.17.0）时加载并正常工作，但使用metpy 0.12.0 或0.12.1 时出现错误（使用 Cartopy 0.18.0b2）。这是我得到的错误：

如果我在我的本地metpy 副本中手动更改cartopy_utils.py，以便USCOUNTIES 和USSTATES 函数调用MetPyMapFeature 使用'10m' 而不是'20m'，那么错误就会消失。这是 cartopy_utils.py 中的“20m”错误吗？Natural Earth 的网页宣传它提供 1:10m、1:50m 和 1:110m 比例的数据。

我的目标是仅在超过感兴趣区域中观察到的最大风速的 2/3 时绘制风倒钩（这是标准米勒复合图的一部分）。

MetPy 示例中为米勒图执行此操作的代码

返回错误

追溯到线

为了排除故障，我打印了 wspd_500：

这很奇怪，因为我不确定为什么所有这些 nan 值不断出现。u_500​​ 和 v_500 似乎包含所有相关值。当我使用不同的方法绘图时：

我在整个领域都遇到了风倒钩。所以我怀疑是否缺少数据。

此外，当我打印 mask_500 时，我得到一个 False。这应该是一组 True 或 False 布尔值，具体取决于风速是高还是低（至少我认为）。

所以我的问题是：我怎样才能实现这个面具只在给定的压力水平下绘制明显的强风？为什么 mpcalc.wind_speed 方法会从完整的 u 和 v 组件数据数组中返回一堆 nan 值？

作为参考，我使用的是通过虹吸管从 UCAR TDS 获得的 HRRR 数据。