構造格子から非構造格子への変換(VTKライブラリ)

  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

 
 
 
 
 
 
 
 
 
 
 
-
-
!
-
|
|
!
|
-
!
|
-
!
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
!
 

#include <vtkSmartPointer.h>
#include <vtkXMLGenericDataObjectReader.h>
#include <vtkXMLUnstructuredGridWriter.h>
#include <vtksys/SystemTools.hxx>
 
#include "vtkCleanUnstructuredGrid.h"
#include "vtkCleanUnstructuredGridCells.h"
 
#include<string>
 
int main(int argc, char* argv[])
{
  // Verify input arguments
  if (argc < 2)
  {
    std::cerr << "Required arguments: InputFilename" << std::endl;
    return EXIT_FAILURE;
  }
 
  // Get the filename from the command line;
  std::string inputFilename = argv[1];
 
  // Get all data from the file
  vtkSmartPointer<vtkXMLGenericDataObjectReader> reader =
    vtkSmartPointer<vtkXMLGenericDataObjectReader>::New();
  reader->SetFileName(inputFilename.c_str());
  reader->Update();
 
  vtkSmartPointer<vtkCleanUnstructuredGrid> cleanFilter =
    vtkSmartPointer<vtkCleanUnstructuredGrid>::New();
  cleanFilter->SetInputConnection(reader->GetOutputPort());
  cleanFilter->Update();
 
  vtkSmartPointer<vtkCleanUnstructuredGridCells> cleanCellsFilter =
    vtkSmartPointer<vtkCleanUnstructuredGridCells>::New();
  cleanCellsFilter->SetInputConnection(cleanFilter->GetOutputPort());
  cleanCellsFilter->Update();
 
  std::string outputFilename = 
    vtksys::SystemTools::GetFilenameWithoutExtension(inputFilename)
                               + ".vtu";
  vtkSmartPointer<vtkXMLUnstructuredGridWriter> writer =
    vtkSmartPointer<vtkXMLUnstructuredGridWriter>::New();
  writer->SetFileName(outputFilename.c_str());
  writer->SetInputConnection(cleanCellsFilter->GetOutputPort());
  writer->Update();
}