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

: T8 K+ p% W% t2 U' ]) t/ a$ a

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

% m8 y: R+ H/ H. I' z

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

2 c4 X+ @& l; l% U! T: J$ ?- S

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

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

import numpy as np

3 W) W! E" A- \+ P8 v5 |

import matplotlib.pyplot as plt

3 H# c8 j/ B( ?4 t! r8 O4 e& q

# 创建随机数

n = 100000

x = np.random.randn(n)

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

fig1 = plt.figure()

# U; G p1 L ~0 Y: X

plt.plot(x,y,.r)

' e! P4 n& V* W) W6 I3 c6 Q7 T! q+ }

plt.xlabel(x)

plt.ylabel(y)

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

: |( d8 F8 g, N+ `9 ]$ j3 _, Y

nbins = 200

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

# H needs to be rotated and flipped

H = np.rot90(H)

H = np.flipud(H)

, y! p' y% Q, n8 Y

# 将zeros mask

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

C6 t. w: s5 I4 n2 Q% j. K

# Plot 2D histogram using pcolor

. }- S* W6 _' @! y# `

fig2 = plt.figure()

plt.pcolormesh(xedges,yedges,Hmasked)

. S0 s5 q: t% G6 K

plt.xlabel(x)

plt.ylabel(y)

& x: J! I+ C9 v% d; ?

cbar = plt.colorbar()

1 W9 {6 A' X; u

cbar.ax.set_ylabel(Counts)

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

1 e; h4 S! H2 B1 l4 D3 G

plt.show()

( s& I$ x7 b1 p0 C# s, r1 J3 j

+ X; z6 z6 L: o2 H* o& x

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

& Q |# U. N6 c2 j

. B* G3 l' `# N. p$ G- Y$ E d9 }

. i) C* E, j- s0 a# A

Example 2 ：双Y轴（Python）

) W4 f* n! y0 F9 J

import csv

import pandas as pd

) Y. ~+ [7 v! m& }

import matplotlib.pyplot as plt

from datetime import datetime

: D) p& s& U- F% a3 n" L

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

time=data[date [AST]]

sal=data[salinity]

# ]" g9 B# k0 q

tem=data[temperature [C]]

print(sal)

# K9 e( w- B1 v3 ?

DAT = []

7 }, J: m% r/ N

for row in time:

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

3 D# J9 w* A3 h2 ?& \( u; U

#create figure

$ d# h, V' D! ~/ n7 B

fig, ax =plt.subplots(1)

8 e+ ?- @# P% l2 i

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

plt.xlabel("Date [AST]")

8 i6 U+ j H! h% x+ L4 L5 i1 h8 j3 \

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

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

- s: g4 B. V8 b1 _

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

# k# d. Q% _3 \& y6 n H

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

; d% M/ @1 P. N

fig.autofmt_xdate(rotation=50)

& x, e4 f; w; e

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

plt.twinx()

/ W# b$ {2 N0 s. C, w

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

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

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

3 v6 [1 v- U2 _5 i& R7 v

#To save your graph

( Q+ S! O5 r* x1 o. X2 N

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

g! C1 o0 s2 r7 B; _5 d

plt.show()

Example 3：拟合曲线（Python）

& v+ P; A/ _5 u9 c6 n _* n

import csv

import numpy as np

' K& ~+ q$ f! @, C( L& y) l

import pandas as pd

1 Z/ o/ h$ r& d8 W3 \; x7 H

from datetime import datetime

) I* H* P/ n7 R! |; b0 {

import matplotlib.pyplot as plt

9 w3 I- F6 k, w7 W, N( O" R

import scipy.signal as signal

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

8 s W5 o2 E9 R3 R6 R( k

time=data[date [AST]]

temp=data[temperature [C]]

! h/ B4 u" h4 k. B

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

. C* V7 V5 T2 [9 v# k9 y( N

DATE,decday = [],[]

for row in time:

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

. }+ `8 R [9 l, n0 u6 i3 P

DATE.append(daterow)

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

# First, design the Buterworth filter

N = 2 # Filter order

Wn = 0.01 # Cutoff frequency

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

! t: J! h: L3 Y7 z

# Second, apply the filter

$ U( G( v! Y4 r$ F; I1 j

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

# Make plots

7 }4 D/ u8 N) z+ M

fig = plt.figure()

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

i3 _& H4 U2 S( V8 t. o

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

) n: i4 \* _% C8 k; S% ]4 |" }

plt.ylabel("Temperature (oC)")

4 M9 O# t. R" ]$ N( V% B

plt.legend([Original,Filtered])

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

% |, x7 E% X+ `: I4 h; h

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

% m6 {' L! B+ r8 n

ax1 = fig.add_subplot(212)

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

plt.ylabel("Temperature (oC)")

1 v0 R/ _1 i. E; _

plt.xlabel("Date")

plt.legend([Residuals])

7 e+ A# b- L1 S2 d2 M- ^( E

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

) `+ Y% Y2 Z; Q! }4 W5 G: q

plt.show()

" I# U& [( _4 E3 B0 g9 f/ l: p( z! `

Example 4：三维地形（Python）

# This import registers the 3D projection

7 {$ A8 C% N! j( V6 y( }

from mpl_toolkits.mplot3d import Axes3D

$ c5 J" @: \* l, U9 o+ B

from matplotlib import cbook

8 a8 d0 L6 I2 k' ~! y, u* T: o

from matplotlib import cm

: L3 d, ^4 ?4 h3 u

from matplotlib.colors import LightSource

/ u5 f" D8 r9 L6 r" A O0 Z

import matplotlib.pyplot as plt

import numpy as np

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

4 N! t1 m; {0 j. u( d

with np.load(filename) as dem:

z = dem[elevation]

nrows, ncols = z.shape

/ v5 g1 Y) G! k

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

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

' J8 ?7 j% t: j$ x4 q

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

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

2 F X5 P) ]5 Q5 f/ B

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

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

- J4 {9 K( I U" S0 R

ls = LightSource(270, 45)

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)

+ p4 F6 S* U% |8 m; h

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

plt.show()

# a2 W" Q! M% x8 S5 U( v

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

3 P6 p$ {) G- @5 o$ P) G! ~ % L2 U6 S, s% i2 a/ Y" w2 p

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

" r! {) `1 H0 B' C( t7 m. f& Z$ K/ ^ % y$ x L$ r% H7 q/ b- p

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

8 h3 p1 }. N0 ]* F7 C9 {8 M* L

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

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

. O! k: [5 }# `) c- k. `- c

+ m" R* z+ `. K7 T2 o5 `9 B( O& {

5 P- f3 d) @$ _# n5 r) N ) i5 y, E D0 m # F# y3 _% L2 Y6 g. B 3 T# z# v& P3 u. |% m