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

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

首先，安装netCDF4和Basemap，Windows下和Linux下会稍微有些不一样，请自行百度。 NetCDF4：​​​ Basemap：​​​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 Title

运行结果：

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。