5个带审图号的中国标准地图

fn='./ETOPO2v2g_f4.nc'
f=addfile(fn)
data=f['z']['15:55','70:140']
# Plot
figure(figsize=[700,400],newfig=False)
ax=axesm(tickfontsize=11,axis=True)
world = shaperead('./country1.shp')
lchina = shaperead('./bou2_4p.shp')
lchina2 = shaperead('./bou2_4l.shp')
geoshow(world,edgecolor='k',size=0.8)
geoshow(lchina,edgecolor='r',size=0.6)
bou1_layer=geoshow(lchina2,edgecolor='r',size=0.6)
levs=arange(-6000,8000,2000)
layer=imshow(data,levs,cmap='MPL_StepSeq_r')
colorbar(layer,ticks=arange(-6000,8000,2000),orientation='vertical',ticklen=0,shrink=1,aspect=25,bold=True,fontsize=12,extendrect=False)
yaxis(tickvisible=True,location='left',tickwidth=2,linewidth=2,ticklength=3) #ticklength刻度线长度，tickwidth刻度线宽度，linewidth边框宽度
yaxis(tickvisible=False,location='right',tickwidth=2,linewidth=2,ticklength=4) #分别调试tick的宽度，边框线宽和tick的长度
xaxis(tickvisible=False,location='top',tickwidth=2,linewidth=2,ticklength=4) 
xaxis(tickvisible=True,location='bottom',tickwidth=2,linewidth=2,ticklength=3)
xlim(70,140)
ylim(15,55)
xticks(arange(70,140.1,10),bold=True,fontsize=12)
yticks(arange(15,55.1,10),bold=True,fontsize=12)
xlabel('Longitude',bold=True,fontsize=12)
ylabel('Latitude',bold=True,fontsize=12)
savefig('./plot179.2.png',dpi=800)

fn='./ETOPO2v2g_f4.nc'
f=addfile(fn)
data=f['z']['15:55','70:140']
# Plot
figure(figsize=[700,400],newfig=False)
ax=axesm(tickfontsize=11,axis=True)
world = shaperead('./country1.shp') 
lchina = shaperead('./bou2_4p.shp')
lchina2 = shaperead('./bou2_4l.shp')
geoshow(lchina,edgecolor='k',size=0.4)
bou1_layer=geoshow(lchina2,edgecolor='k',size=0.4)
levs=arange(0,8000,1000)
layer=imshow(data,levs,cmap='MPL_StepSeq_r')
masklayer(lchina,[layer])
colorbar(layer,ticks=arange(0,8000,1000),orientation='vertical',ticklen=0,extendrect=False,shrink=1,aspect=25,bold=True,fontsize=12)
yaxis(tickvisible=True,location='left',tickwidth=2,linewidth=2,ticklength=3) #ticklength刻度线长度，tickwidth刻度线宽度，linewidth边框宽度
yaxis(tickvisible=False,location='right',tickwidth=2,linewidth=2,ticklength=4) #分别调试tick的宽度，边框线宽和tick的长度
xaxis(tickvisible=False,location='top',tickwidth=2,linewidth=2,ticklength=4) 
xaxis(tickvisible=True,location='bottom',tickwidth=2,linewidth=2,ticklength=3)
xlim(70,140)
ylim(15,55)
xticks(arange(70,140.1,10),bold=True,fontsize=12)
yticks(arange(15,55.1,10),bold=True,fontsize=12)
xlabel('Longitude',bold=True,fontsize=12)
ylabel('Latitude',bold=True,fontsize=12)
ax.scale_bar(0.25,0.24,width=200,linewidth=1,bold=True,fontsize=12,bartype='alternating_bar')
ax.north_arrow(0.4,0.9,width=40,height=40,linewidth=1,bold=True,fontsize=12)
#Add south China Sea
axesm(position=[0.779,0.134,0.15,0.2],axison=False,frameon=True)
geoshow(lchina,edgecolor='k',size=0.2)
bou1_layer=geoshow(lchina2,edgecolor='k',size=0.1)
layer=imshow(data,levs,cmap='MPL_StepSeq_r')
masklayer(lchina,[layer])
xlim(106,123)
ylim(2,23)
savefig('./plot179.1.png',dpi=800)

免费分享5个带审图号的中国矢量地图：

原文链接https://mp.weixin.qq.com/s?__biz=MzUyNzczMTI4Mg==&mid=2247683716&idx=1&sn=d3427f54fff6307ba6b99387c61b78c3&chksm=fa7751b9cd00d8af614d44a54c430c10c7b43e487f2de7fe4d5b4397c79a9988f381f9ff16d2&token=1722637836&lang=zh_CN#rd