if True:
print('With AxisReflection Filter (HTG)')
reflection = vtk.vtkHyperTreeGridAxisReflection()
if isFilter:
reflection.SetInputConnection(htg.GetOutputPort())
else:
reflection.SetInputData(htg)
reflection.SetPlaneToZ()
reflection.SetCenter(0)
isFilter = True
else:
print('No AxisReflection Filter (HTG)')
reflection = htg
# Geometries
geometry = vtk.vtkHyperTreeGridGeometry()
if isFilter:
geometry.SetInputConnection(reflection.GetOutputPort())
else:
geometry.SetInputData(reflection)
print('With Geometry Filter (HTG to NS)')
# Shrink Filter
if True:
print('With Shrink Filter (NS)')
# En 3D, le shrink ne doit pas se faire sur la geometrie car elle ne represente que la peau
shrink = vtk.vtkShrinkFilter()
shrink.SetInputConnection(geometry.GetOutputPort())
shrink.SetShrinkFactor(.8)
else:
print('No Shrink Filter (NS)')
if True:
print('With AxisReflection Filter (HTG)')
reflection = vtk.vtkHyperTreeGridAxisReflection()
if isFilter:
reflection.SetInputConnection(htg.GetOutputPort())
else:
reflection.SetInputData(htg)
reflection.SetPlaneToY()
reflection.SetCenter(0)
isFilter = True
else:
print('No AxisReflection Filter (HTG)')
reflection = htg
# Geometries
geometry = vtk.vtkHyperTreeGridGeometry()
if isFilter:
geometry.SetInputConnection(reflection.GetOutputPort())
else:
geometry.SetInputData(reflection)
print('With Geometry Filter (HTG to NS)')
# Shrink Filter
if True:
print('With Shrink Filter (NS)')
# En 3D, le shrink ne doit pas se faire sur la geometrie car elle ne represente que la peau
shrink = vtk.vtkShrinkFilter()
shrink.SetInputConnection(geometry.GetOutputPort())
shrink.SetShrinkFactor(.8)
else:
print('No Shrink Filter (NS)')
#!/usr/bin/env python
import vtk
from vtk.util.misc import vtkGetDataRoot
VTK_DATA_ROOT = vtkGetDataRoot()
pioreader = vtk.vtkPIOReader()
pioreader.SetFileName("" + str(VTK_DATA_ROOT) + "/Data/PIO/simple.pio")
pioreader.SetCurrentTimeStep(1)
pioreader.Update()
grid = pioreader.GetOutput()
block = grid.GetBlock(0)
geometryFilter = vtk.vtkHyperTreeGridGeometry()
geometryFilter.SetInputData(block)
geometryFilter.Update()
# ---------------------------------------------------------------------
# Rendering
# ---------------------------------------------------------------------
ren = vtk.vtkRenderer()
renWin = vtk.vtkRenderWindow()
renWin.AddRenderer(ren)
iren = vtk.vtkRenderWindowInteractor()
iren.SetRenderWindow(renWin)
lut = vtk.vtkLookupTable()
lut.SetHueRange(0.66, 0)
lut.SetSaturationRange(1.0, 0.25)
