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

, `8 p7 X/ c; J: ^2 o# ~

. a! h7 @4 P- `

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

* `% i! U. f3 J0 o% _! J9 d

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

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

3 \/ G" T6 ^* T4 Z/ @9 B* u

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

import numpy as np

import matplotlib.pyplot as plt

7 y6 `. ^& I0 B; \" ~$ v

# 创建随机数

7 q: Z% M6 s) H/ ]3 o, A$ `; m, r

n = 100000

1 o" C: a ~' |9 h

x = np.random.randn(n)

2 f8 j7 U6 d) _9 U& o, `, q8 N

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

fig1 = plt.figure()

6 H' o& l# v7 m" A# z. w

plt.plot(x,y,.r)

& F, Y5 R$ ]( Z f2 ]

plt.xlabel(x)

plt.ylabel(y)

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

nbins = 200

. d' U& ~% T, u ]2 X, {+ _

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

# H needs to be rotated and flipped

8 G8 f3 h" j% M

H = np.rot90(H)

H = np.flipud(H)

# 将zeros mask

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

: f9 C# X% i0 _# u# q

# Plot 2D histogram using pcolor

1 [0 q5 [7 O& O: `4 G v/ }& h

fig2 = plt.figure()

plt.pcolormesh(xedges,yedges,Hmasked)

* t% Z, A9 M+ I4 O8 g1 S

plt.xlabel(x)

plt.ylabel(y)

* Y- O+ s2 n; j

cbar = plt.colorbar()

cbar.ax.set_ylabel(Counts)

" E; d) V/ \6 }# H2 s ?

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

5 z3 @0 y: }" ]0 W. R/ J; U2 x7 C, b

plt.show()

; A4 F7 |% n9 C0 ?0 p

' v; i6 Q% }$ |8 u% T) g3 t3 X

G9 e6 Y/ S3 R W/ a9 U

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

! U7 K: x" ~* \

% Q/ O9 j+ E6 c* ?

Example 2 ：双Y轴（Python）

) ~# s6 i6 v1 m" y/ @+ N

import csv

8 q3 r/ y0 S% D$ c

import pandas as pd

1 X2 M5 K9 t' b/ s7 N

import matplotlib.pyplot as plt

$ n. z) K- }$ X3 I/ N9 T

from datetime import datetime

+ Q ^/ p6 Q% B3 N

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

time=data[date [AST]]

3 k0 ]2 t6 c |0 A

sal=data[salinity]

tem=data[temperature [C]]

9 S. `. t7 A: ^( v& H6 @

print(sal)

* u0 N0 V# A* u) r- t# g

DAT = []

for row in time:

) A7 R# v; u( u; E, J/ e

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

#create figure

# e9 i( Q+ G0 x0 K

fig, ax =plt.subplots(1)

9 [/ B1 M) V- S

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

plt.xlabel("Date [AST]")

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

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

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

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

# Z' L( F) c2 \# \

fig.autofmt_xdate(rotation=50)

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

* l8 t. r# |! d! v. b4 k$ e" M8 ]

plt.twinx()

- Y/ n% {8 h+ G1 y" W

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

) z& ?3 L# {$ ^9 G

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

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

+ \, z. R+ R" n% T

#To save your graph

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

plt.show()

2 t* i7 C8 [2 R* \

( @- r0 L) r1 X( p9 r

Example 3：拟合曲线（Python）

9 r: B4 i/ F- p5 f6 Q

import csv

, E2 H/ |# y( S, e) Y& J/ A

import numpy as np

import pandas as pd

from datetime import datetime

0 d( z) ^( Q' U% k5 N2 _

import matplotlib.pyplot as plt

import scipy.signal as signal

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

time=data[date [AST]]

4 E+ V3 N/ J2 l+ d0 S2 u$ t

temp=data[temperature [C]]

) j. ?/ P4 O* a; U

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

DATE,decday = [],[]

! S; j% _4 R; c$ R) S# B

for row in time:

: e- r1 d5 f8 e# W2 z

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

DATE.append(daterow)

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

# First, design the Buterworth filter

N = 2 # Filter order

* }- @+ H" ?- ?

Wn = 0.01 # Cutoff frequency

2 w9 v2 ?2 z) n d! V, O

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

0 a8 n9 c& y# ^& k3 D ~+ z* c

# Second, apply the filter

, [, E% E. B. ^/ L; x

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

# Make plots

fig = plt.figure()

% `) }7 o; N# g4 I

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

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

plt.ylabel("Temperature (oC)")

, `' [+ w- g% H+ q

plt.legend([Original,Filtered])

' ]7 j4 l" m- E, q

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

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

ax1 = fig.add_subplot(212)

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

W. D# G1 {4 f! ?

plt.ylabel("Temperature (oC)")

plt.xlabel("Date")

7 `3 K+ A3 e, Y4 L+ v

plt.legend([Residuals])

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

plt.show()

Example 4：三维地形（Python）

4 A7 J3 k$ {, B, n6 \/ J3 A

# This import registers the 3D projection

from mpl_toolkits.mplot3d import Axes3D

from matplotlib import cbook

from matplotlib import cm

. G3 |1 U# }7 ~- B% s+ g

from matplotlib.colors import LightSource

, r2 L/ F8 T" y8 i8 O' W

import matplotlib.pyplot as plt

% I( p+ W0 j* ~( E! a2 U+ M

import numpy as np

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

% {# B1 J/ Q5 r2 o4 Z: J

with np.load(filename) as dem:

z = dem[elevation]

( T# m0 {1 {9 E8 s

nrows, ncols = z.shape

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

k7 m0 d- `& E

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

' L6 W6 S' m4 @, @

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

g& \3 D8 A: n( l

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

. m7 K5 D5 h, c1 S i5 U

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

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

ls = LightSource(270, 45)

9 z- e0 S) Z* X' h

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)

m$ `5 k$ @7 s0 ~

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

plt.show()

; U' Z# D" a5 H: l+ b$ R) {7 K2 Q V, i! C" [9 Q4 u( J

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

# W% t8 h/ i! F$ g, V * t% z S/ T0 u- [

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

q7 h0 j4 F2 h0 C0 s# v8 k* Q

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

' u, ] a2 I1 { ' O1 P+ H: @( x. @

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

2 C2 A$ Y' b# _6 v2 v% Q & @1 w% r2 k' z4 w, z+ N

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

! P, M% S3 e$ Q/ _+ G- R3 C( z 7 m, `2 m3 m r" q

3 j! p3 M( Y: c* @* c: q, N

3 `. K; B( Q. A/ a3 ^2 B2 L

) }- A$ H* v5 T& w9 Q 8 D* q5 k8 i9 R! ]# i- n : Q/ [2 P9 [) e" I