本节提要:尽量符合实际应用要求的刻度与经纬度
基于很多同志询问添加经纬度办法,系统性重编了地图的经纬度添加方式。各种投影中以矩形投影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}#修改经纬度字体大小 gl.ylabel_style=