根据经纬度画区域python_用basemap画气象图

basemap 是Python的一个实现地理信息可视化的库，是Matplotlib的一个附加工具包。GrADS在利用“.nc”格式文件的数据绘图时稍显繁琐，basemap就方便很多。

例如现在有一个1948-2010年的逐月月平均气温数据文件“air.mon.mean.nc”，我们试着用basemap画一张1948-2010年1月气温气候场图。

在读取“.nc”文件的数据时，用到了“netCDF4”库。因此我们先import这个库：

from

接着读取“air.mon.mean.nc”文件中的内容，并赋值给fh：

fh=Dataset(r'D:fortran5air.mon.mean.nc')

如果想要查看“.nc”文件中有哪些变量：

>>> fh.variables.keys()
dict_keys(['lat', 'lon', 'time', 'air'])

上面结果表明，“.nc”文件中有纬度数据’lat’，经度数据’lon’，时间数据’time’和变量’air’的数据（也就是气温数据）。

如果想要查看air变量的精度：

>>> fh.variables['air'].shape
(769, 73, 144)

这表示，air有769个时次的数据，纬度有73个格点，经度有144个格点。

最后，我们把这些数据赋值给几个变量：

lons = fh.variables['lon'][:]
lats = fh.variables['lat'][:]
air = fh.variables['air'][:]

这一步就要用到basemap库了，同时为了显示图形，还要导入

Tags: map