vtkStreamTracer处理非结构网格无输出数据问题求助

• 未配置流线种子源
  vtkStreamTracer需要显式指定流线积分的起始点（种子点），无种子点的情况下不会生成任何流线。你可以根据需求选择不同的种子源生成器，示例配置如下：

# 示例：在网格包围盒中心生成100个随机种子点
seed_source = vtk.vtkPointSource()
seed_source.SetCenter(reader.GetOutput().GetCenter())
seed_source.SetRadius(reader.GetOutput().GetLength()/4) # 种子范围设为网格尺寸的1/4
seed_source.SetNumberOfPoints(100)
seed_source.Update()

# 将种子源关联到流线生成器
streamer.SetSourceConnection(seed_source.GetOutputPort())

如果需要自定义种子位置，也可以用vtkPolyData手动构造种子点集。

• 未指定速度场数组
  vtkStreamTracer默认不会自动选择速度场，需要你显式声明用于积分的矢量数组位置和名称：

# 第一个参数0是固定标识，第二个参数是端口号，第三、四、五个参数分别指定：
# 关联场的类型（点数据/单元数据）、数组关联类型、数组名称
# 示例：使用点数据中名为"velocity"的矢量数组作为速度场
streamer.SetInputArrayToProcess(0, 0, 0, vtk.vtkDataObject.FIELD_ASSOCIATION_POINTS, "velocity")

你需要将上述代码中的"velocity"替换为你的VTK文件中实际存储的速度矢量数组名称，可以先通过以下代码打印所有点数据数组名确认：

point_data = reader.GetOutput().GetPointData()
for i in range(point_data.GetNumberOfArrays()):
    print(point_data.GetArrayName(i))

• 参数不匹配问题
  你当前设置的SetIntegrationDirection(2)对应双向积分，本身不会导致输出为空，但需要确认单位参数是否匹配：SetIntegrationStepUnit(2)指定步长单位为单元长度倍数，如果你设置的MaximumPropagation(500)远小于单步步长、或种子点落在网格域外，也会导致流线生成失败，可以先将种子点设置在网格中心区域排查该问题。

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