问题标签 [metpy]

您好，我收到此错误：ValueError: setting an array element with a sequence。当我运行此代码时：

有问题的行是这样的：

因为我不确定这是如何导致它的，因为有问题的行适用于它所说的不同代码

因为第二个版本中的 p、t 和 dp 也是数组

我正在用绘图中包含的一些复合参数绘制一些观察到的声音，但我对 STP 输出有一个小问题：它总是带有方括号。

我试图检查是否有统一的东西导致了这种情况，但我没有发现任何问题。

想法？

编辑：我很笨，忘记添加代码行

我试图从垂直压力速度估计垂直风速。我找到了metpy.calc.vertical_velocity（欧米茄，压力，温度，混合=0）。

https://unidata.github.io/MetPy/latest/api/generated/metpy.calc.vertical_velocity.html

但是，没有关于它需要什么单元的文档。我给它以 Pa s-1 为单位的 Omega，以 Pa 为单位的压力，以 K 为单位的温度，所有这些都采用 xarray 格式。我希望输出也是一个 xarray 。我收到此错误：“无法从 'pascal * second' 转换为 'dimensionless'”。任何想法为什么它需要无量纲变量以及如何解决这个问题？

我正在尝试为 2016 年密歇根州的龙卷风爆发绘制 level2 雷达数据。当我尝试调用雷达数据文件时，它给了我一个“文件'myfilename'不在注册表中”。但是，当我使用 Metpy 网站示例中使用的确切示例 KTLX 文件时，它可以正常工作。他们住在同一个目录和一切。无论我从 Metpy 站点复制并粘贴确切的代码并更改调用的文件，或者我只是想读取数据，我都会遇到同样的错误。

其次，我注意到示例中使用的文件是 .gz 文件。当我从 NCEI 下载相同的示例文件时，它以 .gz 的形式出现，但 2016 年的文件不是 .gz。为了尝试，我将我的一个文件转换为 .gz ，但仍然没有运气。我应该将所有数据压缩为 .gz 文件还是不需要？

如果您好奇，雷达数据文件是 KGRR20160820_163004_V06 到 KGRR20160820_170654_V06。

感谢您提供任何帮助！

在我的站点上，我们遇到了 MetPy 在尝试调用 surface_based_cape_cin 时返回单位错误的问题

我看到以下错误：

当我检查传入的值 p1、T1 和 Td1 时，它们都有正确的单位（百帕斯卡、度_摄氏度）。

只是为了确保我在调用 surface_based_cape_cin 之前添加了以下内容并检查了结果：

我正在运行以下版本的 MetPy

我不记得在更新到这个版本之前有这个，但我不能确定我看到的问题是在更新之后出现的。

感谢您的任何帮助，您可以提供。

我正在使用 MetPy 生成 SkewT 图。我试图抑制轴标签的自动生成并使用我自己的。生成的自动标签是 hectopascal 和 degree_Celsius。我找不到禁用自动轴标签的方法，因此最终结果是一个看起来凌乱的图。不需要的自动生成的标签用泥红色圈起来：

我尝试通过添加关闭它们

但它不仅删除了标签，还删除了刻度线和所有内容，甚至是情节周围的边框。

我试着变得更有创意

或者

但这些都不起作用。我尝试删除标签的其他方法也删除了刻度值。

我知道我可能遗漏了一些简单的东西，但我似乎找不到它。

这是一个示例代码，它使用 MetPy 中稍微修改的 SkewT 绘图示例演示了这一点：

谢谢

我正在尝试使用metpy为倾斜图上的线条添加标签。我不确定是否可能（在 GitHub 上查看了metpy，似乎他们还没有实现这个功能）。但本质上，我只想为潮湿绝热、干燥绝热和混合比贴上标签。目前，我用来调用和绘制这些函数的代码是：

和

我刚刚安装了 MetPy 第一次使用

在 Mac OS X 10.14 (Mojave) 下。我似乎得到了metpy-0.11.1

然后下面的行生成复制在底部的错误诊断。

TypeError: find_intersections takes 5 parameters, but 3 units were passed如果我改为尝试，我会得到完全相同的错误 ( )

如果我使用导入整个模块，我不会收到错误

我以前从未遇到过任何已建立的 Python 包的导入错误，因此如果能帮助我调试我的设置，我将不胜感激。

----> 1 from metpy.plots import StationPlot 2 中的 TypeError Traceback (last last call last)

~/miniconda3/envs/scipy/lib/python3.6/site-packages/metpy/plots/ init .py in 11 from .ctables import * # noqa: F403 12 from .declarative import * # noqa: F403 ---> 13 from .skewt import * # noqa: F403 14 from .station_plot import * # noqa: F403 15 from .wx_symbols import * # noqa: F403

~/miniconda3/envs/scipy/lib/python3.6/site-packages/metpy/plots/skewt.py in 27 28 from ._util import coloured_line ---> 29 from ..calc import dewpoint, dry_lapse, wet_lapse, vapor_pressure 30 从 ..calc.tools 导入 _delete_masked_points 31 从 ..deprecation 导入metpyDeprecation

~/miniconda3/envs/scipy/lib/python3.6/site-packages/metpy/calc/ init .py in 5 6 from .basic import * # noqa: F403 ----> 7 from .cross_sections import * # noqa : F403 8 from .indices import * # noqa: F403 9 from .kinematics import * # noqa: F403

~/miniconda3/envs/scipy/lib/python3.6/site-packages/metpy/calc/cross_sections.py in 12 13 from .basic import coriolis_parameter ---> 14 from .tools import first_derivative 15 from ..package_tools import Exporter 16 从 ..xarray 导入 check_axis, check_matching_coordinates

~/miniconda3/envs/scipy/lib/python3.6/site-packages/metpy/calc/tools.py 在 104 @preprocess_xarray 105 @units.wraps(('=A', '=B'), ('= A', '=B', '=B')) --> 106 def find_intersections(x, a, b, direction='all', log_x=False): 107 """计算交点的最佳估计。 108

~/miniconda3/envs/scipy/lib/python3.6/site-packages/pint/registry_helpers.py in decorator(func) 246 raise TypeError( 247 "%s 采用 %i 个参数，但通过了 %i 个单位" -- > 248 % (函数名称, count_params, len(args)) 249 ) 250

TypeError: find_intersections 接受 5 个参数，但传递了 3 个单位

我正在尝试按照 Metpy 网站上的示例代码对网格数据集执行横截面分析。

首先，我可以使用 xarray 读取原始 grib2 数据集，并将其转换为 netCDF 格式。

我成功导入了所有必需的包：

然后，我确保可以使用以下命令读取文件：

添加 MetPy 的 CF 解析命令会给我以下 ValueError：

我正在使用 Metpy v0.12.1 和 xarray v0.15.1

而且，有趣的是，当我重新运行相同的代码时，我收到另一条错误消息： NameError: name 'ctables' is not defined

请告知解析数据时出了什么问题。

我想将沿两个坐标之间的横断线的横截面数据保存为 netCDF 文件。

我正在遵循Metpy 的横截面分析中描述的示例代码：

这是添加删除 crs 命令的修改代码：

当我尝试将横截面写入/保存xarray.dataset为 netCDF 文件时，我收到以下错误消息：

我究竟做错了什么？