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

& l8 P5 A9 [4 x. f' ~" `& a

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

6 Q4 V: Z. H" T3 ?& G, T, J" e

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

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

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

import numpy as np

6 c2 o7 I# K! p; f$ p6 B2 ^3 a

import matplotlib.pyplot as plt

# 创建随机数

8 O, M6 f' ~$ ]/ q

n = 100000

* D4 y8 Q/ W3 a8 Q6 t+ Z4 [. I

x = np.random.randn(n)

5 w9 B7 |: f) j" x

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

fig1 = plt.figure()

3 v1 k4 f% _# N- h6 \2 r# Q3 O

plt.plot(x,y,.r)

plt.xlabel(x)

" _/ S7 p& A/ _, F5 ~7 M

plt.ylabel(y)

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

2 o/ L1 J% [9 z. X9 K

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)

# 将zeros mask

2 a+ D4 z L, X5 o9 W5 k/ y

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

2 t; e: E* ~1 q1 a0 G' T+ w+ n7 r

# Plot 2D histogram using pcolor

# f8 ]: c0 d, }$ I

fig2 = plt.figure()

9 V+ I. I5 a2 C; t7 P

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

$ O( N l2 W! M7 c

plt.ylabel(y)

7 x- o% X6 M) w6 G: P. V; l

cbar = plt.colorbar()

cbar.ax.set_ylabel(Counts)

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

& V5 D; x! p$ D, e2 F5 O# Z j) E

plt.show()

: h$ A+ |. C5 } d! C

: k, f |% Z# D

8 ]% o) E( N: b7 n% L: ]& J

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

y, k# T5 R3 E% M

* ~' F& E% W2 t2 g5 ~% w- V

3 R& `* o: `( S8 o$ N

Example 2 ：双Y轴（Python）

- P7 f9 e% H5 a! M) F( o3 _

import csv

; S) Y) s- v0 \" N) U1 Y5 `+ m4 @4 E

import pandas as pd

import matplotlib.pyplot as plt

7 Y0 T" E5 u# j

from datetime import datetime

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

) Q5 t2 B! U* o

time=data[date [AST]]

sal=data[salinity]

) x5 T) P# B2 P# L* v# f, f! q

tem=data[temperature [C]]

print(sal)

0 `9 F# ]% k5 u( x! L

DAT = []

9 V- v/ p4 p% Z4 d3 g

for row in time:

1 W; I* Z q7 p% _* Z7 N

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

#create figure

fig, ax =plt.subplots(1)

% n4 d' M& N: Q% R

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

plt.xlabel("Date [AST]")

% A* W5 I. f6 a( o' K: G

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

0 w6 u e$ v% U: o u4 Y+ U6 V

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

) q0 c! E K H% h3 z* A

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

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

* ]- b& A5 h6 ^+ X, C7 z

fig.autofmt_xdate(rotation=50)

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

8 z5 e" @6 g) {! p' P

plt.twinx()

: F9 w7 @# |) ?9 c7 z2 a5 H9 O. d

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

& {5 N+ E! h, R, n9 N& h4 F# F5 L

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

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

9 e# L' K3 T1 P# Q6 T9 B

#To save your graph

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

9 t3 ?3 x+ \ F" k! f

plt.show()

7 ~5 w1 k# O3 U" Y7 a: ~( O! P9 x

Example 3：拟合曲线（Python）

import csv

Z! a U6 a$ D' j

import numpy as np

7 v3 R; h+ o7 j3 I1 j

import pandas as pd

from datetime import datetime

import matplotlib.pyplot as plt

# G5 r5 {( P. Y2 w" e. f, A4 |

import scipy.signal as signal

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

( J, W( n' {7 n) I6 O- ~; b) Q* `: L

time=data[date [AST]]

( [5 w+ @3 n) W% A7 X1 N

temp=data[temperature [C]]

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

DATE,decday = [],[]

4 p" m: i0 b6 b) _# H

for row in time:

$ l. @! W5 E8 i. D7 P7 D

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

1 @0 J7 r( }# ?% g

DATE.append(daterow)

/ j* \1 Q" {5 E' z* S5 |4 i

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

5 k' m3 P% H) @" v' T

# First, design the Buterworth filter

% x' U3 J; C- T, v

N = 2 # Filter order

* W0 k1 b2 a8 t8 D/ E

Wn = 0.01 # Cutoff frequency

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

# Second, apply the filter

; _1 `8 G* N; W+ {

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

5 ^: \) S" E7 _0 |$ N, k% J

# Make plots

fig = plt.figure()

ax1 = fig.add_subplot(211)

) a- O- S4 x; Z, E+ y

plt.plot(decday,temp, b-)

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

plt.ylabel("Temperature (oC)")

( |3 J" P, O+ H4 V, O$ p

plt.legend([Original,Filtered])

. B7 [( E7 D0 W

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

plt.ylabel("Temperature (oC)")

4 a' d8 E' S# I( Q) N% g3 u

plt.xlabel("Date")

5 P7 ~/ @, |; a; q |4 Y

plt.legend([Residuals])

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

plt.show()

7 Q" M1 m$ o7 ~

Example 4：三维地形（Python）

' P: z1 ~7 o7 N' r& G

# This import registers the 3D projection

) E) k, O$ Z' a+ _6 ^

from mpl_toolkits.mplot3d import Axes3D

# {# b4 D* h# f6 A& [

from matplotlib import cbook

" Y& g: [2 P3 e- F

from matplotlib import cm

5 ]; ~ n' |" K' n5 k

from matplotlib.colors import LightSource

5 j1 ?: @9 L8 l8 ?3 F! R ?

import matplotlib.pyplot as plt

9 K) q; i/ S, K5 k7 A& z; x7 y8 x

import numpy as np

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

with np.load(filename) as dem:

z = dem[elevation]

9 W8 J! c2 a7 x! c: X8 M

nrows, ncols = z.shape

N: |. I3 F2 \ ?

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

% C* k* s. Y* V n1 o) I5 F6 P

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

3 e9 d4 G) W( w2 x( P* g, V1 ?

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

& Y$ |* ]2 r& I/ {0 `! m7 V

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

- K; j1 j6 X, j' \& s- ?) C1 C# c

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,

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

6 z: M" V0 J. _& }5 g/ \

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

2 Y0 G' w0 A2 H; v( E" L2 z

plt.show()

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

. F% S$ D0 l P8 c) d _

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

% v4 F6 u/ n0 k7 K% d % M2 [! J9 w! g u3 S, u

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

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

2 [1 l' K; m3 ]. x/ e3 q& b' r 9 ]: a% V/ Z* X A2 R$ ^

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

2 F5 q7 Y4 O+ n) b: \

6 y2 }2 _9 e! F" i% s% n2 ^

1 A* e6 L8 c/ I* \

3 _9 t4 E( z. J6 {& t t' q0 I& @ : E' {3 B2 O; Y9 } * W: v* H* q% C