Paraview合并非结构化vtk文件后无可视化显示问题咨询

问题根因 #

• 数据类型强制转换错误：vtkUnstructuredGridReader读取vtk文件输出的是vtkUnstructuredGrid（非结构化网格）类型，代码中强制将其浅拷贝到vtkPolyData（多边形数据）对象时，非结构化网格特有的高阶单元、拓扑连接信息会丢失，最终生成的文件结构不符合规范，导致ParaView无法正常渲染。
• 多余的类型转换操作：vtkAppendFilter本身原生支持vtkUnstructuredGrid作为输入，不需要额外转换为vtkPolyData，这步多余操作是问题核心来源。
• 复用同一个reader实例循环读取文件时，未清空上一次读取的缓存，可能导致数据残留异常。

修复后代码 #

from vtk import vtkAppendFilter, vtkUnstructuredGridWriter, vtkUnstructuredGridReader

append = vtkAppendFilter()
filenames = ['file1.vtk', 'file2.vtk']

for file in filenames:
    # 每次读取新建reader实例，避免缓存残留
    reader = vtkUnstructuredGridReader()
    reader.SetFileName(file)
    reader.Update()
    # 直接传入非结构化网格对象，不需要转PolyData
    append.AddInputData(reader.GetOutput())

append.Update()    

writer = vtkUnstructuredGridWriter()
writer.SetFileName('output.vtk')
writer.SetInputData(append.GetOutput())
writer.Write()

额外注意事项 #

• 若两份输入文件的点数据/单元数据属性名称、数据类型不一致，合并后属性会分别存储，在ParaView中需要分别选择渲染，不会自动合并。
• 如果需要合并后保留原始的分块信息，可以直接在ParaView中导入两个文件后使用Group Datasets过滤器，不需要额外编写代码处理。

内容的提问来源于stack exchange，提问作者rasunag27