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

" L* W- ?8 n& s; b* F

- B( L- _: U: K

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

2 s# |$ I6 W+ q( p

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

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

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

6 R6 R. q5 x! ^% J$ q# h

import numpy as np

import matplotlib.pyplot as plt

# 创建随机数

' f4 `' y, p# n, d/ N; }2 t" ^

n = 100000

; |0 |4 P- V5 W6 @. w' `

x = np.random.randn(n)

3 s- J# d. ^5 Y( z, Y% D& [

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

! Z! D! V8 j, d* a

fig1 = plt.figure()

2 I) y# k# a7 m- t$ {/ \9 T4 ^

plt.plot(x,y,.r)

. ]1 O5 ]$ T2 ~3 f" `6 |

plt.xlabel(x)

( L* f8 G7 q4 U" L

plt.ylabel(y)

* S; e2 N8 o) ]3 M

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

; y( ^" ~# A2 Y# u/ s

nbins = 200

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

# H needs to be rotated and flipped

5 [+ x2 {9 p5 l$ o

H = np.rot90(H)

H = np.flipud(H)

7 q6 {7 H% `8 V. w% V- d! V$ M: D7 M

# 将zeros mask

1 |+ R1 L- _1 y7 \% H: {/ e

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

8 {7 T7 W2 {# q S. K

# Plot 2D histogram using pcolor

fig2 = plt.figure()

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

plt.ylabel(y)

cbar = plt.colorbar()

cbar.ax.set_ylabel(Counts)

. }7 @7 R% s3 k2 \

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

. Z4 X& T' q3 j; Z1 @9 `( C- [; [

plt.show()

3 B8 d9 ~9 A6 Y2 ?( ?

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

, D& |; c* R# w1 P% P

7 i S' ^) k" v8 |8 O! V

Example 2 ：双Y轴（Python）

: t* p" M& I& Q' I7 X' [+ l

import csv

) E2 y! Z& p7 X( g0 O4 G

import pandas as pd

import matplotlib.pyplot as plt

from datetime import datetime

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

% A* o3 r2 {8 G0 b% L

time=data[date [AST]]

sal=data[salinity]

tem=data[temperature [C]]

1 m) U+ R8 N+ V. p( l

print(sal)

$ e( ?$ `8 W) _. A! K

DAT = []

for row in time:

# r% w% {4 ? N$ Y* V

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

#create figure

5 V; A/ n2 i7 y r" Z

fig, ax =plt.subplots(1)

$ r$ b4 D2 R- ~% c

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

& y4 n2 n$ I4 z2 e

plt.xlabel("Date [AST]")

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

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

b8 r! [( x5 [; r* C$ `

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

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

fig.autofmt_xdate(rotation=50)

; N/ Z7 h5 ^- A0 I) ~: @/ W' @

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

plt.twinx()

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

3 a( u6 m: K9 u# b

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

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

7 v- v8 }! y5 n6 k- z

#To save your graph

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

plt.show()

- n; Z: ]9 S g9 E8 c- w

! M' S; L/ {- ^

Example 3：拟合曲线（Python）

% ]) a9 x! k, Y3 h& P% v) j% [* H) L

import csv

import numpy as np

3 g/ G" e$ \) i

import pandas as pd

& z5 {. i1 q9 s- E' \6 r# e

from datetime import datetime

+ O& V% g# h c N& W3 p

import matplotlib.pyplot as plt

import scipy.signal as signal

4 R1 g/ X" R2 L1 X/ _* A# D) o& _

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

time=data[date [AST]]

temp=data[temperature [C]]

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

DATE,decday = [],[]

) \3 {: k2 R% q# X

for row in time:

) N, R9 O: _# j: d1 h; ]5 S5 a2 A

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

DATE.append(daterow)

1 S" e. b( p/ h& {8 ?

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

# First, design the Buterworth filter

N = 2 # Filter order

( \7 J' H+ s5 W+ l6 {4 Z5 g

Wn = 0.01 # Cutoff frequency

( k l$ {& ^. t& c+ s Z$ X! ^

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

# Second, apply the filter

4 @/ Z" J0 B: M' [: G7 m

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

, C. k, a6 N2 O; G' P

# Make plots

fig = plt.figure()

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

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

plt.ylabel("Temperature (oC)")

* I4 ^" G+ ]- ?

plt.legend([Original,Filtered])

4 {0 W- S* N+ g' _- w, p7 _

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

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

* e, W' q& `. L1 Y3 i0 a

ax1 = fig.add_subplot(212)

9 h8 N' W" \4 {( U

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

plt.ylabel("Temperature (oC)")

. \, i; |& Z6 U2 K ?! P

plt.xlabel("Date")

plt.legend([Residuals])

+ o2 }/ e) ~' o- e: l& U% o& u

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

plt.show()

+ b1 j$ Q& H- d8 t/ ^( s

Example 4：三维地形（Python）

# This import registers the 3D projection

* ~+ ?' h( e2 d7 k

from mpl_toolkits.mplot3d import Axes3D

- J# W5 k% U" }' C/ S: o. u

from matplotlib import cbook

from matplotlib import cm

& }, a' A# ~" b1 z) `5 r q

from matplotlib.colors import LightSource

& l! c+ N6 H" k4 w7 r" K" x

import matplotlib.pyplot as plt

import numpy as np

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

with np.load(filename) as dem:

z = dem[elevation]

7 Q% L( C5 n4 O, l

nrows, ncols = z.shape

' p1 M. p& }% f9 ]0 X0 W, i. \

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

1 l. n$ `9 R" \! c

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

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

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

2 j4 |6 T$ o+ @2 Z+ |# C, T

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

, ~, i4 ]/ D5 F$ s3 O1 Q+ i

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)

' y1 Q1 D. l+ `# j7 h

surf = ax.plot_surface(x, y, z, rstride=1, cstride=1, facecolors=rgb,

' P4 n; P* d: a0 }5 t( {

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

, g$ n! p$ N/ N. W0 Q/ K) N

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

plt.show()

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

0 a( A1 L5 ?8 F$ B! N 2 C6 _3 `3 f+ x& x1 y- L

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

- ]- c" C0 D% W% g8 \& V

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

- M8 t" k! g) q" G+ i6 ^ - O3 I6 U7 r% m I3 w& _

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

/ X- a6 X; W) J4 h4 M, |

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

# W) Z" J; w! U. h; {' U+ X

% n p. m, s* r& B% W" H

: l$ m, i: r% w b5 U5 @* `, k / K+ z/ [2 D& `# q, M' N3 V9 z5 t 5 E& U3 E& s" X: f/ U 4 G* @: K/ B5 g+ z