快捷導(dǎo)航

如何用Matplotlib 畫三維圖的示例代碼

更新時間：2020年07月28日 09:49:07 作者：鹽味橙汁

這篇文章主要介紹了如何用Matplotlib 畫三維圖的示例代碼，文中通過示例代碼介紹的非常詳細，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧

三維的線圖和散點圖

#繪制三角螺旋線
from mpl_toolkits import mplot3d
%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np

ax = plt.axes(projection='3d')

#三維線的數(shù)據(jù)
zline = np.linspace(0, 15, 1000)
xline = np.sin(zline)
yline = np.cos(zline)
ax.plot3D(xline, yline, zline, 'gray')

# 三維散點的數(shù)據(jù)
zdata = 15 * np.random.random(100)
xdata = np.sin(zdata) + 0.1 * np.random.randn(100)
ydata = np.cos(zdata) + 0.1 * np.random.randn(100)
ax.scatter3D(xdata, ydata, zdata, c=zdata, cmap='Greens')

三維等高線圖

def f(x, y):
 return np.sin(np.sqrt(x ** 2 + y ** 2))
x = np.linspace(-6,6,30)
y = np.linspace(-6,6,30)
X, Y = np.meshgrid(x, y)
Z = f(X,Y)

fig = plt.figure()
ax = plt.axes(projection='3d')
ax.contour3D(X, Y, Z, 50, cmap='binary')
ax.set_xlabel('x')
ax.set_ylabel('y')
ax.set_zlabel('z')
#調(diào)整觀察角度和方位角。這里將俯仰角設(shè)為60度，把方位角調(diào)整為35度
ax.view_init(60, 35)

線框圖和全面圖

全面圖和線框圖相似，只不過線框圖的每一個面都是由多邊形構(gòu)成。只要增加唉一個配色方案來填充這些多邊形，就可以感受到可視化圖形表面的拓撲結(jié)構(gòu)了。

#線框圖
fig =plt.figure()
ax = plt.axes(projection='3d')
ax.plot_wireframe(X, Y, Z, color='c')
ax.set_title('wireframe')

#曲面圖
ax = plt.axes(projection='3d')
ax.plot_surface(X, Y, Z, rstride=1, cstride=1, cmap='viridis', edgecolor='none')
ax.set_title('surface')

#使用極坐標(biāo)可以獲得切片的效果
r = np.linspace(0, 6, 20)
theta = np.linspace(-0.9 * np.pi, 0.8 * np.pi, 40)
r, theta = np.meshgrid(r, theta)
X = r * np.sin(theta)
Y = r * np.cos(theta)
Z = f(X, Y)
ax = plt.axes(projection='3d')
ax.plot_surface(X, Y, Z, rstride=1, cstride=1, cmap='viridis', edgecolor='none')

曲面三角剖分

在某些應(yīng)用場景下，上述這些要求均勻采樣的網(wǎng)格數(shù)據(jù)顯得太過嚴(yán)格且不太容易實現(xiàn)。這時就可以使用三角剖分部分圖形。

theta = 2 * np.pi * np.random.random(1000)
r = 6 * np.random.random(1000)
x = np.ravel(r * np.sin(theta))
y = np.ravel(r * np.cos(theta))
z = f(x, y)

ax = plt.axes(projection='3d')
ax.scatter(x, y, z, c=z, cmap='viridis', linewidth=0.5)

#上圖還有許多地方需要修補，這些工作可以由ax.plot_trisurf函數(shù)幫助我們完成。它首先找到一組所有點都連接起來的三角形，然后用這些三角形創(chuàng)建曲面
ax = plt.axes(projection='3d')
ax.plot_trisurf(x, y, z, cmap='viridis', edgecolor='none')

莫比烏斯帶（應(yīng)用曲面三角剖分）

#繪制莫比烏斯帶
#由于它是一條二維帶，因此需要兩個內(nèi)在維度。theta維度取值范圍是0～2pi，寬度維度w取值范圍是-1～1
theta = np.linspace(0, 2 * np.pi, 30)
w = np.linspace(-0.25, 0.25, 8)
w, theta = np.meshgrid(w, theta)
phi = 0.5 * theta
#x-y平面內(nèi)的半徑
r = 1 + w * np.cos(phi)

x = np.ravel(r * np.cos(theta))
y = np.ravel(r * np.sin(theta))
z = np.ravel(w * np.sin(phi))

#要畫出莫比烏斯帶，還必須保證三角部分是正確的。最好的方法是首先用基本參數(shù)化方法定義三角部分，然后用Matplotlib將
#這個三角剖分映射到莫比烏斯帶的三維空間里
from matplotlib.tri import Triangulation
tri = Triangulation(np.ravel(w), np.ravel(theta))
ax = plt.axes(projection='3d')
ax.plot_trisurf(x, y, z, triangles=tri.triangles, cmap='viridis', linewidth=0.2)
ax.set_xlim(-1, 1);ax.set_ylim(-1,1);ax.set_zlim(-1,1)