一：主要的知识点
1、说明
本文只是教程内容的一小段，因博客字数限制，故进行拆分。主教程链接：vtk教程——逐行解析官网所有Python示例-CSDN博客
2、知识点纪要
本段代码主要涉及的有①vtkRectilinearGridGeometryFilter从非均匀三维网格提取表面集合进行展示
二：代码及注释
import vtkmodules.vtkRenderingCore import vtkmodules.vtkInteractionStyle from vtkmodules.vtkCommonDataModel import vtkRectilinearGrid from vtkmodules.vtkCommonColor import vtkNamedColors from vtkmodules.vtkCommonCore import vtkDoubleArray from vtkmodules.vtkFiltersGeometry import vtkRectilinearGridGeometryFilter from vtkmodules.vtkRenderingCore import ( vtkActor, vtkPolyDataMapper, vtkRenderWindow, vtkRenderWindowInteractor, vtkRenderer ) def main(): colors = vtkNamedColors() x = [-1.22396, -1.17188, -1.11979, -1.06771, -1.01562, -0.963542, -0.911458, -0.859375, -0.807292, -0.755208, -0.703125, -0.651042, -0.598958, -0.546875, -0.494792, -0.442708, -0.390625, -0.338542, -0.286458, -0.234375, -0.182292, -0.130209, -0.078125, -0.026042, 0.0260415, 0.078125, 0.130208, 0.182291, 0.234375, 0.286458, 0.338542, 0.390625, 0.442708, 0.494792, 0.546875, 0.598958, 0.651042, 0.703125, 0.755208, 0.807292, 0.859375, 0.911458, 0.963542, 1.01562, 1.06771, 1.11979, 1.17188] y = [-1.25, -1.17188, -1.09375, -1.01562, -0.9375, -0.859375, -0.78125, -0.703125, -0.625, -0.546875, -0.46875, -0.390625, -0.3125, -0.234375, -0.15625, -0.078125, 0, 0.078125, 0.15625, 0.234375, 0.3125, 0.390625, 0.46875, 0.546875, 0.625, 0.703125, 0.78125, 0.859375, 0.9375, 1.01562, 1.09375, 1.17188, 1.25] z = [0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.75, 0.8, 0.9, 1, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.75, 1.8, 1.9, 2, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.75, 2.8, 2.9, 3, 3.1, 3.2, 3.3, 3.4, 3.5, 3.6, 3.7, 3.75, 3.8, 3.9] print(len(x), len(y), len(z)) xCoords = vtkDoubleArray() for value in x: xCoords.InsertNextValue(value) yCoords = vtkDoubleArray() for i in range(0, len(y)): yCoords.InsertNextValue(y[i]) zCoords = vtkDoubleArray() for i in range(0, len(z)): zCoords.InsertNextValue(z[i]) # 构建vtkRectilinearGrid数据模型 rgrid = vtkRectilinearGrid() rgrid.SetDimensions(len(x), len(y), len(z)) rgrid.SetXCoordinates(xCoords) rgrid.SetYCoordinates(yCoords) rgrid.SetZCoordinates(zCoords) """ vtkRectilinearGridGeometryFilter 主要功能：从 vtkRectilinearGrid（非均匀三维网格）中提取出 表面几何（vtkPolyData） 简单说，它可以让一个三维体数据“露出”一层表皮（例如某个切片面），用于渲染或可视化 也可以整个三维体数据，但是会看不见网格结构 """ plane = vtkRectilinearGridGeometryFilter() plane.SetInputData(rgrid) plane.SetExtent(0, len(x) - 1, 16, 16, 0, len(z) - 1) rgridMapper = vtkPolyDataMapper() rgridMapper.SetInputConnection(plane.GetOutputPort()) wireActor = vtkActor() wireActor.SetMapper(rgridMapper) wireActor.GetProperty().SetColor(colors.GetColor3d("Banana")) wireActor.GetProperty().EdgeVisibilityOn() # Create the usual rendering stuff. renderer = vtkRenderer() renWin = vtkRenderWindow() renWin.AddRenderer(renderer) iren = vtkRenderWindowInteractor() iren.SetRenderWindow(renWin) renderer.AddActor(wireActor) renderer.SetBackground(colors.GetColor3d("Beige")) renderer.ResetCamera() renderer.GetActiveCamera().Elevation(60.0) renderer.GetActiveCamera().Azimuth(30.0) renderer.GetActiveCamera().Zoom(1.0) renWin.SetSize(640, 480) renWin.SetWindowName('RGrid') # Interact with the data. renWin.Render() iren.Start() if __name__ == '__main__': main()
