以下是一段实现这个需求的代码:
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
import pygrib
# 读取grib格式的文件
grbs = pygrib.open('文件路径')
grb = grbs.select(name='T')[0]
# 获取经纬度数据
lats, lons = grb.latlons()
# 创建地图实例
map = Basemap(projection='cyl', llcrnrlat=-90,