【气候软件】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("···")

  
  4 n, d' f' R4 N

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

                               
登录/注册后可看大图

绘制中国区域地图

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

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

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

• 
  + r" w, L! ^0 k+ U3 F
• 
  % ]5 B& E" G4 u( x' s4 p
• 
  1 o9 J  V) ^2 A
• 
  : K" f; H  x3 q* _
• 
  
• 
  
• 
  
• 
  
• 
  % S4 v/ ~' w% a) ?' N
• 
  
• 
  
• 
  
• 
  1 q* |& z+ X7 m# L
• 
  
• 
  ( {4 f6 b" v5 S
• 
  . t; T8 b1 J' J1 o
• 
  
• 
  8 J6 a6 [) {* ^
• 
  5 U+ j/ D! M9 z  ^. }
• 
  , w. j$ z- I3 T8 f1 ~
• 
  
• 
  4 b5 t" c' g, O
• 
  
• 
  
• 
  
  8 P0 [: e( {; @/ b- |8 f

#绘制中国区域基本地图轮廓，包含省界#引库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 d9 V! {1 |" e6 x1 \/ W
6 [5 y- w/ F% t7 F

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

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

• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  - c- y5 T# W  ~% w+ E1 Q
• 
  + k& e3 c, s4 K, U' m# B( c
• 
  3 j% k0 e2 N6 j( g
• 
  0 _4 y1 P6 N# T0 k
• 
  
• 
    Z' r- l' x3 _5 T
• 
  ; Z0 @& |3 O3 P" l! }3 H$ i0 m
• 
  $ q4 X! e, r) y3 ]# c1 [+ _9 G4 q
• 
  
• 
  0 ?9 a' x2 @, a1 \' l
• 
  0 y0 A0 T. S; x- P$ C6 t7 j1 N
  $ p& l9 M) Q+ `5 z0 S

#绘制中国地图，并标记各省会城市名#引库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")
' x* J9 `  Q& m6 {: a

! ]. Z, I7 P' ^4 k2 e- H; p) O- Z# V
7 G6 F' U. m& c7 g( E" l+ ^

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

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

• 
  4 j! h- z5 C. ~
• 
  * T0 r( u# y3 |# h5 r! O8 d* \
• 
  $ o" w0 ~0 f$ B+ d9 _
• 
  2 k. n  N" o4 {9 c1 P' v9 f+ W6 Y) I
• 
  6 J; g/ o6 p5 p; U1 S' Z1 O# `3 r
• 
  
• 
  
• 
  - {* v* m$ C  Y; J
• 
  
• 
  ' R4 j3 @0 i2 D: n" M2 u
• 
  
• 
  
• 
  
• 
  
• 
  & N, @6 W2 b# q: N+ b
• 
  
• 
  
• 
  
• 
  # w( h% s. X) \
• 
  2 E8 d% O; D! F3 }
• 
  
• 
  
• 
  4 W; P6 n1 J1 w/ {! l0 G) N0 @% ?# T
• 
  9 v) r9 ?2 t( O/ W4 Z# z7 U
• 
  
• 
  % S# k  V0 K6 R+ {+ t
• 
  ( ]8 T! L/ Z; A4 c
• 
  
• 
  
• 
  
• 
  6 A8 z: ]' h/ y" k
• 
  & A) ~1 m( p! J8 l5 v
• 
  
• 
  
• 
  $ L) ?- A9 ^$ E. g# Q
• 
  
• 
  
• 
  5 U% r% w  C8 Y! M6 {) ]% f6 e/ K
• 
  
• 
  ) c, a+ g% j' w5 c
• 
  4 ^/ i. y4 i" g7 r0 J7 S
• 
  
• 
  8 ]0 `! ~& p+ f$ }
• 
  
• ​
  , O/ N3 `, g' r, b3 X
• 
  ) l! p9 S$ K% F+ ]2 q$ {
• 
  , c  }5 ]7 d7 U9 [0 V# B4 V) o$ s
• 
  
• 
  # f; I! V. c/ u1 G
• 
  
• 
  
  + J8 v. H+ N" P  ?2 k' y

##绘制有海岸线的中国地图#引库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()
9 R, L0 L' y$ @4 H4 _. j& k) C

+ a" t; L+ m9 a
0 N2 M2 x* W$ v: h) ?' U, \' B% _3 O

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

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

; ^! e& `% I8 N$ K

• 
  8 H' c( B' L7 |
• 
  ! Y4 @, o+ b+ T& a0 G. @3 E
• 
  
• 
  & l, v' m4 Z7 P4 Q  W+ d" K
• 
  . }$ ^! @( v* K* A8 @& X
• 
  4 J( {+ L( O1 d( N
• 
  
• 
  
• 
  
• 
  
• 
  ( c4 ~: I0 K2 ^4 _& y* ^1 ?; N
• 
  
• 
  5 o! R2 n% k3 K- m. a
• 
  
• 
  - l1 z. O6 \, g3 K. Z0 U
• 
  
• 
  
• 
  4 W) E$ j( U0 I
• 
  
• 
  ! O  f! ^8 Y" v
• 
  . u0 P+ n! W8 L  v
• 
  
• 
    T4 f' f! n8 I  H
• 
  - t% e" a; L- Q, E* R
• 
  + p& }; K$ u" U" V# E: r' x, L
• 
  * S- ?  M8 g5 {1 ?9 S
• 
  
• 
  ' c" S# s  ~3 M$ p
• 
  
• 
  $ Z6 T) |$ N( s, |
• 
  ( \9 @0 c  s7 z$ {. b, b3 M* Y. o
• 
  1 f; h, n/ \0 |+ R, m7 j
• 
  , l( c% B# m1 O1 i1 s1 y
• 
  ( g  T2 j# F8 a# _
• 
  / o( Y* Z. V2 D
• 
  & W" X  g6 }- `
• 
  & K' G$ n  ]# Q" u
• 
  ' n9 i- n9 |1 i0 p$ `6 E
• 
  
• 
  6 D  r6 W3 Q+ H/ ^4 Q
• 
  1 v3 E- W3 x& y) {( [5 n! x
  

#例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))
; @& V4 M# @+ M! u* M6 |

5 a) ]- k2 U  ]- w                               
登录/注册后可看大图

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

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

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

proj = ccrs.PlateCarree(central_longitude=130)  

6 O2 L$ Q! }6 x9 B% J

• 
  
• 
  & _" G- n* V% G
• 
  8 T: ?3 r( [$ B2 m8 j* J& p. P
• 
  + d. i8 O# X3 H
• 
  
• 
  
• 
  $ R8 S, e3 d6 `+ J" U
• 
  
• 
  1 c0 n* ]9 U/ e  ^* @, }' X. k
• 
  ! Q( i3 U+ o; I7 w5 z  P5 R
• 
  ' ]1 q# m6 _" v0 X/ N
• 
  
• 
  
• 
  " I8 Q, ]/ _/ O; H' f
• 
  
• 
  " s2 c4 O2 ?; N1 H& F4 Y
• 
  9 m( m8 x" U# K- L
• 
  9 T# |* N8 R% r' S/ e0 h
• 
  
• 
  ) k0 t/ j& ~' W# H: q
• 
  # q- u, [9 ^8 u; q1 L- w
• 
  
• 
  
• 
  
• 
  4 E0 M8 H3 A2 N8 ]" D( k+ C8 R
• 
  + i$ }/ ~! Q, x6 _, Z) y
• 
  ; o' E' Q' P: H
• 
  - [0 }; d( O) h
• 
    U" s* p3 n# a# S  T- s+ T6 I8 L5 [
• 
  
• 
  8 D3 e- c; V6 e- O; d
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  1 c3 {5 \) U/ X5 P- `' u4 c
• 
  
• 
  3 ~3 @% V. W- Q5 C$ |
  

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))
. ?" Q3 i) y: P) U, _/ c% N

* t) h: X6 ?3 j2 K( u% A
' \4 i9 K" j; P

颜色表：

. @  n9 P7 i2 M$ C' W/ m
, G# t, P" |* z. A1 r! c! e" }