並列フィルタ

VTKライブラリーは並列ファイル・フォーマットをサポートしています．複数のファイルに分割されたメッシュの読み込みは，VTKおよびPyVistaでネイティブにサポートされています．

import os

import pyvista as pv
from pyvista import examples

いくつかのファイルに分割された pyvista.UnstructuredGrid を含むサンプルデータセットをダウンロードしてみましょう．

load=False を設定し，ダウンロードしたファイルがあるディレクトリを見て，データセットがどこにダウンロードされたかを調べてみましょう．

filename = examples.download_blood_vessels(load=False)
path = os.path.dirname(filename)
os.listdir(path)

['T0000000500.pvtu', 'T0000000500']

os.listdir(os.path.join(path, "T0000000500"))

['002.vtu', '000.vtu', '001.vtu', '003.vtu']

.pvtu ファイルはディレクトリの横にあるので注意してください．このディレクトリには，メッシュ全体を構成するすべての並列ファイルまたはピースが含まれています． .pvtu ファイルを読み取るだけで，VTKがメッシュをまとめます．

mesh = pv.read(filename)
mesh

Header	Data Arrays
UnstructuredGrid Information N Cells 39353 N Points 48823 X Bounds 5.300e+01, 1.210e+02 Y Bounds 5.000e+01, 9.700e+01 Z Bounds 6.400e+01, 1.820e+02 N Arrays 5	Name Field Type N Comp Min Max node_value Points int32 1 0.000e+00 3.000e+00 simerr_type Points int32 1 0.000e+00 3.000e+00 density Cells float32 1 2.203e-01 5.232e-01 velocity Cells float32 3 -3.607e-01 8.989e-02 shearstress Cells float32 1 6.160e-05 1.726e-02

つなぎ合わせたメッシュをプロットします

mesh.plot(scalars="node_value", categories=True)

mesh.plot(scalars="density")