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

  
  

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

) y1 w. q5 t1 _4 S& s" D" m% R4 D                               
登录/注册后可看大图

绘制中国区域地图

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

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

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

• 
  ) ^0 F' c) M0 L- B  K" n+ {1 F
• 
  
• 
  
• 
  
• 
  / P, Z1 z" W  |& a/ f
• 
  8 B; n! O1 {  b3 i
• 
  0 N$ f3 a% v, ~0 y
• 
  
• 
  0 X, c: D% c0 S1 l/ x+ D0 G
• 
  * g6 G+ Q" X5 k& B/ i0 i
• 
  
• 
  & a/ Q- Z, Z6 W  s
• 
  & ]0 x& b. D' l1 |* G1 s
• 
  
• 
  ) M4 V& C( P# W2 [, x  v: A/ E) n& R
• 
  + u. T. c4 w' |
• 
  
• 
  8 y* p4 `0 h) M, N7 p7 v3 Y/ R
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  7 h/ _. ]. ]3 h& v
  4 _5 {$ Z" C. }1 p1 L4 B# r+ R; q

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

: `9 j0 u" c* _! N' P8 E* e- u

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

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

• 
  
• 
  ! @9 @8 z% f1 H! e9 X
• 
  
• 
  
• 
  - ]8 S  G5 v# l1 e( F/ {; s% l
• 
  
• 
  
• 
  3 ^9 A" p! a/ _) q8 d! F# u
• 
  
• 
  
• 
  
• 
  
• 
  / I% f$ }! A! u
• 
  5 ]/ Q% u, Y. X8 P( ]
• 
  
• 
  
• 
  
• 
  + \1 F3 d& @( W; `6 ?
  

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

/ H) w- v/ S6 C3 S

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

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

) Z2 X% ^  _" ~  A% b7 [  ?: g& a

• 
  5 O4 \# u' }9 Y0 j- l  |! P& b
• 
  
• 
  
• 
  & F& B. s1 ~0 i  N" Y3 r3 W' q
• 
  ' J' f% L* `! k- e6 e/ c
• 
  
• 
  
• 
  4 N& s% {3 e, i" o) H: u  ?
• 
  
• 
  ; S/ q( B4 v6 P$ i
• 
  
• 
    s2 u& ^# s7 w
• 
  ) |6 H; a6 k. t. u2 r
• 
  ! V- ?; G3 u' \% K! N
• 
  
• 
  ; X5 [# K8 T6 _( g& Q
• 
  : ^: K) ^& n( \4 ^& @+ G
• 
  
• 
  2 h3 v7 D4 r1 `
• 
  9 t/ J7 V, l5 L. U
• 
  
• 
  - [; O7 H, \( u8 \$ k
• 
  
• 
  ) Y# g& K9 @! O% _6 p% P! t; p7 B
• 
  5 @$ ^2 r% p4 D9 w- S* E( ^
• 
  + h4 ^  D8 E6 B, D) S
• 
  
• 
  : W" Y( v# u/ O3 b, W
• 
  
• 
  # M. r  O7 q+ k0 L7 ?, \
• 
  
• 
  1 r* D' E! {, n& I  b4 X2 o  j( b6 V
• 
  5 ?- R- U! ], ~( q: H
• 
  
• 
  
• 
  1 u0 b- U! ]5 ?7 n6 p
• 
  
• 
  6 V4 i2 V* v9 j
• 
  3 \! r+ T: E5 r6 X$ z- {; u, O- [. G
• 
  
• 
  
• 
  
• 
  / _" H- R# l8 h8 G* b9 d
• 
  
• ​
  ( y" b  ^( ?2 b6 \" m
• 
  
• 
  0 K2 v# W! u& s
• 
  
• 
  
• 
  ) p& G5 x. `* O
• 
  
  

##绘制有海岸线的中国地图#引库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()
1 `; E, ~4 x, x3 o

! r, Q: E1 D; @& R) e6 i% L& S6 e

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

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

• 
  
• 
  
• 
  & g5 ]6 h1 ~: z& c5 \
• 
  ! K: a5 T5 m8 K* H
• 
  9 A7 X9 I7 h3 ^
• 
  
• 
  4 N" c1 w! H+ v' S' @
• 
  " P8 T6 J: k3 u# {4 B4 _3 O4 P9 G/ |
• 
  ! f  ?) }2 Q) D3 {7 y" ~5 t
• 
    g5 L) I5 F( {
• 
  
• 
  
• 
  6 b9 W# q3 e1 R) e
• 
  ) ?5 m& G1 ]& @! g: `
• 
  4 f( x( f9 o6 C: I# o" Q1 x9 O
• 
  
• 
  / P8 b. N5 i; ], d; N6 l
• 
  ( T# C8 |4 v& J
• 
  1 w/ [, }8 _7 P8 ~( X( |
• 
  
• 
  
• 
  . ^% ^6 E3 T4 u7 Q+ y& b4 Q+ ?3 f
• 
  & a2 ?+ E: D, q/ b3 K- d
• 
    V, }4 m0 M1 ?) s; e) f
• 
  * {, L* p  t# ?) g8 E, s
• 
  , @/ g% K+ d+ O8 M; ^
• 
  ( B$ L: w0 R: v1 o
• 
  
• 
  
• 
  
• 
  
• 
  ! Z$ A( F, Y( c9 R+ t3 s
• 
  
• 
  4 Q" E7 }( {$ `( k; O( b* q; i
• 
  
• 
  
• 
  
• 
  
• 
    K6 b1 Q; h, c
• 
  
• 
  6 w* z! @  Z# X, k6 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))
; Y$ J$ p- _: a* y. j

" ^% C( N4 y" X: W- L) r                               
登录/注册后可看大图

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

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

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

proj = ccrs.PlateCarree(central_longitude=130)  

& p: x$ j6 A  N. ^8 t" \( o

• 
  ! Z1 i' Y7 c1 U, V7 |7 V
• 
  
• 
    C$ n! T3 }4 T7 u5 H. |  t
• 
  1 t2 P9 \, |. U! q9 ]  ]
• 
  
• 
  8 n/ m6 W; ]* M; ~' F
• 
  
• 
  1 d* s# c: M& T: r1 n) |
• 
  4 z' P* W, u& z9 d* T$ B
• 
  ) r9 ?2 P5 n/ A* f7 ?' m) g! q: V. o
• 
  
• 
  4 h$ N8 p/ @  ]( p9 O
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  ) S/ o% C- [" n# F/ v* z" m
• 
  
• 
  
• 
  , F* W8 `6 h. S  ~! S: a
• 
  
• 
  , h; c! H7 n: e- G  O0 r
• 
  - I$ T0 H( Y7 t! \" y; P
• 
  $ F& u* x3 I( J7 ]: y( i
• 
  ( }8 D: x' w+ K, }! `6 u. H9 y0 H
• 
  6 B, U8 z: c0 r! u* E1 {
• 
  
• 
  0 r5 z+ a" R7 m! _( F. ~
• 
  
• 
  $ t3 V9 T) w5 t$ l+ T
• 
  
• 
  
• 
  % }( l3 f9 e8 W" Y/ U* g
• 
  % D/ l2 B# Z' r; ?
• 
  
• 
  
• 
  
• 
  $ L7 E$ Z+ `  E, L% g2 P! a! k6 s+ L
  1 G) }5 K4 A! E! k9 V5 e' X

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

颜色表：