monman53のぶろぐ

いろいろ載せるよ

vtkフォーマットの雛形

数値計算結果などをさくっとParaViewで見たいときに,vtkのレガシーフォーマットで出力するのは楽で良い.

リファレンス https://www.vtk.org/wp-content/uploads/2015/04/file-formats.pdf

頂点が置ければええんじゃ

立方体の頂点を置いてみる.

$ cat out.vtk
# vtk DataFile Version 3.0
title
ASCII
DATASET UNSTRUCTURED_GRID

POINTS 8 float
0.0 0.0 0.0
0.0 0.0 1.0
0.0 1.0 0.0
0.0 1.0 1.0
1.0 0.0 0.0
1.0 0.0 1.0
1.0 1.0 0.0
1.0 1.0 1.0

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

CELL_TYPES 8
1
1
1
1
1
1
1
1

POINTSで頂点を定義する.このままでは実体が表示されないが,CELL_TYPESで点(VTK_VERTEX=1)として指定することで点として表示されるようだ.

f:id:monman53:20181115162053p:plain
立方体の頂点

頂点と枝が表示されればええんじゃ

枝のぶんCELLSを増やして,CELLS_TYPESで線(VTK_LINE=3)を指定する.適当に枝を追加してみる.

$ cat out.vtk
# vtk DataFile Version 3.0
title
ASCII
DATASET UNSTRUCTURED_GRID

POINTS 8 float
0.0 0.0 0.0
0.0 0.0 1.0
0.0 1.0 0.0
0.0 1.0 1.0
1.0 0.0 0.0
1.0 0.0 1.0
1.0 1.0 0.0
1.0 1.0 1.0

CELLS 12 28
1 0
1 1
1 2
1 3
1 4
1 5
1 6
1 7
2 4 5
2 5 7
2 7 6
2 6 4

CELL_TYPES 12
1
1
1
1
1
1
1
1
3
3
3
3

f:id:monman53:20181115162706p:plain
立方体の1面だけ枝を

面も表示したい!

同じ要領でCELL_TYPESをいろいろ指定することでいろいろできる.やはり冒頭のリファレンスを参照されたい.