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

6 s7 ]- C7 c. l6 G2 `' r

6 p. f# L! G8 R8 B& M* b3 E, n

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

' k) U( h5 s, f: M; M5 r; x+ f

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

# k Y7 Z8 h( d) w* k

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

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

: L8 y; t' q0 d6 P& s

import numpy as np

7 {$ e2 j6 e; \

import matplotlib.pyplot as plt

# 创建随机数

n = 100000

x = np.random.randn(n)

* T8 k" R4 b+ `8 o) i

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

) e5 F! W C+ D' A7 H5 x j

fig1 = plt.figure()

4 N3 r* m$ D" b( \% r f" R0 {

plt.plot(x,y,.r)

plt.xlabel(x)

2 q) d3 z/ b/ _% l4 T! }

plt.ylabel(y)

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

+ [$ V, x+ J( \2 L" h

nbins = 200

2 w0 t7 T4 h5 J5 [4 B* T; z

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

# H needs to be rotated and flipped

' Y4 k2 H. _ V7 C4 V# A" f$ P

H = np.rot90(H)

: t, j- w7 a5 _# H

H = np.flipud(H)

# 将zeros mask

# [- X2 B2 S; a& G* W# J. |) |1 u- a

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

' s% c) A0 G, `

# Plot 2D histogram using pcolor

0 T. f$ X; \$ g* z/ N3 r

fig2 = plt.figure()

! H3 N5 ~4 u0 J( t6 L; w

plt.pcolormesh(xedges,yedges,Hmasked)

w7 s y7 k; n/ k5 Q2 \- Q" j

plt.xlabel(x)

' P0 ?* ]6 y* t- m! C/ S) ^" H

plt.ylabel(y)

cbar = plt.colorbar()

( B& @+ q) z4 A% s

cbar.ax.set_ylabel(Counts)

' ]: G9 }/ Z5 ^- D( b* ?

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

" U' l5 k3 f) w4 f

plt.show()

% l" G: ?4 c0 M4 ?+ U, V& T6 w ?2 n& M

) s, p( t2 \) X1 P* Q0 l

L4 {( {8 }' h$ P( j

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

7 n; f9 `% C: m+ s) W

- q6 B7 F; e0 X* `8 g

+ F; r0 K1 _; Z# Y

Example 2 ：双Y轴（Python）

" I3 B p L z6 C5 ^/ P) T

import csv

import pandas as pd

% Z; E7 `% u' ]5 | N/ N" A; W

import matplotlib.pyplot as plt

from datetime import datetime

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

time=data[date [AST]]

sal=data[salinity]

. q. n6 [% U+ E% ~9 }' K! s9 r

tem=data[temperature [C]]

3 j. @% z. ]1 D! y' b0 V( F

print(sal)

DAT = []

1 d! o9 I) J: z

for row in time:

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

#create figure

: S. p$ y, B( _; X# r

fig, ax =plt.subplots(1)

) A7 \" h) r6 v* r

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

E' l* X+ }( u1 {2 D

plt.xlabel("Date [AST]")

4 L9 Q( a& o5 T1 w5 I, m

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

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

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

6 I$ a8 R4 y+ {" \% w! r% ]

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

( W$ }& v8 J2 m3 H

fig.autofmt_xdate(rotation=50)

( L; Y2 V3 ]! J! H

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

: p2 u8 P6 J* T3 Q! G0 K7 E! r

plt.twinx()

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

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

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

/ z# Y% P" A3 X

#To save your graph

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

3 B/ f( o- N G/ F" m2 ~% h

plt.show()

Example 3：拟合曲线（Python）

import csv

6 P- E6 q& e+ d

import numpy as np

import pandas as pd

0 z% Z6 e$ |1 \+ X

from datetime import datetime

import matplotlib.pyplot as plt

import scipy.signal as signal

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

7 r. i* P& _9 p7 @

time=data[date [AST]]

temp=data[temperature [C]]

8 P2 D. P4 r6 L' N

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

! {! N2 n7 q `# R& }) I

DATE,decday = [],[]

for row in time:

$ d8 C4 l& i1 `/ P

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

DATE.append(daterow)

' n0 q% }! X1 c4 n7 }7 ^

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

! M- w0 @- s8 ]# M% D; u# l3 ]

# First, design the Buterworth filter

N = 2 # Filter order

A: X# q6 N+ n

Wn = 0.01 # Cutoff frequency

, [& Q a' S$ I- p" S

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

# Second, apply the filter

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

8 u: b9 q- {! F3 C' N, _8 E4 T

# Make plots

fig = plt.figure()

. B' y( m* `. l7 m

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

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

! o K/ d5 A3 q( R4 s5 p

plt.ylabel("Temperature (oC)")

plt.legend([Original,Filtered])

5 n, w5 n9 i' ]/ n7 @$ z4 p

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

. B0 i# F; N* \* [) D

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

* `. V7 Y) t" h' R9 p

ax1 = fig.add_subplot(212)

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

6 p/ z6 N3 K$ }

plt.ylabel("Temperature (oC)")

8 n2 k/ {1 N" |- b& ~- X2 H

plt.xlabel("Date")

; n( Y) _) J% [7 l! J2 m" Q) `/ w

plt.legend([Residuals])

# x) ~$ C( |* H. @

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

plt.show()

Example 4：三维地形（Python）

# This import registers the 3D projection

$ k" ?) x2 A7 B4 `7 {4 g3 T

from mpl_toolkits.mplot3d import Axes3D

: v8 t5 J$ T4 k

from matplotlib import cbook

7 c$ E3 V0 W' p$ F$ S J

from matplotlib import cm

3 [+ @6 y" U0 g/ s2 U( u+ Y1 w

from matplotlib.colors import LightSource

import matplotlib.pyplot as plt

% L$ a3 f" U# L: w

import numpy as np

; X/ S6 h% _: E9 f! L

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

with np.load(filename) as dem:

; }- T2 P/ m1 n6 x

z = dem[elevation]

3 |& ^5 r9 \3 i" L& p

nrows, ncols = z.shape

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

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

$ ]( v& k! M* ]

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

2 F7 Y9 i0 ]# }$ O: i' |$ M8 U, q. s( W

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

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

" X; t5 u/ @% v

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

4 n: w" f* e4 V

ls = LightSource(270, 45)

5 b5 G) T \/ u$ 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,

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

, S/ j# C7 {# B! o

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

plt.show()

0 K% @4 J2 q+ K& { 0 R2 U& I% n* I: N9 G

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

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

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

5 S$ f( q1 O5 g( U9 a; [: D * ?% i) `. z6 |" R

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

2 }/ B* V' R7 e8 Z% L2 l

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

9 Y, x) \7 R0 n% z

% r- `! t% J0 E

8 E9 `; P+ S- Q' [0 f+ Y$ @, o X4 {; V2 ] Q7 O9 z