Python VTK 代码：存储三角形网格不同时间步的属性数据

以下是使用 VTK 库在 Python 中存储三角形网格不同时刻属性数据的示例代码：

import vtk

# 创建一个三角形网格
points = vtk.vtkPoints()
points.InsertNextPoint(0, 0, 0)
points.InsertNextPoint(1, 0, 0)
points.InsertNextPoint(0, 1, 0)

triangles = vtk.vtkCellArray()
triangle = vtk.vtkTriangle()
triangle.GetPointIds().SetId(0, 0)
triangle.GetPointIds().SetId(1, 1)
triangle.GetPointIds().SetId(2, 2)
triangles.InsertNextCell(triangle)

polydata = vtk.vtkPolyData()
polydata.SetPoints(points)
polydata.SetPolys(triangles)

# 添加属性数据
temperature = vtk.vtkFloatArray()
temperature.SetName('Temperature')
temperature.InsertNextValue(10.0)
temperature.InsertNextValue(20.0)
temperature.InsertNextValue(30.0)
polydata.GetPointData().AddArray(temperature)

# 存储多个时间步的数据
for i in range(3):
    # 修改属性数据
    temperature.SetValue(0, 10.0 + i)
    temperature.SetValue(1, 20.0 + i)
    temperature.SetValue(2, 30.0 + i)

    # 创建一个writer并写入文件
    writer = vtk.vtkXMLPolyDataWriter()
    writer.SetFileName('output_%d.vtp' % i)
    writer.SetInputData(polydata)
    writer.Write()

这个示例代码创建了一个三角形网格，并添加了名为'Temperature'的属性数据。然后，它遍历三个时间步，每个时间步都修改属性数据并将网格写入一个 VTK 格式的文件，文件名中包含了时间步的信息。这样就可以存储三角形网格不同时刻的属性数据了。