python气象绘图速成_分享一个Python绘图例子的网站（气象绘图），建议有兴趣者学习...

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看，没有帐号？立即注册

x

本帖最后由 又是那隻貓 于 2017-5-7 22:38 编辑

分享一个Python绘图例子的网站；计算数据，最终都为了绘制图形。去年也无意中浏览过这个网站，今年再次捡起Python，再一次打开这个绘图指导的网站，有很多py的现成例子，大家有兴趣的可以去浏览和下载！

http://matplotlib.org/basemap/users/examples.html

例如：

from mpl_toolkits.basemap import Basemap, cm# requires netcdf4-python (netcdf4-python.googlecode.com)from netCDF4 import Dataset as NetCDFFile

import numpy as np

import matplotlib.pyplot as plt

# plot rainfall from NWS using special precipitation

# colormap used by the NWS, and included in basemap.

nc = NetCDFFile('../../../examples/nws_precip_conus_20061222.nc')

= nc.variables['amountofprecip']

data = 0.01*prcpvar[:]

latcorners = nc.variables['lat'][:]

loncorners = -nc.variables['lon'][:]

lon_0 = -nc.variables['true_lon'].getValue()

lat_0 = nc.variables['true_lat'].getValue()

# create figure and axes instances

fig = plt.figure(figsize=(8,8))

ax = fig.add_axes([0.1,0.1,0.8,0.8])

# create polar stereographic Basemap instance.

m = Basemap(projection='stere',lon_0=lon_0,lat_0=90.,lat_ts=lat_0,\

llcrnrlat=latcorners[0],urcrnrlat=latcorners[2],\

llcrnrlon=loncorners[0],urcrnrlon=loncorners[2],\

rsphere=6371200.,resolution='l',area_thresh=10000)

# draw coastlines, state and country boundaries, edge of map.

m.drawcoastlines()

m.drawstates()

m.drawcountries()

# draw parallels.

parallels = np.arange(0.,90,10.)

m.drawparallels(parallels,labels=[1,0,0,0],fontsize=10)

# draw meridians

meridians = np.arange(180.,360.,10.)

m.drawmeridians(meridians,labels=[0,0,0,1],fontsize=10)

ny = data.shape[0];

nx = data.shape[1]

lons, lats = m.makegrid(nx, ny) # get lat/lons of ny by nx evenly space grid.

x, y = m(lons, lats) # compute map proj coordinates.

# draw filled contours.

clevs = [0,1,2.5,5,7.5,10,15,20,30,40,50,70,100,150,200,250,300,400,500,600,750]

cs = m.contourf(x,y,data,clevs,cmap=cm.s3pcpn)

# add colorbar.

cbar = m.colorbar(cs,location='bottom',pad="5%")

cbar.set_label('mm')

# add title

plt.title(prcpvar.long_name+' for period ending '+prcpvar.dateofdata)

plt.show()

再附赠几个网站：

Python 常用的安装Module方式汇总http://www.jb51.net/article/113084.htm

利用numpy+matplotlib绘图的基本操作教程 http://www.jb51.net/article/112815.htm

利用matplotlib+numpy绘制多种绘图的方法实例 http://www.jb51.net/article/112805.htm