[英]Finding x and y values corresponding to given z value in 3D numpy object
我在Python( np
)中有一个numpy
对象,我正在用pyplot
绘制图形。 它是一个三维图形,其中x
, y
和z
是包含要绘制图形的值的列表。 下面是我的代码(从这里 ):
data = np.c_[x,y,z]
# regular grid covering the domain of the data
mn = np.min(data, axis=0) - 0.009
mx = np.max(data, axis=0) + 0.009
X,Y = np.meshgrid(np.linspace(mn[0], mx[0]), np.linspace(mn[1], mx[1]))
XX = X.flatten()
YY = Y.flatten()
A = np.c_[np.ones(data.shape[0]), data[:, :2], np.prod(data[:, :2], axis=1), data[:, :2] ** 2]
C, _, _, _ = scipy.linalg.lstsq(A, data[:, 2])
# evaluate it on a grid
Z = np.dot(np.c_[np.ones(XX.shape), XX, YY, XX * YY, XX ** 2, YY ** 2], C).reshape(X.shape)
fig = plt.figure()
ax = fig.gca(projection='3d')
ax.plot_surface(X, Y, Z, rstride=1, cstride=1, alpha=0.2)
ax.scatter(data[:, 0], data[:, 1], data[:, 2], c='r', s=50)
plt.xlabel('X')
plt.ylabel('Y')
ax.set_zlabel('Z')
ax.axis('equal')
ax.axis('tight')
plt.show()
所以在这里,我想得到x
和y
的值,其中z
的值最接近-30
。 如何获得x
和y
值。
您可以尝试以下(另请参阅在numpy数组中查找最近的值 )
value = 30.
zid = (np.abs(Z-value)).argmin()
xval = np.reshape(X, -1)[zid]
yval = np.reshape(Y, -1)[zid]
print xval, yval, np.reshape(Z, -1)[zid]
我在问题中提到的示例( https://gist.github.com/amroamroamro/1db8d69b4b65e8bc66a6 )上检查了它,它似乎有效。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.