【气候软件】Python10：绘制标准的中国地图

( b) i: a" O. L* u

气象数据可视化主要依靠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("···")

  
  

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

/ z' `- `4 F8 V

绘制中国区域地图

前面python9的文章（Python9：绘制中国/世界地图）关于中国区域地图的绘制方法除例3，主要使用的是pyecharts自带的地图，其中中国区域的地图有问题，所以不采用这种方法。下面根据其他博主大佬的文章（python绘图 | 中国地图最正确的使用方式（九段线&南海子图&高分辨率地形图一个都不少））进行适当改进，并采用权威的国家标准地图比对无误的地图文件绘制中国区域地图。以一个具体实例完成：

（1）正确的中国国界线，省界线，九段线的绘制；

（2）全球地形图的添加；

（3）中国标准地图文件的添加。

提前需要下载准备好一些必要文件：

（1）正确的中国行政区划矢量文件；

（2）全球地形图文件（从大佬提供的方式中获取的全球50m分辨率的tif图）；

以上的文件可后台发送“中国全球地图文件”获取。

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

• 
  % X# j7 y1 }8 b/ E# G$ A
• 
  
• 
  
• 
  
• 
  
• 
  " A8 c& x1 H' z% \9 ]* P0 Y
• 
  
• 
  . i* A/ x& `  C+ ?4 u* h
• 
  
• 
  ( u7 u; b$ V+ h& ~3 w
• 
  - N+ h8 @( j; [. O
• 
  7 M5 H  i  [6 y. O
• 
  0 K9 h+ w/ N5 _
• 
  # X- T+ c$ `& x) o" y
• 
  # A4 B1 H* Q6 ~
• 
  
• 
  & j% I8 @! i, o9 u1 ?7 u9 b5 T
• 
    E8 Y5 J0 g! S1 z1 p. G6 }: \
• 
  ; J5 e1 o  ?  ?
• 
  ( Q0 G0 I( O; o8 q5 ~
• 
  ; M3 M  M3 a0 `
• 
  * p5 _' J$ m3 w) m+ F" y+ k8 w
• 
  
• 
  : U- r% R( w+ N
• 
  / _  }4 l( ~, J4 }2 K
• 
  9 d: j* w4 ]2 a+ e. U
• 
  
• 
  9 g! E9 O2 t4 Z
• 
  ; d. h3 X7 W. e' C" H
• 
  ; n1 m3 }, b4 o
• 
  $ I% n) ^5 C6 b" f$ S3 N1 Y
• 
  7 z8 G) R# U* O0 u$ `3 R# d7 W& r% m
• 
  
• 
  + f  a# l7 B; z/ G$ N& W7 e: @# N
  $ J6 l4 x( E& t% J6 ]. J( i

importnumpy asnpimportcartopy.crs asccrsimportcartopy.feature ascfeat   ##添加地图其他地理信息所用fromcartopy.mpl.gridliner importLONGITUDE_FORMATTER, LATITUDE_FORMATTER   #添加经纬度所用fromcartopy.io.shapereader importReaderimportmatplotlib.pyplot aspltimportmatplotlib.ticker asmtickerimportwarningsfrommatplotlib.image importimreaddefcreate_map():extent = [70, 140, 0, 60]shp_path = r'./cn_shp\Province_9/'#可后台回复索取（这里很容易出现问题）# --创建画图空间proj = ccrs.PlateCarree()  # 创建坐标系fig = plt.figure(figsize=(6, 8), dpi=350)  # 创建页面ax = fig.subplots(1, 1, subplot_kw={'projection': proj})  # 创建子图# --设置地图属性reader = Reader(shp_path  +  'Province_9.shp')provinces = cfeat.ShapelyFeature(reader.geometries(), proj, edgecolor='green', facecolor='none')ax.add_feature(provinces, linewidth=0.6)ax.set_extent(extent, crs=proj)ax.set_title('China map', fontsize=5)# --增加高分辨率地形图（需自行下载） #可后台回复索取ax.imshow(imread('./NE1_50M_SR_W.tif'), origin='upper', transform=proj, extent=[-180, 180, -90, 90])returnaxif__name__ == '__main__':warnings.filterwarnings('ignore')ax = create_map()plt.show()

% n' {: D0 W  x9 t: J& H

3 S# n: j6 n6 f8 I% P3 u, c2 A5 V                               
登录/注册后可看大图

上篇文章（Python9：绘制中国/世界地图）的例3：绘制中国区域地图，包含海岸线

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

                               
登录/注册后可看大图

5 J+ b6 E8 Z6 t* `8 t

• 
  
• 
  3 U# `( X9 p& D7 c* H# v# _! f% T+ Z
• 
  
• 
  
• 
  
• 
  5 \+ P& B  v% k4 U- M
• 
  # E6 Z9 Y* a2 P/ W  z' G9 b; z$ z  U
• 
  
• 
  
• 
  
• 
  
• 
  
• 
  , a9 R- V, Q" A5 R& ]
• 
  
• 
  ) R. E. M1 K% X4 d3 N) G0 w
• 
  7 Q8 A. Q) u# k) Y; n
• 
  
• 
  
• 
  
• 
  - L" i$ r0 P7 M3 L: A6 \
• 
  2 H8 [& o( ?, @: k* {8 s
• 
  ( C; a- T) F0 B, u# l
• 
  9 H) }9 `' J+ q2 m9 \% z
• 
  & F0 ^' H! B: v: D; |7 I! {2 w
• 
  
• 
  ( S6 m# r- M7 M% m: q- x2 G$ c
• 
  
• 
  9 ^& N: P+ A: W! k9 g4 U; D
• 
  4 y. [0 b0 f  h4 ~7 G) P4 R: m% K$ ]+ o% a
• 
  ; I; t2 }7 D3 A# f- ~5 F
• 
  
• 
  
• 
  
• 
  
• 
  + C3 Q1 M4 x# x8 i. s4 u
• 
  ' N+ [9 k6 t) N
• 
  
• 
  
• 
  " }# u) |) m' F0 {9 z" h
• 
  + ?  F3 f9 r! U; ^- k  D3 n. F" [; m* Z
• 
  ! n) R2 P% V; T% ~: `: q* s
• 
  * L0 I" t: }0 T8 j2 ?  d6 O1 S8 K  K
• 
  : {" y' q3 D1 T: V1 D
• 
  
• 
  . `- d- ^* S) ^( X
• 
  
• 
  / H" s+ Y& S1 w7 d
• 
  
• 
  
• 
  8 b7 r' J* O8 v
• 
  
  5 R* @2 n; X: v, B

##绘制有海岸线的中国地图#引库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()
+ E  I6 Y" F! @1 W* T8 A* Z! R, e

: k# y. q7 |2 M' V
) G2 I* T2 r4 l% Q

颜色表：