收藏本站 劰载中...网站公告 | 吾爱海洋论坛交流QQ群:835383472

[Python] 【气候软件】Python10:绘制标准的中国地图

[复制链接]
0 U3 P9 Y/ a/ V6 w; ]5 S% E+ @" p

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


    2 J8 _8 \7 q, |2 K

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

! X  J* K* R+ w3 g/ X) M) l

绘制中国区域地图

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

(1)正确的中国国界线,省界线,九段线的绘制;

(2)全球地形图的添加;

(3)中国标准地图文件的添加。

2708a4a2d4c2f8f69c367e9c20ffc43a.png

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

(1)正确的中国行政区划矢量文件;

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

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


/ n# Z. J0 ~3 h1 |% I* c

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


  • 5 q; S% t' v) S8 ^

  • 3 ^  [% u* F' g
  • " ^% ?+ _, F: s- h8 b6 J
  • 2 N' g' {. h3 c1 P' L

  • . h: D* H/ s* k, u7 `0 u; o3 m

  •   t* ?- e3 r6 U. R! G4 t

  • # Z, J  z* z7 }1 \0 q) L' }
  • % D6 q* E- ^1 q3 Q; i
  • " J, ^( P8 G3 i  {

  • 4 c5 e2 g; H. a
  • 2 ^& A) Q& }4 E6 a) L4 v

  • ! X: `( J, g; G8 O+ c
  • 2 Z, R6 o) S5 O6 L% s1 |

  • 6 A) |* ]& `/ ?. }

  • : J! c: x8 Y5 l. }9 }
  • 7 g5 g9 S0 `0 R8 V& i3 e# h" }

  • 6 ?( m+ z8 ~! Y4 x

  • $ ~$ L  ~6 \" Q! F
  • 7 X- U* k- {7 ?+ Q5 z. C

  • 7 I4 e: Z" u  m1 k7 A

  • " l: F, M: m: J) ?( s6 g  q

  • * h2 R+ s  h) K( h) `  f9 u
  • 9 n- t" k( S8 r2 Q0 M1 k

  • 3 k& e1 x1 G1 M! ^

  • 5 q, j( T! |: o2 _# T
  • 1 ^. W3 k( X$ z6 M1 r- V% h

  • 5 A) v9 T: B$ [, e. ?; p5 h& N
  • " e. X2 I& ~3 w" `) B7 e  |: N

  • 1 o" A3 r" y% y$ p

  • 9 C7 q0 {/ b. e  w) }

  • , O& O+ n! o! m4 p" q0 V
  • ! E% |% b- N! x. ?' Y+ b

  • ! |9 y5 U& f: v: _( {8 Q" W
  • & a# C7 E5 J( f% n/ c

    ( f+ l- \3 D9 Z' P' }) L9 o

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()
) h6 e8 N& ^) H& s  `6 c8 K

5 F' t( k, v7 r, z

d52d01fb8ce51133f92387cb24035cb2.png


6 t% F% P4 c6 A. T5 @' V% c5 K* j                               
登录/注册后可看大图

上篇文章(Python9:绘制中国/世界地图)的例3:绘制中国区域地图,包含海岸线

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


2 m0 |* K1 L, X  ~                               
登录/注册后可看大图


5 W1 j- Z, i- C7 Y
  • ' i1 a2 v1 B% L0 s, \0 w' `) T
  •   R6 {  m6 b  S8 b( e- T- \3 ]

  • : `* w8 t( b! u( n. @3 \

  • 2 X- f+ [/ F8 Q

  • 4 ?0 u  `# z9 {9 Q/ L; G

  • 8 `+ ]# s/ t5 m, o/ p

  • ! Q$ o3 L6 k# M- F% L
  • ' ]4 o# C: M5 i$ Q
  • 9 k' B1 C( E* z  B$ |
  • $ x  i$ V1 V; |7 V9 k

  • + Q/ U" K7 t9 F. t, X! D

  • ! H& h" X" ^1 G1 S

  • / d/ c9 k" Y$ V) }
  • 1 g, H5 N* L9 r. D7 V, y

  • 6 T2 [& f, j! V; _

  • " G  L' B) b/ _$ i4 c, [/ h9 W
  • & w' \' Z! i0 Y+ b- R4 Q7 G

  • ' h9 k0 c0 V/ S# D. @' `1 m2 |
  • % S! i5 ~8 P3 s9 c; i, q7 E
  •   S; r: m# ^6 T8 A* z6 F6 G2 f

  • 0 r% _( V3 `" M
  • 2 r3 X$ X- q5 }- e: A

  • $ Q& Y( q7 I/ G2 N8 D
  • : \4 U9 [  ~! z) @  w
  • ; M4 ?# V  p& J2 }

  • & l9 I) Y. D- `) _5 M
  • " @$ k0 i# m" x  {+ t% I
  • ( i% e6 w; a4 H7 c6 Q7 F

  • ; `+ m# F5 L0 W# H: y2 c3 k- L
  • 0 k$ p# q' T% M, t1 ]  I
  • 1 z7 O. H8 {9 z$ Z
  • 1 p) T, I8 W6 q0 B' {7 U
  • " [# N  J+ l( s; c/ W8 Z$ N  G

  • 0 z# [+ ~8 g, H0 r7 |# e* j9 p( V

  • ; ]! m8 }9 R0 G. {4 }) {2 d6 i
  • # p. `9 [. G/ i& i  F0 \/ p, s
  • ) R* P; E+ g. n

  • % j  C4 T) _2 B" s$ Q; o0 ]) \
  • + G0 u/ @7 h/ Z* G! l

  • 0 e. F2 h4 e1 ^* Z

  • # X& t/ m1 b! L, F, m* e8 P5 [

  • ' G6 X2 Z, p3 G
  • : |3 S: Q! ], I

  • $ W5 v; e! P3 A9 \0 B9 w
  • 0 C2 p6 O8 O! w/ V* K
  • / U5 [: ^. H; c/ D. ~

  • / R; C  s, k0 u! D6 o

  • % @; K; T* j3 t# t+ Y* y
  • - d: f3 B/ X% H3 [1 I4 P

  • 5 o+ N6 T. B" j0 _( ?6 g  S

  • 9 L: q: K+ \- [6 z( x; B& G  t" M  G
    - U8 _/ Y& I* G) `2 O. {; 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()
+ \3 C' h4 i: g


* G# I, g! {, a' {+ o
# K& c2 K6 q9 |5 h4 A# P- f) b

颜色表:

706a2a5d8e3967476b47c483a33ce236.png


# O. e. s/ D! ~) N& @/ k3 j1 F& O4 T: O# x; Z$ ^
f5c3384582002193c2c6c354a3d26729.png
回复

举报 使用道具

相关帖子

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
有风
活跃在2022-10-29
快速回复 返回顶部 返回列表