感谢在这个网站上获得的知识。我能写一个简单的脚本,打印出等高线图的坐标(x和y)。在
下面是一个例子:from numpy import *
from pylab import *
# generate a set of random points
npts = 500
phi = random(npts)*2*pi
theta = random(npts)*pi
u = random(npts)
x = u * sin( theta) * cos( phi )
y = u * sin( theta) * sin( phi )
z = u * cos( theta )
bins = linspace(-1,1,100)
H = histogram2d(x,y,bins=(bins,bins))
fig=figure()
ax1 = fig.add_subplot(1,1,1)
levels = [1,2,3]
CS = ax1.contour(H[0].T,levels=levels)
show()
# print contours coordinates
for kk in range(len(levels)):
print CS.allsegs[kk][0][:,0],CS.allsegs[kk][0][:,1]
我也想做同样的事,但在3D里,我很挣扎。在
第一步是存储数据,我可以使用:
^{pr2}$
但是现在,我怎样才能得到三维轮廓的坐标呢?我甚至不能想出一个简单的算法。在