如何将vtk遗留文件转换为vtu格式

时间:2014-10-02 10:06:08

标签: vtk paraview

我想将旧的.vtk文件转换为二进制文件,最好是.vtu文件,因为我使用的是非结构化网格。 为此,我改编了http://www.vtk.org/Wiki/VTK/Examples/Cxx/IO/ConvertFile

中的ConvertFile-Example
#include <string>
#include <vtkSmartPointer.h>
#include <vtkGenericDataObjectReader.h>
#include <vtkVersion.h>
#include <vtkXMLUnstructuredGridWriter.h>
#include <vtkUnstructuredGrid.h>

int main(int argc, char *argv[])
{
  if(argc < 3)
    {
    std::cerr << "Required arguments: input.vtk output.vtu" << std::endl;
    return EXIT_FAILURE;
    }

  std::string inputFileName = argv[1];
  std::string outputFileName = argv[2];

  vtkSmartPointer<vtkGenericDataObjectReader> reader = vtkSmartPointer<vtkGenericDataObjectReader>::New();
  reader->SetFileName(inputFileName.c_str());
  reader->Update();

  vtkSmartPointer<vtkXMLUnstructuredGridWriter> writer = vtkSmartPointer<vtkXMLUnstructuredGridWriter>::New();
  writer->SetFileName(outputFileName.c_str());
  writer->SetInputConnection(reader->GetOutputPort());
  writer->Update();

  return EXIT_SUCCESS;
}

但是当我使用它来转换我的遗留文件时,我在第一集之后丢失了所有Cell Data。在我的遗留文件的最小示例中,Scal_1位于.vtu文件中,但Scal_2不存在。

# vtk DataFile Version 3.1
Lattice Boltzmann data
ASCII
DATASET UNSTRUCTURED_GRID
POINTS 9 INT 
0 0 0 1 0 0 2 0 0 
0 1 0 1 1 0 2 1 0 
0 2 0 1 2 0 2 2 0 

CELLS 4 20
4 0 1 3 4 
4 1 2 4 5 
4 3 4 6 7 
4 4 5 7 8 

CELL_TYPES 4
8 8 8 8 

CELL_DATA 4
SCALARS Scal_1 DOUBLE
LOOKUP_TABLE default
1 2 1 0 

SCALARS Scal_2 DOUBLE
LOOKUP_TABLE default
1 3 2 1

我还是vtk的新手。我应该使用其他读者或作者吗?或者是完全错误的?

2 个答案:

答案 0 :(得分:3)

这里的问题是,您选择的阅读器会因为包含2个单元格数据的输入文件都标记为标量而感到困惑。因此,阅读器只输出一个单元数据阵列。我的建议是使用ParaView,特别是pvpython可执行文件来转换文件。相应的Python代码如下所示:

from paraview.simple import *
r = LegacyVTKReader( FileNames=['input.vtk'] )
w = XMLUnstructuredGridWriter()
w.FileName = 'output.vtu'
w.UpdatePipeline()

答案 1 :(得分:0)

您可以只使用meshio(属于我的项目)。使用

安装
pip3 install meshio

然后运行

meshio-convert in.vtk out.vtu