我想使用Mayavi以与本页第三个图形完全相同的方式绘制三维等高线图(氢电子云模型):
http://www.sethanil.com/python-for-reseach/5
我有一组数据点,我使用自己的模型创建,我想使用它.数据点存储在一个多维的numpy数组中,如下所示:
XYZV = [[1, 2, 3, 4],
[6, 7, 8, 9],
...
[4, 5, 6, 7]]
Run Code Online (Sandbox Code Playgroud)
数据点不是在XYZ空间中均匀分布,也不是以任何特定顺序存储.我认为这个例子使用了一个meshgrid来生成数据点 - 我已经看了这个但完全不理解它.任何帮助将非常感激?
H http://www.sethanil.com/_/rsrc/1267943775903/python-for-reseach/5/Article5-fig3.png
我正在尝试matplotlib.tri.Triangulation用于为 matplotlibs 生成三角形plot_trisurf。我想指定三角形而不是让 Delaunay 三角剖分matplotlib.tri.Triangulation使用,因为它不适用于某些情况,例如 xz 或 yz 平面中的三角形。我不确定自己指定三角形是否会解决问题,但我似乎是一件值得尝试的事情。
问题是三角剖分需要一个 (n,3) 数组,其中 n 是三角形的数量。引用 matplotlib.org 上的页面“对于每个三角形,构成三角形的三个点的索引,以逆时针方式排序。” https://matplotlib.org/api/tri_api.html#matplotlib.tri.Triangulation 。我无法辨别如何以正确的形式创建数组,那就是我需要帮助。我很感激任何帮助。
到目前为止,我已经尝试了一些东西,但这是我最后一次尝试的样子:
import matplotlib.pyplot as plt
import numpy as np
import matplotlib.tri as mtri
fig = plt.figure()
ax = fig.gca(projection='3d')
x1=0
x2=1
x3=1
x4=0
y1=0
y2=0
y3=2
y4=2
x=[]
y=[]
x.append(x1)
x.append(x2)
x.append(x3)
x.append(x4)
y.append(y1)
y.append(y2)
y.append(y3)
y.append(y4)
z=np.zeros(8)
triang = mtri.Triangulation(x, y, triangles=[[[x1,y1],[x2,y2],[x3,y3]],[[x3,y3],[x4,y4],[x2,y2]]])
ax.plot_trisurf(triang, z, linewidth=0.2, antialiased=True)
ax.view_init(45,-90)
ax.set_xlabel("X")
ax.set_ylabel("Y")
ax.set_zlabel("Z")
ax.set_aspect("equal")
fig.set_size_inches(8,8)
plt.show()
Run Code Online (Sandbox Code Playgroud)