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

: `; B! E: | R. @

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

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

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

% I* }9 X! F* w3 l/ T. \

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

import numpy as np

import matplotlib.pyplot as plt

! @4 Y% M& |4 N$ o$ b' a, X& u

# 创建随机数

; |, i$ }6 _2 |

n = 100000

x = np.random.randn(n)

7 l$ I% W, Y# G) G, s/ }5 m6 R

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

* E$ O, g6 D& D

fig1 = plt.figure()

( ^) t# {- P0 C2 E8 x/ q6 }

plt.plot(x,y,.r)

plt.xlabel(x)

plt.ylabel(y)

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

7 l" Z3 r; I5 h8 o0 ]6 C8 K) h

nbins = 200

, ]6 O# q3 t7 k7 Z8 j# q1 I- w

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

# H needs to be rotated and flipped

6 I. Y& Z+ {* G9 Z* y! O: K( y1 t

H = np.rot90(H)

H = np.flipud(H)

# 将zeros mask

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

# Plot 2D histogram using pcolor

+ U7 I- X9 |3 A/ y. }

fig2 = plt.figure()

plt.pcolormesh(xedges,yedges,Hmasked)

+ j( K/ V+ E4 ]7 @

plt.xlabel(x)

plt.ylabel(y)

3 o( n/ ?0 ?) {9 m. Z9 \) }, C

cbar = plt.colorbar()

cbar.ax.set_ylabel(Counts)

, J5 V0 s3 b" y+ p8 G! m0 P$ b

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

plt.show()

" K# X0 @# @* z* r; N% ^# r

: ~9 W* J% X: @

1 V1 Q+ e" G; N

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

1 x1 [3 M1 x) N% q, F' }& h% W

Example 2 ：双Y轴（Python）

import csv

/ b: `4 `% }2 g0 r

import pandas as pd

import matplotlib.pyplot as plt

from datetime import datetime

4 g* w( D7 e; p0 v6 {0 e

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

time=data[date [AST]]

& A. @* J7 t Z$ l+ Q4 Y# U1 X' |. B

sal=data[salinity]

tem=data[temperature [C]]

* ]0 Q! m: T/ o8 W* g

print(sal)

DAT = []

# i ]: @/ r3 {: y8 P

for row in time:

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

7 b+ L5 k5 I! y# Y& y; Y ?1 c% G

#create figure

fig, ax =plt.subplots(1)

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

plt.xlabel("Date [AST]")

3 K2 S8 R6 O5 `4 g* Q9 K3 n ?

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

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

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

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

1 {8 g4 r5 I* `+ r/ v

fig.autofmt_xdate(rotation=50)

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

plt.twinx()

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

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

% V! D2 a. d7 r' m/ T

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

. h7 h+ X$ I) ^* U. M k4 d

#To save your graph

- t. l7 o. d2 \3 c

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

' n& n' `% @2 x2 Y/ R

plt.show()

* G5 R# r! o. k; @! s$ L. k: B

9 c3 G0 b8 C: } ~7 C8 u

Example 3：拟合曲线（Python）

import csv

import numpy as np

5 F6 o' v3 J& G- F

import pandas as pd

$ `/ r8 {9 M: c" R7 i: D5 i2 b/ U

from datetime import datetime

3 T- c: T" [' G/ b, }" ?

import matplotlib.pyplot as plt

import scipy.signal as signal

4 g, e. q) H( y% u

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

7 W' U3 t) X9 ?/ d5 ?. ^: e

time=data[date [AST]]

temp=data[temperature [C]]

$ n6 }7 ?' J6 I

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

9 V3 E: _% }* t9 u" p$ [

DATE,decday = [],[]

for row in time:

/ q9 N* o3 ^3 B. a

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

DATE.append(daterow)

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

# First, design the Buterworth filter

0 C, f- M* d& H L

N = 2 # Filter order

6 J6 o% h1 ^/ d

Wn = 0.01 # Cutoff frequency

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

# Second, apply the filter

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

# Make plots

, A0 d9 V" } I1 S

fig = plt.figure()

ax1 = fig.add_subplot(211)

6 ]& K1 w" U, i5 P8 W/ ^

plt.plot(decday,temp, b-)

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

plt.ylabel("Temperature (oC)")

& G2 U: Z; }% H5 f" {) }+ s. l

plt.legend([Original,Filtered])

; ]( p, p. z/ r0 y, @

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

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

ax1 = fig.add_subplot(212)

! k4 N# r# ?; L6 r

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

! z% I$ A# C9 }. C/ |, U

plt.ylabel("Temperature (oC)")

9 r2 t% e: D- o& A5 W# S

plt.xlabel("Date")

plt.legend([Residuals])

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

plt.show()

* W3 T+ y& i" \0 J1 B

Example 4：三维地形（Python）

0 G* c* C; v0 @

# This import registers the 3D projection

from mpl_toolkits.mplot3d import Axes3D

from matplotlib import cbook

from matplotlib import cm

from matplotlib.colors import LightSource

- z, g) g1 O$ c

import matplotlib.pyplot as plt

import numpy as np

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

2 R+ i- C% F, H3 K0 I4 e+ H3 F9 u; @" A

with np.load(filename) as dem:

2 ]1 Q( S) \. I8 g3 m

z = dem[elevation]

4 d8 ~2 X2 ^7 D+ a: Y

nrows, ncols = z.shape

U9 \) h2 Y+ ~3 s r7 ^

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

; l H6 @" V6 C

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

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

- ^3 Q8 C; M% R" k

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

( {' K# X B+ e* x+ U

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

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

8 S! T( X5 J Q& U) K

ls = LightSource(270, 45)

0 r% I5 T8 {( x4 ^+ C7 ]# H

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)

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

plt.show()

4 w6 ?) O q" b. N0 h

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

1 H6 z2 a6 v X: v

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

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

7 Z* M# U' s: s- Q5 T, y0 a / b( g$ S6 _9 G% j' e9 K2 _, A# t$ m

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

u! t. a$ ? T( I & G! N; @) t# {3 E

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

2 f0 b" J/ H9 J5 { O - z4 o+ b# \* ^! K: }0 M 9 |* F: ?: f- k $ f7 B: y" n6 L! W7 i5 U