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

' N; I' L; N& A7 o6 B

1 ^* \! I7 l1 j4 v4 \- E a3 F* z

: z* c1 ?2 _0 O, c, o: I f8 p

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

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

5 w& W, e* V; B$ ^

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

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

% o4 H. L9 N9 ~* m6 E) W$ B

import numpy as np

import matplotlib.pyplot as plt

1 P% j( N6 K k% p; t+ {

# 创建随机数

n = 100000

9 S" t3 L7 S* |9 E ~" L, U

x = np.random.randn(n)

6 g$ m9 s! Z5 l8 P) m* {

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

% \+ _0 v0 z) N$ y( @

fig1 = plt.figure()

9 V8 ?0 s8 ?- E/ C! N$ C. `

plt.plot(x,y,.r)

% ~/ p$ b& D, c% d8 l

plt.xlabel(x)

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

H = np.rot90(H)

H = np.flipud(H)

% c* Z0 n& A8 S2 p/ Y4 E

# 将zeros mask

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

' Y" n" \( x( O4 t! h

# Plot 2D histogram using pcolor

. }+ x- M7 {5 D8 ]5 H+ E1 @, j% R8 P

fig2 = plt.figure()

- \3 V) E, i* M" O

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

. J5 l( _1 ?, x- {

plt.ylabel(y)

5 Z5 ]) B: K( o8 T( v# E

cbar = plt.colorbar()

4 h; q6 K. e! \; _) ~

cbar.ax.set_ylabel(Counts)

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

plt.show()

: p$ r0 l$ M D: P3 R- P: f

9 Z' k1 U: [ X5 Z! P# ]- R

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

$ T4 r# I$ [7 C. ~

$ }6 c) g x* B; t% t. o

- T( `+ A3 U& ?' w

Example 2 ：双Y轴（Python）

5 |+ P( z, e8 C z; X5 z. W

import csv

import pandas as pd

- J6 z- I9 z/ G

import matplotlib.pyplot as plt

from datetime import datetime

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

, c- o" [! G+ B4 O$ u1 s) u

time=data[date [AST]]

" j0 ?# G" c) Q

sal=data[salinity]

tem=data[temperature [C]]

( l. ^! H6 I5 [3 p; ^1 _ y

print(sal)

/ x9 J% O( @ U3 c: k

DAT = []

for row in time:

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

#create figure

9 H$ O) c. g; B( c( s

fig, ax =plt.subplots(1)

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

plt.xlabel("Date [AST]")

' ]3 Z$ n$ T, ^, n* b1 A C$ L

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

/ m+ ^: \9 H, [8 L# l( v+ R8 x

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

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.

plt.twinx()

$ i( Z9 C( u+ A" p% k

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

4 G$ c- R- l- r4 N

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)

plt.show()

3 M3 h6 w+ L6 U. m3 _' L

Example 3：拟合曲线（Python）

import csv

import numpy as np

import pandas as pd

! w( s; r; M5 ]0 T# z# K

from datetime import datetime

1 J; R7 g' w) s. i. j

import matplotlib.pyplot as plt

import scipy.signal as signal

+ p3 ~. L9 X! [, W m% ?* L

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

time=data[date [AST]]

temp=data[temperature [C]]

9 b4 g) k5 w& }4 Q$ w/ l1 ~

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

DATE,decday = [],[]

$ W+ ]1 g8 r0 A2 a& N

for row in time:

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

4 @* C4 R4 P+ X& A

DATE.append(daterow)

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

$ G) c9 R# e$ O

# First, design the Buterworth filter

) s- O6 ^/ U0 d

N = 2 # Filter order

Wn = 0.01 # Cutoff frequency

5 D" H. n3 E/ R7 ]6 g6 L! k( i

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

5 S, g, f( T7 g! x6 w

# Second, apply the filter

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

* l" f, N! N" x: ]8 J( S

# Make plots

' S# V) M9 u) ~- w& I" M% G! _

fig = plt.figure()

7 t1 g3 |7 T" i V

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

. }% A, n; n$ \

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

plt.ylabel("Temperature (oC)")

7 z4 a0 m3 v# `8 g: W% E# Q

plt.legend([Original,Filtered])

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

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

ax1 = fig.add_subplot(212)

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

plt.ylabel("Temperature (oC)")

plt.xlabel("Date")

plt.legend([Residuals])

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

$ b z' y4 D# L, i3 C3 j6 v

plt.show()

& | k5 `+ |. Y8 {( |

Example 4：三维地形（Python）

" G1 w1 @2 k2 W$ C& R

# 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

import matplotlib.pyplot as plt

import numpy as np

# v0 x2 B1 p2 M

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

& C9 r' O8 L4 V0 _

with np.load(filename) as dem:

z = dem[elevation]

+ T$ y- |3 J3 ^. ]% _9 L5 V

nrows, ncols = z.shape

# R0 K7 J5 }" N& g7 j

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

9 U7 l6 m& S* n8 T5 B

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

8 g% y n' N; E6 M9 P

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

: K# _" w2 A ?1 @

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

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)

) B' s4 O* A7 v; `& i

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

plt.show()

( x3 v; k: k# Y" n6 r4 B8 }. v6 R6 D4 g

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

1 Q7 W+ Z; e/ `/ s( ~

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

) a7 V. ]" L; R) @; j3 r9 ~

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

$ [3 K* H9 I, z7 K4 v; V9 y L9 L

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

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

1 p9 d6 u$ K" ^8 {

: f( |$ |: S# ^4 Y# x ~7 ^

; ~+ l# ^! o) T7 m/ S0 ? 0 r6 P7 K0 h0 L. P 5 W- X8 @; J+ C; y1 K, }" F" h: P