我使用虹吸请求了一个netCDF子集,并形成了一个查询来检索边界框中的变量:
from siphon.catalog import TDSCatalog
cat = TDSCatalog("https://thredds.ucar.edu/thredds/catalog/grib/NCEP/GFS/Global_onedeg/catalog.xml?dataset=grib/NCEP/GFS/Global_onedeg/Best")
ncss = cat.datasets[0].subset()
query = ncss.query()
query.variables("Absolute_vorticity_isobaric")
query.lonlat_box(north=34., south=33., west=-102., east=-101.)
query.accept("netcdf4")
我正在寻求一种可靠、简洁的方法来获取变量的坐标值,特别是时间和垂直水平。一种可行但不切实际的方法是请求并使用整个数据集。你知道吗
import xarray as xr
query.all_times()
data = ncss.get_data(query)
datastore = xr.backends.NetCDF4DataStore(data)
ds = xr.open_dataset(datastore).metpy.parse_cf()
xarray.DataArray
对于作为xarray.DataArray
的数据集的每个变量,调用ds.VARIABLE.metpy.DIMENSION
使MetPy自动返回适当的坐标变量(不管它叫什么,例如lat
、lon
、time
、time1
、altitude_above_msl
、isobaric3
、height_above_ground1
),其中DIMENSION
是time
、vertical
、x
和y
之一。你知道吗
在本例中,ds.Absolute_vorticity_isobaric.metpy.time
返回ds.time
,ds.Absolute_vorticity_isobaric.metpy.vertical
返回ds.isobaric2
。将.values
添加到调用只返回numpy.ndarray
,其中包含我一直试图获取的值。因此,调用ds.Absolute_vorticity_isobaric.metpy.time.values
会产生以下结果(在下面被截断):
array(['2019-11-17T00:00:00.000000000', '2019-11-17T03:00:00.000000000',
'2019-11-17T06:00:00.000000000', ..., '2020-01-02T06:00:00.000000000',
'2020-01-02T09:00:00.000000000', '2020-01-02T12:00:00.000000000'],
dtype='datetime64[ns]')
调用ds.Absolute_vorticity_isobaric.metpy.time.values
和ds.Absolute_vorticity_isobaric.metpy.vertical.values
将只返回NumPy数组,这就是我要寻找的。
虽然上面的代码确实满足了我的要求,但只运行一个变量就花了将近一分半钟,而且(我假设)对UCAR服务器不必要地征税。有没有任何方法可以在不加载所有数据本身的巨大开销的情况下获得上述输出?你知道吗
如果您关心原始方法的性能,并且只希望提取时间和垂直坐标,我建议使用OPENDAP来访问您的数据,而不是NCSS。这将首先简单地获取元数据,然后缓慢地加载您请求的数据(在您的情况下是时间和垂直坐标)。使用MetPy v0.11或更新版本,使用您感兴趣的TDS目录的示例脚本如下所示:
在我的系统上运行大约需要半秒钟。你知道吗
如果MetPy早于v0.11,则在打开数据集时需要使用
.metpy.parse_cf()
,如下所示:相关问题 更多 >
编程相关推荐