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

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

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

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

& B, p% z' [9 J

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

- w! [( @! Y/ ~

import numpy as np

8 V+ v6 v3 u% e3 @6 O6 m' Z5 d

import matplotlib.pyplot as plt

: f" a$ g j) w+ N( R. e

# 创建随机数

n = 100000

x = np.random.randn(n)

6 t! r1 K" g& `. U. J8 V

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

fig1 = plt.figure()

plt.plot(x,y,.r)

plt.xlabel(x)

plt.ylabel(y)

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

nbins = 200

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

# H needs to be rotated and flipped

7 ?% I! s) v- A

H = np.rot90(H)

H = np.flipud(H)

6 k, \5 }/ p( a; A& N9 }) L

# 将zeros mask

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

- ]( B6 x9 p1 Z t# f0 `

# Plot 2D histogram using pcolor

" v, h; V9 g' N+ T9 s4 d

fig2 = plt.figure()

plt.pcolormesh(xedges,yedges,Hmasked)

! P: {/ ], Y4 I. s( t

plt.xlabel(x)

: u1 V0 C8 h1 o1 v/ t0 h; J

plt.ylabel(y)

cbar = plt.colorbar()

3 X: k& x- M% \; n2 Q

cbar.ax.set_ylabel(Counts)

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

/ j3 @5 d3 p! K! ~* U/ _

plt.show()

: N- @$ T1 l, Z9 E

# K8 p! ~. }, C/ p) a

: N) |) N$ f6 O ^

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

" g" T% N4 ^9 |+ J) L" X+ T3 O3 P- t2 v

7 {2 \ c7 U! M1 Q

* m" M8 ]6 ^# Y; f5 V; X, U

Example 2 ：双Y轴（Python）

% U+ c$ H# Z6 Z9 V) c/ Q

import csv

. A6 m6 q$ n" g

import pandas as pd

import matplotlib.pyplot as plt

from datetime import datetime

|/ }6 {3 H4 B) I q% ~- \

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

time=data[date [AST]]

sal=data[salinity]

2 R Q0 E M9 M2 m

tem=data[temperature [C]]

- Y3 f* j; D/ l+ _

print(sal)

: g, {; \; D! Q: v: D

DAT = []

. r2 i1 h6 @- T. {. E/ c6 u

for row in time:

1 @( v% _; I3 v- N- r/ h w7 M8 A7 R

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

#create figure

( Y; r% p$ N' r9 [6 |

fig, ax =plt.subplots(1)

: ~+ m# D+ U5 l1 W5 u# \" D

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

6 ?, ~- r, B7 B: i5 A; L

plt.xlabel("Date [AST]")

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

& p8 F4 x2 |' P% r; G+ ~0 s5 v7 w+ R

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

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

9 d! ?) K/ X) z

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

: s6 I% v- P4 g6 f5 W- v5 }% S

fig.autofmt_xdate(rotation=50)

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

plt.twinx()

- M( ^. X: z- s9 ^( ^

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

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

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

3 {& i/ ^8 I) W( B0 O. Y( Q# P

#To save your graph

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

( S, P+ S) C3 A5 V7 I

plt.show()

( J# u- e8 l9 K) y$ c" n: o6 R- q3 v

Example 3：拟合曲线（Python）

K* x/ K1 N) F- u* g5 m

import csv

' c' R5 b; L8 F- q% ]* }2 S, N5 Z

import numpy as np

import pandas as pd

from datetime import datetime

0 E# O: J9 ~3 X- O; ~# H

import matplotlib.pyplot as plt

import scipy.signal as signal

, X. b' Q0 \5 b8 d+ N+ u" p

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

time=data[date [AST]]

% F- c( M; r% u+ ?( [, j

temp=data[temperature [C]]

: k7 |; C- ]* y2 ?% X

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

DATE,decday = [],[]

, y. u: h- u, } l2 N3 q

for row in time:

3 C$ }! V8 @3 H& ^* _& q8 j) \

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

DATE.append(daterow)

! N- v/ i- R5 h

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

& r4 q/ g' d/ J

# First, design the Buterworth filter

" f& m7 e5 `2 O. n5 u

N = 2 # Filter order

Wn = 0.01 # Cutoff frequency

- [, \( M& {" B9 U

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

# X! Z) o, J7 ?7 P9 ^0 T! b8 M

# Second, apply the filter

, C6 \3 t; B+ k0 U" L

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

4 v0 X3 J N! k

# Make plots

0 {; u! C( E! W

fig = plt.figure()

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

) ?+ N+ `8 }* m7 }" o

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

3 P' N: [' L0 {/ j% `8 S0 K! |

plt.ylabel("Temperature (oC)")

) F* Z& g2 `2 @. n. ^

plt.legend([Original,Filtered])

4 n( K" g1 x. |7 ~0 s- `! A

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

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

& ~, r7 V& M* H2 @$ H

ax1 = fig.add_subplot(212)

! E* [7 M, p0 O9 f4 ^ V( T

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

plt.ylabel("Temperature (oC)")

% V( a& k8 `3 i5 ~

plt.xlabel("Date")

plt.legend([Residuals])

0 O, ~3 ]* e2 a% y" I) F! x

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

; H( C% m# V2 o! [. u" A8 m( t1 `' D

plt.show()

+ y; ^2 c. }; I$ p) e) r- m( U' s

0 G. ^ o' Z; a# \) q F/ ]

Example 4：三维地形（Python）

# This import registers the 3D projection

from mpl_toolkits.mplot3d import Axes3D

from matplotlib import cbook

4 U/ L- _. ~. G, m

from matplotlib import cm

* p! C+ t0 {$ q$ q1 `8 `+ V

from matplotlib.colors import LightSource

import matplotlib.pyplot as plt

( N& d2 R {0 w8 h, a7 F& g/ W

import numpy as np

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)

8 {1 O* h' v5 Q& \( k* B& @

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

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

" H0 _1 i' @$ c

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

3 W# Z8 o1 B' J/ V+ \8 E% B

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,

* S* H+ [# @ h F) O0 A

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

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

, S7 p8 B1 j' s% u, Z* m1 c9 ~

plt.show()

) w H! K9 L( ]4 S! g) [! {

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

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

. e' M: f% z Q2 i$ c+ r# P ; v P7 o I3 K

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

1 _' \3 }" k d- s' v 0 G! q5 t9 \( `, E2 ]) k

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

, E* g. P Y) s+ G0 M

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

$ ?+ |6 y0 F6 ^& P9 r* ?2 M ! x* J) V4 p2 S; H0 ^

2 Y& C- W. d9 U! D, `. x

4 M7 f2 g- U& z' K8 Z3 b