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

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

[复制链接]

0 y% h! f  x: N

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


    9 L$ A4 e2 d% H! E  E

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


) [) ?' G4 R  P

绘制中国区域地图

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

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

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

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

2708a4a2d4c2f8f69c367e9c20ffc43a.png

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

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

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

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

( u# I0 e+ E8 t  p8 A

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


  • ' o. y2 Q1 @3 k5 j; \) D; C( d" D
  • & x1 @) ^3 I0 t+ [- L! \

  • 5 M* |  V/ t9 C
  • , S* D: `( P8 y+ D- w2 o

  • ! _9 L( b+ |6 R" M) b. X

  • ' M5 G$ x' M' Q$ L8 e& k

  • ! }! _9 ?& }$ J" |, C5 t2 R4 l3 K

  • ) C/ W$ k  [1 u$ N0 W0 L

  • # U8 ^8 \  m! D! L8 Z: O6 _
  • - R$ S, E$ \$ p. v* a

  • : ^6 D, K( f5 c7 I/ E" Q/ D, f, A
  • ; N: r4 n8 Q. G7 g% c+ u9 i, \
  • 7 ^9 R$ P, l- e8 N1 o% x5 z9 C

  • % l( U8 {$ q) H4 X# p9 H# J

  • 6 z/ b4 |2 B) f" x
  • 2 \7 h+ ?# O& ~/ D2 o- _0 t
  • 1 ~; ^2 l0 n9 _+ x1 j5 q
  • 6 U0 o' O/ X$ T' A

  • , d- T$ \: \; o/ v

  • % G. x0 f( \3 R$ B
  • ; x1 j5 o& Y  p& p  b2 o7 A
  • # w9 i$ V3 V. c+ Z$ ]. g
  • ) e) J0 Q( p8 y! M' E3 ]) a
  • 5 y, H6 [9 X' g* |

  • ' S0 O9 [0 E; b3 f/ p. u

  • ' X; u: h. l: j0 ]. B! O" u5 H0 y$ {
  • 4 s, E; r+ w$ Z8 Q! C

  • ; Z$ N' Q% x/ \4 z
  • / a: C- j* Y3 }! N  n
  • 6 V& ^5 @6 b9 s- Q9 b* C. E
  • 9 Z9 K9 P+ }$ F( M
  • . V2 p9 q+ j2 b# ?: y( J4 m
  • : v3 }' d5 u" g: U" M6 [. ]. h
  • % {3 s( h! W  Z+ @
    9 e0 s$ n4 E: c  z% Y6 k# v& H

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()/ c" Q- y- B# E9 q


; S- Q. x7 {( C1 t- q2 \' ?5 F# w

d52d01fb8ce51133f92387cb24035cb2.png

' I9 U5 m: w6 v7 }6 Z5 A
                               
登录/注册后可看大图

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

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

: {/ c9 P3 w7 t3 N& D+ r
                               
登录/注册后可看大图


& E* _  d$ M2 b' |$ z' D+ O
  • 3 s- c0 x, A  ~) a+ \
  • + H; M( T* }8 }8 h0 I0 b1 @" x1 H4 q

  • : l; V  \& W  i7 F$ u6 |

  • 1 Y! A: C' i- W" x
  • * c. M% {& {) x

  • ; p: Y+ S4 @( I9 x: ?9 W! O
  • ( |2 Z, D7 O; c  {) @& T

  • - W" s+ z% m: x0 G$ |; K$ L! f

  • ' |) D! N8 t2 ^, z& z

  • 0 V2 J; f! V$ D+ R; M/ Q, v0 b  H$ O9 v
  • % n! I9 Q( S5 K5 ]
  • ) d- v( \2 o; J+ [  F9 p
  • ; c' E8 k* m/ W+ N& {8 ]

  • ( W7 _/ c! h4 l, r  ^, _
  • 9 W! g; w0 {' H, D! B
  • 7 H: e' ~: w) k8 U8 U7 W* U" R

  •   [+ y& q. \0 k# w" T2 e* v% v( ?7 @

  •   e- P2 \3 j) [8 {

  • 6 H* A6 [; X( ?! u8 o
  • # z( y. v- X0 `- B

  • * f  m2 ]5 b  S: j: Y7 y% J9 y  P) d

  • + _' ~7 u; h4 K' [& C0 g# A
  • , a; b$ A+ _& `. O

  • ) D3 ^9 q% Y/ x* r2 w' Z: g( s

  • 8 o" U0 s2 I, H5 U( F3 P

  • $ \& W, D' x( a+ F% f1 ]3 L7 B

  • $ n' K# G% k+ k0 G

  • ' _- y5 g: t0 N8 l) y

  • * G6 |8 X, K+ B& P- W# r; V

  • : G7 b& l4 `- y4 L8 b7 Y6 C' T
  • 2 `- {6 p3 q3 z' ?  P7 L
  • 3 H# m* Y" E8 u. v) v2 I4 Q/ R# c
  • $ ^& V/ }: b" Y5 a, }% B5 Q8 H

  • + A0 g( u+ z# s9 A: ?/ x' H

  • # T4 E# o' b8 v" s, ?- w

  • * U. e0 U% |6 p6 ?$ G& A9 ]# h# @

  • 7 a! ~" O- \/ H& B8 M
  • / |3 [3 ~# F0 `2 o

  • , p% N+ J& s" r% |/ `' S8 s

  • 2 y( [% U: S3 f, ^. K
  • . W6 ^7 a3 l7 ^4 S5 r
  • % b3 H# p  `! J+ X

  • ) |8 q; ^- l% V2 p- T3 W% t& {; c

  • 8 {3 B3 a2 R9 ]$ p7 R8 B$ l
  • 9 a0 @! w/ d  K1 g" u
  • - D- X/ `  y; F8 v% ~0 i1 n, y9 V& f

  • 5 C7 q* C+ i6 Z$ B5 m6 Y
  • / h2 @1 b, \7 a1 v
  • 9 Q; P* k& u! u
  • $ @- \9 a* ]2 I: V

  • ( t* _' i/ r  E& `, |
    : d' Y! F! a. V% i1 p# Y

##绘制有海岸线的中国地图#引库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()
6 U; b% t( e/ Z$ [9 n

/ t% Q$ @- y' q; Y3 Y! W4 [6 X

  g8 a* L) b" c- n; e% ]( L: S4 x

颜色表:

706a2a5d8e3967476b47c483a33ce236.png

9 l: i% g7 D# x7 K& J# G

- _8 u/ h6 H! y8 h1 v1 P1 t
f5c3384582002193c2c6c354a3d26729.png
回复

举报 使用道具

相关帖子

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