这篇博客展示了如何使用Python的netCDF4库读取netCDF文件，并通过mpl_toolkits.basemap库绘制美国降水量的地图。代码详细解释了如何获取并处理数据，以及如何设置地图的各种元素如海岸线、州界、县界、经纬度线等。通过四个不同的时间段，演示了如何展示不同季节的降水量分布。
#coding:utf-8
from mpl_toolkits.basemap import Basemap, cm
from netCDF4 import Dataset as NetCDFFile
import numpy as np
import matplotlib.pyplot as plt
def readNetCDFFile(str):
"Read data from netCDFFile"
nc = NetCDFFile(str)
prcpvar = nc.variables['amountofprecip']
data = 0.01*prcpvar[:]
latcorners = nc.variables['lat'][:]
loncorners = -nc.variables['lon'][:]
lat_0 = nc.variables['true_lat'].getValue()
lon_0 = -nc.variables['true_lon'].getValue()
return prcpvar,data,latcorners,loncorners,lat_0,lon_0
def plotBasemap(str):
"plot the basemap of the USA"
prcpvar,data,latcorners,loncorners,lat_0,lon_0 = readNetCDFFile(str)
m = Basemap(projection='stere', lon_0=lon_0, lat_0=lat_0,\
最低0.47元/天 解锁文章
扫一扫
4275
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
