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

  
  4 s1 b/ Z: N" _% y0 R& [. T9 Y

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

                               
登录/注册后可看大图

绘制中国区域地图

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

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

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

• 
  
• 
  
• 
  
• 
  
• 
  
• 
  5 |% U" J5 u/ s( A' h) q) |
• 
  
• 
  * O9 a, Q6 t+ h# R2 @# [
• 
  
• 
  7 I- U# r# y1 F
• 
  0 S* A# Y  @0 }  {2 u6 o% K
• 
  3 q) D7 b. @& \& W2 Y: j
• 
  
• 
  5 q0 q( l' [' f' A! U! B
• 
  , J; o( D5 B" |6 M/ Q+ _# O  ~2 ]
• 
  0 v# [: {3 @( x. i3 Z3 [
• 
  ; n) C: L8 e& C5 R3 \) P8 L. o5 |& q
• 
  
• 
  
• 
  2 Q' x8 {) G, B) M2 n! S8 a  m
• 
  
• 
  + H# g( K& D$ J! ~; p1 n1 w& a
• 
  % o* [* |" k0 L1 m" @# v
• 
  - M1 o8 F- c. Z" {
• 
  
  

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

/ n& ?2 e, g, z* }- V

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

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

' E, ^0 T# {/ v0 k9 @

• 
  
• 
  , w. ?! V2 L6 X. \' L8 Z
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  6 Q; D3 ]+ u2 V# H! ]
• 
  
• 
  
  

#绘制中国地图，并标记各省会城市名#引库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")
; o1 W4 z# R6 R  u/ Q

& V) [# o5 ^8 |& T. _1 Q! C! Q

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

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

+ i  t- \+ h$ e) H- n% O2 j2 X

• 
  % H* [, W7 w" l, W
• 
  1 o0 k8 z+ Y: p$ y* D
• 
  
• 
  ) E1 E) R% b. ^& M3 m
• 
  8 U4 ^! n( f% S; M- h
• 
  
• 
  
• 
  / W2 B3 M4 @( r2 y9 N
• 
  
• 
  % [1 y! I' S- T- i) s) F" `  \7 q/ H
• 
  
• 
  
• 
  
• 
  
• 
  0 N! c  d% g* A  F
• 
  
• 
  
• 
  5 ]& Q% {% x8 k
• 
  9 c9 N; ^* ]( \+ A: Q: E9 i
• 
  : X) Z. |; s3 r* A$ f* k- q
• 
  % y0 \  O! l- g; Z4 h! V
• 
  + z+ k5 {$ N$ }
• 
  
• 
  ) y+ a; i% f# g8 Z
• 
  
• 
  
• 
  
• 
  1 I7 Q! ^* Y. T4 _% ^  N
• 
  " A8 M- U7 b2 M1 C
• 
  
• 
  
• 
  ! d  S* ?6 H1 n' u+ \
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  ( W% P$ X4 @. u0 q% o, d! c" u
• 
  9 ^4 @; j& B0 N6 y' m
• 
  
• 
  ; \: ]9 B4 S, H/ z  i5 I" ?. g
• 
  
• 
  
• 
  
• ​
  
• 
  
• 
  $ z6 V) Y5 _, n; n0 v
• 
  
• 
  
• 
  , E2 v- H: u; r4 B
• 
  
  ; r3 s- ]2 V7 ?2 A5 v+ k1 W

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

$ I+ n  Y& _2 b6 u# Y

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

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

- t4 r4 r# y; T1 T6 N

• 
  " U. w# F- ~) y* E+ v8 m
• 
  
• 
  ! e- x: x; l3 J. F  u+ X) J! h
• 
  & ]" T. D( T4 d( L" e
• 
  
• 
  2 R# x/ W: _& g  e3 P% `
• 
  
• 
  
• 
  , y/ s3 {6 y. b. Y8 ~5 \0 x
• 
  3 F! D! f$ S- U4 W
• 
  
• 
  & g3 k4 M* z9 L
• 
  5 u1 R0 N; M" h& y: ?
• 
  
• 
  
• 
  ( K! v1 q& Q% w/ L2 Q
• 
  9 M3 f  \! f! P- ^
• 
  * n( o0 _) o3 [' A  \, S/ D
• 
  
• 
  * u5 s; L2 I+ P' m% j. S9 g
• 
  ; l+ ?; a% M7 i$ M$ }' a
• 
  5 e6 _2 G0 E8 j- R4 W( N: m
• 
  % U# ?  n6 _' z9 @+ k
• 
  2 v8 {3 s1 t' ?7 {' E' a" a8 }; @  {
• 
  
• 
  1 j( i% v1 U1 g& s% x
• 
  - r) O% ?! _  z3 Z) o
• 
  
• 
  1 ?# U& a8 y" o5 u
• 
  
• 
  / `' Z/ ^1 S0 u; ^
• 
  & L( h& T4 i6 c0 t9 u( O
• 
  
• 
  
• 
  6 M, _4 u2 l- f$ F
• 
  
• 
  
• 
  4 i. ^* i1 u: N5 K. E
• 
  
• 
    R% f0 S2 n( W* b: t" p
• 
  2 X  g* N2 {! o! R: M- P
  

#例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))
; ?1 e8 k  [9 n( h4 P6 A9 k

7 h. v+ d! r$ I+ z

                               
登录/注册后可看大图

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

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

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

proj = ccrs.PlateCarree(central_longitude=130)  

• 
  * a7 E5 v/ {; U" {
• 
  % y& D! Q& ^2 q$ `2 P/ {
• 
  ( b6 e8 @8 ?1 z$ Z- |; y: A$ B7 F
• 
  $ ^- d" \4 i9 K" U! w
• 
  . c: S2 O" h8 K5 j
• 
  
• 
  1 V* z- _/ F7 u. @+ ~% l
• 
  ' D6 @: P; H& L8 z/ o% y# T! D
• 
  
• 
  3 I! q" ?/ [; d, K7 q
• 
  
• 
  
• 
  
• 
  - _3 t( b- \5 [
• 
  ! l* Z: t' E( i
• 
  
• 
  
• 
  1 s. A4 E+ L, p7 O
• 
  
• 
  * v9 L: l7 [5 t, m/ I& u" x7 a9 F
• 
  9 g# [+ e8 h1 h# b  c. k! G
• 
  , @5 o- Q+ U: x  s) f$ x
• 
  7 m7 t% n( L2 Z- _: b
• 
  
• 
  4 @; K* P! j+ n
• 
  6 `1 V4 W& S  H, F) Y$ k
• 
  
• 
  * G# m8 z4 I' e8 b9 K' d% o$ s
• 
  6 }; H- Z' U$ n$ a! m" [0 h* \7 h/ m
• 
  9 v, N" P) @- ~* N
• 
  
• 
  * J  U3 O5 C0 f# g* y
• 
  
• 
  
• 
  0 l- {* r  y+ X& N  O  P
• 
  6 o3 b" _% t6 x% _
• 
  : Z: W+ I( g5 L# R# x2 W) e; m1 z7 \
• 
  
• 
  
  

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))
0 S2 \# w) y. l+ l

* \* n1 i8 R% e* S) C

颜色表：

+ e. K2 a' [* I9 W4 h2 M