用python和GDAL 读取GRIB数据

from osgeo import gdal
import numpy as np

path = './data/201912/anl_surf125.2019120100'
dataset = gdal.Open(path)
number_bands = dataset.RasterCount

for i in range(1, number_bands):

    band = dataset.GetRasterBand(i)
    data_array = band.ReadAsArray()
    print(data_array,data_array.shape)