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

  
  / J7 h: N4 R6 M, @4 d7 {6 C$ S

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

                               
登录/注册后可看大图

绘制中国区域地图

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

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

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

• 
  - x, ~. k7 s# Y5 C! F8 z
• 
  
• 
  
• 
  
• 
  % W7 ?- `2 y: p: H# h* R
• 
  
• 
  
• 
  
• 
  
• 
    t7 N5 r% ^, q" V/ {4 x$ X" a
• 
  
• 
  : t5 \9 Y: ?" j7 _! m
• 
  8 q' x7 P6 T7 I1 w
• 
  
• 
  
• 
  
• 
  
• 
  2 ?0 ^3 {( v1 K5 T+ [9 N
• 
  8 E( n+ p5 r5 E2 |% w
• 
  
• 
  ) }6 A# Q. m/ c+ j& p) k& _
• 
  , O$ V) S0 Y& |- ?
• 
  2 G) j% F: Q8 J9 ^. ^. \
• 
  
• 
  
  9 S6 g2 F3 q( P; q; M3 y3 A. D

#绘制中国区域基本地图轮廓，包含省界#引库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")
) c6 i# a( Z; u, `3 x' e

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

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

: @( c" F- b+ H- j4 Q, r1 n' |9 C9 _

• 
  , B2 z! {( N* |/ F9 O
• 
  / t( k' l1 V6 T0 x
• 
  * {* b, D- ?) _
• 
  
• 
  
• 
  0 _- `4 V! g6 J' Y* ]* L) w' |
• 
  
• 
  : Z$ W# d: x4 E" v4 u2 Y7 A# d
• 
  $ Q3 `; J( J- C+ j# y& N4 ]3 W
• 
  9 P7 F. s9 l3 p/ R* N
• 
  # B( O( V$ t2 J
• 
    s4 X! r) a7 ?- u, ~. P
• 
  * {. }7 \. s1 `/ a3 H# b& I
• 
  5 t; ^( _  i: m/ U  P+ h
• 
  5 y$ [5 l  j! s" D6 _, T
• 
  " a0 Q- i$ M3 R" e: }) O" y( B5 S
• 
  
• 
  
  

#绘制中国地图，并标记各省会城市名#引库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")
5 }1 ~0 v" i9 C& {

0 l: S( [* y1 w/ t7 c) s

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

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

• 
  5 ]; D1 U* K: u- X$ T1 l" F& {
• 
  
• 
  8 L! O1 @0 h: W
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  1 C. |! l% Q" o: k1 a) q
• 
  ' |% _& n% ~1 a: Y
• 
  0 n" o" W! D2 V8 D" s. }
• 
  - y- C5 T+ ?3 Z) D6 z' e1 Q# P/ c
• 
  ( l9 M5 k# [" F% C
• 
  + d' o0 y; C5 L6 m. I  y
• 
  
• 
  
• 
  
• 
  
• 
  " r3 R* B2 ]# l5 c4 A1 F% }
• 
  + t1 S& t4 Y1 |9 C; N. A3 M/ \
• 
  + s% p& W' G( I( i* n% W& M1 a
• 
  
• 
  
• 
  $ N, w: [: Y! h
• 
  : D' S1 R) K6 U2 h6 _, B2 N- f% W
• 
  1 V0 C/ V1 W% {$ i- ]
• 
  4 @% m2 C4 a: B- t% p# Y, X2 ~2 d8 Y
• 
  
• 
  ' I" f8 Q6 F3 F3 P; C- W! S
• 
  
• 
  ( U* N' c- W* g! N: B* {  f: z3 j
• 
  ) w6 |, _9 V! `9 T
• 
  
• 
  
• 
  
• 
  / j9 G2 N! {5 P9 S! j$ P2 f# a6 c) V& ^
• 
  
• 
  
• 
  
• 
  9 S5 t1 s( p. R% {
• 
  & C8 E; X* Z9 U3 _4 g/ T
• 
  7 \) S2 t' d; m% v
• 
  
• 
  
• ​
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
  5 ~# B$ t) y& G9 u

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

% }$ Z: p; l! Z8 W% J3 k% [- X2 `
" `! f' L3 J) r7 m. l0 p% \9 f

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

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

• 
  ! M" c0 D$ }1 @
• 
    v2 T, p6 R: x8 V2 W- y
• 
    S, |: e" \- `! o9 @% M( s0 Q
• 
  2 u* |( E5 S3 }( P4 C* X/ X9 g
• 
  4 @* V. L( E; Z9 G8 X3 f! [
• 
  6 e# |; G- _( F4 \
• 
  
• 
  
• 
  8 V7 m4 A+ Q% S
• 
  # g$ Z8 ^# w4 [9 Z6 [3 K; e
• 
  , ?& }% F/ h* @6 ^/ Y6 ]/ [
• 
  % V& x  ~" m& L& s
• 
  
• 
  
• 
  
• 
  5 O' f/ Z( C7 Z( n1 B. e, b
• 
  
• 
  
• 
  # T7 v! A) }+ |5 x" \6 {% J
• 
  1 @3 p& T$ r0 i0 k, T8 `
• 
  
• 
  
• 
  
• 
  ( K2 |4 x3 h8 ?1 D
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  1 c- ^& c  l  @3 l9 c
• 
  1 J7 [! o' b2 r# _+ R8 D! o7 D3 k
• 
  
• 
  
• 
  
• 
  
• 
  / }) @% {. r; J/ i" I
• 
  + V7 U( _# y- k9 C1 `- _9 F9 B
• 
  
• 
  8 f, F9 C$ q. W- U* ]- X
• 
  5 P9 q; ]$ Z. z; ~$ R
• 
  
  - P0 G/ w9 ?- r

#例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))
" t5 z5 T8 P6 C/ ?1 z

. F1 y  x0 T( y( j- Z8 \# ~% D

0 k4 \; E: R3 |* v                               
登录/注册后可看大图

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

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

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

proj = ccrs.PlateCarree(central_longitude=130)  

! d) o3 t" {4 D% i1 s( c

• 
  . c* n5 {& r' h  S& z' U
• 
  5 u5 u% j* c0 P) q( d* Z( z
• 
  
• 
  , [7 {3 c( ~, \4 {) s
• 
  & V3 J" A& O* X1 H: w4 h
• 
  5 A* n* a& E0 F. t8 p  ^, y
• 
  
• 
  
• 
  # V; W2 i. `) ]2 A0 E$ o  h8 W
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  ) g5 J" w5 E  c/ s' Z5 Z+ q
• 
  * ]" |6 ]: {' U9 V/ ^  ]
• 
  ; F9 T6 z3 [5 j# `$ M- |1 L  h
• 
  
• 
  
• 
  
• 
  
• 
  5 O& a/ K8 G1 N' m
• 
  
• 
  $ ], Y6 O7 j0 h* E( P' ?
• 
  $ Y7 {! j# e6 h5 J$ m* j9 x1 B
• 
  
• 
  & M5 o- U) L. C
• 
  
• 
  
• 
  
• 
  & J2 X3 Z2 k3 }/ I% r8 f
• 
  : S1 V1 B5 Q, o4 F- J2 |
• 
  ) F# Q0 O, e; Z1 i- X' M. p
• 
  & r/ N+ T3 L0 @3 B
• 
  6 G6 @. w' \) O3 ^7 H1 j
  

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))
/ z, a, w8 \( O3 X, C- Y% g" `: S) `

8 c$ s2 Z" ~+ K1 g1 A0 O

颜色表：

. W" f  A7 p6 [  L' H5 G