python GDAL读取grib2

安装

1. https://www.lfd.uci.edu/~gohlke/pythonlibs/选择合适的gdal版本下载,我下的是GDAL-3.0.2-cp37-cp37m-win_amd64.whl
2. 开始pip 安装

> pip install D:\Python\whl\GDAL-3.0.2-cp37-cp37m-win_amd64.whl
Looking in indexes: http://mirrors.aliyun.com/pypi/simple/
Processing d:\python\whl\gdal-3.0.2-cp37-cp37m-win_amd64.whl
Installing collected packages: GDAL
Successfully installed GDAL-3.0.2

安装成功提示！

python程序运行

from osgeo import gdal

ds = gdal.Open('1.grib2') #读取文件
print(ds.RasterXSize,ds.RasterYSize,ds.RasterCount) #获取经度、纬度和变量数
print('projection is:',ds.GetProjection()) #获取投影信息
geotransform = ds.GetGeoTransform() #获取参数
print(geotransform) 
dt=ds.GetRasterBand(1) #获取第1个通道的值,第几个通道代表什么变量需要wgrib2来获知。gdal没有读取变量名的方法
data=dt.ReadAsArray()  #赋值给data

缺点： gdal无法直接读取变量名，需要借助wgrib2查看变量名所在通道。