MetPy：xarray 数据集中横截面的投影信息？

Question

metpy.cross_section() 的输入是“必须附加投影信息”的 xarray 数据集或 DataArray。

要使用 MetPy 横截面，我需要编写自己的 xarray 数据集，因为我的带有元数据的 NWP 数据无法自动读入 xarray 数据集。但是，我不清楚 MetPy 需要哪些元数据才能理解我的投影。

我意识到是metpy.parse_cf() 为MetPy 提供了投影信息，但我在metpy.parse_cf() 上找不到任何信息，即哪些变量是必要的。目前，metpy.parse_cf() 对我来说是一个黑盒子。

Answer 1

横截面函数所需的“附加投影信息”具有crs 坐标和axis 坐标属性，MetPy 在使用ds.metpy.parse_cf() 时添加到数据集（如the xarray with MetPy tutorial 中所述）。

MetPy 根据CF Conventions 使用数据集中的元数据创建crs 坐标和axis 坐标属性。如果您的 GRIB-to-xarray 数据集读取器包含 Lambert Conformal 网格上数据的此信息，则应正确创建 crs 坐标（您可以通过确保它读取“投影：lambert_conformal_conic”来检查这一点）。但是，如果您的数据集阅读器不包含此元数据，您可能需要自己添加。