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

9 n% H: W# \6 a

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

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

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

i5 e( v" c( C

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

: b: r& }/ {8 }! { t' A

import numpy as np

. A) u6 h! J! V* ]& @8 O7 Y$ K3 h

import matplotlib.pyplot as plt

# 创建随机数

n = 100000

4 b& b" Z1 V l, ]5 ?, W% r% e

x = np.random.randn(n)

+ @8 P9 {# k9 J

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

fig1 = plt.figure()

plt.plot(x,y,.r)

plt.xlabel(x)

plt.ylabel(y)

! T# u: J! J2 }% T

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

J2 U) y* b2 {) ^# e/ X

nbins = 200

% r% {5 E( x' @

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

0 u4 [8 E, g; [ n, {; O+ A

# H needs to be rotated and flipped

! H/ m- [5 U8 r' O3 t

H = np.rot90(H)

' E" } F' a* H5 v) {! a* H1 u

H = np.flipud(H)

# 将zeros mask

6 R3 J- d H: f5 v" J

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

1 D: r9 K# D. ^9 C* \

# Plot 2D histogram using pcolor

fig2 = plt.figure()

# V1 C" W; [; V1 v

plt.pcolormesh(xedges,yedges,Hmasked)

& e7 p1 a- X; @9 n6 S

plt.xlabel(x)

plt.ylabel(y)

cbar = plt.colorbar()

cbar.ax.set_ylabel(Counts)

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

plt.show()

4 {( x* o, V" m- y

; V: \: v2 M3 D1 v" J! x, d6 v

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

' a+ n' o2 Z( K8 X9 i& t u

. R% L9 ` _% b7 V7 |" C

Example 2 ：双Y轴（Python）

import csv

% [( |, T. Q8 G! P: b4 M+ K! z

import pandas as pd

, } H& j' q$ X8 }3 R8 ~. Z1 {

import matplotlib.pyplot as plt

from datetime import datetime

4 T* v a9 V8 t) L& T

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

time=data[date [AST]]

sal=data[salinity]

tem=data[temperature [C]]

print(sal)

" x, L' E! ]' T, _* w

DAT = []

for row in time:

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

. Q6 A P) V: O* ^' Z4 A

#create figure

0 d- i R! O! m1 r! ?

fig, ax =plt.subplots(1)

& D! Z3 M- x" z# N& C

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

+ O, z, J$ f: G3 G' B2 i+ w2 G

plt.xlabel("Date [AST]")

3 r$ L5 S4 \ Z# C

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

- @/ o Y, N1 x" q5 ?! ?$ t

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

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

4 S# S" }2 n0 I" l) x

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

. x4 \$ I! Z* y2 m0 a/ R4 ?. M

fig.autofmt_xdate(rotation=50)

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

2 O- H" I x5 p/ `4 a! M

plt.twinx()

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

0 v8 n, @; j0 J/ w

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

8 |+ X/ {3 x7 s- Y3 ^0 W2 G) y

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

# ]4 l0 e- c3 J7 z, S# |5 K, b$ g

#To save your graph

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

+ i0 }0 v5 e. j( L

plt.show()

/ a: ^8 `2 V( B% b, m

Example 3：拟合曲线（Python）

import csv

& y; |4 H3 Z3 L V( B _ _" v

import numpy as np

) K: h# `/ `* \

import pandas as pd

from datetime import datetime

% Z( F2 Q+ p, A6 A* l9 O

import matplotlib.pyplot as plt

0 k4 X0 ?! n4 W( `& Z% w

import scipy.signal as signal

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

* A' M6 i6 c0 K/ c5 d+ }) I

time=data[date [AST]]

2 Z* {7 L) }9 |, j8 @

temp=data[temperature [C]]

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

! H8 R$ @2 S8 d' Q s! J4 X" O

DATE,decday = [],[]

: a+ W8 J0 W4 f1 `5 W" M3 T1 Q

for row in time:

( L/ B; U, y3 V0 Q$ \3 w6 s, e

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

$ O3 C. d, d) c2 o' I

DATE.append(daterow)

$ W8 Y$ |0 O/ _' ?/ g( L

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

' S: a5 T. ?2 Y' z% C7 z1 _

# First, design the Buterworth filter

N = 2 # Filter order

Wn = 0.01 # Cutoff frequency

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

4 Z' W+ ]) b0 X% ~7 j8 x# y V

# Second, apply the filter

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

# Make plots

" i* ^3 d* y7 M( X5 T

fig = plt.figure()

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

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

5 k( O' g4 t$ G; `8 A. ]1 m! k

plt.ylabel("Temperature (oC)")

plt.legend([Original,Filtered])

" S/ R1 _+ n* j) }: K$ _% [( o0 d8 K

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

2 `; Y# l- n) {4 b6 z4 g* L2 V

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

ax1 = fig.add_subplot(212)

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

plt.ylabel("Temperature (oC)")

1 K! N& o1 s! K$ L! U. y

plt.xlabel("Date")

plt.legend([Residuals])

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

plt.show()

Example 4：三维地形（Python）

" l [ b4 x3 `4 L8 Z3 h* u% m

# This import registers the 3D projection

from mpl_toolkits.mplot3d import Axes3D

from matplotlib import cbook

, m& \, x% ^9 C( `' b) l3 [. m2 Z! Z' [

from matplotlib import cm

from matplotlib.colors import LightSource

import matplotlib.pyplot as plt

import numpy as np

0 x& o; ?" v2 A9 t6 C, Z# f" a

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)

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

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

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

0 v, k: d/ c$ n4 L! h6 Y

ls = LightSource(270, 45)

+ Z a6 ]8 o. B9 ?/ r

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

0 i5 [0 w. P5 z4 t# m% V

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

- I$ }- m6 h6 u, P! l& N# @3 k( `5 N

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

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

7 @; w! U% ?8 J5 u

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

# d2 A+ }6 T$ U0 S2 X

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

3 M3 W* Z6 P6 y2 v

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

/ U$ t+ C' I; X, f, z" h" |

$ O" @( }& @- _& o2 I& L2 H" D T % l: v4 w. c+ R. W8 o$ w + B: a1 g1 s% v x8 G* J2 W