Попытка открыть серию файлов netCDF с помощью OpenDAP

Question

Я хочу открыть все данные за 1950-2005 гг., Используя xarray и open_mfdataset. 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)

Однако, я не могу понять, что мой список интерпретируется как переменные годы в источнике.

Есть идеи?

Заранее спасибо.

EDIT: 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'

Не правильно ли я ввел путь?

Answer 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)]

См. https://pyformat.info/ для получения дополнительной информации о форматировании строки в Python.