求助：如何用Paraview绘制64×32矩形网格中的列数据点？

我完全懂你这种 frustration——ParaView的逻辑和gnuplot那种直接映射列的方式确实不一样，尤其是处理这种有隐含网格结构的数据时。下面是我试过的可行方法，一步步来：

1. 正确导入数据文件 #

首先要确保ParaView正确识别你的列数据和空行分隔：

• 打开ParaView，通过File > Open选择你的dat.25文件。
• 在导入对话框中，设置Field Delimiters为空格+换行，并勾选Treat consecutive delimiters as one（这样空行会被正确识别为分隔符，不会导致数据错位）。
• 确认每一列被识别为单独的点数据数组（比如默认命名为Field 0、Field 1，你可以后续重命名为Col1、Col2等方便识别）。
• 点击Apply完成导入。

2. 生成网格的X/Y坐标 #

你的数据是64（X方向）×32（Y方向）的网格，每个点的X/Y坐标是隐含的，需要手动生成：

• 选中导入的数据源，点击Filters > Alphabetical > Programmable Filter。
• 在弹出的编辑器中，输入以下Python代码来生成X/Y坐标数组：

  import numpy as np
  
  # 定义网格尺寸
  x_count = 64
  y_count = 32
  total_points = x_count * y_count
  
  # 生成X坐标：0~63循环32次
  x_coords = np.tile(np.arange(x_count), y_count)
  # 生成Y坐标：0~31每个值重复64次
  y_coords = np.repeat(np.arange(y_count), x_count)
  
  # 将X/Y添加为点数据
  output.PointData.append(x_coords, "X")
  output.PointData.append(y_coords, "Y")
  

• 点击Apply，现在你的数据源里就有了每个点对应的X、Y坐标，以及原文件的各列Z值。

3. 绘制网格线（模拟gnuplot的wi li效果） #

这里有两种方式实现，按需选择：

方法一：快速生成网格线（适合大部分场景） #

这种方式用内置过滤器组合快速实现：

• 选中Programmable Filter的输出，点击Filters > Alphabetical > Table To Points。
  • 在参数面板中，设置X Column为X，Y Column为Y，Z Column选择你要绘制的列（比如Col1）。
  • 点击Apply将表格数据转为三维点云。
• 接着点击Filters > Alphabetical > Delaunay2D，这个过滤器会把点云自动连成2D网格面。点击Apply。
• 最后点击Filters > Alphabetical > Extract Edges，提取网格的边线。点击Apply后，你就能看到和gnuplot类似的网格线效果了。

方法二：精准控制线条（自定义X/Y方向线） #

如果需要更精准地控制只绘制X或Y方向的线条，可以用另一个Programmable Filter直接生成线单元：

• 选中带X/Y坐标的数据源，再次添加Programmable Filter。
• 输入以下代码（记得把Col1替换成你要绘制的目标列）：

  import vtk
  import numpy as np
  
  input_data = inputs[0]
  x_count = 64
  y_count = 32
  target_col = "Col1"  # 替换为你要绘制的列名
  
  # 初始化输出的PolyData
  output = self.GetOutput()
  points = vtk.vtkPoints()
  lines = vtk.vtkCellArray()
  
  # 先添加所有点到输出
  for idx in range(input_data.GetNumberOfPoints()):
      x = input_data.PointData["X"][idx]
      y = input_data.PointData["Y"][idx]
      z = input_data.PointData[target_col][idx]
      points.InsertNextPoint(x, y, z)
  
  # 生成X方向的线条（每行内的点依次连接）
  for y_idx in range(y_count):
      for x_idx in range(x_count - 1):
          line = vtk.vtkLine()
          line.GetPointIds().SetId(0, y_idx * x_count + x_idx)
          line.GetPointIds().SetId(1, y_idx * x_count + x_idx + 1)
          lines.InsertNextCell(line)
  
  # 生成Y方向的线条（每列内的点依次连接）
  for x_idx in range(x_count):
      for y_idx in range(y_count - 1):
          line = vtk.vtkLine()
          line.GetPointIds().SetId(0, y_idx * x_count + x_idx)
          line.GetPointIds().SetId(1, (y_idx + 1) * x_count + x_idx)
          lines.InsertNextCell(line)
  
  # 赋值给输出
  output.SetPoints(points)
  output.SetLines(lines)
  

• 点击Apply后，视图中会直接显示完整的网格线，效果和gnuplot的splot命令完全一致。

4. 多列绘制（可选） #

如果需要同时绘制多列数据，比如Col1和Col2，只需要复制上述流程：

• 分别为每列创建独立的Table To Points + Delaunay2D + Extract Edges（或自定义Programmable Filter）流水线。
• 在视图中选中不同的流水线输出，通过Properties面板调整颜色、线型等属性来区分不同列的网格。

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