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

0 d, g! O0 k2 p& B0 A' C2 N% j+ x. W3 R

5 O% k$ ~& D0 s1 t( Y4 Q8 R$ h7 [0 b

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

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

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

. ~; w: k9 x Y6 [

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

import numpy as np

import matplotlib.pyplot as plt

# 创建随机数

n = 100000

x = np.random.randn(n)

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

fig1 = plt.figure()

& w" @1 }1 N# T( T$ C

plt.plot(x,y,.r)

$ l4 @% a+ n1 ?3 N( B

plt.xlabel(x)

1 A4 V. C2 @7 v- H. c% _4 I$ `

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 ^- l2 f; u5 V# Z

H = np.rot90(H)

& d5 I& c# _4 }

H = np.flipud(H)

7 X* s" H) U: s3 b

# 将zeros mask

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

# Plot 2D histogram using pcolor

fig2 = plt.figure()

% P+ F0 @& q+ {8 F

plt.pcolormesh(xedges,yedges,Hmasked)

4 c2 s3 v* ?% `: I* @5 U

plt.xlabel(x)

2 b; h; z9 l, [) H

plt.ylabel(y)

9 l7 `. y6 m- X+ B

cbar = plt.colorbar()

* E+ w: e/ t. C6 Q* E- R0 ?) P

cbar.ax.set_ylabel(Counts)

* [% @, U1 Z, K/ R- E1 t9 |; x

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

plt.show()

) z) {3 O. n2 }8 U$ ^# G

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

2 |& L, f0 a$ O; w, Z& D

- S+ k5 m6 o) u* w. Q

Example 2 ：双Y轴（Python）

import csv

import pandas as pd

$ ^( r! [5 q! q1 S

import matplotlib.pyplot as plt

4 t: I8 Z' e( H+ o

from datetime import datetime

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

/ y- B3 U6 T( o' Y3 s/ g& W

time=data[date [AST]]

sal=data[salinity]

tem=data[temperature [C]]

print(sal)

DAT = []

& D3 d% }! e( v2 P4 @1 \/ q

for row in time:

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

: [$ ~+ k* Q' ` x; A: @9 m

#create figure

) k# l6 E8 }0 d

fig, ax =plt.subplots(1)

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

( x5 A. t) x8 R: a: \# x+ d

plt.xlabel("Date [AST]")

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

D; |. H1 d# y/ t+ Q

fig.autofmt_xdate(rotation=50)

- F( ]" B; A# J. C! I5 H9 h

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

plt.twinx()

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

1 H- f- p2 i& L/ g& V8 d ~

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

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

#To save your graph

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

, c) q. _$ J( c8 T, R" e& w6 N

plt.show()

, C* `5 u- o. \6 C8 d3 h7 h

Example 3：拟合曲线（Python）

3 Y2 k# e) C& a9 f/ ?7 `

import csv

import numpy as np

* w6 {7 }2 r# h6 f X) J, @

import pandas as pd

from datetime import datetime

import matplotlib.pyplot as plt

2 _8 H6 Z( T, w' B- n

import scipy.signal as signal

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

time=data[date [AST]]

temp=data[temperature [C]]

$ S% w& P# N, m" o

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

8 n8 z) E; F& b" \

DATE,decday = [],[]

' ~1 |- V: X5 c4 `% X

for row in time:

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

N = 2 # Filter order

Wn = 0.01 # Cutoff frequency

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

; W' Z% _5 ?9 w7 g8 q9 E

# Second, apply the filter

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

# Make plots

fig = plt.figure()

. B$ `8 @9 |6 @9 j- X1 I* o

ax1 = fig.add_subplot(211)

& P6 Q; q( B+ e+ K/ J, Z- E

plt.plot(decday,temp, b-)

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

8 m5 M$ a; \$ Q; ~8 ~2 c

plt.ylabel("Temperature (oC)")

plt.legend([Original,Filtered])

, T: |1 n6 d$ O5 I& Q' s: `

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

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

; K9 ]1 b* I$ r1 o

ax1 = fig.add_subplot(212)

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

plt.ylabel("Temperature (oC)")

9 i; t4 @: p+ O

plt.xlabel("Date")

plt.legend([Residuals])

7 d* C F9 K+ o

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

plt.show()

/ D' [ ^3 B4 m* v

Example 4：三维地形（Python）

# This import registers the 3D projection

9 Y) K; X7 c& X; S# W

from mpl_toolkits.mplot3d import Axes3D

from matplotlib import cbook

& C' v$ d8 y }; F, m+ Q% J

from matplotlib import cm

from matplotlib.colors import LightSource

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:

. R' C# S% C7 L9 u& `+ Q9 W0 L

z = dem[elevation]

nrows, ncols = z.shape

8 a t2 q& T6 `: { b

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

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

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

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

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

6 z% o9 Q. S' I. V# f8 m/ z

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

ls = LightSource(270, 45)

6 c' E/ H s8 ^7 P* ]" ~

rgb = ls.shade(z, cmap=cm.gist_earth, vert_exag=0.1, blend_mode=soft)

( x& M- u# ?2 ~* U! z+ a

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

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

- E" D9 d5 \3 f# B4 S' K, q

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

9 S& u4 k: N8 ^9 G

plt.show()

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

. Y e. }7 A4 N z

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

' @. o) O n9 l/ Z: `) m 4 d! S- w1 h# B1 z7 q7 x

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

; \* l4 v- V( d

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

2 ? Y5 V/ J9 H * m& o, h9 {, Q

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

8 ?( t# p3 I x/ F# P

/ R+ O& \9 y6 I/ h; U1 M% D) X' c

8 R* a) T5 J. w9 @2 A( d4 O

y P4 H# }! l: C/ X + o7 c% K' z6 F' T3 U- l