% ^* V! n2 C: T, p) K3 l8 z
气象数据可视化主要依靠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("···")
' |" B6 P1 E4 J c$ y
只要按照以上六步,基本绘图没有问题!!!
+ ~/ L! P& D8 R, J# y绘制中国区域地图 前面python9的文章(Python9:绘制中国/世界地图)关于中国区域地图的绘制方法除例3,主要使用的是pyecharts自带的地图,其中中国区域的地图有问题,所以不采用这种方法。下面根据其他博主大佬的文章(python绘图 | 中国地图最正确的使用方式(九段线&南海子图&高分辨率地形图一个都不少))进行适当改进,并采用权威的国家标准地图比对无误的地图文件绘制中国区域地图。以一个具体实例完成: (1)正确的中国国界线,省界线,九段线的绘制; (2)全球地形图的添加; (3)中国标准地图文件的添加。
提前需要下载准备好一些必要文件: (1)正确的中国行政区划矢量文件; (2)全球地形图文件(从大佬提供的方式中获取的全球50m分辨率的tif图); 以上的文件可后台发送“中国全球地图文件”获取。
# \1 x2 A9 F) U1 ^' ]0 M(如果没有下载相应库包,请先在anaconda或cmd命令行去下载!!!代码的详细说明已备注在注释中,不再额外赘述。)
: W& j) q: K% K" U
$ u0 m, ^, L" c ^0 U- F
! M% I3 L8 X% a# I: x- 5 S# t8 ~9 l- d' N ^( c
- 5 a$ {3 z* q. Y! }
- * R6 t* t: z" g' f5 ^
9 G8 H5 Y( N9 Z7 N+ }" L9 N( t. R- & {, g$ R, W8 s: [# I+ o4 v( @
- R1 @6 g5 @, g/ r C
3 o2 c" A) U$ ~, N4 E- # o3 B1 D" T9 x5 b
. I) `1 _/ [" T/ W6 z) b
$ x& y0 G% D. j. M
: [5 W% l! c5 d; z& R6 w& l3 z- k& b/ h4 i4 y( ]; z3 d
- 4 Q) Y+ h, e) v6 z+ N2 N" r% z3 F
- & A7 h+ F9 \4 }; H: n: M
u) }: }5 S @# O& {- ' H9 h5 M# ]9 C0 E4 ~
- + n0 ?$ S) m6 {" q: K- I+ H
- 2 B) t3 M2 P* T, n5 M* u
- 0 A _: a4 P% D" j+ k1 P7 D
& s( _' B! P; a1 b/ X5 V- - P' m2 }' Y* [/ P7 \
- / g# r& D9 O/ F# p
* W! r. }& O$ H" L2 p
5 v8 Y0 v8 i1 A" B. D X. W' _- 7 e. d9 V' `: N8 t0 P5 v6 k
. s2 B1 Q1 `+ y9 ^" J$ J- B/ B+ A h- - c: s! ?' K4 \$ Y6 k n: ~
9 o9 t _" t- F* `" x1 N
7 f5 B: A1 L4 [; O# M X- + v9 _# t8 \/ w* w$ ^+ r
( Z( @6 h! B; ]2 D- L3 j) L6 j1 [' K. s& [8 M
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()/ D$ I$ F$ X/ F
% _8 I, V- i% h* I
上篇文章(Python9:绘制中国/世界地图)的例3:绘制中国区域地图,包含海岸线 (如果没有下载相应库包,请先在anaconda或cmd命令行去下载!!!代码的详细说明已备注在注释中,不再额外赘述。) ) U5 e1 G! [' A2 ^3 M
- 4 y) j' T, Z) h6 \4 `$ q) G. y
& U( e3 {/ ]1 p0 p/ i! l& g# C& D
* ], R; P8 p0 U2 \* a) H1 }- : }/ F, r2 l4 ?" s$ `; M
- 3 h$ }% s/ Y6 F3 w! p9 ?1 k' X
- 1 A r, k& g9 v" y' G
- & s5 E, a5 ^3 }. M. H2 I
7 Y$ V. l* v6 d$ j+ A/ b+ D% c- 3 e9 C$ J5 [* _; C" b$ n
+ t) W9 q! m2 c/ L3 g
% T' R y5 F% Y& q- / g+ g/ B. Y8 {: Q7 X& y
- - i" n) Y/ _9 q, k; l1 b; T6 y3 E
) j$ D/ N8 r- S
9 ~% ~/ ?( j+ F$ \; ]- s
2 g8 I8 J7 E0 f7 ^& j% {7 Y# o6 O- + J$ K' g9 g5 k( s- r
: @: v2 D: o" u- r
3 u/ H% l# k( U- I+ \- 6 t6 c. m+ {' D& b- _- r% x. Y
- ) o* p [0 m* \# G t1 u
- 6 E/ ^; g+ V) X Z; k
7 V5 q' `6 z* f+ g9 P) h$ ^- ) |$ w5 h: S# Q& V
- 9 ^( _- o) H0 l7 x; {
3 \4 |/ z4 y/ H6 `1 S$ n) e! {' D- * J; c. h) `- @. B% h
% b$ O* V2 c/ V/ x1 B. G" _- ! r5 ?& l) Q" h% _) V3 W) h6 m( t
0 l# t3 r3 g) R3 c! r2 V2 a/ T- ) w* X2 e9 M( r
- 5 ~' K3 O# H& R9 h; l
9 u& ~; \. Y4 s. b; h0 f
, O) `; B6 `6 g/ t; Q$ f) o- 9 B! [9 c* d) @3 P6 U
- + k; e4 Y+ ~, R/ x: @' K
' D# Y+ ]. O2 J
' n. f( D4 m- Z$ v( b% K- 5 j" c! V w5 ~& ?: m/ b4 u5 O" X
- * ^% L# ^. R% l) C( G% K
/ \ |* W% I. T- ' ?0 q4 ^; Q: s! h* }
- 3 G! h$ B* X4 N
$ Q5 O3 x1 g2 ?# B- 0 v; N# a, S# w+ W+ R% p& {
' L7 H4 {0 g0 h S5 Y/ i- % |* ]2 {+ d* j# D
- - P0 f, a* s5 J# L
- : E" o( a6 G. {# g- i+ _! W8 d
- ' f- W: x l% Q! ^6 R% P
- $ w) O3 }) [) Y! U/ G$ H, ]& K/ {/ o U
, c6 W) g8 L9 Y4 ^$ u
##绘制有海岸线的中国地图#引库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()
1 E' G2 K! ]8 ]) ?
1 U) `/ M3 ]- b. l' b! T) O. `" z% L P! h, z
颜色表:
. h4 S& w" j' ?
# z/ P0 g$ I7 z0 T2 K3 t, w
|