如何在VTK中制作简单的2D等高线图?

时间:2011-07-27 08:52:44

标签: c++ visual-c++ data-visualization vtk

我正在填写VTK示例,以了解如何制作我创建的数据的简单等高线图。假设我在71x71域中解析了一个结构化矩形网格中的简单PDE,其中节点之间的距离相等。

我已经检查了filledContours的示例,它告诉我们如何可视化来自* .vtp文件的数据。只要我能写出.vtp文件,我就可以了。问题是我不知道如何创建适合轮廓绘图的.vtp数据文件。我设法从嵌入式演员(如球体或圆锥体)创建.vtp文件,但我不知道如何使用我的2D数据阵列。

2 个答案:

答案 0 :(得分:1)

由于你需要一个等距网格,我会使用基于vtk xml的ImageData格式(* .vti)。

来自http://www.cacr.caltech.edu/~slombey/asci/vtk/vtk_formats.simple.html

“ ImageData - 每个ImageData片段在数据集的整个范围内指定其范围。点和单元格...... 由范围,起源和间距隐含地描述。请注意,原点和间距在所有部分中都是常量,因此它们被指定为ImageData XML元素的属性,如下所示。

<VTKFile type=" ImageData" ...> 
    <ImageData WholeExtent=" x1 x2 y1 y2 z1 z2" 
        Origin=" x0 y0 z0" Spacing=" dx dy dz"> 
        <Piece Extent=" x1 x2 y1 y2 z1 z2"> 
            <PointData>...</ PointData> 
            <CellData>...</ CellData> 
        </ Piece> 
    </ ImageData> 
</ VTKFile> 

- 链接信息的结尾“ 请注意,只有x0 y0 z0和dx dy dz是真实的,WholeExtent和PieceExtent是指像素索引。

2D map with 10x10 pixels visualized on Paraview此示例将显示一个10x10像素的地图,其温度从左下角到右上角。值与每个单元格相关联。您可以将此格式调整为2D数据。文件内容(注意我只使用CellData):

example.vti:

<VTKFile type="ImageData" version="0.1" byte_order="LittleEndian">
    <ImageData WholeExtent=" 0 10 0 10 0 1" Origin=" 0 0 0" Spacing=" 1 1 0">

        <Piece Extent=" 0 10 0 10 0 1">

          <CellData Scalars="scalars">
              <DataArray type="Float32" Name="Temperature[C]" format="ascii">
                  1       2       3       4       5       6       7       8       9       10
                  11      12      13      14      15      16      17      18      19      20
                  21      22      23      24      25      26      27      28      29      30
                  31      32      33      34      35      36      37      38      39      40
                  41      42      43      44      45      46      47      48      49      50
                  51      52      53      54      55      56      57      58      59      60
                  61      62      63      64      65      66      67      68      69      70
                  71      72      73      74      75      76      77      78      79      80
                  81      82      83      84      85      86      87      88      89      90
                  91      92      93      94      95      96      97      98      99      100
              </DataArray>
          </CellData>

        </Piece>

    </ImageData>
</VTKFile>

答案 1 :(得分:0)

一种简单的方法是以CSV格式输出数据,每个网格方格都有一个新行,即

....
grid_idx_i, grid_idx_j, grid_idx_k, val_ijk
....

然后,如果您将此加载到paraview中,您可以应用过滤器TableToPoints以在paraview可以处理的形式中获取它。然后应用过滤器Delaunay2D将其从点转换为单元格。完成此操作后,Contour过滤器应该可以正常工作。

如果数据以CSV格式存储为直接网格,但是我不确定是否可能以某种方式加载数据。

编辑:对不起我隐含地假设你试图在paraview中可视化数据。是这种情况吗?