分享Python读取NC文件并画阴影图的源代码

victory_19

这段时间在家园里学习Python，从安装到画图，受益匪浅，学习了很多老师的Python代码，这里分享一下我学习之后画的nc阴影图~谢谢家园~代码如下：
#!/usr/bin/python
#coding=utf-8
# -*- coding: gbk -*-
#coding=gbk
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np
import netCDF4
from netCDF4 import Dataset

#绘制基础地图
map = Basemap(projection='cyl',llcrnrlat=15,urcrnrlat=55,llcrnrlon=70,urcrnrlon=140,resolution='l')
map.drawmapboundary()   
map.drawstates()        
map.drawcoastlines()   
map.drawcountries()     
lons,lats=map.makegrid(71,41)
lats=lats[::-1]
x,y=map(lons,lats)
map.drawparallels(np.arange(15.,56.,10.),labels=[1,0,0,0],fontsize=12)
map.drawmeridians(np.arange(75.,141.,15.),labels=[0,0,0,1],fontsize=12)
plt.title('China basemap u-shaded')

nc_obj=Dataset('E:\\uwnd.2015.nc')

#读取数据
LEVEL=(nc_obj.variables['level'][1])
TIME=(nc_obj.variables['time'][1])
LAT=(nc_obj.variables['lat'][:])
LON=(nc_obj.variables['lon'][:])
UWIND=(nc_obj.variables['uwnd'][1,1,:,:])

#将原始数据变成网格数据形式
X,Y = np.meshgrid(LON,LAT)

#绘制等值线和阴影图
curve=map.contour(X,Y,UWIND)
#这里是使用了NCL的colorbar
cmap_color=plt.cm.get_cmap('RdYlBu_r')      
shade=map.contourf(X,Y,UWIND,cmap=cmap_color)
cbar=map.colorbar()
cbar.ax.tick_params(labelsize=12)
plt.savefig('u-shaded.png')

效果图如下：

                               
登录/注册后可看大图

xiaofeixiazyh

不错不错，鼓励一下

小其其格

不错不错，鼓励一下

分享是最高级的学习方式！

zoeyyuanyuan

谢谢楼主分享~~

贫道敬孔

很通俗易懂了，多谢楼主分享

stand33

感谢楼主分享，提供学习参考，非常感谢！

dongfangxing

提供学习参考，分享和介绍！

subtropical

谢谢楼主分享，楼主我试用你的程序时候出错了，请问怎么解决？

珮瑶瑶

楼主您好，TIME=(nc_obj.variables['time'][1])，这里把1改成别的时间，怎么图形没有变化呢？请指导，谢谢

珮瑶瑶

小其其格 发表于 2019-5-10 11:52
不错不错，鼓励一下

分享是最高级的学习方式！

小其其格您好，TIME=(nc_obj.variables['time'][1])，这里把1改成别的时间，怎么图形没有变化呢？请指导，谢谢