並列フィルタ

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

from __future__ import annotations

from pathlib import Path

import pyvista as pv
from pyvista import examples

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

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

filename = examples.download_blood_vessels(load=False)
path = Path(filename).parent
list(path.iterdir())

[PosixPath('/home/runner/.cache/pyvista_3/pvtu_blood_vessels/blood_vessels.zip.unzip/blood_vessels/T0000000500.pvtu'), PosixPath('/home/runner/.cache/pyvista_3/pvtu_blood_vessels/blood_vessels.zip.unzip/blood_vessels/T0000000500')]

list((path / 'T0000000500').iterdir())

[PosixPath('/home/runner/.cache/pyvista_3/pvtu_blood_vessels/blood_vessels.zip.unzip/blood_vessels/T0000000500/002.vtu'), PosixPath('/home/runner/.cache/pyvista_3/pvtu_blood_vessels/blood_vessels.zip.unzip/blood_vessels/T0000000500/001.vtu'), PosixPath('/home/runner/.cache/pyvista_3/pvtu_blood_vessels/blood_vessels.zip.unzip/blood_vessels/T0000000500/000.vtu'), PosixPath('/home/runner/.cache/pyvista_3/pvtu_blood_vessels/blood_vessels.zip.unzip/blood_vessels/T0000000500/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