尝试使用 OpenDAP 打开一系列 netCDF 文件
Trying to open a series of netCDF files with using OpenDAP
我想使用 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'
我是不是输入的路径不对?
所有文件都命名为 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/。
我想使用 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'
我是不是输入的路径不对?
所有文件都命名为 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/。