使用Python的netCDF4和matplotlib.basemap包进行气象数据的可视化

来源：互联网发布：钻石篮球联赛官网数据编辑：程序博客网时间：2024/05/17 06:58

首先，安装netCDF4和Basemap，Windows下和Linux下会稍微有些不一样，请自行百度。
NetCDF4：https://github.com/Unidata/netcdf4-python
Basemap：http://matplotlib.org/basemap/users/installing.html

下面是一个示例代码：

from netCDF4 import Datasetimport numpy as npimport matplotlib.pyplot as pltfrom mpl_toolkits.basemap import Basemapmeteo_file = "/home/theone/Data/GreatKhingan/MERRA/MERRA2_400.inst1_2d_lfo_Nx.20131201.nc4"fh = Dataset(meteo_file, mode='r')# 获取每个变量的值lons = fh.variables['lon'][:]lats = fh.variables['lat'][:]tlml = fh.variables['TLML'][:]tlml_units = fh.variables['TLML'].units# 经纬度平均值lon_0 = lons.mean()lat_0 = lats.mean()m = Basemap(lat_0=lat_0, lon_0=lon_0)lon, lat = np.meshgrid(lons, lats)xi, yi = m(lon, lat)# Plot Data# 这里我的tlml数据是24小时的，我这里只绘制第1小时的（tlml_0）tlml_0 = tlml[0:1:, ::, ::]cs = m.pcolor(xi, yi, np.squeeze(tlml_0))# Add Grid Lines# 绘制经纬线m.drawparallels(np.arange(-90., 91., 20.), labels=[1,0,0,0], fontsize=10)m.drawmeridians(np.arange(-180., 181., 40.), labels=[0,0,0,1], fontsize=10)# Add Coastlines, States, and Country Boundariesm.drawcoastlines()m.drawstates()m.drawcountries()# Add Colorbarcbar = m.colorbar(cs, location='bottom', pad="10%")cbar.set_label(tlml_units)# Add Titleplt.title('Surface Air Temperature')plt.show()fh.close()

运行结果：
Basemap可视化

0 0