間引き

メッシュを縮小する

from __future__ import annotations

import numpy as np

import pyvista as pv
from pyvista import examples

mesh = examples.download_face()

# Define a camera position that shows this mesh properly
cpos = [(0.4, -0.07, -0.31), (0.05, -0.13, -0.06), (-0.1, 1, 0.08)]
dargs = dict(show_edges=True, color=True)

# Preview the mesh
mesh.plot(cpos=cpos, **dargs)

スタティックシーン

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

ここで，削減目標を定義し， pyvista.PolyDataFilters.decimate() フィルタと pyvista.PolyDataFilters.decimate_pro() フィルタを比較してみましょう．

target_reduction = 0.7
print(f'Reducing {target_reduction * 100.0} percent out of the original mesh')

Reducing 70.0 percent out of the original mesh

decimated = mesh.decimate(target_reduction)

decimated.plot(cpos=cpos, **dargs)

スタティックシーン

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

pro_decimated = mesh.decimate_pro(target_reduction, preserve_topology=True)

pro_decimated.plot(cpos=cpos, **dargs)

スタティックシーン

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

並べて比較:

pl = pv.Plotter(shape=(1, 3))
pl.add_mesh(mesh, **dargs)
pl.add_text('Input mesh', font_size=24)
pl.camera_position = cpos
pl.reset_camera()
pl.subplot(0, 1)
pl.add_mesh(decimated, **dargs)
pl.add_text('Decimated mesh', font_size=24)
pl.camera_position = cpos
pl.reset_camera()
pl.subplot(0, 2)
pl.add_mesh(pro_decimated, **dargs)
pl.add_text('Pro Decimated mesh', font_size=24)
pl.camera_position = cpos
pl.reset_camera()
pl.link_views()
pl.show()

ポリラインメッシュを分解する

かなり遅いスパイラルポリラインメッシュを生成します。

n_points = 100
n_rotations = 5
phi = np.linspace(0, 2 * np.pi * n_rotations, n_points)
ratio = 1.1
r = ratio ** (phi)

points = np.zeros((n_points, 3))
points[:, 0] = r * np.cos(phi)
points[:, 1] = r * np.sin(phi)

spiral = pv.PolyData(points, lines=np.append([n_points], np.arange(n_points)))

将来使用するために、再利用可能なプロット関数を構します。

def compare_decimation(spiral, decimated):
    pl = pv.Plotter()
    pl.add_mesh(spiral, line_width=5, color='r', label='Original')
    pl.add_mesh(decimated, line_width=3, color='k', label='Decimated')
    pl.view_xy()
    pl.add_legend(face='line', size=(0.25, 0.25))

pyvista.PolyDataFilters.decimate_polyline() フィルタを使用して、50%を目標にデシメートします。

decimated = spiral.decimate_polyline(0.5)
print(f'Original # of points:  {spiral.n_points}')
print(f'Decimated # of points: {decimated.n_points}')

Original # of points:  100
Decimated # of points: 50

この程度ならデシメーションは問題なさそうです。

compare_decimation(spiral, decimated)

80％という大きな削減レベルを使用すると、表現がより粗くなります。

decimated = spiral.decimate_polyline(0.8)
print(f'Original # of points:  {spiral.n_points}')
print(f'Decimated # of points: {decimated.n_points}')

Original # of points:  100
Decimated # of points: 20

渦巻きの内側の構造は完全に失われています。

compare_decimation(spiral, decimated)

素早く変化するフィーチャーのエラーを避けるには、 maximum_error パラメータを使用します。 単位はバウンディングボックスの最大長の何分の一かです。 これは、達成される削減のレベルを制限することに注意してください。

decimated = spiral.decimate_polyline(0.8, maximum_error=0.5)
print(f'Original # of points:  {spiral.n_points}')
print(f'Decimated # of points: {decimated.n_points}')

Original # of points:  100
Decimated # of points: 45

The structure of the inner part of the spiral is captured adequately.

compare_decimation(spiral, decimated)

Tags: filter