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

) @3 _; G% b: @& r; D! W4 n v8 \% D

- t% Y& P) E; U) J& Z

0 D1 ?8 t& q/ n5 {

1 ]8 u1 n! j% X9 _; c- M' ~

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

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

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

& J' O5 C7 g4 d5 M8 n0 `- Z9 C

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

import numpy as np

import matplotlib.pyplot as plt

# 创建随机数

$ l6 K% A7 ]$ o% m/ Q0 @% L, X/ H

n = 100000

/ r6 I* W; I+ z- a5 B. J! d$ |1 a+ `

x = np.random.randn(n)

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

) p! ^7 F) X" ]$ K, V; r( l

fig1 = plt.figure()

; y3 }) ]$ |% ~3 c% w% u9 v( h

plt.plot(x,y,.r)

plt.xlabel(x)

0 z2 ?3 r. O/ P- `7 u L- r

plt.ylabel(y)

& X2 H' I+ [$ P5 E$ m9 i

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

nbins = 200

1 E9 M3 z6 E( ] u! }2 x3 U

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

# H needs to be rotated and flipped

" e) @- u! A4 j# w0 ^2 i; x n

H = np.rot90(H)

+ ^2 F" q4 U Y9 |4 o0 z

H = np.flipud(H)

' B/ z) Z0 D( c6 r) ]

# 将zeros mask

1 I* v) y, K* ]8 k/ ~$ r

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

: r* b! T- Y4 m

# Plot 2D histogram using pcolor

5 N; G" e3 z2 G9 I

fig2 = plt.figure()

% s4 `: [ Z, M

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

5 v' x! l1 ^6 ]5 D6 F

plt.ylabel(y)

cbar = plt.colorbar()

cbar.ax.set_ylabel(Counts)

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

plt.show()

1 G2 d9 h: O! |- m) I

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

+ A7 t/ U! n" [4 Y3 H$ P& F- y

2 K# s* C" i) [0 }; p! ^. M( H

Example 2 ：双Y轴（Python）

X# N6 J. n( I( X5 j2 l6 s; g

import csv

import pandas as pd

import matplotlib.pyplot as plt

l7 L6 i. Q* s8 X3 Q) U1 u

from datetime import datetime

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

time=data[date [AST]]

, a' \) M2 Q/ e1 F" K

sal=data[salinity]

tem=data[temperature [C]]

print(sal)

2 e- N7 j- `* F& i( j0 E: U

DAT = []

for row in time:

0 v7 g# P9 d; T; s

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

; k! m& k3 h5 n$ ~

#create figure

fig, ax =plt.subplots(1)

) X% u) d! C0 s" X$ V7 G8 s( m+ u

# 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)

H# O" Z1 j( f9 c% m9 p; o

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.

/ p4 k: {) z9 w2 s

plt.twinx()

7 z, |2 J: d7 J2 r2 p5 }# J, B( H

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

6 m. a( b/ q+ Z

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

/ q `. r" u# W! } b3 \

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

#To save your graph

2 N! i+ ]- G+ w

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

; {( \, l( _/ ^# \# V

plt.show()

Example 3：拟合曲线（Python）

+ y3 a" n0 @' G

import csv

4 f& F- O# ~/ a

import numpy as np

import pandas as pd

from datetime import datetime

import matplotlib.pyplot as plt

import scipy.signal as signal

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

3 M" R% G, ~5 m3 @* `; V9 q

time=data[date [AST]]

temp=data[temperature [C]]

* `! p& |' @4 p+ J; d

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

DATE,decday = [],[]

for row in time:

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

9 O0 B9 z% Q, y9 i) k% [5 D

DATE.append(daterow)

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

# First, design the Buterworth filter

( ]8 g3 y; P9 m; P: k$ M- f& Q

N = 2 # Filter order

- a! i, y( b% T9 n& `1 ?% S

Wn = 0.01 # Cutoff frequency

7 g6 O5 \3 v7 S7 K% [

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

# Second, apply the filter

7 K$ n }7 X$ X# W4 U0 _0 C% m$ t

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

! b$ D& N& e) n& Y

# Make plots

" \6 W) r+ N+ K8 L0 h

fig = plt.figure()

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

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

plt.ylabel("Temperature (oC)")

& G/ q5 z5 \7 k2 K+ t# @- h

plt.legend([Original,Filtered])

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

; A% L0 J* k- W' a

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

ax1 = fig.add_subplot(212)

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

$ y4 L2 C5 j T1 v! k

plt.ylabel("Temperature (oC)")

plt.xlabel("Date")

, |2 J' F* D4 q' v

plt.legend([Residuals])

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

plt.show()

6 T$ G7 m; |, G G! E# V

8 H' H# T+ i, \; ?0 A7 C6 \+ s* k

Example 4：三维地形（Python）

# This import registers the 3D projection

from mpl_toolkits.mplot3d import Axes3D

from matplotlib import cbook

from matplotlib import cm

from matplotlib.colors import LightSource

7 ?& b! w" I/ A0 G% u

import matplotlib.pyplot as plt

. d- D9 S' P/ l

import numpy as np

4 Z) m, E3 G0 x0 g

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

1 {) x$ k' _! E( I @

with np.load(filename) as dem:

z = dem[elevation]

% G" L( j% s3 z% P! C

nrows, ncols = z.shape

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]

/ W5 P9 L. @" H- ^; e, X0 h8 S

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

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

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,

' s8 g. t5 t6 a, I/ }# x

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

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

plt.show()

/ O! h e- i: R + l# L8 g$ ]2 u1 M4 p

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

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

9 ^3 v& i5 R2 w2 x9 D' h+ n . a# U' V7 R( U8 @. M% C* i

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

% l' o" c8 A$ `( y 6 b8 d2 M- C7 W& b2 ~3 i

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

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

0 W$ Q }4 d: d" G- i% I4 T

, g$ D$ |8 e1 L. A9 w