问题描述
我正在使用 netCDF 文件绘制地图。但是,我是新手,进展很小,我所做的只是学习了一些教程。在地图上叠加颜色时我遇到了一些困难。这是我的代码(带有下载的 file):
from mpl_toolkits.basemap import Basemap,cm
from netCDF4 import Dataset
import numpy as np
import matplotlib.pyplot as plt
data = Dataset('filename','r')
lat = data.variables['lat'][:]
lon = data.variables['lon'][:]
time = data.variables['time'][:]
lwe_thickness = data.variables['lwe_thickness'][:]
data.close()
mp = Basemap(projection = 'mill',llcrnrlon=lon.min(),urcrnrlon=lon.max(),llcrnrlat=lat.min(),urcrnrlat=lat.max(),resolution = 'c')
x,y = np.meshgrid(lon,lat)
x,y = mp(lon,lat)
mp.drawcoastlines()
mp.drawstates()
mp.drawcountries()
mp.drawmapboundary()
parallels = np.arange(0.,90,30.)
mp.drawparallels(parallels,labels=[1,0])
meridians = np.arange(-180.,180.,30.)
mp.drawmeridians(meridians,labels=[0,1])
cmesh = mp.pcolormesh(x,y,lwe_thickness,shading='flat',cmap=plt.cm.jet,latlon=True)
cbar = mp.colorbar(cmesh,location='right')
plt.show()
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)