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

  
  ' ~3 c: f8 [# B( q: F

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

4 H$ M( Z& y; I0 z                               
登录/注册后可看大图

绘制中国区域地图

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

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

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

• 
  + Z# Z/ g9 M7 F' h4 \# E" {7 ^
• 
  
• 
  & A% |2 K0 F' l
• 
  3 p) Z' Z5 `. Z
• 
  
• 
  / N% J$ y' I* F+ W8 x4 g, c2 K+ s
• 
  
• 
  
• 
  0 |6 }7 d$ ], G2 c; v) N: p
• 
  0 D4 }2 z! U2 G, v5 @. k" ?' ~
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  9 [) g2 p, \0 V
• 
  / v, Q" W$ ^8 I" _4 o' ]" U
• 
  
• 
  
• 
  
• 
  / a8 r0 @2 ^! j( u
• 
  , C& }/ [) n: R& ]/ v% F& ?
• 
  
• 
  4 @$ p# k! U, I" u
  + s+ y! O( i- q9 M$ I6 {5 m9 A

#绘制中国区域基本地图轮廓，包含省界#引库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")
; b1 R: }- s+ Q3 p8 c' w% K( G( b; d

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

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

$ v! B; [/ Z% O. T6 Q2 k+ g0 J

• 
  
• 
  
• 
  + {; Z$ w, K0 ?3 M% i' K
• 
  ( s  g" G% M; \) b* [& e7 E
• 
  ( Q, c0 t) ~  ^" w
• 
  
• 
  
• 
  1 W3 J9 N2 b( ^1 g8 W  e
• 
  
• 
  
• 
  ) h6 B0 `$ J0 c* i4 C( C
• 
  5 B4 z) r) j* m/ v8 v
• 
    s- p9 j3 t8 x  C
• 
  + [4 |" a* P5 i* Y  B' L
• 
  
• 
  * r! c$ p' Q; U0 z+ [
• 
  
• 
  - I1 B7 L0 d1 P) y
  + q' e9 b7 r: d9 i5 z( o

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

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

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

• 
  
• 
  4 e. V* {0 k; N+ S( `3 H
• 
  * Z( C- {1 C4 V% F6 i4 c
• 
  
• 
  
• 
  
• 
  ) k( J, M" T, |0 Q" r/ ]
• 
  / W* K1 n1 X- x
• 
  9 L- T0 r3 Y% S- d/ w9 e
• 
  
• 
  * G/ {2 F# C2 S5 s) H
• 
  % \% i& g8 o" [: b
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  , O, f( v4 M  |/ M
• 
  
• 
  
• 
  1 o3 I9 x3 u8 C9 Z# A
• 
  
• 
  ) h' P6 z' w2 H7 I3 J0 k; T
• 
  
• 
  ' m+ M( S' V9 R; k. T2 p! f# D+ s
• 
  0 [- Q' H# s  o7 P# ^' v2 K
• 
  ! P7 Y9 a3 C7 t) g
• 
  ) O0 a; T4 R0 c" W- w2 X* T
• 
  + G9 g$ T1 Q1 D$ W% \" y7 z
• 
  
• 
  
• 
  
• 
  ' B0 `  ?  M2 v; a$ r; T
• 
  " `" u) |. r4 v
• 
  
• 
  
• 
  
• 
    z: E2 Y/ z* z2 K
• 
  
• 
  
• 
  % w5 u0 `5 [; A6 r! U& g6 f
• 
  - T) D1 K1 s  ]- a3 e& o! Z5 |
• 
  0 \6 m, y% ~  N1 [0 r2 T
• 
  ( e2 `4 C$ E/ V% H
• ​
  5 W/ a6 ?: z$ v0 {
• 
  
• 
  4 D7 f! S4 m" g; Y$ _" o
• 
  8 ], ?; U2 z  c, P' N
• 
  
• 
  ' I8 G! G% C- i$ R& Z$ N# q, {
• 
  + p  P2 ?6 m: N* n
  

##绘制有海岸线的中国地图#引库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()
- G* B5 p" Y# `6 N5 R

. @' u' \4 `, h5 v3 H
( P8 I" [0 @1 W& u

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

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

3 `9 N4 V" X: r) M& q

• 
  
• 
  
• 
  ( Y% {4 y& \7 V# S2 M
• 
  
• 
  2 c0 u3 V8 d+ ]8 \
• 
  
• 
  % G* H& i- I* q7 \* l/ c- r) I9 {, _
• 
  + M8 @% j/ j' c5 r% t
• 
  
• 
  
• 
  6 t/ W- A* T, z7 J
• 
  
• 
  
• 
  / \5 M0 J8 ~, F% R8 Q) P1 b
• 
  
• 
  ! R; u5 u- g$ o) o
• 
  
• 
  9 d1 \7 X. ^; J9 l
• 
  
• 
  
• 
  
• 
  ) u$ k7 h, `4 }0 G" ~
• 
  % c$ g- H; H: E! e( G
• 
  
• 
  6 Z7 h$ ]7 A: Y9 x, S, t5 T/ ]
• 
  
• 
  9 K0 l; b  D5 T  f" }2 D
• 
  5 J+ ]( I$ y: U6 U9 K& l3 t; R
• 
  
• 
  
• 
  . @1 l( G# {3 ?% O& B; ~4 u5 W9 c# d, g
• 
  
• 
  
• 
  4 U# d5 s* F! B5 }, T2 {& N
• 
  8 K. Y  [# }% U% V& C/ ?7 Q$ u: J5 p
• 
  + B( H, G! x" c  m% D$ {/ l
• 
  
• 
  ; Q. \; T2 F4 `% F5 d$ c
• 
  / }5 Q3 |  L) T; l+ ?9 [
• 
  # x( \- w! ^/ J' L! k, l- V
• 
  
  

#例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))
- M# N) v2 j4 P6 p7 l" m

  E2 s* G/ r/ a7 F: I9 F8 T

                               
登录/注册后可看大图

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

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

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

proj = ccrs.PlateCarree(central_longitude=130)  

• 
  
• 
    |0 s$ G1 i3 I2 I* ]( A
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  % x, G& H% O$ {  Z7 z  ^
• 
  * Y5 E/ C. ^2 N1 O) y. t% d
• 
  
• 
  
• 
  
• 
  " |1 j5 {" S- |- t2 W0 T1 s! f
• 
  
• 
  8 N) p3 U. l9 `. N" [1 c6 {
• 
  
• 
  ( w! g) B# i2 ]. z2 c
• 
  9 E( P* @9 }2 ]5 h' ^! t. g. ~
• 
  
• 
  ; r- T. ]- P" ]
• 
  / o* ^/ Q! w3 k
• 
  
• 
  / ]% z' j9 ^# R% C% \
• 
  ; ~' F- [. @& o
• 
  
• 
  
• 
  
• 
  & ?: p' J( Q/ k( c0 E+ t9 @/ V
• 
  8 J* K9 }7 b, ]$ G! n& [
• 
  7 |3 x! U: Z) c& i3 m4 X, q
• 
  0 U0 c$ ~) N1 q' H. D; U
• 
  # y) ~4 l2 s) y/ K5 o# b
• 
  8 O, Y# o. i2 p! Z3 o5 s% f8 y
• 
  
• 
  
• 
  
• 
  $ I. j$ F6 R8 Q& S6 b
• 
  . e' v. E4 O# o9 r( V+ w
• 
  
  

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))
; f- P# |* g  L" m

# R+ W5 n) h& Z! {

颜色表：

9 H1 M5 A& \1 y& u8 D