- H# J7 R3 n6 l3 i( q" I; G+ f气象数据可视化主要依靠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("···") 1 E+ q5 p' P) D4 K
只要按照以上六步,基本绘图没有问题!!!
0 q, t3 P" ^, |) Y7 {/ R绘制中国区域地图 前面python9的文章(Python9:绘制中国/世界地图)关于中国区域地图的绘制方法除例3,主要使用的是pyecharts自带的地图,其中中国区域的地图有问题,所以不采用这种方法。下面根据其他博主大佬的文章(python绘图 | 中国地图最正确的使用方式(九段线&南海子图&高分辨率地形图一个都不少))进行适当改进,并采用权威的国家标准地图比对无误的地图文件绘制中国区域地图。以一个具体实例完成: (1)正确的中国国界线,省界线,九段线的绘制; (2)全球地形图的添加; (3)中国标准地图文件的添加。
提前需要下载准备好一些必要文件: (1)正确的中国行政区划矢量文件; (2)全球地形图文件(从大佬提供的方式中获取的全球50m分辨率的tif图); 以上的文件可后台发送“中国全球地图文件”获取。 4 E$ h- s1 i( N/ E) {
(如果没有下载相应库包,请先在anaconda或cmd命令行去下载!!!代码的详细说明已备注在注释中,不再额外赘述。) - - ?7 Y5 T2 h' y1 ~/ y5 A
- 8 R# h1 V& P: V* K) f) c
$ B4 c: G: Q$ ^- h2 U
- P+ \- b# n+ g# O6 R- [$ K" _6 l; z( @) i9 o' F; f
- 9 C$ R0 j" H3 s5 O* F6 A1 N6 U/ d* }
- - |0 T5 @' j3 N# {0 I
' z3 L5 F& U) S4 c' ~+ T* M
, ~! P ?5 w" S: E2 M
+ m) Y3 U9 t, Y
& ]' W: w1 P: e) S r- # Y3 U$ X- C$ O9 ?# n0 e7 _
- % D7 N+ m8 F1 O+ b! |- p
- 6 f9 M3 R5 t& ]; q
# Y( f; r( D& w0 k$ l5 b
) P) X. K% U0 Q- 3 N8 a6 s9 \3 c" Y8 {
- & v- \. {+ r: ?% k# |1 c
3 I2 g6 U0 _4 y, s- # f4 A% v7 p' U
- # E6 _2 h: o r7 s3 q
; D/ |5 S2 [1 H& y1 q5 i6 j" S
8 J) ?% x! }: o$ ?
) J4 T( i0 h% Y# L c+ j) p1 {, S
( N) V' Z4 p+ [: ^% `. `0 r
" S. u9 y* H1 j" ~
5 N6 s9 `- Q. K1 Q% w, r. X# p
/ W6 R# c. U; |* F% h( M- * m) m8 V) E! ?# }' i4 Y; S( T
- $ `* C# n9 s: c9 R% `5 z+ x
- 7 s8 |& z2 T% V9 E# Z6 Q" D
- - U7 J( s3 s' u% \0 G
) L0 s( e8 A$ ]) P1 S: r
9 w6 {2 \% H0 d! f k1 N0 C" E
& @8 I7 T" { J ]4 `
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()2 Q$ n- y! e4 F
- s2 M/ Z6 `; ~ P) e4 ^: j
上篇文章(Python9:绘制中国/世界地图)的例3:绘制中国区域地图,包含海岸线 (如果没有下载相应库包,请先在anaconda或cmd命令行去下载!!!代码的详细说明已备注在注释中,不再额外赘述。) " o6 g3 R* E& e! [. J; U7 Y
- # U/ f" G/ M' i7 k: A: f
- ' b3 g# I7 j) j t3 L9 ]
( l, o* H, H* @- - I; B! }# I9 E0 v m4 d& ^
- ) D* x3 \1 D3 T' p/ ^: ?
- ! ?$ w3 B# S+ D, L* o+ k
- 8 F6 q- ]# S P/ c: Z
( M) q) P z/ |- u
8 D6 U( S5 @0 S) \2 C* }
# g5 L) E4 G5 @8 }- ; E3 o# T; n( {) @
- 9 C6 z/ p3 `2 F* s' _6 i8 b, p- _
$ m" ?- T' u6 {/ f) o/ M- / l' B- W7 p& z# @# J
- 3 D: l5 D* f5 c0 q, [% r# C
4 c1 r2 G0 q) B& Q5 Q- 7 Q/ g; g F2 ^+ o4 ]; M
- * [& a7 Q& d8 q1 L9 w
- 5 Q, f! D: v% w
' N1 W5 o$ ~ Z% b; F2 j: v
2 x: G1 @! Z8 L0 M- {7 @
/ e8 Q, A2 T. `1 d; n1 J- " U1 o9 w7 V/ ?6 k- ?3 X. M
- 1 U/ w4 e2 j0 o/ y
! W2 o- H, ]4 j2 K- B3 h
2 {) b0 Q, @* \/ @! V9 N; ?- 7 R' @* n: Z# l
- 0 [* a) L& r3 ~" x3 k
- ( q8 n+ A) B9 }1 A
3 Y+ U6 B+ X1 d A' [, K/ Q* A. W/ ?, b
# m& A6 [) ]; {$ U+ ], |+ }8 X: l6 A- % X o" r5 E+ v+ R) a( {, \
- * G* @& Y" E3 t) _
- 4 {" x! F! r( p3 x0 s4 I7 p
- , A- P+ i- k1 M5 k/ ?" G0 b
- , r1 B1 \, \! P0 z( B! q% T# ~
% K- N( C6 B8 O7 H% S- . \9 X0 S# a, v7 P3 L% u
/ l, k# i% B9 {9 I
5 K& @8 }! q& L/ J. t+ T" e& B+ i
+ R+ Q3 `+ x P& k3 E- 9 Q. c% C+ j3 c: U7 O @
- $ v' Q: `" R! f& x9 s5 M
- & q7 i# A2 k5 Z" l
- ' O @6 O* v6 m" w. K) ^- [& r% d
) Z L. Y6 y, K5 X( Y: a6 M- ( K, b- u4 g( q) ` Q. R! Y2 k4 y
- + y; {3 z/ P, ]
- - N/ ~# J: ]% h( O, j" ]
- 6 B/ Z- Z+ N" e7 R: U: z6 V
8 d7 b4 x8 z! E7 J3 g: M7 ~7 E
' L. @" u( p. R1 M" S
##绘制有海岸线的中国地图#引库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()
4 T6 f& C+ D2 m
5 w' R$ a4 N8 R/ j, u' }( ~, ~ Q; s! t
颜色表:
; D ~9 g. j" a; \8 \
5 q5 d) _, k: f |