matplotlib - 从一个矩形阵列的高度的3d表面

时间:2012-08-01 19:51:42

标签: python matplotlib 3d hdf5 geometry-surface

我试图在matplotlib中绘制一些HDF数据。使用h5py导入它们后,数据以数组的形式存储,如下所示:

array([[151, 176, 178],
       [121, 137, 130],
       [120, 125, 126])

在这种情况下,x和y值只是数组字段的索引,而z值是特定字段的值。在(x,y,z)形式中,它看起来像:

(1,1,151)
(2,1,176)
(3,1,178)
(1,2,121)
...

等等。

有没有一种简单的方法可以从这种数据中获得表面图?我知道我可以通过遍历整个数组将其更改为(x,y,z)元组,但也许不需要它?

1 个答案:

答案 0 :(得分:40)

如果您需要三维曲面图,则必须先创建meshgrid。你可以尝试:

from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
import numpy as np

X = np.arange(1, 10)
Y = np.arange(1, 10)
X, Y = np.meshgrid(X, Y)
R = np.sqrt(X**2 + Y**2)
Z = np.sin(R)
fig = plt.figure()
ax = fig.gca(projection='3d')
surf = ax.plot_surface(X, Y, Z, rstride=1, cstride=1, cmap='hot', linewidth=0, antialiased=False)
ax.set_zlim(-1.01, 1.01)

fig.colorbar(surf, shrink=0.5, aspect=5)
plt.show()

会生成, enter image description here

但是,如果唯一相关信息在z值中,则只需使用imshow即可。这里,z值由它们的颜色表示。你可以通过以下方式实现这一目标:

im = plt.imshow(Z, cmap='hot')
plt.colorbar(im, orientation='horizontal')
plt.show()

哪个会给,

enter image description here