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

  
  

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

                               
登录/注册后可看大图

绘制中国区域地图

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

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

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

• 
  
• 
  # {( `; k% _* h- j9 Q
• 
  " Q/ D' f% M1 P( d, u
• 
  ; q4 X1 ^: E/ c/ h
• 
  
• 
  ' F: [! a' B" R% }7 E
• 
  0 g1 i) V" T+ c3 H7 S- x
• 
  
• 
  
• 
  
• 
  
• 
  3 G0 Y' F. a: T
• 
  8 h+ u3 q' {6 z8 Q: D
• 
  1 h/ Y' G$ `- w. e( F/ \# c
• 
  
• 
  
• 
  
• 
  * G& h, V+ S4 J4 k# S6 t% E
• 
    g2 I3 ?, r" d* W, N! i( l/ f
• 
  
• 
  
• 
  ; F1 Z3 d, J% e0 }# Q6 Y& q1 \9 t
• 
  
• 
  5 i6 B! O$ `# x% }: G
• 
  % M1 M. Y2 N7 P2 i  y
  

#绘制中国区域基本地图轮廓，包含省界#引库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")
) d2 b% }$ a% ^5 t& D

+ L% t- y0 j3 [8 k

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

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

% ?0 Y; u% j1 w7 C

• 
  
• 
  ) s# b' z8 T8 [2 W; a
• 
  9 k+ u4 ~3 z5 h( O$ J
• 
  
• 
  
• 
  
• 
  " v+ g1 X4 ^( q  S2 U
• 
  
• 
  
• 
  
• 
  
• 
  8 `" z' M5 t0 }  u1 W0 F" I
• 
  6 P: ^. h3 u1 m! h$ I5 r. p4 x
• 
  
• 
  4 d' l7 Q7 F7 M6 J
• 
  1 O9 B) x! B" i: F" @$ b5 _
• 
  $ C: n& v5 X. ~$ g: F4 Q
• 
  
  

#绘制中国地图，并标记各省会城市名#引库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")

1 @- c% v1 a* H8 Z5 C) V

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

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

  V- B- a" Q/ i0 D$ u" q4 M

• 
  + y# D; Z* g9 t5 s1 ]- R
• 
  
• 
  7 P. C; `/ V1 t& n0 i" n) j
• 
  
• 
  
• 
  
• 
  + @% t" K' ^7 v3 N  V0 t
• 
  
• 
  
• 
  
• 
  
• 
  , A, O  l1 E) D! u" _- D
• 
  3 [$ m# V4 P! n/ R; O
• 
  3 e# I2 J6 R, |+ Y" v
• 
  
• 
  & b4 F+ S6 |$ M/ I0 p' T0 c
• 
  
• 
    u" O) ^- C6 o9 d
• 
  
• 
  1 L. H" Q$ q' H& o2 v3 T# ~: b6 f7 {
• 
  # ^* f8 Q/ ~' `
• 
  ) r5 R& z9 _, |: g" M
• 
  0 D& u3 u' g+ W
• 
  7 U$ f- O1 J9 v$ C! W& ]$ p
• 
  & L6 ]' z7 R! F( u& L
• 
  
• 
  : B+ J& h# J& c: ?
• 
  
• 
  
• 
  
• 
  : {/ W5 S( y. x. @5 v3 R
• 
  % a" W9 W: X6 b7 }7 P6 t' j
• 
  " P6 Q" E( f1 w
• 
  " z1 U6 N* [; Q' Z- I0 z* O7 W' N
• 
  
• 
  - V# D* C  y6 B. P" H
• 
  9 d* e5 q( r" b( G' _
• 
  
• 
  5 ^/ S& U# y; t/ i) Y- B2 m9 M
• 
  
• 
  
• 
  / C! }  S) t' u/ Y7 D7 q
• 
  
• 
  3 q5 e, [  b0 A; B1 M% S
• ​
  
• 
  & f( F/ N- Q4 j! y
• 
  9 y- x9 z" ^7 T6 ~6 c2 q. M
• 
  
• 
  6 b* p$ Z' A  K: V) E
• 
  
• 
  
  

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

3 k2 O2 ?8 p1 F) M$ o  E

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

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

• 
  
• 
  
• 
  , R5 K. h2 f* }! F3 B+ P1 Q" V) j4 ^. }
• 
  
• 
  
• 
  1 I9 }1 Y* F0 S6 v3 v; D6 D
• 
  ) l) h$ D& S0 k) o. B: K! F$ _
• 
  4 }& q" J4 w! S& Y( Z+ t' y# A
• 
  8 a; v- h+ F; o
• 
  
• 
  ! C/ }* Q5 ]$ b
• 
  
• 
  ( F$ U" P0 \8 S, I
• 
  
• 
  % N! g; I9 X; ~0 f# q' B, S
• 
  : g# k8 @" d' s: D
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  ' W9 ^3 f) q' s5 v5 K( s
• 
  
• 
  
• 
  
• 
  $ T/ u; G! u' B  S! D1 L
• 
  
• 
  
• 
  
• 
  
• 
  1 M& O, K( n! E$ N$ q8 h$ M- s$ Y
• 
  
• 
  
• 
  
• 
  1 O, t( [$ s7 y# A; G4 n; M
• 
  
• 
  # n4 b! U' o% Q3 P1 U7 v
• 
  # f; K2 {$ d: |
• 
  ' Z% i6 s0 X8 P
• 
  
• 
  6 ?( J/ M. B% `; J
  

#例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))
1 A  L& \- F# c2 B& ?- i

                               
登录/注册后可看大图

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

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

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

proj = ccrs.PlateCarree(central_longitude=130)  

2 T, I! s0 t+ @

• 
  1 [# i1 q7 s$ T2 C- t% m
• 
  
• 
  
• 
  
• 
  
• 
  # c2 u7 V. ^/ I
• 
  
• 
  
• 
  / w' A, R. c% Z9 _& o. U2 k
• 
  . @9 r7 l% B" ]/ a3 ~
• 
  
• 
  $ i9 [* s0 N7 [9 t0 h* F5 o
• 
  
• 
  
• 
  
• 
  * j0 n' F' t2 f! v3 k5 a/ e
• 
  
• 
  $ x9 T9 v' F, m8 y  }3 y9 [% I0 w1 u  B
• 
  6 Q( H! V- c% l- m( }
• 
    L- T+ c# ~9 ?1 g" i
• 
  
• 
  : ~1 T5 I3 v; c' S" o, J5 Q# ^
• 
  - i& d; I* t- v. d. t3 d/ h' i
• 
  4 O$ {0 o* K3 X, F
• 
  
• 
  
• 
  : z% X8 s" _$ j- s- D/ X# _
• 
  
• 
  1 R. `+ G4 P, {+ `( o, T9 l# S
• 
  
• 
  ; Z; R) `3 m/ h6 n
• 
  6 _& @' q# y* Z
• 
  
• 
  1 d( T1 m) N$ Y' e
• 
  
• 
  
• 
  9 p/ @5 `7 d/ _1 M
• 
  . W, h  Z7 p9 u& v3 k% b& {
• 
  
  ( R4 ?2 v' V" h9 ]: P

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))
* p  E6 ~  ?# k* A" m5 k

" ~6 p2 ^/ Y3 g9 r1 T) i5 F1 t
; I- l( J# L* q5 m, X

颜色表：

8 k+ k$ r- H4 T* L