可视化由六面体组成的结构

Question

我使用matplotlib编写了一个Python脚本，用于显示和旋转多维数据集以及在多维数据集中显示对称性的平面或轴。从下面的插图中可以看出，我成功了一点。在那里你可以看到七个七个可能的六个立方体的三维立方体之一。我的脚本使用matplotlib，它似乎做得很好：matplotlib允许我旋转多面体并在任何想要的角度查看它。大多数多面体存在一个大问题，如下面的第二张图所示。在我开始旋转多面体之后的一点点，matplotlib显示部分位于其他平面之后的平面，因此部分不可见，不能绘制或仅部分绘制。

我一直在论坛和谷歌上搜索，但无济于事。有点打击表明我应该使用mayavi而不是matplotlib。所以我对mayavi进行了广泛的研究。我花了几周时间试图找出如何让mayavi继续前行。对docs.enthought的点击最初看起来很有希望，但虽然mayavi显然适合目的，并且在可视化对象方面非常出色，但我找不到可理解的文档。我可以在vtk或tvtk上使用真正的程序员指南。有很多文档，但主要是设计师，而不是程序员。如果不可用（？）我也会对带有六面体或不规则网格的脚本示例感到满意，这些脚本在python的canopy（1.6.2）实现中有效。

Answer 1

我已经修改了vtk文档（http://www.vtk.org/gitweb?p=VTK.git;a=blob_plain;f=Examples/DataManipulation/Python/BuildUGrid.py）

中的示例

此示例可以转换为可以从点坐标数组和元素数组创建几何的函数。

我没有使用任何vtk指南，我通常会参考这里找到的python示例：http://www.vtk.org/Wiki/VTK/Examples/Python

import vtk

# Initialize the vtkPoints variable and set the number of points
points = vtk.vtkPoints()
points.SetNumberOfPoints(8)

# Add points to the variable, with the point number first, then the x, y, z coordinates.
# For demonstration purposes, I started numbering the ponts at 10 (normally they would start at 0).
points.InsertPoint(0, 0, 0, 0)
points.InsertPoint(1, 1, 0, 0)
points.InsertPoint(2, 1, 1, 0)
points.InsertPoint(3, 0, 1, 0)
points.InsertPoint(4, 0, 0, 1)
points.InsertPoint(5, 1, 0, 1)
points.InsertPoint(6, 1, 1, 1)
points.InsertPoint(7, 0, 1, 1)

points.InsertPoint(8, 0, 0, 1.1)
points.InsertPoint(9, 1, 0, 1.1)
points.InsertPoint(10, 1, 1, 1.1)
points.InsertPoint(11, 0, 1, 1.1)
points.InsertPoint(12, 0, 0, 2)
points.InsertPoint(13, 1, 0, 2)
points.InsertPoint(14, 1, 1, 2)
points.InsertPoint(15, 0, 1, 2.5)

# Define the hexahedron, then set the point Ids of the hexahedron cell/element.
# From the documentation: points (0,1,2,3) is the base of the hexahedron which, using the right hand rule, forms a
# quadrilaterial whose normal points in the direction of the opposite face (4,5,6,7)
aHexahedron1 = vtk.vtkHexahedron()
aHexahedron1.GetPointIds().SetId(0, 0) # Cell point 0 corresponds to point 0 which was defined above
aHexahedron1.GetPointIds().SetId(1, 1)
aHexahedron1.GetPointIds().SetId(2, 2)
aHexahedron1.GetPointIds().SetId(3, 3)
aHexahedron1.GetPointIds().SetId(4, 4)
aHexahedron1.GetPointIds().SetId(5, 5)
aHexahedron1.GetPointIds().SetId(6, 6)
aHexahedron1.GetPointIds().SetId(7, 7)

# Define a second hexahedron
aHexahedron2 = vtk.vtkHexahedron()
aHexahedron2.GetPointIds().SetId(0, 8) # Cell point 0 corresponds to point 8 which was defined above
aHexahedron2.GetPointIds().SetId(1, 9)
aHexahedron2.GetPointIds().SetId(2, 10)
aHexahedron2.GetPointIds().SetId(3, 11)
aHexahedron2.GetPointIds().SetId(4, 12)
aHexahedron2.GetPointIds().SetId(5, 13)
aHexahedron2.GetPointIds().SetId(6, 14)
aHexahedron2.GetPointIds().SetId(7, 15)

# Define an unstructured grid.
aHexahedronGrid = vtk.vtkUnstructuredGrid()

# Add the hexahedron to the unstructured grid
# Note: this operation defines the point ids, and not the actual point coordinates
aHexahedronGrid.InsertNextCell(aHexahedron1.GetCellType(), aHexahedron1.GetPointIds())
aHexahedronGrid.InsertNextCell(aHexahedron2.GetCellType(), aHexahedron2.GetPointIds())

# Set the points which includes the coordinates. The point ids defined in the line above correspond to the point ids
# that were defined earlier (i.e. points.InsertPoint(10, 0, 0, 0))
aHexahedronGrid.SetPoints(points)

# Now we have defined one hexahedron, and added it an unstructured grid.
# We could create more hexahedrons, and add them to the same unstructured grid.

# To view the unstructured grid, we need to define a mapper and set the unstructured grid as the input
aHexahedronMapper = vtk.vtkDataSetMapper()
aHexahedronMapper.SetInputData(aHexahedronGrid)

# Define an actor, and set the mapper as the input
aHexahedronActor = vtk.vtkActor()
aHexahedronActor.SetMapper(aHexahedronMapper)

# Create the usual rendering stuff.
ren = vtk.vtkRenderer()
renWin = vtk.vtkRenderWindow()
renWin.AddRenderer(ren)
iren = vtk.vtkRenderWindowInteractor()
iren.SetRenderWindow(renWin)
iren.SetInteractorStyle(vtk.vtkInteractorStyleTrackballCamera()) # Change the rotation type from the default to 'trackball'
ren.SetBackground(.1, .2, .4)

# Add the actor to the renderer to actually view the geometry
ren.AddActor(aHexahedronActor)

# Render the scene and start interaction.
iren.Initialize()
renWin.Render()
iren.Start()