我刚开始在 python 中使用 mayavi,对 3d 图形编程了解不多。我想绘制一个凹面/物体,但当我尝试时,结果是这样连接的:
你看图片看不出来,但它是由平行平面组成的(也不是我想要的)
当我尝试绘制具有这种形状的结构时:
所以这是我为获得要点和情节所做的工作:
#extract the black in the image
x,y,z = np.where(quest == 0)
#make z1,z2 the right size(y*0) and value(+12)
z1 = y*0
z2 = y*0+12
x,y,z1,z2= np.atleast_2d(x,y,z1,z2)
#duplicate the shape at z1 and z2
x = np.concatenate((x,x))
y = np.concatenate((y,y))
z = np.concatenate((z1,z2))
#plot the mesh
mlab.mesh(x,y,z)
mlab.show()
我注意到数据点的顺序对于使用网格很重要,因为它会从每个连续的点画一条线。
有没有更简单的方法来组织我的观点?
还有没有办法让形状填充点之间的区域?
感谢所有帮助,谢谢
后续问题:
是否有更好的 python 3d 库可以满足此目的?
我正在考虑用 vtk 做这件事,但是同样的问题是像这样切掉问号的凹面部分:
更长的代码是:
#create the points
points = vtk.vtkPoints()
for i, j in zip(x,y):
points.InsertNextPoint(i,j,0)
#create a polygon
polygon = vtk.vtkPolygon()
polygon.GetPointIds().SetNumberOfIds(len(x))
for i in range(len(x)):
polygon.GetPointIds().SetId(i,i)
#add it to a cell array
polygons = vtk.vtkCellArray()
polygons.InsertNextCell(polygon)
#make the poly data
polyData = vtk.vtkPolyData()
polyData.SetPoints(points)
polyData.SetPolys(polygons)
#set up the window
mapper = vtk.vtkPolyDataMapper()
mapper.SetInput(polyData
actor = vtk.vtkActor()
actor.SetMapper(mapper)
renderer = vtk.vtkRenderer()
rWindow = vtk.vtkRenderWindow()
rWindow.AddRenderer(renderer)
rWindowInteractor = vtk.vtkRenderWindowInteractor()
rWindowInteractor.SetRenderWindow(rWindow)
renderer.AddActor(actor)
renderer.SetBackground(.5,.3,.31)
rWindow.Render()
rWindowInteractor.Start()
最佳答案
这篇文章已有一年多了,但它可能对其他人仍然有用。
我遇到了同样的问题,我使用 vtkTriangleFilter
解决了它.此过滤器将输入多边形转换为三角形,更易于 OpenGL 理解。
这是一个快速脚本,我在其中创建了一个凹形(“看起来像”一头大象),并提供了一个 bool 值来在开始时过滤或不过滤该形状。
#!/usr/bin/env python
from vtk import *
# switch to False if you don't want to filter the shape
FILTER = True
# x and y values of points to create a concave shape
x = [-2.25032, -3.04773, -3.06125, -2.69633, -2.52063, -2.16923, -1.72322, -1.47994, -1.85837, -1.04745, -0.68253, 1.25018, 1.81783, 3.20992, 2.15571, 1.99353, 2.72336, 3.35859, 3.57483, 3.60186, 3.69647, 4.18303, 4.29115, 3.87217, 2.87203, 2.07462, 1.43939, 1.10151, 0.912292, 0.466283, 1.27721, 0.723076, -1.22315]
y = [-1.3371, -0.671689, 0.094162, 0.081607, -0.4457, -0.960452, -0.797238, 0.0313873, 1.43754, 1.76397, -0.00627746, -0.295041, 1.42498, 1.36221, -0.39548, -1.21155, -1.57564, -1.57564, -1.03578, -0.47081, 0.244821, 0.169492, -0.822348, -1.98996, -2.7307, -2.96924, -3.20778, -2.89391, -2.56748, -2.12806, -1.93974, -1.65097, -1.60075]
# creation of the source
positions = zip(x, y)
triFilter = vtkTriangleFilter()
pts = vtkPoints()
polygon = vtkPolygon()
polygon.GetPointIds().SetNumberOfIds(len(positions))
for k, l in enumerate(positions):
pts.InsertNextPoint(float(l[0]), float(l[1]), 1)
polygon.GetPointIds().SetId(k, k)
polygons = vtkCellArray()
polygons.InsertNextCell(polygon)
polygonPolyData = vtkPolyData()
polygonPolyData.SetPoints(pts)
polygonPolyData.SetPolys(polygons)
# let's filter the polydata
triFilter.SetInputData( polygonPolyData)
triFilter.Update()
polygonPolyDataFiltered = triFilter.GetOutput()
# map it
mapper = vtkPolyDataMapper()
if VTK_MAJOR_VERSION <= 5:
if FILTER:
mapper.SetInput(polygonPolyDataFiltered)
else:
mapper.SetInput(polygonPolyData)
else:
if FILTER:
mapper.SetInputData(polygonPolyDataFiltered)
else:
mapper.SetInputData(polygonPolyData)
# render it
actor = vtkActor()
actor.SetMapper(mapper)
ren = vtkRenderer()
ren.AddActor(actor)
renwin = vtkRenderWindow()
renwin.AddRenderer(ren)
iren = vtkRenderWindowInteractor()
iren.SetRenderWindow(renwin)
iren.SetInteractorStyle(vtkInteractorStyleTrackballCamera())
iren.Initialize()
iren.Start()
关于python - 用 mayavi 或 vtk 绘制凹形物体,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/15103489/