读取压缩的二进制文件（.grib2.bz2）

Question

我已经从这个列表中下载了一个文件https://opendata.dwd.de/weather/nwp/icon-eu/grib/03/t_2m/（实际文件名每天都会变化），它们是 bz2 压缩的。

我可以使用例如读取解压缩的文件

import xarray as xr

# cfgrib + dependencies are also required

grib1 = xr.open_dataset("icon-eu_europe_regular-lat-lon_single-level_2020101212_001_ASHFL_S.grib2", engine='cfgrib')

但是，我想读入压缩文件。

我尝试过类似的事情

with bz2.open("icon-eu_europe_regular-lat-lon_single-level_2020101818_002_ASWDIFD_S.grib2.bz2", "rb") as f:

    xr.open_dataset(f, engine='cfgrib')

但这不起作用。

我正在寻找任何以编程方式读取压缩文件的方法。

Answer 1

我在处理数值天气预报数据时遇到了同样的问题。我在这里所做的是下载文件并将其保存为二进制对象（例如使用urlopen或requests）。将此对象传递给以下函数：import bz2, shutilfrom io import BytesIOfrom pathlib import Pathdef bunzip_store(file: BytesIO, local_intermediate_file: Path):    with bz2.BZ2File(file) as fr, local_intermediate_file.open(mode="wb") as fw:        shutil.copyfileobj(fr, fw)解压后的文件将存储在local_intermediate_file. 现在您应该能够打开该文件。