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

  
  ! x+ H0 P. {; {0 V8 |

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

& _1 h2 K9 q: ^, H. _6 O- W$ W                               
登录/注册后可看大图

绘制中国区域地图

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

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

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

• 
  ! j$ S) Z- O% ]/ e5 @* r4 B: f/ s6 e
• 
  2 \, Y, @) \3 E5 i6 `/ e/ a
• 
  
• 
  
• 
  + s9 I6 u- U4 ?( Z' [% k5 g
• 
  : |: |5 T' w  E0 j# i1 ?7 f
• 
  
• 
  
• 
  
• 
  : {/ A& s7 P' s8 Q
• 
  - u; D' ^! x$ j2 o6 |" C# I
• 
    S4 x: l! J- A& h+ }; c3 y- q
• 
  / q% r- _; H- n7 X# A
• 
  . ?$ }5 h* d# x; W4 y3 y3 j. G0 y
• 
  
• 
  ! j6 h0 a% E7 z! y* N, W, S
• 
  2 y  F  r; a% _1 _' k
• 
  
• 
  
• 
  : a) n% z$ u" G4 b& ]- z6 u9 J
• 
  
• 
  
• 
  
• 
  
• 
  3 p' `$ t/ w  Z* @* i0 Q
  

#绘制中国区域基本地图轮廓，包含省界#引库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")

/ q. S, t% s0 ^/ ]! k- U% j

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

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

4 ?' |# C) m# E! Q. ?. l' B

• 
  6 i4 Q7 I( }( ]$ }, j
• 
  
• 
  ) p3 V1 F8 J4 P' w3 B! a
• 
  6 N$ u! J; u  r6 T: k1 P
• 
  
• 
  4 t3 I& a+ f6 n! |* F, {8 f$ Z
• 
  3 t4 c8 q2 T2 i! q5 R7 ?
• 
  * d- t* e' T8 c2 H$ V5 W: T9 w
• 
  % z$ z4 P( C$ e3 H) t3 x% f
• 
  
• 
  
• 
  
• 
  
• 
  + d6 A- e, ~+ L9 A$ {
• 
  , W5 H5 @( G8 f- P9 ~6 ?: w4 n
• 
  
• 
  
• 
  
  . ~$ e( @, p1 D* a& P+ e

#绘制中国地图，并标记各省会城市名#引库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")
" C' r" x6 \4 l- n  Z

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

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

) a+ h/ Y# e6 y/ k/ z" c

• 
  
• 
  
• 
  * \" {1 |+ e; n% C' A' z
• 
  ! ^+ ^9 N$ g, q* C& h
• 
  2 s  [9 U" r* U% z+ \
• 
  
• 
  & E& ~6 h4 r0 y6 r" \: q
• 
  
• 
  ' x7 m$ Z2 u! F2 S! b" |) E
• 
  3 p8 _! E$ |1 L6 V
• 
  % f$ b7 f8 `" t# L4 C
• 
  
• 
  
• 
  
• 
  
• 
  5 D6 u9 k) U9 c8 }3 m1 A
• 
  
• 
  6 O2 f, t8 k$ o! {+ r
• 
  # w6 w5 Y4 t; a4 k
• 
  3 O" Y1 a" X( J$ `- P  Q
• 
  
• 
  
• 
  
• 
  
• 
  6 l+ [! m9 h. ]% n4 L$ Q  ]1 M
• 
  
• 
  
• 
  7 T: u- K0 h6 `/ u# n% e; q
• 
  , V5 x) o+ I+ I8 k* x
• 
  " m! ?7 ^  ^$ q7 }  c  P7 @2 G
• 
  
• 
  
• 
  8 [1 V0 i( r* I  B8 ~
• 
  
• 
  $ K6 N2 L2 M6 T+ q
• 
  
• 
  3 |9 n; T* S: Y7 ?* A, X
• 
  8 E7 M; H, v$ j+ v) o4 [
• 
  : p( l; ?  E% R1 I9 D
• 
  " m' w5 g# i( Y
• 
  * k4 ]; K' g4 l
• 
  
• 
  # ?. h0 o8 p7 e8 Q: g9 P% [9 G, o
• 
  
• ​
  
• 
  
• 
  
• 
  ! U/ N2 B+ m& R; ?/ I% w* n9 ^5 i
• 
  
• 
  
• 
  
  

##绘制有海岸线的中国地图#引库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()

5 n8 {. Y. [! B6 X( V1 q

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

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

4 m2 Q: a! h/ I' Z' I* q

• 
  . R6 e) Z1 r" [: n9 g0 Z8 b
• 
  
• 
  7 h& S; @, m; M6 ]4 b; c4 A
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  . m- r! e+ o# B3 C) v
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  ( h3 h4 v$ W" a" @
• 
  ' _+ U8 H, x( W$ B% @2 {- q
• 
    o3 x0 u' P3 X# t  a2 s
• 
  4 q* i- _: s+ H2 F" `
• 
  ) n! x; P# {2 X
• 
  * ~: q- x% i* \: z) k: q
• 
  
• 
  
• 
  
• 
  
• 
  % k2 \  u- A3 Z) i: I& q8 A8 o
• 
  
• 
  # o4 O( _- w0 q3 E: s% w' H  _
• 
  
• 
  5 F3 X9 ?$ p4 Q( b8 n7 {  @
• 
  
• 
  & Q5 X( S8 ]! t
• 
  8 m( ~3 O0 g1 ]* k! ^4 Y0 ?
• 
  
• 
  
• 
  
• 
  % G2 w7 @  Z0 B
• 
  
• 
  : o8 h. r- f9 M7 q8 `
• 
  ! f  l% ^: Z$ 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))
" p5 ]& Y. u7 i. s/ I

. a+ ~9 B+ G) a                               
登录/注册后可看大图

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

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

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

proj = ccrs.PlateCarree(central_longitude=130)  

3 W4 K7 e; A& e# K

• 
  9 t9 b1 j6 a3 y' A4 m( Q9 u* ~$ J
• 
  - l( A5 I9 |. a% d5 Z
• 
  
• 
  9 N/ m# Z* \0 i8 ^( m1 n
• 
  
• 
  " _( B) P9 w5 Y' G1 x$ a
• 
  5 }; C1 e+ M" H4 a. Y' [. C2 |
• 
  " I% U6 |- D1 E, {% V
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  2 a+ J2 R, l% G, l; g
• 
  
• 
  4 x( q$ u$ o5 W2 q; R
• 
  
• 
  ; S9 I4 m8 R) l( b, l- h
• 
  ' M: h' O- Z. A- a
• 
  1 {, J6 T0 \+ f! H. J
• 
  * p1 \2 X# ^# _8 A/ c0 l  B
• 
  
• 
  ( \& P. Y8 S7 Q' R
• 
  
• 
  ' y: f9 e8 s- I, @- _8 X3 J1 C
• 
  
• 
  % ]$ t) @: z4 }. U& D) v8 K4 b& S1 A" o
• 
  
• 
  * s8 R& Y1 V0 N! c! T" h4 Y3 E
• 
  
• 
  . _  Q" S- m$ |* n: ^
• 
  . j$ k  v, O. W
• 
  
• 
  7 J; [- t( j$ ?. \$ e8 e
• 
  % p+ i; {! v9 Y( R% c; i4 @  m" K
• 
  * |5 \' d  }, b1 E- }' o; _4 \! S0 u
• 
  
  ; v* L$ d+ m* \

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))
2 t/ N- H, @9 E0 K, w. S. g

6 S/ G/ Z' J& L- A  V- b. [

颜色表：

$ d, Z3 a$ {6 w9 q% w, i0 w6 `
. ?2 ^# ]! c6 d8 _