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

[Python] 【气候软件】Python8:绘制面积图和填充图

[复制链接]

气象数据可视化主要依靠matplotlib做绘图,其他库包为辅助,包括数据处理的,地理信息处理的等等。


. Q) H4 t2 N5 _' |  y

绘图主要有六步(“六部曲”)(不是绝对的方法,视具体情况使用):

  • 引入库包: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 Q5 b) |: L: M( A; k" m

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

3 v9 @2 y/ E! u. i7 ?8 m

绘制面积图和填色图

2b2b8cffc67a22fcc69646ca599a3645.png

(有的时候需要画这种交差的填充图,能直观地表达差值的变化和含义。)

014dfdb505f25bf97af4e4013898632f.png

下面直接上代码(关于代码的解释已在注释中详细说明,应该能理解):

  • ( K) J: K; R# ~# ~4 |. D
  • . P* L' k' H  R& E1 U7 G

  • 5 i0 P2 |0 X7 a3 b. k& M4 G! b4 o
  • ( ^- Q) W# O" m9 h( N% n) C0 ~

  • 1 D- e3 o$ v0 `: ~7 V8 [$ K
  • & C6 w- w" W/ {3 F" S7 Y

  • 0 t0 m1 h# t3 z
  • 5 T% t7 `" F0 R" x; ?  P3 V& _

  • . H$ v& l/ a' _( X/ o2 H0 X( w1 B- F

  • ! ]# w& U3 Y& m* Z' O$ n

  • ; J8 S$ k  i. {% x$ ^4 v

  • $ i5 c. ?" h1 d5 d* [+ \' n
  • / d; A: T% X% `8 [3 {
  • 1 h: F8 f, J+ s# g

  • 5 L3 E8 u: C& g

  • 7 I9 l, c8 U# ^, Q: c- u/ J
  • ! P9 J# V; ?7 ~& i( k
  • ( W/ A0 }* D( ~- o% i$ N" q
  • ' K% K0 D. O) y3 k- Z; |
  • - ?3 g$ M" X1 c6 z, h, V# Q
  • 8 f3 f, ^& c  X" K7 A
  • , @2 _7 |& N/ o/ l6 [  U
  • 3 e4 I: J1 y# |2 A; Q" T

  • 0 U; G  N( k+ A

  • 2 z- b7 P% g# s. H" ~

  • 3 N) {8 l2 B/ e! k, Z( f. y
  • $ Y$ e0 m8 A- ]. V

  • 4 j* ?" W' L3 r" A- [* q6 {

  • / @1 Y* k  E' Z2 t
  • / X: `; v. {8 ~$ `' D5 Z; t

  • ! @& q) k, T( c

  • ' p0 y, B9 m2 n% I$ D0 M
  • 4 i" w/ I0 R* Z1 T

  • . q' m, T$ Y3 t8 I4 x  s( k; `. W2 i
  • 1 Z/ j, X1 U: K6 ]+ i

  • 9 k) `% O4 \8 z6 B$ \0 w0 I0 @! A' S3 d
  • 2 ~. X$ V# J/ v2 f

  • 2 ^0 Y6 ^# t  d
  • 5 x) s/ l) a) s& n% E: C2 m7 s/ d

  •   D/ S5 z0 G# E
  • 5 B* v. g& R8 ?$ M  D' [  M$ Y4 o2 U) b

  •   [9 @9 W. ~) I" h' [4 g

  • 6 j- l( F* c' g0 Q4 h: D
  • 2 ~3 U, W# ~5 @- K/ _% A% K) s5 K

    7 ~" E; O. ~  x" ]

#引库importnumpy asnpimportmatplotlib.pyplot aspltimportpandas aspd#中文及负号处理plt.rcParams['font.sans-serif'] = 'Microsoft YaHei'plt.rcParams['axes.unicode_minus'] = False#画布fig = plt.figure(figsize=(15,6), dpi=200)ax1, ax2 = fig.subplots(1, 2)#数据x = np.arange(0.0, 2, 0.01)y1 = np.sin(2*np.pi*x)y2 = 1.2*np.sin(4*np.pi*x)#绘图#子区域1绘制限定区域填充图labels = ["Fibonacci ", "Evens"] ax1.plot(x, y1, ls='-', lw='2', c='k', label='$y1$')ax1.plot(x, y2, ls='--', lw='2', c='r', label='$y2$')ax1.legend(loc='upper right', fontsize=10)   #图例位置#ax1.set_xticks([])  #取消横轴刻度#ax1.set_yticks([])  #取消纵轴刻度ax1.set_title('fill between where')  #图标题ax1.set_xlabel('x')   #横坐标标题ax1.set_ylabel('y')   #纵坐标标题#填充图还支持限定条件下的区域填充ax1.fill_between(x,y1=y1,y2=y2, where=(y1>=y2), color='skyblue')  #对y1>=y2的区域填色ax1.fill_between(x,y1=y1,y2=y2, where=(y2>=y1), color='yellow')   #对y2>=y1的区域填色#添加标注。xy:标注箭头想要指示的点,xytext:描述信息的坐标ax1.annotate('y1>y2', xy=(0.3, 0), xytext=(0, -0.5), arrowprops=dict(facecolor='black', shrink=0.01))ax1.annotate('y2>y1', xy=(1.7, 0), xytext=(1.8, 0.5), arrowprops=dict(facecolor='black', shrink=0.01))#子区域2绘制堆积面积图ax2.stackplot(x, y1, labels=labels, color='yellow') #堆积面积图ax2.set_title('area picture')  #图标题ax2.set_xlabel('x')   #横坐标标题ax2.set_ylabel('y')   #纵坐标标题#出图plt.show()


0 u7 q7 f6 L$ U. ~
回复

举报 使用道具

相关帖子

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