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

  
  

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

. `' @# t! j; [" i8 V5 |6 w- P                               
登录/注册后可看大图

绘制中国区域地图

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

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

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

• 
  
• 
  
• 
  8 U1 ]5 |) D6 V  V7 k/ \
• 
  - B" t9 k( I" M) U* a5 I/ m
• 
  
• 
  
• 
  - Z, Z3 ^6 Z* _$ E. _* s
• 
  9 f1 z$ F+ Y9 X0 U$ `# V- S
• 
  
• 
  
• 
  3 [. R2 i+ l( N* F# O3 s$ r5 x
• 
  
• 
  / O0 B* V) W2 h' n* B
• 
  , E5 T/ i0 S* x# U
• 
  : B4 g; d# T3 ^
• 
  3 G& V% u3 ?2 T) b
• 
  : y+ S& C4 }+ H
• 
  8 ?+ |- E# a1 l* j7 \. ^
• 
  
• 
  # k9 d' v# n( B# F
• 
  
• 
  
• 
  0 @& g" c6 d# r) i' R. t
• 
  5 Y' {/ B5 L* p5 H" l
• 
  , {$ b* Y4 i1 B
  

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

2 N( O. g, |7 f: @. w4 e& |0 T( ~
* V9 ^) Q# j+ F  k

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

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

+ I  c5 y1 c1 z. S6 @6 D3 S. g

• 
  ( y% g% m' ~& y$ O# ^9 U: C* F
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  / [, r$ V0 p! ]; s" e- A3 u+ ]
• 
  
• 
  9 X8 n& {# I) u) I/ [
• 
  # U* _$ N4 B$ j* H
• 
  
• 
  : ?/ ]8 l( G0 `5 A- b
• 
  2 K- V4 l# [- i. ~6 v  `4 f: @
• 
  : ~, P9 Y) G6 n/ D3 r- \8 V
• 
  ; ?3 H3 I8 T: F
• 
  7 j' C+ u, l' y5 K
• 
  % K8 S; U9 a3 B! Y
• 
  
  3 _3 R8 `( |7 b' `. f* Y  }

#绘制中国地图，并标记各省会城市名#引库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")
* D0 b9 `8 j& @

+ P3 ], W' D9 L3 J

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

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

• 
  
• 
  
• 
  
• 
  5 u( `9 ~$ b: b# Z# v9 l8 M- V
• 
  4 J1 q) I: ~; v& w6 a. [
• 
    D- l6 Q  F  ]# ?- c
• 
  $ P! W. Q! K  m- }- _% U
• 
  7 l: A9 C0 q% r+ ~
• 
  : d4 e, Y7 w3 Y4 q8 S
• 
  
• 
  
• 
  # v( N& A$ @# A9 }7 `" a+ H
• 
  
• 
  
• 
  " x/ v# v4 E: @+ B1 ?
• 
  - H; w  a' S# d- o1 q# u( F- ^
• 
  - p2 J; ?  ?6 d  G
• 
  
• 
  , s+ I5 C8 T5 u: m9 `
• 
  & j% \/ c6 B% b$ g3 c' t
• 
  
• 
  
• 
  
• 
  
• 
  % m# t. L3 L& r* R5 \. ?$ v
• 
  1 l* N* Z" t2 Z9 J0 Q5 u" t
• 
  
• 
  7 L- z! I+ ~6 K0 z
• 
  
• 
  & f$ Z  W, M  d* f' X
• 
  
• 
  
• 
  
• 
  + F: T, Q6 S$ e3 `( |
• 
  
• 
  ' b3 Z9 J2 V' {7 P  e% h/ f
• 
  
• 
  & H: H1 R/ y1 X" g  k  k
• 
  % ]  X5 G2 g2 N, E1 e3 T6 l
• 
  
• 
  - ]: v% Z2 |. x$ H' F
• 
  
• 
  3 _" }0 L; f8 S: e7 `9 z
• 
  
• ​
  
• 
  4 F9 v: k+ B! u7 k1 H
• 
  $ J- i* |9 ^9 c# h  U9 `
• 
  ' R: n0 a: u( `! n; g0 x; C8 d* U- p, Z
• 
  7 h: B- m0 p4 V$ p! ^; Z5 _
• 
  2 r/ a3 \" g. i: e
• 
  % w" m6 `, ?" \: r3 x
  

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

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

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

7 A% ]) R0 d) v

• 
  ) {& `0 f6 |: O( j/ V1 A' H
• 
  
• 
    e) b0 ^: A+ S' y2 @& v
• 
  
• 
  + \  B- ?6 y* j4 v3 }7 W
• 
  
• 
  
• 
  
• 
  3 U2 Y5 s7 R/ T
• 
  . ^5 v8 y3 y/ |  L8 e
• 
  4 ^5 S# t$ j, d! g9 J! ^2 s- T  @
• 
  - C& C- d2 z0 C: N* \- F$ l2 K
• 
  1 @# Y( k) e3 `1 e, @# I
• 
  4 x+ J8 l  p- s) B9 m8 Q
• 
  
• 
  1 j& I  O. I6 k' y, c
• 
  
• 
  8 L# T- _, |6 w" t' n
• 
  " e$ x& h6 c8 |2 ^+ F* Y' k
• 
  
• 
  
• 
  - p, p8 }( D: q! ]
• 
  
• 
  
• 
  : `+ y) J$ E' E" Y& t
• 
  
• 
  ) y$ j# N/ m$ u) x- x/ Z
• 
  
• 
  1 G6 |! s9 X# H$ r  A
• 
  
• 
  ; C( y' O8 w. R0 q
• 
  
• 
  
• 
  9 ?: w* ~! ~  g& O7 I, B1 x
• 
  
• 
  " A$ C4 T0 c4 Y5 Y5 e. Q. m
• 
  " L0 N* C4 ~- n% i; M
• 
  
• 
  2 Y3 m9 l2 v' f- v0 H; p9 m
• 
  0 e, ]2 m8 K  X) o. i
• 
  8 ~8 w% l/ x0 U; W. t$ {0 U
  

#例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))
. s" M$ R! t, I3 f6 [; F' W. k. L! w

& d/ L% o0 T; T. _; V& J                               
登录/注册后可看大图

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

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

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

proj = ccrs.PlateCarree(central_longitude=130)  

$ Q! K' J) u$ _9 H  g

• 
  8 o3 V+ q3 O8 w3 ^3 \
• 
    k8 q( `# j3 \5 N$ @/ E* t2 g
• 
  ( B- R- g1 S3 \; R* f8 F
• 
  , z. m$ l* m- @- v  m
• 
  
• 
  8 a" i  c2 a  p( u# _
• 
  
• 
  , j! @; h  s* p, N
• 
  3 ?2 A. H# }: f- R* N
• 
  
• 
  4 s4 w: r, M* D+ a
• 
  
• 
  # E. |4 F7 _7 ~3 W$ |
• 
  
• 
  ( O% P; _2 Y% k% C# z! X& P; z
• 
  
• 
  9 b# ]  ?( V. @/ |
• 
  
• 
  . u6 L% C* |, z; p* e+ `
• 
  2 }; `, l% \, G( H: x! H. g
• 
  * n; q* {; W1 O
• 
  : r( _) J! i: b/ z
• 
  / R) k: g& N3 [6 J
• 
  4 j$ r" Z, n: S0 Z& z( o5 q4 Z2 r- T
• 
  
• 
  
• 
  0 e: F  g5 r" x' b  _' X( V
• 
  2 M& [8 Z& \+ h- p* a8 w* m1 f& V
• 
  ( J" E/ S8 L4 C
• 
  7 `- v$ |- n2 N. S
• 
  4 B$ L' A. |( i6 s! }& _' y$ U
• 
  9 g7 l3 F" o/ W0 n2 L! ~6 L8 m# y
• 
  5 w+ G- o4 o& o+ C0 w8 l
• 
  
• 
  2 h0 T! c8 Y$ \  o2 Z
• 
  
• 
  
• 
  
• 
  
  

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))
. I. r3 n/ L4 k( @

" e: m# Y3 _# I+ ^' ]; p

颜色表：

" L- m7 x2 s' o; U2 t# Y