並列フィルタ

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

from __future__ import annotations

import os
from pathlib import Path

import pyvista as pv
from pyvista import examples

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

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

filename = examples.download_blood_vessels(load=False)
path = str(Path(filename).parent)
os.listdir(path)

['T0000000500.pvtu', 'T0000000500']

os.listdir(str(Path(path) / 'T0000000500'))

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

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

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')

スタティックシーン

インタラクティブなシーン

Tags: load