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

8 t4 p. W. z$ {

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

# |( [+ P9 h8 x* O0 w

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

" E$ z$ k& [4 w4 g

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

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

1 ?2 @" E$ T# Z) L, G7 q# q( `

import numpy as np

0 Q( v* n: d: I+ X

import matplotlib.pyplot as plt

: X2 ]0 [8 O2 Z/ [% q* I3 |

# 创建随机数

5 h. U6 u& k! ~# a5 k9 e9 ?

n = 100000

x = np.random.randn(n)

0 D9 x, H( ~- ^0 q) f

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

& L# R+ a6 ~ L7 Q+ A

fig1 = plt.figure()

plt.plot(x,y,.r)

7 z2 j9 Z7 E& J6 {, G/ c1 m8 y {

plt.xlabel(x)

plt.ylabel(y)

8 }$ h5 ~% S0 T% z. a

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

$ y; W/ M% S: Q4 P

nbins = 200

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

- _' v: }% H# M' T; v$ U9 I a+ d

# H needs to be rotated and flipped

0 K* p' @8 _3 d% ]

H = np.rot90(H)

& F4 X! k8 b' f: r7 W/ q

H = np.flipud(H)

# 将zeros mask

: P% p) h; |4 z$ a. ~

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

7 m3 x' k2 V6 e( X

# Plot 2D histogram using pcolor

. q7 D# g( f6 f" Q& Z

fig2 = plt.figure()

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

- [+ T* L* O! ?" w& f9 l

plt.ylabel(y)

0 g3 G# j1 ~: J0 r9 i* g0 B5 T

cbar = plt.colorbar()

) q- n% e! j: S% R1 q" e

cbar.ax.set_ylabel(Counts)

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

0 l6 t" e0 E0 b" e9 k+ C

plt.show()

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

2 {1 w& `, e( S! V0 Y

* p: u2 [# U' S% M9 Y

Example 2 ：双Y轴（Python）

) f S Q) C" o: i! l

import csv

7 |7 n1 x6 T$ t3 q

import pandas as pd

( z# P+ E3 C c% X u7 @3 z

import matplotlib.pyplot as plt

0 L2 g) B6 g. J

from datetime import datetime

5 S1 c) l2 p3 N" ?3 [. |. z5 t6 E

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

+ Y' j8 H. F7 i' ~3 r) E

time=data[date [AST]]

sal=data[salinity]

tem=data[temperature [C]]

# i* o9 ?, w( d7 l0 l

print(sal)

DAT = []

$ m$ J" X# G/ G, G+ b0 t9 S

for row in time:

; X) a* W% l5 @

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

; ]1 y8 u: N* k( m

#create figure

fig, ax =plt.subplots(1)

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

$ _( I! x$ @5 @- f

plt.xlabel("Date [AST]")

9 ~; p& d1 x: \4 T: k K

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

, L0 z! ?' ^5 v) h' f/ @

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

2 [5 k/ C6 E l/ t I6 N! h

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

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

fig.autofmt_xdate(rotation=50)

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

4 S& S! K4 v8 G* @! t( V

plt.twinx()

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

B- k) Y# G e U8 p8 {

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

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

8 u0 ?/ V' [/ C* L+ S+ X, M: K S

#To save your graph

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

plt.show()

) h* ^2 X& d9 F" W6 z. S; [* n

Example 3：拟合曲线（Python）

" {- Y5 m' z: B5 z/ w7 e% _

import csv

' U( c0 j h3 ?' l) r

import numpy as np

2 f% |4 ^0 B# K! A% J+ W% g7 Y4 K

import pandas as pd

from datetime import datetime

0 e% P7 |5 }$ ]4 v

import matplotlib.pyplot as plt

import scipy.signal as signal

. e+ d$ s5 _) A

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 = [],[]

1 e/ d7 C/ M. r2 N

for row in time:

6 l" k* K0 C/ `6 v: Q

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

DATE.append(daterow)

4 w( } M0 r5 G0 j

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

# First, design the Buterworth filter

: d' R" h1 g9 O' k7 q( }1 U

N = 2 # Filter order

Wn = 0.01 # Cutoff frequency

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

# Second, apply the filter

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

* X9 J8 u8 G: _0 A) ^

# Make plots

5 ~) b6 {7 [3 w' q3 `8 R/ E: Q

fig = plt.figure()

9 i3 I c0 }* l! \% B

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

$ |1 @" ] M# q1 D$ ]

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

plt.ylabel("Temperature (oC)")

, x" m% u: D1 R) c1 ~3 \, j0 k

plt.legend([Original,Filtered])

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

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

$ J0 ?$ g) q3 K/ Y& e

ax1 = fig.add_subplot(212)

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

4 ~* a# | W2 B' L$ D0 a9 e

plt.ylabel("Temperature (oC)")

plt.xlabel("Date")

plt.legend([Residuals])

. f8 A& z! O! w) F" O

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

/ I; e1 F8 z/ A0 L L& [5 D

plt.show()

7 S) }. `, J" ^9 o g" x, W

Example 4：三维地形（Python）

# This import registers the 3D projection

from mpl_toolkits.mplot3d import Axes3D

from matplotlib import cbook

! N( T8 N1 ^# {4 J4 b

from matplotlib import cm

& W* z [. ?- {" F$ }! Z1 E

from matplotlib.colors import LightSource

import matplotlib.pyplot as plt

7 t4 Z$ E3 }! H8 [' a

import numpy as np

V ]' P( F' R

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)

' P2 ^( q+ D7 k2 O' j7 q

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

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

2 h) R! k# Q! h h: j) O, v3 ?

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

6 Z) o s0 P1 o+ W

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

7 D! G+ A% s. F( {5 l

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,

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

9 W: h% d8 U# n

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

plt.show()

7 O7 x+ {0 t+ F, v3 K

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

! ~; ~$ {! T6 d( A( V- K

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

! h2 f6 h( J" k( h3 k 3 ?% H- x0 m" J6 l% `: u: p

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

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

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

: u) ]3 ^: c7 a9 w

. r4 @+ B9 p+ [# q2 @ ! }% ~, L$ D4 W5 W" U0 g( P 9 G* _( o( N7 `# K4 D , [7 {; O, k& I" \7 W2 C* K $ ?) `* m4 }) U