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

[复制链接]
' o/ A6 ~' I; M

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

    : e' |5 h- U- _4 P3 j

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

! U+ k7 l8 _2 L3 K2 h+ O

绘制中国区域地图

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

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

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

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

2708a4a2d4c2f8f69c367e9c20ffc43a.png

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

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

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

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


- @3 t4 Q0 d! i5 K. e8 p

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


  • % v+ I  Y. q: C* ~% j
  • , Y7 x; [( h9 }9 w- W& _# f+ x

  • ) u5 X" I. x( D$ Y9 {% c0 c& w
  • ( S/ L) C- L3 i$ N8 A" y5 ~

  • ! y* Y6 V4 O0 ]' m
  • 0 W: \, g/ m! Q) \6 _0 ^  e/ h5 I

  • 2 Y% V$ K/ Y3 R- U- ^

  • . a- r: }1 u. O! o

  • ; U/ l/ _! V2 z  n) M: p
  • 2 T5 p+ [! l* Q# v8 [
  • # D+ F" T0 ~) l* b6 t  i
  • ! ?! Q+ [6 X( U4 e

  • * `3 O1 R" T" z8 t$ d8 u4 d  a
  • 6 t8 q2 O) M$ f, k

  • 0 C. h6 P7 J5 j
  • , ^! x+ |1 ~8 K+ k; v6 d8 v# C

  • - R, s5 m- a% U# K; f% S8 y5 d

  • 3 s' n/ y6 O# p) @, o9 |

  • ( h& n# Y! Y  ]
  • ( Y* B: `. J# k8 J5 W/ o* h

  • ) T& Y) s* \3 q- m4 R7 J: k/ \0 W

  • : g; T5 w6 D& }0 R* d5 V

  • 7 `4 X/ M# @; H2 t* N- L9 y

  • : V: H0 Y: V* B

  • ) s: d7 U8 ~. L! F" d

  • 7 z) G+ T8 h& w& Q. j5 H
  • / p8 G- T- @6 G6 J! A! u! `* |$ x- ^% e

  • 4 O: y, p) y* V; z: m1 s

  • : _  P1 r6 U' J( m' g1 V2 h8 E

  • / \" h8 }7 Y5 w! }" V$ f+ p

  • ' L  W0 O( k- Z& C/ B" o

  • ! u) J$ ?4 Z$ O0 ~( h+ K3 {8 y

  •   s5 j- p1 }5 G* u

  • 1 f$ j3 o' p7 R& r
    * f, p5 c5 h! v- C+ G4 _5 g

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()5 x0 |+ k5 l- Y# d0 r  S


. z/ i5 _3 {7 r7 e# n. Y

d52d01fb8ce51133f92387cb24035cb2.png

! p3 z/ V  `- i( R) o; C* ~( [
                               
登录/注册后可看大图

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

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


5 h+ A+ \6 `( e2 {  i  S- D- ?0 G$ g1 Q                               
登录/注册后可看大图

, b6 s8 M* T$ _4 S( V9 L; n

  • , K+ _* l5 N( a9 J) c
  • * e. A0 N' b' ?: D5 \

  • 9 a6 r/ V: }" {3 U+ K1 |- C2 g
  • , n* Z, f# X' Z+ G5 o1 l% W9 p0 H
  • 5 l+ z$ m- Y; s8 L

  • & w+ A) P' c3 ~2 i8 @6 W
  • 8 p. S* r9 l; p; Z( Y
  • / a% ~( f* a- A$ |+ t) I. T

  • 5 N- b  i# p" \3 ]& A8 v6 w
  • % }; Z7 X8 g0 b5 o) E- M

  • - E+ p* Z& J! A2 l7 w

  • ) i* [, o& v# Q7 W' o
  • * e" r/ s! O2 s- `
  • 3 r% b* ^0 c! }* K
  • 6 j6 @$ {+ L' Q, Q

  • " C  n+ ~# U( f, X% B: M
  • 0 H2 D& n' ~7 l# _3 k- H. G

  • ) d; Q8 ?, @# _

  • 4 C3 U+ b+ J# y0 q, Q
  • - F5 i6 T" D5 A9 D9 X9 z

  • ! l& E) y% }  j. x" t

  • 7 p+ J! `% I8 [+ K3 [) i

  • 2 f1 s( q" b8 n7 h* z
  • . x5 L0 Q7 A. }

  •   J1 @4 w# G1 M8 D+ x

  • 3 M; C0 _  p  Q1 A% E
  • % P3 b0 h+ S) C: R' m: y+ Q
  • $ W0 n0 V" b8 e# ^
  • ) ?7 K" y9 ^5 E& f7 H
  • 7 D/ ?: B3 U; c/ Z0 G& _
  • $ v' \3 w* Y3 K7 [: {

  • 4 B& V9 V! p2 b. A, [/ K2 ?+ W' V
  • " R: h. h/ V5 {0 \
  • % C6 n  \2 r: t$ f! f5 ?1 S

  • - Q3 C  |4 H4 b. `5 C1 n" _6 n

  • , v+ h) j5 z" n3 O& n' I
  • 4 d$ A8 ~: ^1 y+ q* G7 C
  • : w' W4 c1 C2 W' e* X! W

  • 2 u7 L. l% Q6 @# e. h7 ~# ]

  • - \4 Y9 i! E4 v0 C3 k" P6 b6 u" b/ [

  • ' s1 n$ {# U. H
  • ( q) m9 ?; b: K+ q
  • 5 y. l0 z9 ^& D) l% A9 u2 n/ P

  • $ U" p4 {3 ~2 M0 `5 V

  • ! p4 z/ z8 G: D. u2 n7 ~0 ^

  • 5 [3 a$ J6 K" U& G1 z0 b& Z

  • ; }$ v5 c5 w0 `: d/ I5 `2 _. d
  • . c* H3 q! d) _, C7 `' n
  • , x- @" M! m& d6 T

  • , q3 `2 H$ v2 N; j* m

  • : b" c! r' ~( J
    : [) C& O) r" B7 _

##绘制有海岸线的中国地图#引库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()
! q( h$ Y; S$ m$ H; N. U


2 _* G' R6 B5 c( I7 z& O, q; r( S) V3 v* F. j3 n

颜色表:

706a2a5d8e3967476b47c483a33ce236.png


, O7 N; o% l- d# ]* p
5 C1 [" w: y( s+ e4 e! B
f5c3384582002193c2c6c354a3d26729.png

相关帖子

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