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

) R0 V i- r: r/ T9 f

; u) c$ ]. e2 Q, n4 Z3 D

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

6 ^, j5 C1 X3 `9 g k* d* {+ M

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

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

5 |" P' G& [. Y/ K" P0 j

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

; B) [# h S3 c% y

import numpy as np

import matplotlib.pyplot as plt

# 创建随机数

n = 100000

x = np.random.randn(n)

' h+ B# @! p8 M

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

/ d/ e& j1 E$ X! Z! }1 s

fig1 = plt.figure()

5 J& x$ p* v0 d3 x6 S/ N3 I1 x

plt.plot(x,y,.r)

: ]& y+ ~3 _- Z' @. o- G

plt.xlabel(x)

0 C, S/ X, n4 u* [: B. [

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

; o3 l& g8 ~4 u/ H8 X

H = np.rot90(H)

H = np.flipud(H)

# 将zeros mask

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

# Plot 2D histogram using pcolor

$ X I9 E3 F9 b0 F# { r# h4 F

fig2 = plt.figure()

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

plt.ylabel(y)

cbar = plt.colorbar()

% ?/ M( Y' D4 X

cbar.ax.set_ylabel(Counts)

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

plt.show()

" L u6 o6 t: f# s5 B- x. ^% ] L

3 v5 z- W% M/ u2 o) J- M

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

Example 2 ：双Y轴（Python）

# N. J; m J: T& D

import csv

8 i+ G F/ x5 n

import pandas as pd

import matplotlib.pyplot as plt

: z' K& k0 W4 j9 d

from datetime import datetime

4 n- c% ], [" M: L8 c3 ?. ?

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

4 s/ x6 D _9 S& i1 w) ]6 {4 ]

time=data[date [AST]]

sal=data[salinity]

tem=data[temperature [C]]

print(sal)

DAT = []

" _/ ~; {: u" G: a4 r1 O# a

for row in time:

0 z# a( w0 o- o: a

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

#create figure

( D1 ?. K7 `: Y3 S

fig, ax =plt.subplots(1)

; Z/ a7 S. X# R- c' B p9 j4 |$ i

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

) F6 H5 f M7 v5 n

plt.xlabel("Date [AST]")

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

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

" {3 ?" M1 X2 D# @" j

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

: h4 `* _( M2 Y5 n# b2 U. y

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

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

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

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

9 q' v6 `( U1 N* B! O2 F

#To save your graph

7 A& g* `$ b3 T& D" m* a

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

plt.show()

" s: D+ Z7 L; [3 j- @ G

Example 3：拟合曲线（Python）

import csv

9 `# o+ X: q! j1 Z. @' f0 _6 b

import numpy as np

import pandas as pd

from datetime import datetime

import matplotlib.pyplot as plt

import scipy.signal as signal

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

% v# l$ F+ ]* u9 s' I, b; ~% M

time=data[date [AST]]

temp=data[temperature [C]]

0 D& p0 a3 _, c! B+ j3 f5 w1 i

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

! E+ @# ^) ^: K$ o

DATE,decday = [],[]

# ? N1 i# H$ R7 C2 c# B

for row in time:

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

DATE.append(daterow)

7 } @; v5 r* \0 Z

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

- `8 A( j( L) j' \- [0 L

# First, design the Buterworth filter

7 r) o+ |* f$ {1 T- S# _! T

N = 2 # Filter order

: C, e- g. E; S! P# @( K

Wn = 0.01 # Cutoff frequency

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

# Second, apply the filter

0 b( R& C2 J: W' }9 N4 K3 V

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

' p/ [' W1 z2 z$ j* p, k

# Make plots

fig = plt.figure()

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

% p- y* g% l0 g7 _

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

plt.ylabel("Temperature (oC)")

# M( `" o8 g+ \) K$ U

plt.legend([Original,Filtered])

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

1 m0 P! h/ f' f* w5 {

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

8 E! A, v) p$ n6 q/ C2 [+ |: s

ax1 = fig.add_subplot(212)

) P9 y: E: J3 P9 ~$ _

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

- h% d# r5 r8 A$ F5 o

plt.ylabel("Temperature (oC)")

& l t( f- u2 X+ ^# i* R

plt.xlabel("Date")

plt.legend([Residuals])

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

plt.show()

; D- f, e7 n9 o$ J2 ~1 C B( K

Example 4：三维地形（Python）

# This import registers the 3D projection

/ y# G. x7 w* Y$ L3 V, g

from mpl_toolkits.mplot3d import Axes3D

from matplotlib import cbook

from matplotlib import cm

1 R, U% Y9 W' ~6 \- U# E D+ i

from matplotlib.colors import LightSource

import matplotlib.pyplot as plt

. D/ X; R2 _( @# Z& r

import numpy as np

$ b* V1 E. z+ U! y+ N

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)

' h$ B! c4 k( G1 s

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

6 r8 | l7 ]" |, Q6 L, R

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

$ |9 d4 P- u0 }; a. T

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

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

ls = LightSource(270, 45)

; r) G! f3 v) j$ ^- n

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

E0 }6 }1 I r2 R5 b1 ]& g! H5 X

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

7 K# v& l/ A) g0 q* s

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

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

- q; e% A! q( @% y

plt.show()

: w W7 w- m9 A+ z/ c

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

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

; C9 w" y3 L+ m0 j

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

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

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

7 {/ H, N# K0 Y7 f

$ T' D1 s' R5 V* A/ b8 o

3 S7 w B. p' f/ B" Y. N , D7 n7 ?1 B* `% ]1 o, ]% I9 W