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

  
  8 x. @3 [5 L/ ^2 I" f; ~( R3 c9 s

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

% ?7 n: U) V/ K1 `" B                               
登录/注册后可看大图

绘制中国区域地图

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

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

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

• 
  , e9 ^2 U; |, @/ e% P( R2 K: h
• 
  
• 
  : C7 e+ W7 i7 H( x. ~# w
• 
  ' r# E0 c6 c8 z
• 
  
• 
    v" B/ V+ G( O: h0 g& u7 e8 O9 ?
• 
  
• 
  , H' o" }4 ]% t
• 
  
• 
  . \0 K, q' K3 W* Z  }0 o
• 
  ! R# }; o* K, Z$ s( Z( |8 n9 C( k9 T
• 
  . d4 d& e$ Y* U3 D
• 
  
• 
  
• 
  6 B- `0 V; Y4 @" `3 t
• 
  
• 
  
• 
  
• 
  2 @! k6 N+ e& `( X7 J
• 
  ( F( G- n. }2 ], v* g/ }
• 
  3 T9 R$ A# o: W+ c* d) B( g4 x( N! }
• 
  * |3 G/ k  n% F# l7 n
• 
  8 N5 q& B7 w( l8 R. P! S
• 
  
• 
  
  ! [. c/ @. L: X: g/ z6 x6 F; C3 ]# o7 \

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

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

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

5 ]8 z% y7 Q5 J$ q% f5 i

• 
  
• 
  
• 
  6 Z9 z5 \& T* l; ?1 J
• 
  % t5 K" x5 r5 h! y% ^
• 
  
• 
  8 g- I. k' m  y' H! r  r. z
• 
  + U) o4 y. _4 z+ M6 @4 _9 i( X8 ~- ~
• 
  
• 
  
• 
  
• 
  ' A" }  `: O( G4 l
• 
  4 B0 }% H4 _2 M" j. U6 x- v# s
• 
  6 p% f8 N8 v( E: |% a, m: W3 D7 c; s
• 
  
• 
  
• 
  
• 
  * G% h- P/ C" c' c- R
• 
  
  7 n" O( G( Q4 h$ W5 n* Q& E0 c

#绘制中国地图，并标记各省会城市名#引库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")
& m0 c; o! Y0 l

7 a; a2 r0 `" @8 `6 j

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

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

• 
  
• 
  ! Z4 C% Z9 T1 U. d8 S( U& [& E
• 
  
• 
  " A) i, P  o+ r: Z6 X- e9 ^$ k
• 
  
• 
  
• 
  
• 
  ) Y4 @. _, \" i/ K
• 
  ( Y; p5 S$ B  r! \  G5 Q
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  3 E# C9 \; x7 P/ |7 o2 i
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  " R' h/ R  F# |3 \
• 
  ' E6 i% _3 E' J  S4 L+ ?) P
• 
  : j2 a" u" @( V' l! e2 @4 b
• 
  
• 
  * ^" d4 J& A, a
• 
  
• 
  
• 
  - c% ?1 ^. F0 d. L6 z- z
• 
    k) E9 e4 I8 p0 I+ E
• 
  & A% c' \9 e* k, {$ {$ Y3 M
• 
  * ^  y* h3 `! a# T
• 
  # y% f8 r# h+ Q8 b1 t. [) J5 N
• 
  # Y9 I# h3 f! B) Y0 H# j1 y' H
• 
  7 R& U# [" w$ x7 M' l# o  l
• 
  
• 
  ( i9 o. s$ ^& o! I( \9 m; U
• 
  + j; \" U9 M0 k: U! s8 J' S  W
• 
  ; h/ l0 g0 T7 E  L* o
• 
  . J* t% G5 {; S( a* z
• 
  + I. r4 x5 D; l+ M! s
• 
  
• 
  
• ​
  2 K+ P! M+ Y% a& a% D1 i& x; j4 K& p
• 
  
• 
  % g' ]6 ?# Y0 R) F
• 
  
• 
  
• 
  % C6 g& x5 P- t; V, H$ ~6 I/ L3 m
• 
  ( C+ s4 t. v: D; `6 v$ K  {. t
  

##绘制有海岸线的中国地图#引库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()
7 {  l; ]% O0 L% _( [+ |& r7 S9 F

" I# o5 z7 L, I$ l. S7 X/ ?

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

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

• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  " h  z% t1 r7 W4 B0 O
• 
  
• 
  8 B" j3 i/ Q7 m0 i
• 
  
• 
  
• 
  
• 
  " X. z! p$ b1 z* h6 e
• 
  2 A5 J: `1 R; r) t9 P0 z
• 
  
• 
  : {! d# b* |7 r
• 
  
• 
  ) `# X* A6 y7 Y
• 
  
• 
  
• 
  
• 
  
• 
  : X0 q/ S- L/ f% r* S" H1 c5 U
• 
  1 Y! x; d, V* O& y% o
• 
  
• 
  4 P& s  _$ P1 W4 ?
• 
  
• 
  3 a  N& `6 d7 U% [( \* G
• 
  
• 
  
• 
  7 k4 s" K; U) S; C1 p/ c2 n
• 
  % I2 ?, u5 v% q8 \
• 
  0 N+ C* I6 b% C2 S* c5 M
• 
  
• 
  2 c# }# a  N; Q1 K$ q/ a0 @
• 
  
• 
  6 y3 Q( y, O8 O. N
• 
  
• 
  6 b, C& c6 I) j* j& y" P8 y
• 
  2 z' s6 v: K) V/ y. B! ]
• 
  
  : I# c4 J: j! H7 A+ g# u  o$ K

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

. ?$ p# k& |' m9 F( ?4 O                               
登录/注册后可看大图

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

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

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

proj = ccrs.PlateCarree(central_longitude=130)  

• 
  
• 
  
• 
  
• 
  + L( Z  q) b+ |" e$ P- q" A
• 
  
• 
  
• 
  & T7 ?0 g$ \1 S' e! J3 @
• 
  " T1 J0 `$ t* Y! p. y9 J/ F5 r+ [3 M% E
• 
  
• 
  ( F5 P! j- t, i2 p0 s" {+ q$ l+ W
• 
  " }$ D8 @; [& b3 _4 K! a2 @( J
• 
  
• 
  
• 
  ! }: b, l, n8 E0 w1 O- }3 B2 q
• 
  5 D9 x9 r$ f1 ]$ F8 A
• 
  
• 
  
• 
  1 N/ d0 O6 L; _, \- d6 M
• 
  
• 
  ' R" ?) U! Z. k( ^: d
• 
  & M: Y! G- R8 f6 X2 q+ X- a
• 
  6 L! h7 s/ U# E7 }% U
• 
  
• 
  : k" G4 r) P* Z1 Z, L: z- |/ f& y
• 
  * H; S3 R7 _$ q) v0 G
• 
  , g8 n* k1 Y& U& \7 u4 l) Q
• 
  
• 
  ( w4 e: g0 G7 i9 ^2 d
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  ' ^+ n0 y8 g; j$ a0 h4 W
• 
  3 c0 q6 z9 p) `
• 
  # A' _" u. z! l* c" \
• 
  
• 
  
  3 y1 ^0 a/ ^# {! t. ~' 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))
) A* N/ L3 d1 R5 Y8 ^! K2 I

# r4 y; @0 ?8 r$ }' T! Z8 i
, h. e1 N1 ~9 W- H3 J6 Z

颜色表：

' I' W' I: w% M' E$ c7 `, b5 o
+ n% I* O2 @+ x% V5 P( z