【气候软件】Python9：绘制中国/世界地图

气象数据可视化主要依靠matplotlib做绘图，其他库包为辅助，包括数据处理的，地理信息处理的等等。

绘图主要有六步（“六部曲”）（不是绝对的方法，视具体情况使用）:

• 引入库包：import matplotlib.pyplot as plt

• 设定画布：fig=plt.figure()

• 导入数据：之前已有介绍导入nc文件格式数据（【气候软件】Python读取气象数据 NetCDF文件（***.nc））和导入txt文本格式数据（【气候软件】Python2：读取TXT文本格式的数据）。

• 线图命令：plt.plot(x,y,lw=,ls=，c=,alpha=)

• 出图：plt.show()

• 存图：fig.savefig("···")

  
  

只要按照以上六步，基本绘图没有问题！！！

0 P0 b. J2 X& \7 Y2 q; C                               
登录/注册后可看大图

绘制中国区域地图

下面直接上代码（关于代码的解释已在注释中详细说明，应该能理解）

例1：绘制中国区域基本地图轮廓，包含省界

（如果没有下载相应库包，请先在anaconda或cmd命令行去下载！！！代码的详细说明已备注在注释中，不再额外赘述。）

• 
  , _! \3 r) t% g$ I, c. q4 S
• 
  
• 
  . u6 L) Z- t( B1 T6 _2 S
• 
  . h1 p9 U- F8 r
• 
  
• 
  - T; J4 J, L2 G# C+ g' d. f
• 
  
• 
  
• 
  
• 
  + T# y/ m" X3 Y2 X  ]' t
• 
  
• 
  , l8 q8 R; R, x7 g8 \1 Z
• 
  
• 
  
• 
  
• 
  ; w2 O$ q3 {( L8 Z' a
• 
  
• 
  / G& Y) h, w0 w
• 
  
• 
  ; L& b( ~, u: u% Z4 B/ Y0 U0 \" _
• 
  / {/ W4 q. B+ W' J0 G
• 
  
• 
  8 K0 V9 X8 g/ a0 V! `  M
• 
  , c6 ^) [: f8 g! d
• 
  0 A. E" S1 n' T
  

#绘制中国区域基本地图轮廓，包含省界#引库frompyecharts.charts importGeofrompyecharts.render importmake_snapshotfromsnapshot_phantomjs importsnapshot  #下载库包snapshot-phantomjs（中间不是下划线，是短横线）# ->Geo 是函数注解，表示该函数返回值为Geo对象defgeo_effectscatter()-> Geo:# 以下为链式调用方法声明对象c = (Geo()# 添加底部地图.add_schema(maptype="china"))returnc# 生成对象c = geo_effectscatter()# 渲染地图c.render()# 生成图片make_snapshot(snapshot, c.render(), "map0.png")
6 {4 r& i2 Z) v5 X

4 P6 p" Z5 a) Y7 m0 C

例2：绘制中国区域基本地图轮廓，并标记各省会城市名

（如果没有下载相应库包，请先在anaconda或cmd命令行去下载！！！代码的详细说明已备注在注释中，不再额外赘述。）

• 
  , Q+ a9 A' S% R! r
• 
  
• 
  3 G- z. Q9 x; {7 C, j% E
• 
  2 U$ T2 f9 E: v( @( e* `5 H
• 
  
• 
  ! \+ S/ {. o+ b
• 
  7 x  {, F6 ]" I0 n8 j2 i
• 
  ) k" z% H9 e) a6 R# `. Z8 D
• 
  $ L& L3 Y% {7 L4 o
• 
  # w2 U% j" [+ e' c0 _+ ~: [
• 
  % G( D6 c* G+ S
• 
  
• 
  
• 
  ) x6 u2 b7 I1 ~$ N
• 
  - P- d) r' W+ v  _& {, M5 p
• 
  " h# D! {% [3 [9 ?' E$ B
• 
  
• 
  1 M/ y4 ]* T2 a$ e
  ( p- E9 W. |) L* j" k; z! {

#绘制中国地图，并标记各省会城市名#引库frompyecharts importoptions asoptsfrompyecharts.charts importMapfrompyecharts.faker importFakerfrompyecharts.render importmake_snapshotfromsnapshot_phantomjs importsnapshotc = (Map().add("中国基本地图", [list(z) forz inzip(Faker.provinces, Faker.values())], "china").set_global_opts(title_opts=opts.TitleOpts(title="Map-中国地图")).render("map_china.html"))# 生成图片make_snapshot(snapshot, c,  "map.png")
+ B$ T* x8 n5 R5 w- x! L

9 X1 p6 R) c* o4 K4 W* j

例3：绘制中国区域地图，包含海岸线

（如果没有下载相应库包，请先在anaconda或cmd命令行去下载！！！代码的详细说明已备注在注释中，不再额外赘述。）

+ ~8 D; ~# o8 D# q/ |  `2 i9 f

• 
  
• 
  $ A, @9 e! R4 H6 B$ a2 v6 w
• 
  
• 
  
• 
  ! D) K  x9 y3 I$ \, K
• 
  : t" U0 H& @0 n1 A# R
• 
  
• 
  
• 
  
• 
  
• 
  - @$ j1 x( l5 U7 I3 Y  k
• 
  3 a& D  C9 u' y
• 
  . `$ ?! E0 E7 ^, D
• 
  
• 
  / a& N# i, a! m  d5 i
• 
  
• 
  
• 
  ) h/ x; S& \" [7 b
• 
  
• 
  
• 
  + k: ^! n$ Z$ F3 V
• 
  
• 
  3 i5 G5 o1 A+ _
• 
    C9 ^. R0 f  ?% ~$ \  Z8 S
• 
  
• 
  * b7 q2 t9 X, f: _( M& |
• 
  0 f, S% n- ~2 f& k! j& m# ~
• 
  
• 
  
• 
  
• 
  " k* ?* l- {- z) ^, u% r- E$ q
• 
  
• 
  
• 
  % y1 ~5 c* l& p* h5 Q4 Z8 X5 U
• 
  
• 
  3 o4 I" O4 Z7 o
• 
  
• 
  
• 
  ! a7 B) `3 V2 R) C5 u
• 
  
• 
  
• 
  + |) K' L4 w) Q2 U
• 
  
• 
  
• ​
  $ O# r! r( X8 Y- s$ l+ `
• 
  & J$ d; s, X3 U! Q3 C  N) T2 u: S
• 
  . N) z) L' d8 q* \3 p, Y6 J# T
• 
  
• 
  
• 
  
• 
  
  

##绘制有海岸线的中国地图#引库importpandas aspdimportnumpy asnpimportcartopy.crs asccrsimportcartopy.feature ascfeaturefromcartopy.mpl.gridliner importLONGITUDE_FORMATTER, LATITUDE_FORMATTERimportcartopy.mpl.ticker asctickerimportcartopy.io.shapereader asshpreaderimportmatplotlib.pyplot asplt#数据data = pd.read_csv("meteo.txt", skiprows=1, sep='\s+', header=None, names=['station','lon','lat','tem', 'pre', 'sun'])     print(data)#建立画布fig2 = plt.figure(figsize = (15,15))proj = ccrs.PlateCarree(central_longitude = 115) #设置一个圆柱投影坐标，中心经度115°Eleftlon, rightlon, lowerlat, upperlat = (70,140,15,55)#建立子区域#ax1 = fig.add_axes([left, bottom, width, height])f2_ax1 = fig2.add_axes([0.2, 0.2, 0.6, 0.6], projection = proj)#在画布的绝对坐标建立子图f2_ax1.set_extent([leftlon, rightlon, lowerlat, upperlat], crs=ccrs.PlateCarree())#海岸线，50m精度f2_ax1.add_feature(cfeature.COASTLINE.with_scale('50m'))#以下6条语句是定义地理坐标标签格式f2_ax1.set_xticks(np.arange(leftlon,rightlon+10,10), crs=ccrs.PlateCarree())f2_ax1.set_yticks(np.arange(lowerlat,upperlat+10,10), crs=ccrs.PlateCarree())lon_formatter = cticker.LongitudeFormatter()lat_formatter = cticker.LatitudeFormatter()f2_ax1.xaxis.set_major_formatter(lon_formatter)f2_ax1.yaxis.set_major_formatter(lat_formatter)f2_ax1.set_title('Chinamap_coastline', loc='center', fontsize =15)  #图标题名#读取shp文件china = shpreader.Reader('cn_province.shp').geometries()#绘制中国国界省界九段线等等f2_ax1.add_geometries(china, ccrs.PlateCarree(),facecolor='none', edgecolor='black',zorder = 1)#添加南海，实际上就是新建一个子图覆盖在之前子图的右下角f2_ax2 = fig2.add_axes([0.8-0.6/7+0.01, 0.25, 0.08, 0.13], projection = proj)f2_ax2.set_extent([105, 125, 0, 25], crs=ccrs.PlateCarree())f2_ax2.add_feature(cfeature.COASTLINE.with_scale('50m'))china = shpreader.Reader('cn_province.shp').geometries()f2_ax2.add_geometries(china, ccrs.PlateCarree(),facecolor='none', edgecolor='black',zorder = 1)#出图plt.show()
8 ~! t/ P9 k/ B3 g  M) X' S

例4：绘制世界区域地图，包含海岸线（中心经线一般为本初子午线）

（如果没有下载相应库包，请先在anaconda或cmd命令行去下载！！！代码的详细说明已备注在注释中，不再额外赘述。）

$ X9 d) D2 R9 [

• 
  
• 
  
• 
  : S& v9 _! B+ w% {1 o
• 
  2 v  e1 M3 B5 ?; X# Z  b+ y. U# s
• 
  
• 
  1 B; ^# t, |' n7 `
• 
  # A' F! V; K# P2 [0 b5 Y' B: ?
• 
  " a4 h% d$ E- S: `& p
• 
  / T! W) b" c7 A# p. I
• 
  
• 
  
• 
  
• 
  
• 
  : o8 ^5 J5 t( }) ~' v2 ~
• 
  $ `/ p# X! h; Z2 T0 W3 g) j
• 
  : R# J% a* T! y* g" \  |5 P" [
• 
  - e5 o% ]0 e% U; Z5 U# Q
• 
  
• 
  , x# x' e& g4 D4 t' Q; @
• 
  3 @3 t! n) P* e, C' L' w
• 
  
• 
  
• 
  1 A7 ]6 O' P/ k8 z& A: h2 D9 f2 K
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  8 }2 i- s7 X! g. r. |# Y
• 
  . f6 [8 ~1 m4 J, a9 u$ u
• 
  # y- e! `4 b: [$ W# H' |% W* L
• 
  
• 
  
• 
  1 C' o: i% S5 A& K# l
• 
  
• 
  $ k0 K# ^& Q& d' C( I1 B7 U0 ?
• 
  
• 
  % O& |3 X* \( r" K
• 
  
• 
  
  

#例1：绘制粗糙的全球世界图#引库importnumpy asnpimportmatplotlib.pyplot aspltimportcartopy.crs asccrsimportcartopy.feature ascfeature  #添加地图其他地理信息所用fromcartopy.mpl.gridliner importLONGITUDE_FORMATTER, LATITUDE_FORMATTER   ##添加经纬度所用importmatplotlib.ticker asmticker   ##添加经纬度所用#中文及负号处理 plt.rcParams['font.sans-serif'] = ['SimHei']  # 中文字体设置-黑体plt.rcParams['axes.unicode_minus'] = False# 解决保存图像是负号'-'显示为方块的问题#画布fig = plt.figure(figsize=(4, 4), dpi=200)#子图proj = ccrs.PlateCarree()  #默认制图时，中心经线一般为本初子午线ax = fig.subplots(1, 1, subplot_kw={'projection': proj})   #一行一列子区域#默认海岸线ax.coastlines()#添加其他地理信息特征ax.add_feature(cfeature.LAND, color='limegreen')  ####添加陆地######ax.add_feature(cfeature.COASTLINE,lw=0.1)  #####添加海岸线#########ax.add_feature(cfeature.RIVERS,lw=0.25)  #####添加河流######ax.add_feature(cfeature.LAKES)  ######添加湖泊#####ax.add_feature(cfeature.OCEAN, color='skyblue')  ######添加海洋########ax.set_title('worldmap', loc='center', fontsize =15)  #图标题名extent=[-180,180,-90,90]##经纬度范围gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True, linewidth=0.2, color='k', alpha=0.5, linestyle='--')gl.xlabels_top = False##关闭上侧坐标显示gl.ylabels_right = False##关闭右侧坐标显示gl.xformatter = LONGITUDE_FORMATTER ##坐标刻度转换为经纬度样式gl.yformatter = LATITUDE_FORMATTER gl.xlocator = mticker.FixedLocator(np.arange(extent[0], extent[1], 30))gl.ylocator = mticker.FixedLocator(np.arange(extent[2], extent[3], 30))#更改刻度字体大小gl.xlabel_style={'size':3.5}gl.ylabel_style={'size':3.5}gl.xlocator = mticker.FixedLocator(np.arange(extent[0], extent[1]+10, 30))gl.ylocator = mticker.FixedLocator(np.arange(extent[2], extent[3]+10, 30))
- B+ E" Q9 [7 x" L! p$ w+ {

/ k9 g4 M" b3 i5 o# ?* s

# P' K: x8 b. H/ c                               
登录/注册后可看大图

例5：绘制世界区域地图，包含海岸线（中国移动至地图中心）

（如果没有下载相应库包，请先在anaconda或cmd命令行去下载！！！代码的详细说明已备注在注释中，不再额外赘述。）

与例4不同的是只需要改动：

proj = ccrs.PlateCarree(central_longitude=130)  

+ v/ f6 c- k/ R0 G- Y: G. B

• 
  
• 
  
• 
  
• 
  
• 
  . L  W; t1 X9 S% f
• 
  
• 
  
• 
  & }$ X" E) X) S: q. O
• 
  
• 
  + q* \5 X8 n  E9 u3 [# O) A
• 
  
• 
  3 ]0 n" ]+ F$ d2 Q, x
• 
  " J! u3 p+ j$ P, y( F4 E
• 
  
• 
  
• 
  # m  k  E9 U4 S- s2 k
• 
  
• 
  2 E2 e4 s5 m" |% [
• 
  
• 
  $ X( w, a! s  R  ^' G$ n
• 
  
• 
  + N0 E/ @9 u5 r
• 
  5 N+ I4 _' D' M# u9 n. p
• 
  % }2 ?+ v- B9 f" p* x
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  / p: N# j) N1 X
• 
  ; o; u3 p) q( O
• 
  
• 
  
• 
  
• 
  
• 
  : c5 N" |& _/ o0 j+ |/ C7 s
• 
  
• 
  ' V" }. `. {9 k! v+ I! C
• 
  ' ]4 r& p- U+ B: N+ ?& S# ^
  

importnumpy asnpimportmatplotlib.pyplot aspltimportcartopy.crs asccrsimportcartopy.feature ascfeature  #添加地图其他地理信息所用fromcartopy.mpl.gridliner importLONGITUDE_FORMATTER, LATITUDE_FORMATTER   ##添加经纬度所用importmatplotlib.ticker asmticker   ##添加经纬度所用#中文及负号处理 plt.rcParams['font.sans-serif'] = ['SimHei']  # 中文字体设置-黑体plt.rcParams['axes.unicode_minus'] = False# 解决保存图像是负号'-'显示为方块的问题#画布fig = plt.figure(figsize=(4, 4), dpi=200)#子图proj = ccrs.PlateCarree(central_longitude=130)  #中心经线为中国ax = fig.subplots(1, 1, subplot_kw={'projection': proj})   #一行一列子区域#默认海岸线ax.coastlines()#添加其他地理信息特征ax.add_feature(cfeature.LAND, color='limegreen')  ####添加陆地######ax.add_feature(cfeature.COASTLINE,lw=0.1)  #####添加海岸线#########ax.add_feature(cfeature.RIVERS,lw=0.25)  #####添加河流######ax.add_feature(cfeature.LAKES)  ######添加湖泊#####ax.add_feature(cfeature.OCEAN, color='skyblue')  ######添加海洋########ax.set_title('worldmap', loc='center', fontsize =15)  #图标题名extent=[-180,180,-90,90]##经纬度范围gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True, linewidth=0.2, color='k', alpha=0.5, linestyle='--')gl.xlabels_top = False##关闭上侧坐标显示gl.ylabels_right = False##关闭右侧坐标显示gl.xformatter = LONGITUDE_FORMATTER ##坐标刻度转换为经纬度样式gl.yformatter = LATITUDE_FORMATTER gl.xlocator = mticker.FixedLocator(np.arange(extent[0], extent[1], 30))gl.ylocator = mticker.FixedLocator(np.arange(extent[2], extent[3], 30))#更改刻度字体大小gl.xlabel_style={'size':3.5}gl.ylabel_style={'size':3.5}gl.xlocator = mticker.FixedLocator(np.arange(extent[0], extent[1]+10, 30))gl.ylocator = mticker.FixedLocator(np.arange(extent[2], extent[3]+10, 30))

颜色表：

% H" x$ r) i* X
( R: F- z7 c8 f' z& Q