0

我想使用 xarray 和 open_mfdataset 打开 1950-2005 年的所有数据。 https://www.esrl.noaa.gov/psd/thredds/catalog/Datasets/ncep.reanalysis/surface/catalog.html

这是我到目前为止所做的：

来源= https://www.esrl.noaa.gov/psd/thredds/catalog/Datasets/ncep.reanalysis/surface/air.sig995.years.nc

    files = [source for years in range(1950,2005,1)]
    ds=xr.open_mfdataset(files)
    print(ds)

但是，我似乎无法将我的列表解释为读入源中的变量年份。

有任何想法吗？

先感谢您。

编辑： path = 'https://www.esrl.noaa.gov/psd/thredds/catalog/Datasets/ncep.reanalysis/surface' files = ['{0}/air.sig995.{1:04d}.nc'.format(path, years) for years in range(1950,2005,1)] print(files) nc = netCDF4.MFDataset(files)

这是我正在使用的代码。当我尝试打开这些文件时，出现错误：

OSError: [Errno -90] NetCDF: file not found: b'https://www.esrl.noaa.gov/psd/thredds/catalog/Datasets/ncep.reanalysis/surface/air.sig995.1948.nc'

我没有正确输入路径吗？

4

1 回答 1

1

所有文件都命名为air.sig995.YYYY.nc，因此您需要以下内容：

files = ['air.sig995.{0:04d}.nc'.format(years) for years in range(1950,2005,1)]

产生：

In [2]: files
Out[2]: 
['air.sig995.1950.nc',
 'air.sig995.1951.nc',
 'air.sig995.1952.nc',
 'air.sig995.1953.nc',
 .....

您还可以在此处轻松包含（远程）路径（如果需要）：

path = '/some/file/path'
files = ['{0}/air.sig995.{1:04d}.nc'.format(path, years) for years in range(1950,2005,1)]

有关 Python 中字符串格式的更多信息，请参阅https://pyformat.info/。

于 2019-03-25T13:02:54.113 回答