* K3 S- K, t5 Y4 B& X
气象数据可视化主要依靠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("···")
+ q6 x4 P5 M7 E4 D0 t4 E( k
只要按照以上六步，基本绘图没有问题！！！
" m' Z- l( M0 X绘制中国区域地图 前面python9的文章（Python9：绘制中国/世界地图）关于中国区域地图的绘制方法除例3，主要使用的是pyecharts自带的地图，其中中国区域的地图有问题，所以不采用这种方法。下面根据其他博主大佬的文章（python绘图 | 中国地图最正确的使用方式（九段线&南海子图&高分辨率地形图一个都不少））进行适当改进，并采用权威的国家标准地图比对无误的地图文件绘制中国区域地图。以一个具体实例完成： （1）正确的中国国界线，省界线，九段线的绘制； （2）全球地形图的添加； （3）中国标准地图文件的添加。
提前需要下载准备好一些必要文件： （1）正确的中国行政区划矢量文件； （2）全球地形图文件（从大佬提供的方式中获取的全球50m分辨率的tif图）； 以上的文件可后台发送“中国全球地图文件”获取。
. Y9 J$ j `2 m& F. Q（如果没有下载相应库包，请先在anaconda或cmd命令行去下载！！！代码的详细说明已备注在注释中，不再额外赘述。） - 3 A$ _/ N& W- ]9 r6 i! X, U
- R: f: K O# w- a4 J2 P: x1 k" p- 7 U' v+ |5 y; u1 O- |
- o" J2 V2 }0 G; |- 9 F0 G0 D% B7 c. {. v$ s2 \. I
; {* n. u' |& p- $ Y6 G P* p2 [) E2 P! u
- : H; @8 {& \( ~
- 2 i+ L) l6 S3 s% e5 y4 S
- ) M) A. v6 d4 j- v3 P
C) x6 X6 e: m, Q- k- 0 Q& Q9 J% _! z6 c6 G
- 8 w6 l2 }- Z) f; E( f& K5 {- ?
3 x! M8 _5 L0 I3 C
# m: b' L: D: t [
& v$ p1 w0 K' E6 X/ r P- o- v
: K3 n; E( l3 G4 r9 U
( u& f. u' ^& Y0 `' d0 ~- * m1 B1 Q; O6 L9 q$ A; H
- 4 v, j$ z4 j( j! q; a
3 z9 p) Y. w" \ u9 z& ^& G ~
7 _, I0 ?+ X) Q0 j" u# a: ^) @
, m+ n5 M# x* A( l, B5 U- # t5 l$ W9 v$ Q; h, }% S
& s G% Y! R9 D+ a2 A6 q* D1 S, t& A- 6 W x8 v; M) q2 U' B1 N
- _3 ~5 ?3 M( X' t- D0 O- ) {0 |1 S# b$ `. }: V
- 9 l% r) m3 k+ v% h: t1 \
- / \3 r, a6 M% c. o
3 @4 x2 C% \- r; V! I- 1 E+ y' A' n- s, e
5 K" t. y" N( K5 R- a) x7 a4 ^- 4 [! }- j9 z- A% q. E" D
, g5 q6 |3 u, s5 z$ p% v3 K! w
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 u7 g) U% u5 g0 A% V6 v7 d8 o
2 E: f. T/ W, b
上篇文章（Python9：绘制中国/世界地图）的例3：绘制中国区域地图，包含海岸线 （如果没有下载相应库包，请先在anaconda或cmd命令行去下载！！！代码的详细说明已备注在注释中，不再额外赘述。）
6 d3 t7 I& A( D* V* R0 T
$ U2 u4 C' c( m/ l: X
" @5 t: w* k, E+ H, y
) Q. Y( C( w+ k/ i* `' T
. ` X7 t3 M: Y% Y6 a9 r% H- " g- M* x* r. x* _: U8 f3 _
' y0 M6 a' s. f. o' {- ( d; a% y8 T$ A
- 7 T7 t7 Y1 W. Z8 ]/ C* c
- 3 {' u8 E8 Y$ l" o# g! |% r. K
- : H: z/ k' g2 }. w8 D
- $ @4 y3 G0 Z S2 B! W
- , c! |4 Q; T. h
]! D/ {, e, i; _- % h( z/ Z, r% n8 y* M+ C" y+ O& e
- * V' D; H5 p* b. g
- # ?0 n' w; G1 z- [
- * n j' e2 J# M5 Y, T* q0 ]7 U
- ' c1 I$ y7 E% e! G
- 3 ?" |4 t ^, S, ~9 q
- 9 v7 f \) Z3 b2 b% q0 r# x: O/ M
% v' w' S% Z9 b7 j% s( Q" m
7 H& u' U$ @0 t8 o
* V8 }: K9 s9 D7 h4 g
! |5 ?8 o4 `* Z* j
8 v4 F3 @1 f5 q; h( Q% Y/ ~. |
' Z0 K# f) F9 A: U( f' i! U
( k3 F- n" `) ?! Z7 ]+ @! c7 Y/ `
: a$ C5 u& T3 O2 i2 {) A- h2 u: l- ( s+ ?$ ~; }2 Y8 D0 U: q
- & v1 g2 m$ Y' K# o8 k9 ~2 [* p
- 3 H- g, b8 f) J6 x: H. `* ~% j& y/ v
. u- ]1 i) k3 H9 k5 m
4 u, [$ F5 S! H }- Z
) \: }0 q f" L9 j. k
9 a- r7 h, b3 h/ c7 z# ^) U% d- 7 F' h2 ^% q9 K, c1 A
- ( t" y2 {* u7 c$ G: i
- ; Y- z$ ? {) l3 A
/ w1 g t# o( H: Z3 f- 1 D! T8 k6 `. v
- $ S2 P# `$ p" z* x
- 0 n O3 {3 G6 y$ ]0 s( K! [. U
- 9 t- r* ~% [) D# a i
- $ W' B) Y5 _( y# H8 U! H
+ c2 h2 Z! n: D2 |' r1 q" n- 2 m5 ?- W m: n. [
: f" k3 K$ P# T3 X4 D, S5 f
2 |8 K' ]' @$ q: g3 {. f& G
9 \% h# U# t; t% |6 Q3 L% J% {( U- % |! u4 f q z8 T6 P% Z) R/ K
- 2 w- H H$ n1 H
9 h: B- S# b$ c7 v: r% _9 H; k
##绘制有海岸线的中国地图#引库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()
. @8 h% L) |1 O" j+ a
: Z W' q1 _7 F! ~4 ?: `9 J3 Q5 T. x. T0 m7 F8 s/ Q+ X
颜色表：
% U2 k& U ~6 Z7 Q5 O3 q5 u
+ ^8 _1 O& e' K3 ?6 i# u |
) i5 v( w% C( M! y& u) p
$ c+ n F" i3 y, p% [/ y
