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)