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

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

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

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

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

import numpy as np

+ N( o5 G6 h% q) U, i% \& g

import matplotlib.pyplot as plt

# 创建随机数

2 v+ k9 f7 V/ p' F$ _3 _+ M8 N( v4 v# b

n = 100000

; ^3 }2 a' K3 }# \5 V* ~0 J% m

x = np.random.randn(n)

+ g+ S' c% V3 G: j$ O

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

fig1 = plt.figure()

plt.plot(x,y,.r)

plt.xlabel(x)

/ D1 y# Q9 Q5 Q- @" j

plt.ylabel(y)

1 g2 T+ m/ y s

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

nbins = 200

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

8 d; a& }. j% r9 \ n/ a

# H needs to be rotated and flipped

^) M( ^/ y6 T# k

H = np.rot90(H)

, e3 a! t& ?3 {6 Y4 F

H = np.flipud(H)

# 将zeros mask

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

# Plot 2D histogram using pcolor

fig2 = plt.figure()

& `- n" N6 g7 @6 C+ t6 p

plt.pcolormesh(xedges,yedges,Hmasked)

plt.xlabel(x)

plt.ylabel(y)

# e, W" w9 {& l' P& p4 F% h

cbar = plt.colorbar()

3 T9 _" v: @* }7 \

cbar.ax.set_ylabel(Counts)

. ?9 Y Q, ^$ C- ^3 i5 Q

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

plt.show()

0 ?8 _* D- ^( s3 @5 W

; Y' v1 F) R$ U5 W4 J! _

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

* z& [2 J9 G( t: u+ @

Example 2 ：双Y轴（Python）

* t4 v4 Y0 N8 m. j6 ?2 m

import csv

import pandas as pd

import matplotlib.pyplot as plt

+ m) n+ `$ r8 |' s, _' z" C

from datetime import datetime

1 r/ E i; h7 m% Z6 T

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

time=data[date [AST]]

sal=data[salinity]

tem=data[temperature [C]]

. `- N5 T+ P: r

print(sal)

DAT = []

for row in time:

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

/ U* N& ~8 [* `

#create figure

fig, ax =plt.subplots(1)

8 i; F$ _/ A$ o6 F' \+ v. U

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

) h5 p" K1 d5 Y" ^

plt.xlabel("Date [AST]")

" ^9 l1 s- D( C" t5 \# y9 ^

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

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.

4 v- I. W- [9 r0 I

plt.twinx()

5 O# M, Y9 _0 q, G/ k3 [+ R

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

( `5 ~( ~& d+ J9 N5 h4 M6 F

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

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

" z" ?7 N4 X+ ^) `" C+ h

#To save your graph

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

0 ?7 S# ]5 P9 Z( I& S2 k

plt.show()

- e! U7 J2 g# F& L5 o. \ ?5 g& ^

2 ]$ R0 e2 l. L* \$ e

Example 3：拟合曲线（Python）

import csv

& l( S* q7 S- x$ Z

import numpy as np

import pandas as pd

( y; H! m- e6 g k" I7 R# G

from datetime import datetime

import matplotlib.pyplot as plt

import scipy.signal as signal

2 x( h) B m) }5 J

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

0 i6 x) {- ?$ N% N- L

time=data[date [AST]]

temp=data[temperature [C]]

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

DATE,decday = [],[]

W; J* X- j' _" ~! k4 G4 E2 B( q4 M) K

for row in time:

' T( X! P+ J4 Q; |9 ^& T, m) s

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

DATE.append(daterow)

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

# First, design the Buterworth filter

+ [# ]3 {! v! ~; M; q* e' b

N = 2 # Filter order

Wn = 0.01 # Cutoff frequency

! [2 @ _) l4 a5 f6 `7 ~

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

# Second, apply the filter

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

# Make plots

9 [8 E% L+ x8 _8 A$ U

fig = plt.figure()

ax1 = fig.add_subplot(211)

plt.plot(decday,temp, b-)

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

plt.ylabel("Temperature (oC)")

4 Q3 r1 x" g1 F: I. w& ~

plt.legend([Original,Filtered])

4 O m+ h6 P" K; P/ w* x% o

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

9 e/ x: f# @# }4 [, B

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

plt.show()

Example 4：三维地形（Python）

/ e$ b% a) W9 g9 ~: g! n$ z! j

# 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

; z( Y- Z' G: Y+ e" |

import matplotlib.pyplot as plt

import numpy as np

& n2 k/ E6 `' `9 C0 B. J' A) q

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

8 o4 O. Y# }1 M! |( g$ {( ~

with np.load(filename) as dem:

' D* _1 s- E- ?9 t- V

z = dem[elevation]

nrows, ncols = z.shape

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

- w* N4 x4 t3 u# M7 V! I% j2 P/ ?

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]

, f- r, l8 ~7 O: a$ W: N

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)

: P0 P! ~8 {* h/ ^5 G' j

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

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

7 U3 E- m& Y; v, @( ~: _

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

plt.show()

. w" N0 E5 T: b * i5 Z. t& @3 Z

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

3 I6 C- t6 Y5 n6 k

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

* q D4 s1 v; P8 C6 P. C T2 c6 J$ d* M

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

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

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

P: ]! E- w& N2 v

% [9 U2 [) g; K# i5 K6 d

* U1 @. S6 G' D, A% d0 o. s! c0 u , }- C6 z! Z Z- V5 i" H8 U - R8 N, `0 Q# O1 I0 v I ?