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

/ l) F/ r% ~9 B K( E$ `

+ d! D) j p7 |/ H; v/ j

5 S. z# ]( g) ~% {2 K

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

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

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

6 {6 m6 Y Q8 y" G$ M a& f+ z

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

: e! ]+ x6 y9 p1 Q# w, D

import numpy as np

import matplotlib.pyplot as plt

; w( {" J' q$ P2 l5 X

# 创建随机数

n = 100000

2 |& V5 }6 T" `

x = np.random.randn(n)

$ O6 ~( @* `! r/ D7 z! T- B% }

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

# A9 |/ K' F& U

fig1 = plt.figure()

plt.plot(x,y,.r)

plt.xlabel(x)

plt.ylabel(y)

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

% ~+ b8 R5 v1 W: C I! B; [$ V8 O+ E

nbins = 200

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

* \* g3 L. h0 C9 a/ Z$ I W( k

# H needs to be rotated and flipped

1 q2 ^5 V4 O8 w8 }* X

H = np.rot90(H)

* i& A3 u% n; x; s( _

H = np.flipud(H)

# 将zeros mask

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

: M- P1 C" y7 \# i" i ^

# Plot 2D histogram using pcolor

G1 z; M4 d! f, Q8 T* v

fig2 = plt.figure()

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

: m0 s3 _; n: n- \) t

plt.ylabel(y)

cbar = plt.colorbar()

cbar.ax.set_ylabel(Counts)

& r/ D5 j$ p! u {2 k

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

plt.show()

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

; s3 U" i' g# l# g S

% p/ [& s; W8 B2 q% o; g8 z

+ V) l2 f u* f: b& r: b5 x

Example 2 ：双Y轴（Python）

2 X8 m0 G% U. `4 l; C

import csv

1 _5 B- f! V( X L0 O! V

import pandas as pd

' N x/ R7 Z7 Y

import matplotlib.pyplot as plt

from datetime import datetime

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

3 K) @" k% A8 C9 T3 T3 T; Z0 |

time=data[date [AST]]

" I! p2 z: x) u

sal=data[salinity]

( p( T% H) l+ o% ?

tem=data[temperature [C]]

print(sal)

DAT = []

for row in time:

) K) c( Q! x4 C0 l* M# r+ |4 X9 t

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

; V) l/ j) O9 }) z; \* K

#create figure

fig, ax =plt.subplots(1)

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

F4 U$ t& w8 l0 ^7 [5 o

plt.xlabel("Date [AST]")

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

% ]7 v. L9 D; G+ t. T6 ] B) y

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

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

* k. D$ q. @0 }5 D- _

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

) s" G0 q. X# v, z0 A0 I4 d2 X

fig.autofmt_xdate(rotation=50)

, }2 P6 q1 ^( t+ ]: y, `

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

plt.twinx()

- E. s! }! U) _3 l9 I

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

# |& b" z" p3 _( \" c1 }: p

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

5 t' l2 F$ L' O& ~

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

- _ ~2 _- x+ e9 p* }5 u; O

#To save your graph

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

plt.show()

Example 3：拟合曲线（Python）

: e: e2 d8 t5 \+ ?" ?

import csv

import numpy as np

import pandas as pd

1 |1 O) v* ?- {

from datetime import datetime

import matplotlib.pyplot as plt

import scipy.signal as signal

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

time=data[date [AST]]

8 ]: g9 R' A4 w3 V

temp=data[temperature [C]]

5 U4 U- y4 I) R) E4 o7 _9 W# Y

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

( o: |9 p; D' V' W+ Q

DATE.append(daterow)

. }+ ]9 q5 l9 Y

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

# First, design the Buterworth filter

6 m/ Q6 z7 N8 n# d

N = 2 # Filter order

7 |* A7 ^9 M1 h& F6 e

Wn = 0.01 # Cutoff frequency

0 C- v7 G/ z: q7 S: X

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

# Second, apply the filter

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

2 S* o4 g2 \0 A$ c6 r

# Make plots

! T/ A% W8 b( r8 A1 N5 f

fig = plt.figure()

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

$ j# L; H2 }" o/ e: k. u

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

/ e+ Z) C6 m4 m* J" M

plt.ylabel("Temperature (oC)")

0 o" q& b0 w6 Q& Z

plt.legend([Original,Filtered])

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

. c9 K5 N9 P8 `" Q

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

ax1 = fig.add_subplot(212)

/ x6 G& W2 ^' W$ G* g% C/ N9 x, j w6 G

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

5 ]( N5 V H4 d( K& [3 |

plt.ylabel("Temperature (oC)")

5 F: b: {" ?, B% M9 U9 b

plt.xlabel("Date")

plt.legend([Residuals])

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

+ v3 r5 W+ r5 K- T) [* T0 _

plt.show()

, B3 I' B* p- u* ?7 t2 o2 @

- A4 ^$ T/ @! ^+ \" G+ e9 r

Example 4：三维地形（Python）

# This import registers the 3D projection

3 a# u# d: x9 Z2 U0 ~: T: l

from mpl_toolkits.mplot3d import Axes3D

from matplotlib import cbook

from matplotlib import cm

$ f a4 @' A! Q6 a( ?- @

from matplotlib.colors import LightSource

! B+ n3 L' B) @6 u ~! M! b

import matplotlib.pyplot as plt

import numpy as np

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

- k5 W2 R2 I9 ~$ U6 f2 H/ [3 _

with np.load(filename) as dem:

z = dem[elevation]

) m+ m% L4 e5 e- ]% Z

nrows, ncols = z.shape

! z: i! p+ e* A4 n+ d

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

; U7 \4 H$ Z T _. l8 J

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

8 _. i) i6 a" g. \8 p' j

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

' L: R$ ~9 J) E- I; W

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

' R& \2 N1 i% w6 o* U2 S

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

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

9 C6 O- j% Y3 D2 h# Q3 `

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,

- C3 g% l2 s2 _& |* B

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

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

plt.show()

& x9 l. `" J: k& c5 A) B / u' K1 A* r: ~

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

7 w. r2 ^2 O- Y9 I* Y3 e, S

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

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

* @, ?# W3 l( l4 Q

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

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

6 w7 E7 N( K" D* W8 O5 U3 v) k8 x2 p % I4 j. _6 u3 O9 ]0 B9 r# P

' e' x- _7 B& t+ n8 V

" O* X) j5 \6 e3 ^' J