Как читать файл gcip netcdf в python?

У меня есть работающая программа python, которая читает в большом количестве больших файлов netCDF, используя команду Dataset из модуля netCDF4. Вот фрагмент соответствующих частей:

from netCDF4 import Dataset import glob infile_root = 'start_of_file_name_' for infile in sorted(glob.iglob(infile_root + '*')): ncin = Dataset(infile,'r') ncin.close() 

Я хочу изменить это, чтобы читать в файлах netCDF, которые gzipped. После создания самих файлов были созданы gzip; они не внутренне сжаты (т.е. файлы * .nc.gz). Если бы я читал текстовые файлы gzipped, команда была бы:

 from netCDF4 import Dataset import glob import gzip infile_root = 'start_of_file_name_' for infile in sorted(glob.iglob(infile_root + '*.gz')): f = gzip.open(infile, 'rb') file_content = f.read() f.close() 

После того, как он отправился на поиски в течение полчаса и прочитал документацию netCDF4, единственный способ сделать это для файлов netCDF:

 from netCDF4 import Dataset import glob import os infile_root = 'start_of_file_name_' for infile in sorted(glob.iglob(infile_root + '*.gz')): os.system('gzip -d ' + infile) ncin = Dataset(infile[:-3],'r') ncin.close() os.system('gzip ' + infile[:-3]) 

Возможно ли напрямую читать файлы gzip с помощью команды Dataset? Или иначе не называть gzip через os?

• Как gzip при загрузке в s3 с помощью boto
• Python Gzip – добавление к файлу на лету
• Ответ Gzip в Flask / Tornado
• Python распаковывает gzip chunk-by-chunk
• чтение gsipped csv файла в python 3

 import gzip import os import shutil import tempfile import netCDF4 def open_netcdf(fname): if fname.endswith(".gz"): infile = gzip.open(fname, 'rb') tmp = tempfile.NamedTemporaryFile(delete=False) shutil.copyfileobj(infile, tmp) infile.close() tmp.close() data = netCDF4.Dataset(tmp.name) os.unlink(tmp.name) else: data = netCDF4.Dataset(fname) return data