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

+ H/ k' O) f6 U: w

" z3 ?( ?) r- k; s x

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

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

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

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

. U- S1 U: M9 ^% E9 z: {; ]% u

import numpy as np

: N5 V, t& d4 `! \$ V

import matplotlib.pyplot as plt

# 创建随机数

n = 100000

7 r v/ }2 y# m

x = np.random.randn(n)

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

fig1 = plt.figure()

plt.plot(x,y,.r)

plt.xlabel(x)

plt.ylabel(y)

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

! h* n1 `5 J: p9 [$ q6 ~

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)

# 将zeros mask

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

, _# l' s" B0 I* A5 c" P0 ^) M! V

# Plot 2D histogram using pcolor

fig2 = plt.figure()

3 R5 S% Z5 k+ j0 S5 f

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

plt.ylabel(y)

cbar = plt.colorbar()

cbar.ax.set_ylabel(Counts)

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

plt.show()

& @; {8 w6 c7 [$ o8 _, x# ^/ @. n6 J8 C

! u' w2 ]" Z% K+ i

) i$ w2 n) P) B6 H% ]" D

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

9 e" \" i4 @1 d6 M# y

Example 2 ：双Y轴（Python）

4 d- h! M! C) a$ C1 E5 a1 a& s

import csv

# }2 N+ ]8 y" S+ P0 b

import pandas as pd

8 b& _- Q5 T! \

import matplotlib.pyplot as plt

from datetime import datetime

$ Z% N, r8 H+ G

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

time=data[date [AST]]

sal=data[salinity]

tem=data[temperature [C]]

" L4 ` W! |9 Z5 b, O) m. E

print(sal)

4 A: c8 V8 R& D5 k( \% b4 k/ s

DAT = []

for row in time:

9 h; g* ~ Q D# C& g- d6 f1 D

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

#create figure

fig, ax =plt.subplots(1)

$ U8 Y+ x/ M2 p7 W+ |+ w7 Q

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

plt.xlabel("Date [AST]")

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

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

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

- q$ I& ^3 B( K

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

. ?! r/ L# w i# x

fig.autofmt_xdate(rotation=50)

, H7 O1 I0 Z2 c5 E- L5 M/ z

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

_, T* Z. p* S% X& v" X& |# u

plt.twinx()

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

$ T2 N5 [5 J# r w+ E: x# e

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

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

8 J' [2 X; ~. `, P: b) m( y9 @

#To save your graph

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

$ w0 `5 ^3 O! L& }( B8 w

plt.show()

, e3 }/ H9 r( T+ j$ |

3 T* B# a' G! j7 h/ b

Example 3：拟合曲线（Python）

8 x/ f5 P: s6 M$ h7 t. X% P

import csv

import numpy as np

import pandas as pd

from datetime import datetime

: G5 a; d! V: f- ?2 w

import matplotlib.pyplot as plt

import scipy.signal as signal

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

+ L; }' E" P! @0 ?

time=data[date [AST]]

temp=data[temperature [C]]

1 q" K' s5 {6 N$ C; N k( l6 S

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

DATE,decday = [],[]

for row in time:

* c- m6 M4 H8 K( I

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

|/ h+ ]5 H+ N3 P: b4 D

DATE.append(daterow)

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

# First, design the Buterworth filter

N = 2 # Filter order

2 \4 h3 f+ p D

Wn = 0.01 # Cutoff frequency

) k7 i( C8 [7 H9 f

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

5 x1 p# X6 j) _$ ?) n3 u W* m

# Second, apply the filter

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

& |/ o X+ V9 K/ \# n

# Make plots

8 v. l* R! M5 f y, I5 E, f

fig = plt.figure()

8 X& I. u# m Q( g) M3 w

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

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

plt.ylabel("Temperature (oC)")

plt.legend([Original,Filtered])

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

?& t2 J0 s! e( A# l6 R" m

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

ax1 = fig.add_subplot(212)

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

" g$ n! [! |0 w" b+ M4 o& V

plt.ylabel("Temperature (oC)")

* y+ X0 }( t: W

plt.xlabel("Date")

plt.legend([Residuals])

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

plt.show()

Example 4：三维地形（Python）

# This import registers the 3D projection

from mpl_toolkits.mplot3d import Axes3D

from matplotlib import cbook

from matplotlib import cm

, g/ x8 v3 C* \8 G

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:

6 D8 A1 M5 K) x% ?6 c" Y% W5 E

z = dem[elevation]

nrows, ncols = z.shape

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

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

! l- k; ^: n7 I* ^ C

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

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

, P3 i$ c' p- f

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

; O3 E2 l4 b0 V

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)

; Q4 ~# z/ S! j, m

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

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

+ u8 N# n" n) D/ F

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

plt.show()

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

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

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

* R- \8 x( j5 g% k! k$ k$ z! I

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

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

5 A l0 n+ n) ^ ]& T

" j8 V/ O# C/ d& ^5 `

* K& Z9 f2 B- E/ Z" B * \" h* \* [5 u ; {- M+ N5 W" o( m" J$ k9 Q. W 5 { s! l! V1 ^0 a