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

3 u; O( p; `) B0 h3 P

% X1 [+ F* y5 d N! }

$ Q- l% P' e- d; o8 C" N8 e1 Y

, d F" T+ ^( O4 S7 ?

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

/ p8 O/ Q' _0 Y9 R# Y d8 Z

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

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

( D; I% D7 ~0 O5 ]8 q% d( K0 [! i

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

* o! K n/ p* ?$ [' ~

import numpy as np

6 f# T3 ]$ {) S+ F6 [0 S2 A. `- j

import matplotlib.pyplot as plt

' C2 y$ E( Y, I* G

# 创建随机数

3 j7 w; S2 w% D4 b& X, n

n = 100000

0 X v; u3 D% l1 m" A

x = np.random.randn(n)

& n( c$ j s/ }4 f7 P2 H

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

fig1 = plt.figure()

plt.plot(x,y,.r)

. y: c8 p) {5 M4 h p

plt.xlabel(x)

plt.ylabel(y)

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

nbins = 200

: m- m% D& X' I: O' E

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

# H needs to be rotated and flipped

H = np.rot90(H)

% n/ c. l6 u3 u) `

H = np.flipud(H)

* a3 Q$ y" H# d+ {& i; d6 b: f

# 将zeros mask

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

# Plot 2D histogram using pcolor

. k% @6 Z" E' y: S; @

fig2 = plt.figure()

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

plt.ylabel(y)

cbar = plt.colorbar()

cbar.ax.set_ylabel(Counts)

1 v, ~: N4 t/ k8 c, d( V3 z

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

plt.show()

@# H8 { |7 c

" Q- H) u2 M, ]. ^1 w

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

5 ] R/ `6 V- q g

# i* X1 b7 {* L

6 y6 v0 u/ M$ u+ E7 g7 K9 z

, _( D5 }" O% A7 p7 J& P0 ~/ n/ J

Example 2 ：双Y轴（Python）

import csv

import pandas as pd

8 U8 Y4 \) h, H4 h

import matplotlib.pyplot as plt

from datetime import datetime

) T; G! M: d9 V1 K4 d

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

time=data[date [AST]]

sal=data[salinity]

+ \) Q$ O( D4 p

tem=data[temperature [C]]

print(sal)

DAT = []

for row in time:

. D2 x) S# E- s& M1 `5 B5 \6 k8 _% O M

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

" b/ R3 J) C. I: P! n" v6 K' s

#create figure

( A" @! E8 M' j- u' G

fig, ax =plt.subplots(1)

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

& m( o! j; U( L+ I$ h1 I3 L. k

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

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

fig.autofmt_xdate(rotation=50)

* V) }( i X5 N8 Y

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

plt.twinx()

% o2 s* J" J3 P. `

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

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

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

. Q6 L7 i( J- p9 c

#To save your graph

4 Q) o7 R3 A- l9 {# l

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

plt.show()

+ y" K0 @/ ]# R$ o6 A d" C

Example 3：拟合曲线（Python）

import csv

: {% i, V2 }1 W% }3 b9 E8 G* \

import numpy as np

3 ?( g* s' c* X" L! D

import pandas as pd

from datetime import datetime

E3 s$ i8 _! C( `' H5 ^1 `

import matplotlib.pyplot as plt

import scipy.signal as signal

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

! J1 o+ h3 {- j$ ]/ a- d( Q

time=data[date [AST]]

temp=data[temperature [C]]

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

DATE,decday = [],[]

+ K$ M) U+ t0 Q# S7 U3 X

for row in time:

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

% F/ h# {; E j* `# p/ y

DATE.append(daterow)

+ `. Y$ v4 X2 H

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

' m! d$ T* w' Y* h

# First, design the Buterworth filter

# D4 ^, F# j2 Z

N = 2 # Filter order

Wn = 0.01 # Cutoff frequency

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

# Second, apply the filter

$ c( f5 j7 T# h/ _& W/ m$ O$ M

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

+ d& K: N% U9 I* m

# Make plots

fig = plt.figure()

ax1 = fig.add_subplot(211)

( A' I0 m" q4 @) T3 {- M. q$ z

plt.plot(decday,temp, b-)

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

plt.ylabel("Temperature (oC)")

plt.legend([Original,Filtered])

- t& J7 j: J9 h

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

* S1 d3 o# D! `6 u+ H

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

ax1 = fig.add_subplot(212)

3 E0 d( Y; ]& D- K# I

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

% y s F9 a* C2 ~

plt.ylabel("Temperature (oC)")

plt.xlabel("Date")

$ a% Y; y6 g: e9 H

plt.legend([Residuals])

& ? y2 D. J( ?/ R' `

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

: Z9 I I/ L' t) T. D! T

plt.show()

7 Q9 h1 L e1 Y" |: A7 Y# U( _7 v

Example 4：三维地形（Python）

5 Z# w+ b3 I/ g( a) u( V+ ?" V

# This import registers the 3D projection

# u" b$ y8 }- j0 J

from mpl_toolkits.mplot3d import Axes3D

! S6 \) q6 l; |2 x7 a1 u

from matplotlib import cbook

from matplotlib import cm

from matplotlib.colors import LightSource

import matplotlib.pyplot as plt

import numpy as np

* w8 H2 o! Q( T7 S) G5 {6 q

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

with np.load(filename) as dem:

z = dem[elevation]

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]

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

$ }+ y2 G! J( s, f

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

ls = LightSource(270, 45)

9 u4 Q+ @. c l0 d% _- |2 E

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,

0 O" L0 S0 Z1 E* k M& W

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

; A& N- t' j, b* k# i4 P k8 Z7 g

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

plt.show()

+ q/ Y* E& N* |+ t d, K1 A ! b& i9 v' p5 m& y r/ h6 M1 H

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

2 F3 f3 g7 E* d' J: E$ p

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

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

. p1 f3 e o& h1 s ~, _" T" U5 q

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

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

, i; p$ ]* ~6 I0 _( \ B

* ]: A: p3 s7 w* K) }# D- f9 ^+ }