很多教程里往大地图上插入南海小地图是通过 fig.add_axes([left, bottom, width,
height]) 方法添加的。为了控制小地图的大小、让它位于大地图的右下角并与边框对齐，需要手动
调整 [left, bottom, width, height] 这四个参数的值。这样不仅非常繁琐，并且在画组图时更是会
让人无从下手。
本文将介绍的代码可以自动设定小地图的位置，不必再反复试错。实现上先用 get_position
方法获取大地图在画布上的位置作为作为参考。考虑到 GeoAxes 的 aspect_ratio 固定为 1，小地
图在经纬度范围决定好之前形状是未知的，所以等到小地图的经纬度范围设定好后再等比例缩小其宽
高，最后利用 set_position 方法移动到大地图的右下角。
这个方法的缺点有两个：一是在调用前需要提前设定好大地图和小地图的经纬度范围，这有点
反直觉；二是仅适用于方形边框的地图。如果读者有更好的方法还请多多交流。
import matplotlib.pyplot as pltimport numpy as npimport
matplotlib.transforms as mtransformsfrom matplotlib import
rcParamsimport cartopy.crs as ccrsimport cartopy.feature as
cfeaturefrom cartopy.io.shapereader import Readerfrom
cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter#
设置字
体
.config={'font.family':'Times New Roman','font.size':
10,'mathtext.fontset':'stix'}rcParams.update(config)def
adjust_sub_axes(ax_main, ax_sub, shrink): ''' 将 ax_sub 调整到
ax_main 的右下角. shrink 指定缩小倍数. 使用前需要设定好两个 Axes 的经纬度范
围. ''' bbox_main = ax_main.get_position() bbox_sub =
ax_sub.get_position() #
使
shrink=1
时
ax_main
与
ax_sub
等宽或等高
.
if bbox_sub.width > bbox_sub.height: ratio = bbox_main.width /
bbox_sub.width * shrink else: ratio = bbox_main.height /
bbox_sub.height * shrink wnew = bbox_sub.width * ratio hnew =
bbox_sub.height * ratio bbox_new =
mtransforms.Bbox.from_extents(bbox_main.x1 - wnew,
bbox_main.y0,bbox_main.x1, bbox_main.y0 + hnew)
ax_sub.set_position(bbox_new)#
地图范围
.extents_main = [75, 135, 5,
55]extents_sub = [105, 125, 0, 25]lonmin, lonmax, latmin, latmax =
extents_mainxticks = np.arange(lonmin, lonmax + 1, 10)yticks =
np.arange(latmin, latmax + 1, 10)#
分别初始化主图和子图
.proj =
ccrs.PlateCarree()fig = plt.figure(figsize=(10, 8))subplot_kw =
{'projection': proj}axes_main = fig.subplots(3, 3,
subplot_kw=subplot_kw)axes_sub = fig.subplots(3, 3,
subplot_kw=subplot_kw)#
读取
shp
文件
.geom_river =
list(Reader('F:/Rpython/lp27/data/river1.shp').geometries())geom_china1
=
list(Reader('F:/Rpython/lp27/data/china1.shp').geometries())geom_china2
=
list(Reader('F:/Rpython/lp27/data/china2.shp').geometries())geom_lakes
= list(Reader('F:/Rpython/lp27/data/ne_50m_lakes.shp').geometries())for
ax_main, ax_sub in zip(axes_main.flat, axes_sub.flat): #
设置坐标轴范
围
. ax_main.set_extent(extents_main, crs=proj)
ax_sub.set_extent(extents_sub, crs=proj) adjust_sub_axes(ax_main,
ax_sub, shrink=0.3) #
设置坐标轴刻度
. ax_main.set_xticks(xticks,
crs=proj) ax_main.set_yticks(yticks, crs=proj)
ax_main.xaxis.set_major_formatter(LongitudeFormatter())
ax_main.yaxis.set_major_formatter(LatitudeFormatter()) #
添加
shp
文件
.
for ax in [ax_main, ax_sub]: ax.stock_img()
ax.coastlines(resolution='10m', lw=0.2)
ax.add_feature(cfeature.LAKES, alpha=0.5)
ax.add_geometries(geom_river, proj, fc='none', ec='b', lw=0.2)
