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

7 D/ q% F5 I; ^7 [# o9 i5 y1 n

5 [7 Y. M0 G9 R7 X+ f6 C" y) Y: X

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

5 @# D* t2 \3 J) V6 B* [4 \

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

2 X3 w. W6 B B1 }+ ]6 ~. h: L

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

5 h S" i: j6 ^

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

U& H. ]" g$ z! [" z" ]

import numpy as np

1 T# |- h) `8 \: \7 H5 s

import matplotlib.pyplot as plt

$ l# j% q# Z) @5 V |

# 创建随机数

n = 100000

) K0 l( S" F) m$ D ]) y

x = np.random.randn(n)

& r$ y; I$ m+ H) d" ]5 u7 {

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

" P0 W: F2 f* ^9 u0 N

fig1 = plt.figure()

plt.plot(x,y,.r)

plt.xlabel(x)

( ]+ w7 L& M. U1 ?3 Q% f T

plt.ylabel(y)

* L- Z3 I p! x6 A7 J& ?1 ?& u: ?

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

! s; F: H- p8 A- @

nbins = 200

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

; x8 \) I* _3 [& Y+ [1 u

# H needs to be rotated and flipped

" |9 R% T3 O/ S8 R4 b) L8 I

H = np.rot90(H)

H = np.flipud(H)

$ w/ r, L# U& b( {8 ^

# 将zeros mask

/ \- Q) x) F% v

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

# Plot 2D histogram using pcolor

) M2 ~* {2 H) e/ z `

fig2 = plt.figure()

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

plt.ylabel(y)

cbar = plt.colorbar()

+ r- C; k" \. N, C3 ~& c' F0 w

cbar.ax.set_ylabel(Counts)

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

9 |" M- W- p; z! Z

plt.show()

6 F# M. g& O6 {& W' `( R

5 Y) J# b7 i! F9 u5 K

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

Example 2 ：双Y轴（Python）

9 g* [: G" g( }$ I3 e: B

import csv

import pandas as pd

import matplotlib.pyplot as plt

( f7 N3 a/ o8 K0 ?

from datetime import datetime

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

! d. a: d# ~* n

time=data[date [AST]]

sal=data[salinity]

tem=data[temperature [C]]

print(sal)

DAT = []

for row in time:

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

. f+ h- P5 ~2 O- Q( J% p

#create figure

& q# F" h' P1 e$ v

fig, ax =plt.subplots(1)

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

' }0 @% }; q# G9 @( D m. M

plt.xlabel("Date [AST]")

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

3 w9 v B! g2 y: ?

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

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

8 H. \+ Y. {; g6 d6 d. }5 R

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

fig.autofmt_xdate(rotation=50)

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

! ~1 {& {& V1 L# P) Y

plt.twinx()

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

9 W$ E# E" z# i+ B

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

, W) I0 q$ n. g/ T Z9 M4 x6 V

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

#To save your graph

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

1 ^3 W* D! d$ h" p" o+ d

plt.show()

9 d _2 [0 g6 t6 p1 b

* u# f4 W( l6 J- h" x

Example 3：拟合曲线（Python）

[% o# u' v, ~, }: ]0 L

import csv

! ~" |- Y9 C2 N& K1 g: Y/ b

import numpy as np

: s8 \: ~/ n1 A# g2 l/ v

import pandas as pd

C; ~. q R0 [! ]- |0 Y

from datetime import datetime

) X: C$ t% a' \8 f

import matplotlib.pyplot as plt

import scipy.signal as signal

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

; c/ Y) N+ N' v" {- A* z+ c

time=data[date [AST]]

temp=data[temperature [C]]

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

) M0 z4 E/ m O% ]+ \

DATE,decday = [],[]

for row in time:

5 k& ~3 ]" R, P+ g' g# a

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

7 A- {) |' a! [

DATE.append(daterow)

9 F1 ~8 V0 w; ?7 R+ K: V

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

% G8 X! W5 ^! ]. c+ O) A

# First, design the Buterworth filter

N = 2 # Filter order

Wn = 0.01 # Cutoff frequency

5 K# |; H: z7 o1 j5 f

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

e3 ?# @, h, [2 C Z" D

# Second, apply the filter

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

* P/ `( I6 e% S

# Make plots

G2 R( O9 W' q2 Y

fig = plt.figure()

. T" X$ _* F0 q

ax1 = fig.add_subplot(211)

9 \9 o8 h/ H4 _- m, A0 }0 P

plt.plot(decday,temp, b-)

8 X2 ?, e" X- N* j

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

plt.ylabel("Temperature (oC)")

plt.legend([Original,Filtered])

* {' R, }# g9 f+ O8 g" d! L

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

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

; w6 d1 L9 }1 |! a) t( t- U, @

ax1 = fig.add_subplot(212)

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

plt.ylabel("Temperature (oC)")

5 S7 f+ o2 W2 Z

plt.xlabel("Date")

2 ~$ A! b6 Q o q' v

plt.legend([Residuals])

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

plt.show()

; t) h& m+ X8 [- k

0 f- ], R4 _! E+ W/ b

Example 4：三维地形（Python）

# This import registers the 3D projection

. b1 ?6 V+ L8 a% t, U

from mpl_toolkits.mplot3d import Axes3D

& W5 B/ g6 k( k

from matplotlib import cbook

6 f$ C, {; l0 G' H/ H& i2 U

from matplotlib import cm

3 E+ C1 U1 A$ D+ q' }' x j$ H

from matplotlib.colors import LightSource

import matplotlib.pyplot as plt

/ F" Z( J& P2 e- F3 v# Y4 O) ?/ O

import numpy as np

+ q8 V# e* n- g. g$ U: m

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

+ W5 D0 _' b5 K6 j% q

with np.load(filename) as dem:

. ^# L/ ?" y8 @' U1 L& d, A; E

z = dem[elevation]

0 p: ~ j) A) R1 |& I+ k

nrows, ncols = z.shape

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

2 ]* A1 H1 F# s5 ?

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

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

6 K: n5 G; s( m. k9 ?' I* C! n

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

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

) E, i, v( f, N1 i H2 w, H

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

, A! N: F" S' y

ls = LightSource(270, 45)

; _5 e5 `- K2 W- Q1 e& c

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,

8 R: v9 K& O) |- _+ ?5 C1 X( O

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

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

4 m1 i8 C7 C. n& D

plt.show()

4 m! F3 T) c2 g* F9 c

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

, }' q- n6 O3 m

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

& D6 Y) d) z3 z7 X* s9 f 9 t9 A& I. J' _3 T

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

5 T0 b. O2 u, F$ l 6 ^) y0 K* A' V" B0 t

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

$ r0 o0 p0 m% m: X/ E ; L- g. h6 u4 I" r# B6 s8 q

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

- N7 e. w3 b' l 2 d( y6 }$ A% q1 q

9 D/ |6 o# m3 m/ x3 r, G

. s! _6 t, m ~! [/ v1 ]