Python 使用metpy和虹吸获取数据时未找到NetCDF属性

我试图以NetCDF格式绘制一些气象数据，通过Unidata虹吸包访问

我已经导入了MetPy文档建议的相关库

import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt
from netCDF4 import num2date
import numpy as np
import xarray as xr
from siphon.catalog import TDSCatalog
from datetime import datetime

import metpy.calc as mpcalc
from metpy.units import units

我已经根据虹吸文档构造了一个数据查询

best_gfs = TDSCatalog('http://thredds.ucar.edu/thredds/catalog/grib/NCEP/GFS/Global_0p25deg/catalog.xml?dataset=grib/NCEP/GFS/Global_0p25deg/Best')

best_ds = best_gfs.datasets[0]
ncss = best_ds.subset()
query = ncss.query()
query.lonlat_box(north=55, south=20, east=-60, west=-90).time(datetime.utcnow())
query.accept('netcdf4')
query.variables('Vertical_velocity_pressure_isobaric','Relative_humidity_isobaric','Temperature_isobaric','u-component_of_wind_isobaric','v-component_of_wind_isobaric','Geopotential_height_isobaric')

data = ncss.get_data(query)

不幸的是，当我尝试使用Metpy文档中的代码解析数据集时

data = data.metpy.parse_cf()

我收到一个错误：“AttributeError:NetCDF:未找到属性”

---

在试图解决这个问题时，我遇到了似乎有相同问题的问题，但解决方案建议——将我的metpy更新到最新版本——对我不起作用。我使用Conda更新了metpy，但遇到了与更新前相同的问题。关于如何解决这个问题还有其他想法吗

现在虹吸中的以下代码

data = ncss.get_data(query)

将从

netcdf4 python

返回一个

Dataset

对象。您还需要一个额外的步骤将其交给xarray，这将使MetPy的

parse\u cf

可用：

from xarray.backends import NetCDF4DataStore
ds = xr.open_dataset(NetCDF4DataStore(data))
data = ds.metpy.parse_cf()

---

现在，虹吸中的以下代码

data = ncss.get_data(query)

将从

netcdf4 python

返回一个

Dataset

对象。您还需要一个额外的步骤将其交给xarray，这将使MetPy的

parse\u cf

可用：

from xarray.backends import NetCDF4DataStore
ds = xr.open_dataset(NetCDF4DataStore(data))
data = ds.metpy.parse_cf()