python matplotlib cartopy学习笔记-练习2

JOJOweibo

参考资料：
1.调整字体，colorbar, colormap, 刻度 https://www.cnblogs.com/qqhfeng/p/5567539.html
2. cartopy 官网例子 http://scitools.org.uk/cartopy/d ... xes_grid_basic.html
3. matplotlib 上面关于colorbar的例子 https://matplotlib.org/examples/api/colorbar_basics.html
https://matplotlib.org/examples/api/colorbar_basics.html
4. http://bbs.06climate.com/forum.p ... 3601&extra=page%3D1
注意这帖子里面的第二步程序会自动完成，地图数据自动导入，无需手动设置
地图数据文件路径：.../.local/share/cartopy/shapefiles/natural_earth/physical

其他说明请看代码中的注释，不明白的可以发邮件联系，欢迎一起交流学习

1. #==============================================================================
   
2. 
   
3. #JOJO, IAP, Beiing, Email:mtjsummer@163.com
   
4. #2018-02-07
   
5. 
   
6. #练习2目标：
   
7. #  1.使用函数形式调用cartopy
   
8. #  2.绘制2个子图
   
9. #  3.调整colorbar
   
10. #  4.叠加显著性检验
    
11. #  5.添加标题
    
12. 
    
13. #==============================================================================
    
14. 
    
15. import matplotlib.pyplot as plt
    
16. import matplotlib as mpl
    
17. import cartopy.crs as ccrs
    
18. import cartopy.feature as cfeature
    
19. from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter
    
20. from netCDF4 import Dataset
    
21. import numpy as np
    
22. 
    
23. #==============================================================================
    
24. 
    
25. #读入数据
    
26. f = Dataset('...slp.nc')
    
27. #print(f)
    
28. 
    
29. slp1    = f.variables['slp1'][0,:,:]
    
30. slp2 = f.variables['slp2'][0,:,:]
    
31. lats         = f.variables['lat'][:]
    
32. lons         = f.variables['lon'][:]
    
33. 
    
34. #==============================================================================
    
35. 
    
36. #创建图纸,以及两个子图
    
37. fig = plt.figure(figsize=(8,6))
    
38. ax  = fig.add_subplot(211, projection=ccrs.PlateCarree(central_longitude=180))
    
39. #定义投影方式
    
40. projection = ccrs.PlateCarree()
    
41. 
    
42. #==============================================================================
    
43. 
    
44. #函数形式，调用cartopy，绘制全球地图
    
45. def make_map(ax):
    
46.     projection = ccrs.PlateCarree()
    
47.     ax.set_global()
    
48.     ax.coastlines(linewidth=0.5)
    
49.     '''标注坐标轴'''
    
50.     ax.set_xticks(np.linspace(-180, 180, 5), crs=projection)
    
51.     ax.set_yticks(np.linspace(-90, 90,5), crs=projection)
    
52.     '''zero_direction_label=True 有度的标识，False则去掉'''
    
53.     lon_formatter = LongitudeFormatter(zero_direction_label=True)
    
54.     lat_formatter = LatitudeFormatter()
    
55.     ax.xaxis.set_major_formatter(lon_formatter)
    
56.     ax.yaxis.set_major_formatter(lat_formatter)
    
57.     '''添加网格线'''
    
58.     #gl = ax.gridlines()
    
59.     #ax.grid()
    
60.     return ax
    
61. 
    
62. #==============================================================================
    
63. 
    
64. #设置colorbar
    
65. # Set the colormap and norm to correspond to the data for which the colorbar will be used
    
66. # 指定colorbar为bwr
    
67. cmap = mpl.cm.RdBu_r
    
68. # 设定每个图的colormap和colorbar所表示范围是一样的，归一化
    
69. norm = mpl.colors.Normalize(vmin=-2.0, vmax=2.0)
    
70. # 设置levels
    
71. levels = mpl.ticker.MaxNLocator(nbin=16).tick_values(-2.0, 2.0)
    
72. 
    
73. #==============================================================================
    
74. 
    
75. # 开始绘图
    
76. 
    
77. # 在ax上画地图
    
78. ax = make_map(ax)
    
79. # 填充颜色
    
80. p  = ax.contourf(lons, lats, slp1, levels=levels, cmap=cmap, transform=projection)
    
81. # 画等值线
    
82. pl = ax.contour(lons, lats, slp1, levels=levels, linewidths=0.5, colors='k', transform=projection)  
    
83. # 设置title
    
84. ax.set_title('Reg_SLP1')
    
85. # 添加colorbar
    
86. fig.colorbar(p,ax=ax, extend='both')
    
87. 
    
88. # 在ax1上画图
    
89. ax1  = fig.add_subplot(212, projection=ccrs.PlateCarree(central_longitude=180))
    
90. # 在ax1上画地图
    
91. ax1 = make_map(ax1)
    
92. # 填充颜色
    
93. p1  = ax1.contourf(lons, lats, slp2, levels=levels, cmap=cmap, transform=projection)
    
94. # 画等值线
    
95. pl1 = ax1.contour(lons, lats, slp2, levels=levels, linewidths=0.5, colors='k', transform=projection)  
    
96. # 设置title
    
97. ax1.set_title('Reg_SLP2')
    
98. # 添加colorbar
    
99. fig.colorbar(p1,ax=ax1, extend='both')
    
100. 
     
101. #==============================================================================
     
102. 
     
103. # 保存图片
     
104. # adjust spacing between subplots so ax and ax1 title and ticks labels don't overlay
     
105. fig.tight_layout()
     
106. plt.savefig("./fig/06carotpy.multimap.png")
     
107. 
     

复制代码

随缘

{:eb502:}{:eb502:}

weinihou

{:eb502:}{:eb502:}

一只little_muer

{:eb502:}{:eb502:}{:eb502:}{:eb502:}

叶臻臻

{:eb502:}{:eb502:}{:eb502:}

13top

先分享，后下载，再学习，谢谢楼主{:5_213:}

Byby

楼主我想请问 第69行的norm作用是什么 因为看到后面没有使用它

1521288006

great  i also need  it

弓长铭

感谢！

燃烧的洋葱

学习一下cartopy~~~