Python气象绘图教程（二十一）—Cartopy_10

本节提要：尽量符合实际应用要求的刻度与经纬度

---

 

---

基于很多同志询问添加经纬度办法，系统性重编了地图的经纬度添加方式。各种投影中以矩形投影PlateCarree最为方便，可以套用matplotlib.mticker的形式。在最新的0.18版本的cartopy中，虽然还不完善，但是终于能直接绘制兰勃脱下的标签了。墨卡托在官网上有示例。

一、PlateCarree投影下的两种投影方式

（1）gridlines方式

这个方式是cartopy自带的，也只能在cartopy中使用。需要引入cartopy的经纬度格式部件，然后正常添加。这种方式属于对全部投影的普适办法，但是在矩形投影上就不如matplotlib的方式。

•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  
•  

import matplotlib.pyplot as plt                                          import cartopy.crs as ccrs                                               import cartopy.feature as cfimport numpy as npimport matplotlib.ticker as mtickerfrom cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTERplt.rcParams['font.sans-serif']=['SimHei']fig=plt.figure(figsize=(2,2),dpi=400)ax=fig.add_axes([0,0,1,1],projection=ccrs.PlateCarree(central_longitude=120))ax.add_feature(cf.LAND.with_scale('110m'))ax.add_feature(cf.OCEAN.with_scale('110m'))ax.add_feature(cf.COASTLINE.with_scale('110m'),lw=0.4)ax.add_feature(cf.RIVERS.with_scale('110m'),lw=0.4)################################################################gl=ax.gridlines(draw_labels=True,linestyle=":",linewidth=0.3,color='k')gl.top_labels=False #关闭上部经纬标签                                  gl.right_labels=Falsegl.xformatter = LONGITUDE_FORMATTER  #使横坐标转化为经纬度格式            gl.yformatter = LATITUDE_FORMATTER                                        gl.xlocator=mticker.FixedLocator(np.arange(-180,180,30))      gl.ylocator=mticker.FixedLocator(np.arange(-90,90,30)) gl.xlabel_style={'size':3}#修改经纬度字体大小