python讀取NC文件畫風速等值線圖

来源：互联网发布：mac怎么免费下载软件编辑：程序博客网时间：2024/06/06 17:34

關於用basemap畫圖，可以去看我的另一篇博客，那篇很全面，這篇主要學習怎麼讀取nc等文件

有些部分沒有調好，比如等值線的間隔、顏色等，不過會繼續學習的。現在先做一下筆記。

#coding=utf-8__author__ = 'kong'import matplotlib.pyplot as pltimport numpy as npfrom mpl_toolkits.basemap import Basemap, cmimport netCDF4 as ncfilename='C:\\Users\\kong\\PycharmProjects\\replacencl\\meow'data=nc.Dataset(filename)for i in data.variables.keys():  print(i)lat= data.variables['XLAT_U'][:].squeeze()lon= data.variables['XLONG_U'][:].squeeze()u=data.variables['U'][:].squeeze()idate=data.variables['Times']m=Basemap(resolution='l',area_thresh=10000,projection='cyl',llcrnrlon=50,urcrnrlon=150,llcrnrlat=0,urcrnrlat=70)fig1 = plt.figure()ax = fig1.add_axes([0.1,0.1,0.8,0.8])#clevls_hgt=[-80,-60,-40,-20,0,20,40,60,80]x, y = m(lon, lat)#parallels = np.arange(-90.,90,30.)#meridians = np.arange(0.,360.,60.)CS2 = m.contour(x,y,u[0,:,:],origin='lower',cmap=plt.get_cmap('Blues'),extend='both')CS2.clabel()m.drawcoastlines(linewidth=0.2)plt.title('U',size=20)#plt.xlim(50,150)#plt.ylim(0,70)plt.show()

真是丑到一個境界ORZ

改了，好看多了

filename='C:\\Users\\kong\\PycharmProjects\\replacencl\\meow'data=nc.Dataset(filename)for i in data.variables.keys():  print(i)lat= data.variables['XLAT_U'][:].squeeze()lon= data.variables['XLONG_U'][:].squeeze()u=data.variables['U'][:].squeeze()idate=data.variables['Times']m=Basemap(resolution='l',area_thresh=10000,projection='cyl',llcrnrlon=50,urcrnrlon=150,llcrnrlat=0,urcrnrlat=70)fig1 = plt.figure()ax = fig1.add_axes([0.1,0.1,0.8,0.8])#clevls_hgt=[-80,-60,-40,-20,0,20,40,60,80]x, y = m(lon, lat)m.drawparallels(np.arange(0,70,10.),labels=[1,1,0,0],fontsize=15)m.drawmeridians(np.arange(50,150,20),labels=[0,0,0,1],fontsize=15)#curve=m.contour(x,y,u[0,:,:],colors='k')CS2=m.contourf(x,y,u[0,:,:])m.colorbar(CS2)m.drawcoastlines(linewidth=0.2)plt.title('U',size=20)plt.show()

0 0