如何读取VTK文件中的Field Data？已掌握PolyData读取方法

下面针对你提供的VTK文件，给出两种常用的读取Field Data的实现方式：

一、使用VTK原生Python库 #

如果你习惯用VTK原生API，读取PolyData后可以通过以下步骤提取Cell Data中的Field Data：

import vtk

# 读取PolyData文件
reader = vtk.vtkPolyDataReader()
reader.SetFileName("sampleSurface.vtk")  # 替换为你的文件名
reader.Update()
polydata = reader.GetOutput()

# 获取单元数据集合
cell_data = polydata.GetCellData()
# 获取Field Data对象
field_data = cell_data.GetFieldData()

# 遍历所有Field Data数组
for idx in range(field_data.GetNumberOfArrays()):
    array = field_data.GetArray(idx)
    print(f"数组名称: {array.GetName()}")
    print(f"分量数: {array.GetNumberOfComponents()}")
    print(f"元素总数: {array.GetNumberOfTuples()}")
    
    # 读取单个元素示例（以第0个单元为例）
    if array.GetNumberOfComponents() == 1:
        print(f"第0个单元的值: {array.GetTuple1(0)}\n")

# 直接通过名称获取指定数组（比如你的文件中的"p"数组）
p_array = field_data.GetArray("p")
if p_array:
    # 获取第10个单元的p值
    target_value = p_array.GetTuple1(10)
    print(f"第10个单元的p值: {target_value}")

二、使用PyVista简化操作 #

PyVista封装了VTK的复杂接口，代码更简洁：

import pyvista as pv

# 读取PolyData
mesh = pv.read("sampleSurface.vtk")

# 直接通过cell_data字典获取Field Data中的"p"数组
p_values = mesh.cell_data["p"]
print(f"p数组长度: {len(p_values)}")
print(f"第0个单元的p值: {p_values[0]}")
print(f"第10个单元的p值: {p_values[10]}")

对你提供的VTK文件结构说明 #

你的文件中CELL_DATA 3600表示后续数据对应3600个单元；FIELD attributes 2说明包含2个字段（你提供的片段中只显示了p字段）；p 1 3600 float表示该字段名为p，是单分量数组，包含3600个浮点型元素，对应每个单元的数值。

内容的提问来源于stack exchange，提问作者Muhammed Ibraaheem Noufal