【气候软件】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命令行去下载！！！代码的详细说明已备注在注释中，不再额外赘述。）

• 
  
• 
  3 @) A& r! {. B% W7 I
• 
  
• 
  
• 
  
• 
  
• 
  ) P! l4 V* g8 J' q
• 
  
• 
  
• 
  * j1 R# B6 C  G+ \; F
• 
  
• 
  6 I  A/ q/ E) L
• 
  
• 
  
• 
  ) P3 l: b- q0 x  ]6 t3 f
• 
  3 ^1 X  F; D1 B& c' n
• 
  0 i; r7 G/ r/ ?8 l
• 
  ! Z. \; [5 `+ a! G  n
• 
  ! g+ }! x" ]& P
• 
  
• 
  
• 
  * \9 c$ ^  k1 y
• 
  
• 
  ) _- e" s  V0 Q) d8 `9 F2 W& ]
• 
  
  

#绘制中国区域基本地图轮廓，包含省界#引库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")
  N8 Q+ B; y$ x9 Q/ W; ^" J3 r$ p

' x# X4 z" B7 B, B+ G) C

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

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

• 
  * v6 D( \! _. ^  i
• 
  6 u. K- S; f' t- g% j# X
• 
  
• 
  , ^0 p9 D% L! a; d3 _! Z5 [8 O
• 
  ) _; R( d3 t9 e0 Q
• 
  * a: V( u/ }. q% q* A
• 
  
• 
  6 e! |$ C: B5 m9 H$ N! A. A  c3 B
• 
  4 s; ^  ^- x6 y5 y0 o; y" E
• 
  
• 
  
• 
  ! M8 B) s3 Z! ?/ Y7 E
• 
  + Q; _3 T, {( r% F
• 
  3 j! A: J7 n2 u
• 
  ' ]/ W+ K" T0 h
• 
  
• 
  
• 
  
  # ~8 M6 K" J* n, {  U( \2 K. 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")

  t* _( K2 J: j) \' o0 u! L

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

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

$ S8 k; W8 I* k1 t. L" H

• 
  
• 
  
• 
  , e$ V0 p: h+ M) K4 _
• 
  
• 
  
• 
  
• 
  8 X  v5 U5 P7 p
• 
  ; u; ^1 o3 B& l8 |
• 
  : u# L3 D* E; v0 {( Y1 I; j' ^
• 
  4 K0 x6 \2 x4 J) W6 y0 O, F9 `* X* N
• 
  
• 
  
• 
  # R% d( i$ `$ E5 @: C
• 
  : v) q: p2 e1 q  t* b
• 
  
• 
  5 V6 s, R1 o8 y0 X$ k+ w. x
• 
  
• 
  " X  b# w; K1 W4 h' `) Y- f
• 
  
• 
  ( e  i0 X9 p$ {5 ?# S8 M: o* k+ K
• 
  
• 
  1 C2 i7 ^! |1 r- n! k: ~/ i
• 
  8 w  _/ c2 I0 A# n& S. x
• 
  ; S. @' Y9 a8 m5 w5 c; z! m( {
• 
  
• 
  
• 
  & @9 E2 s3 D0 W% w
• 
  ( h1 H# w  X5 d9 ~* L% p7 J% q+ j$ _2 q
• 
  : L8 y3 l- Z- g! f/ d
• 
  
• 
  
• 
  ) ]% q5 p* \$ K& e
• 
  
• 
  
• 
  
• 
  4 a$ q+ o; Z. q' h
• 
  1 d7 T. B0 T5 f2 x% D2 q
• 
  
• 
  # Z2 I% D& X# N# T9 X8 a
• 
  
• 
  
• 
  
• 
  
• 
  
• ​
  : A  @( }& ~: d6 a& h( f5 X
• 
  
• 
  
• 
  
• 
  
• 
  0 [( D- Y1 u. }5 T' E" P! d
• 
  
  ! D& }8 |  m+ u3 P

##绘制有海岸线的中国地图#引库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命令行去下载！！！代码的详细说明已备注在注释中，不再额外赘述。）

• 
  / U! ]9 n) |. t: ?" w1 w$ ]
• 
  ! @; _5 N: P0 q6 Y+ \+ s( ]
• 
  9 j- o% u1 o# H  ]9 X) Z2 w( l$ c
• 
  & r" ~$ j* P: U* W9 q% w& ^+ A
• 
  
• 
  
• 
  3 M. t, x/ Y% U
• 
  
• 
  # l# W" d5 B' _# i
• 
  : o$ r7 s- C9 k- v; b9 o
• 
  7 L' q& F$ u2 ~: K5 H
• 
  
• 
  9 S. F# L! C- |/ G
• 
  
• 
  . c2 ?5 K# N8 w. Z# j' z
• 
  
• 
  
• 
  4 _7 g/ w, ~2 S1 J, s8 Z# I: }4 a; W
• 
  . {- b. \9 }+ q" Q- M: [* a: _
• 
  
• 
  9 o+ i' g: h' |8 t% B
• 
  
• 
  ' ]4 d! N) l( [5 l7 F, c) U
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  % F* r6 w, U' W& p5 _4 j- w8 s. d
• 
  
• 
  
• 
  . a, _) p! ]$ n  q
• 
  * O7 l+ k# ^; @# s: @0 L
• 
  
• 
  
• 
    o6 S! e8 A- Y
• 
  
• 
  
• 
  # K/ C" P% b$ w1 x; \
• 
  + w6 \1 z7 R5 a" G, F
  

#例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))

" }; r5 t1 p& i" G9 h2 {5 A+ e6 a

                               
登录/注册后可看大图

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

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

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

proj = ccrs.PlateCarree(central_longitude=130)  

• 
  
• 
  % p. T% z3 O( [% m3 H, }, `. E
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  7 k" _& ^: h# G9 f
• 
  
• 
  ! l6 f* G6 n& [( h' C! Z2 l) f
• 
  
• 
  
• 
  5 z1 |$ T! l! W2 J9 {$ h
• 
  9 n1 P& ^% H; B0 o) h  o) L; r# e" U
• 
  
• 
  + `7 v% Y/ C/ R4 y* P
• 
  
• 
  ) b8 ~5 e4 j' A/ `9 ?
• 
  
• 
  
• 
  
• 
  $ o$ C6 ]" j" M$ P- r8 Q. x  J
• 
  
• 
  & i( h* ~& ~' D
• 
  
• 
  + A& y) Q8 v. M
• 
  
• 
  
• 
  
• 
  " O% W% @9 \( r; N
• 
  
• 
  
• 
  / d% x* X( \% E: |7 X4 k
• 
  ! e$ U6 T$ M* B/ x
• 
  : Q( q( F" [9 ?4 y5 q* G; F
• 
  5 u& \7 h$ ]* i) w
• 
  
• 
  $ N1 k$ u8 v; ]$ i
• 
  + F* f( M7 q" E2 e% E8 _! A  v! p" {
  ! ~6 M* x- n" Y& c; B9 S2 R

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))
$ ?/ s9 i% p5 E+ b" s

: Q4 J6 l) _% p

颜色表：

( l# g) w3 X/ x
6 [3 h& X) B. E& y8 z9 y6 @. }, D  ^