标签:python matplotlib matlab mplot3d
这是3D绘图代码的MATLAB版本:
编辑:
这是当前的代码:
plt.figure(2)
fig_b = Axes3D(fig2)
xx2 = np.arange(0, L+h_grid*L, h_grid*L)
yy2 = np.arange(-b, b+h_grid*b, h_grid*b)
X, Y = np.meshgrid(xx2, yy2)
W = np.zeros((41,21), float)
mx = len(xx2)*len(yy2)
X = np.reshape(X, (1, mx))
Y = np.reshape(Y, (1, mx))
W = np.reshape(W, (1, mx))
for j in range(0, mx):
W[0][j] = np.sin(np.pi*X[0][j]/L)
surf = fig_b.plot_surface(X, Y, W, rstride=1, cstride=1, cmap=cm.jet, linewidth=0, antialiased=False) # This is the line number 168
plt.xlabel('x')
plt.ylabel('y')
这是我收到的错误消息:
Traceback (most recent call last):
File "nonhomog.py", line 247, in <module>
main()
File "nonhomog.py", line 245, in main
nonhomog(nu)
File "nonhomog.py", line 168, in nonhomog
surf = fig_b.plot_surface(X, Y, W, rstride=1, cstride=1, cmap=cm.jet, linewidth=0, antialiased=False)
File "/usr/lib/pymodules/python2.6/mpl_toolkits/mplot3d/axes3d.py", line 618, in plot_surface
polyc = art3d.Poly3DCollection(polys, *args, **kwargs)
File "/usr/lib/pymodules/python2.6/mpl_toolkits/mplot3d/art3d.py", line 290, in __init__
PolyCollection.__init__(self, verts, *args, **kwargs)
File "/usr/lib/pymodules/python2.6/matplotlib/collections.py", line 668, in __init__
self.set_verts(verts, closed)
File "/usr/lib/pymodules/python2.6/mpl_toolkits/mplot3d/art3d.py", line 312, in set_verts
self.get_vector(verts)
File "/usr/lib/pymodules/python2.6/mpl_toolkits/mplot3d/art3d.py", line 305, in get_vector
xs, ys, zs = zip(*points)
ValueError: need more than 0 values to unpack
解决方法:
在为X和Y设置网格网格后,您需要为Z值设置网格.
我目前在我的代码中执行此操作的方式是:
# [ (x1, y1, z1), (x2, y2, z2), ... (xN, yN, zN) ]
all_vals = ...
# (x1, x2, ... xN) , (y1, y2, ... yN) , (z1, z2, ... zN)
all_xvals, all_yvals, all_zvals = zip(*all_vals)
fig = plt.figure()
ax = Axes3D(fig)
X, Y = np.meshgrid(xvals, yvals)
# This is the part you want:
Z1 = np.zeros(X.shape, float)
for (x, y, z) in all_vals:
x = find_in_sorted_list(x, xvals)
y = find_in_sorted_list(y, yvals)
Z1[y,x] = z
surf = ax.plot_surface(X, Y, Z1, rstride=1, cstride=1, cmap=cm.jet,
linewidth=0, antialiased=False)
plt.xlabel('Blur standard deviation')
plt.ylabel('JPEG quality')
ax.w_zaxis.set_major_locator(LinearLocator(10))
ax.w_zaxis.set_major_formatter(FormatStrFormatter('%.03f'))
fig.colorbar(surf, shrink=0.5, aspect=5)
plt.show()
这给了我一个如下图:
我已将其保存为文件,但是当您调用plt.show()时,您将获得一个交互式窗口,您可以在其中将视点更改为您想要的任何内容.
标签:python,matplotlib,matlab,mplot3d 来源: https://codeday.me/bug/20190610/1211801.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。