【气候软件】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 e3 l. T, R. m

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

                               
登录/注册后可看大图

绘制中国区域地图

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

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

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

• 
  
• 
  % n5 H. \. o' y. i
• 
  - Z. V+ F4 x% E; ?' e
• 
  
• 
  5 }( z& D; B! ^' K! a
• 
  
• 
  
• 
  
• 
  : q7 O6 d; t2 J5 P6 E* H2 \  o6 R
• 
  
• 
  % F4 L6 z9 K& L0 |& l, a: x
• 
  8 \: r) Z8 _& e- X. D' e
• 
  
• 
  
• 
  6 s0 P5 s% r* u) E" P0 [2 d2 s! A
• 
  
• 
  
• 
  
• 
  * P2 w$ B4 U; n- [9 c5 z) |5 K
• 
  % H1 c1 a& J( ?3 D$ u
• 
  
• 
  
• 
  1 h9 m, \- y+ ?0 K4 R$ j
• 
  
• 
  , N& v/ |4 ~- o9 m$ J
  9 u! O. y4 {+ E; ^  R/ I$ 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")

3 K( Y. U! x) p6 \

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

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

1 H1 v2 s0 r* D) ^0 W/ |, n

• 
  % R0 X9 U/ T# k8 z5 ^! [
• 
  9 K% \% X2 c" D, Z
• 
  6 @- E2 O% n! ?" [
• 
  
• 
  * @$ B! U. e! C$ q: `- [
• 
  : y4 y% z% c# e% v6 r2 @# f" p, E
• 
  ! G; `* i  V; _3 m! G5 ?
• 
  , X* u3 S: ~, f: f) k3 J
• 
  
• 
  ( A, {- G2 o- G6 M
• 
  3 s1 Q5 y3 i9 _! }- l( F
• 
  " m5 o! J/ m7 W# {) i
• 
  * c' e' V% z# Z. ]* @
• 
  
• 
  9 T6 Z" g/ `. c7 S& v# _
• 
  8 ]) O' L# A. Y9 [+ J
• 
  : j% R1 ^1 {. `
• 
  
  

#绘制中国地图，并标记各省会城市名#引库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")
, `" `( q. u8 {( c1 e, D6 K5 x2 `

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

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

• 
  7 F. Z- v- S" r. y2 O& E" a/ p
• 
  
• 
  1 a. e/ [/ d$ O& ^
• 
  & N+ J" ]% s4 V6 c7 D/ x
• 
  + x# X6 F5 O8 [+ H
• 
  2 u9 V- i2 v. O+ t/ I2 P
• 
  ! w6 M; J0 N! t( m" F
• 
  
• 
  
• 
  * j' u: A9 G- i; P# @/ o
• 
  
• 
  + n1 C5 f% V# I
• 
  7 c5 b/ u/ k( L. ?5 i7 Q
• 
  
• 
  * ]& \- I8 i* g$ N8 ~9 C
• 
  ( B: w* v0 F. M0 E+ ~. @% \7 z
• 
  % [+ d0 P6 F. l
• 
  
• 
  $ O0 x3 l5 o" Z3 @# e5 z
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  , H% o4 y2 l- {& V. N) \! L6 C# {
• 
  7 l3 Z9 j( k- j) _# ?0 V8 b
• 
  8 M, W; R6 o9 v$ I3 P
• 
  & u9 c$ N( G3 O  M9 @; K
• 
  $ b$ `' {$ s4 r6 V
• 
  
• 
  
• 
  
• 
  * }& n/ F0 V6 [. g' g. h
• 
  
• 
  
• 
  
• 
  ; t" m( {  y5 T: u& n" t
• 
  % F, u- h* P: Z. K$ X( }) A
• 
  3 A% F0 V8 m8 l$ d/ x0 R
• 
  
• 
  
• 
  
• 
  , f. M+ V' e7 I  o0 Y& n9 ^  M/ Z7 l& l
• 
  
• ​
  
• 
  ! \, x& |+ x8 [; X: x' m
• 
  0 w7 l. `7 G. V
• 
  
• 
  4 Y2 ?# f; o) x8 R: o
• 
  
• 
  
  ) i7 W0 p" m7 Z7 Y- f

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

+ S. _5 I4 D3 ~+ L$ k/ }

• 
  
• 
  # P6 {6 G& z* h# q
• 
  
• 
  6 k# ?9 s( [: I5 y2 o4 E1 `" \0 S
• 
  / F4 q. s7 X9 Z! s/ }
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  , e  E) T+ B9 C9 b5 P: N" r% C' @
• 
  2 a3 t. B+ M1 B) R1 }6 b
• 
  4 R' n) L' ?# d
• 
  : U* f! I; D2 ~9 j; l
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  2 H2 q2 |/ o8 J5 a- ~6 q1 m; y
• 
  4 _; ?6 Z  F3 ~! o  y- P; m* E
• 
  - n& I; N" R; H2 j8 x
• 
  
• 
  
• 
  6 n4 R2 ~* O" W
• 
  - Q' P5 g+ E; m* g- N) {
• 
  
• 
  ; B5 q0 f  b/ O. a3 o
• 
  9 i! }* a; D2 z/ \4 {: Z
• 
  
• 
  
• 
  
• 
  
• 
  5 S* W2 H+ Y% G7 U0 V; @% ^
• 
  1 P, u" q' e6 \, v( `
• 
  2 _- `/ W! W" _0 w+ S  i
• 
  
• 
  
• 
  
• 
  
• 
  
  

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

! w" b. R  g  D9 q/ K                               
登录/注册后可看大图

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

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

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

proj = ccrs.PlateCarree(central_longitude=130)  

• 
  ) M8 }5 G; p) N' a4 R
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  9 p" s1 h0 G* n# f  s3 |6 Y
• 
  
• 
  
• 
  # V% `) B4 ~' l# b1 R- K
• 
  
• 
  
• 
  
• 
  ; o9 @3 r$ `( y/ |
• 
  4 |+ P, r; q9 b1 |6 C" p% r
• 
  6 R$ @: h2 h0 X
• 
  1 R/ e% U; U/ ^6 z4 W8 S& I
• 
  
• 
  
• 
  
• 
  
• 
  4 m+ A. U% ~# k2 e: {9 b' H0 ]/ r
• 
  
• 
  
• 
  3 R+ E2 n4 G9 R) Q$ ~+ i
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  7 e# E2 h. a7 b. L( x+ o) A8 t
• 
  
• 
  8 q- {( Z1 b/ H0 o5 @- k; A
• 
  . x4 ]1 O( B6 `0 Z; j5 ^" g! S
• 
  
• 
  
  

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))
2 R* w" ~* M: C

颜色表：

' h* v- S) r" c* s4 [8 R