数据可视化之美 -- 以Matlab、Python为工具

9 T" e( c4 z3 A4 K" v+ S! a2 y2 f: G

在我们科研、工作中，将数据完美展现出来尤为重要。

9 _; _7 M/ u9 K; _0 k

数据可视化是以数据为视角，探索世界。我们真正想要的是 — 数据视觉，以数据为工具，以可视化为手段，目的是描述真实，探索世界。

下面介绍一些数据可视化的作品（包含部分代码），主要是地学领域，可迁移至其他学科。

Example 1 ：散点图、密度图（Python）

' m! U; v# \6 _; x

import numpy as np

6 X O- A# c- Q& [- [

import matplotlib.pyplot as plt

# 创建随机数

n = 100000

- `' ~* i; i6 l/ A& }

x = np.random.randn(n)

" `4 R) M, e( ]

y = (1.5 * x) + np.random.randn(n)

d7 R5 N4 i0 C8 P

fig1 = plt.figure()

- P5 P4 Z% f5 Z- j" M# W

plt.plot(x,y,.r)

. e2 `5 b8 c; `, Q4 \

plt.xlabel(x)

plt.ylabel(y)

plt.savefig(2D_1V1.png,dpi=600)

nbins = 200

H, xedges, yedges = np.histogram2d(x,y,bins=nbins)

/ t5 Q/ W" i- \6 y5 z- H& W

# H needs to be rotated and flipped

( X v+ c) |9 J) ~

H = np.rot90(H)

4 x% R. V/ J# v0 ?4 N! l2 q- o# F

H = np.flipud(H)

: n; d* g+ ~) a; ?6 D8 \5 @

# 将zeros mask

+ K" ?3 O. G |7 m/ [$ H) G/ F

Hmasked = np.ma.masked_where(H==0,H)

# Plot 2D histogram using pcolor

fig2 = plt.figure()

) l4 S% D( z- P/ G$ f, B

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

plt.ylabel(y)

3 M( i3 F, r9 p' z

cbar = plt.colorbar()

cbar.ax.set_ylabel(Counts)

plt.savefig(2D_2V1.png,dpi=600)

plt.show()

" i. m2 n a% b* W) w% U

3 L4 H) v. x' y: u% [, T

& x' t9 ^7 b. ^; R$ J

打开凤凰新闻，查看更多高清图片

, P; ]. M# i2 i% h& Z9 {& i

; f2 h4 b+ k e4 W9 h

Example 2 ：双Y轴（Python）

import csv

& c9 U- ^( r. a$ r. I

import pandas as pd

import matplotlib.pyplot as plt

( z3 D9 P( Q% B0 w$ A0 E5 G6 X

from datetime import datetime

5 c# X& i) Q- U- H; a. s6 [

data=pd.read_csv(LOBO0010-2020112014010.tsv,sep=)

) h V5 D) ?" W' y! F. j

time=data[date [AST]]

$ }/ A% q$ c# G3 B/ U# Z

sal=data[salinity]

tem=data[temperature [C]]

2 Y! z8 I. ?6 t( {

print(sal)

. f8 ~8 [4 }; k, V

DAT = []

- n: ?& W) `. a( o/ b

for row in time:

5 V7 P; ^8 }7 |1 ~, l+ {

DAT.append(datetime.strptime(row,"%Y-%m-%d %H:%M:%S"))

#create figure

fig, ax =plt.subplots(1)

# Plot y1 vs x in blue on the left vertical axis.

plt.xlabel("Date [AST]")

% P$ T* Y5 ]0 W# `3 F; A

plt.ylabel("Temperature [C]", color="b")

plt.tick_params(axis="y", labelcolor="b")

plt.plot(DAT, tem, "b-", linewidth=1)

, l) A ^/ I4 d* Q+ z! b

plt.title("Temperature and Salinity from LOBO (Halifax, Canada)")

# }; k! y% y$ Q, d

fig.autofmt_xdate(rotation=50)

# Plot y2 vs x in red on the right vertical axis.

plt.twinx()

plt.ylabel("Salinity", color="r")

plt.tick_params(axis="y", labelcolor="r")

plt.plot(DAT, sal, "r-", linewidth=1)

7 d0 y9 O9 N7 T- @' ~

#To save your graph

plt.savefig(saltandtemp_V1.png ,bbox_inches=tight)

- ^3 ~1 [5 O- h

plt.show()

4 E' y: V" {$ J- ^4 \& q! w0 P

Example 3：拟合曲线（Python）

import csv

import numpy as np

, \& ?% C9 }7 B

import pandas as pd

from datetime import datetime

import matplotlib.pyplot as plt

8 X5 n0 H% ^* ]1 M$ e% N

import scipy.signal as signal

# ^; c5 h1 z" h

data=pd.read_csv(LOBO0010-20201122130720.tsv,sep=)

time=data[date [AST]]

temp=data[temperature [C]]

( n# I' C* c0 Y1 T

datestart = datetime.strptime(time[1],"%Y-%m-%d %H:%M:%S")

2 B4 F& n3 i: N* r

DATE,decday = [],[]

: B0 n: H6 h: \* m7 |$ O# E& Z

for row in time:

daterow = datetime.strptime(row,"%Y-%m-%d %H:%M:%S")

5 O5 o0 N. x- G% q2 Q# r4 O

DATE.append(daterow)

decday.append((daterow-datestart).total_seconds()/(3600*24))

# First, design the Buterworth filter

* M& W% k. J% z# o$ Y& o

N = 2 # Filter order

! M. d& V& D9 C. M' ]

Wn = 0.01 # Cutoff frequency

B, A = signal.butter(N, Wn, output=ba)

$ F8 e2 a9 h" v- e" P

# Second, apply the filter

* }) X* Q0 u' H& ]: k/ [

tempf = signal.filtfilt(B,A, temp)

& n7 j& N" r( T' R c }

# Make plots

. H1 V$ I6 `2 Z, C* V0 p

fig = plt.figure()

ax1 = fig.add_subplot(211)

' Y1 a: K, B* ^' J

plt.plot(decday,temp, b-)

plt.plot(decday,tempf, r-,linewidth=2)

plt.ylabel("Temperature (oC)")

" S- i i$ a% o( x: q: x& B' ~+ k/ Q

plt.legend([Original,Filtered])

: x! M9 h$ D* e

plt.title("Temperature from LOBO (Halifax, Canada)")

ax1.axes.get_xaxis().set_visible(False)

6 V+ t0 M0 G' V

ax1 = fig.add_subplot(212)

plt.plot(decday,temp-tempf, b-)

- P2 Z; b, ~3 r$ j/ L

plt.ylabel("Temperature (oC)")

plt.xlabel("Date")

, z5 x1 k/ \$ p- j9 T! w- I3 p& Z# |

plt.legend([Residuals])

/ J- p2 V) ]0 T1 M- o

plt.savefig(tem_signal_filtering_plot.png, bbox_inches=tight)

plt.show()

! ?$ [* F& v" R2 g- B

Example 4：三维地形（Python）

# This import registers the 3D projection

4 n p8 q7 P N" Z8 H6 l

from mpl_toolkits.mplot3d import Axes3D

from matplotlib import cbook

from matplotlib import cm

from matplotlib.colors import LightSource

9 q' R5 y0 I+ m6 |7 h

import matplotlib.pyplot as plt

import numpy as np

filename = cbook.get_sample_data(jacksboro_fault_dem.npz, asfileobj=False)

8 D [$ m2 m1 u8 m

with np.load(filename) as dem:

! `, B8 i" p7 C1 }

z = dem[elevation]

nrows, ncols = z.shape

0 T8 C1 S5 I: b8 M g' _

x = np.linspace(dem[xmin], dem[xmax], ncols)

y = np.linspace(dem[ymin], dem[ymax], nrows)

x, y = np.meshgrid(x, y)

region = np.s_[5:50, 5:50]

2 F! w4 Z! ~2 t5 r& L: V' b' H

x, y, z = x[region], y[region], z[region]

fig, ax = plt.subplots(subplot_kw=dict(projection=3d))

ls = LightSource(270, 45)

- d& J* f$ `/ @: u, N4 }, r

rgb = ls.shade(z, cmap=cm.gist_earth, vert_exag=0.1, blend_mode=soft)

surf = ax.plot_surface(x, y, z, rstride=1, cstride=1, facecolors=rgb,

linewidth=0, antialiased=False, shade=False)

plt.savefig(example4.png,dpi=600, bbox_inches=tight)

( ]4 [4 u; L- g. B% `

plt.show()

]& X: m+ a* r

Example 5：三维地形，包含投影（Python）

Example 6：切片，多维数据同时展现（Python）

3 h/ V {$ `$ X* L# [, K * M, h$ M# @ j; N, \

Example 7：SSH GIF 动图展现（Matlab）

Example 8：Glider GIF 动图展现（Python）

- |# @5 F$ g9 N+ q2 V+ z: E# n

Example 9：涡度追踪 GIF 动图展现

" [! z+ r, P1 }4 h4 ^4 s

; {* K2 u. u- [: @0 q3 [/ a2 y " f' T/ b) k" Q: C, f