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