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

以下是一个示例代码，它生成一个简单的三角形网格并存储不同时刻的属性数据：

import vtk

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

triangles = vtk.vtkCellArray()
triangles.InsertNextCell(3)
triangles.InsertCellPoint(0)
triangles.InsertCellPoint(1)
triangles.InsertCellPoint(2)

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

# 创建属性数据
temperature = vtk.vtkDoubleArray()
temperature.SetName('Temperature')
temperature.InsertNextValue(10.0)
temperature.InsertNextValue(20.0)
temperature.InsertNextValue(30.0)

pressure = vtk.vtkDoubleArray()
pressure.SetName('Pressure')
pressure.InsertNextValue(1000.0)
pressure.InsertNextValue(2000.0)
pressure.InsertNextValue(3000.0)

# 将属性数据附加到网格
polydata.GetPointData().AddArray(temperature)
polydata.GetPointData().AddArray(pressure)

# 存储不同时刻的属性数据
writer = vtk.vtkXMLPolyDataWriter()
writer.SetFileName('mesh.vtp')
writer.SetInputData(polydata)

for i in range(3):
    # 修改属性数据
    temperature.SetValue(i, temperature.GetValue(i) + i * 10)
    pressure.SetValue(i, pressure.GetValue(i) + i * 1000)
    
    # 存储网格和属性数据
    writer.UpdateTimeStep(i)
    writer.Write()

这段代码创建了一个简单的三角形网格，并在每个时间步长上修改温度和压力属性数据。然后，它使用VTK的XML格式存储三角形网格和属性数据。在保存每个时间步长时，它使用UpdateTimeStep()方法来指定当前时间步长，以便VTK可以将不同时刻的数据存储在同一个文件中。