1#!/usr/bin/env python 2import vtk 3from vtk.util.misc import vtkGetDataRoot 4VTK_DATA_ROOT = vtkGetDataRoot() 5 6# create a rendering window and renderer 7ren1 = vtk.vtkRenderer() 8renWin = vtk.vtkRenderWindow() 9renWin.AddRenderer(ren1) 10renWin.StereoCapableWindowOn() 11iren = vtk.vtkRenderWindowInteractor() 12iren.SetRenderWindow(renWin) 13reader = vtk.vtkGenericEnSightReader() 14# Make sure all algorithms use the composite data pipeline 15cdp = vtk.vtkCompositeDataPipeline() 16reader.SetDefaultExecutivePrototype(cdp) 17reader.SetCaseFileName("" + str(VTK_DATA_ROOT) + "/Data/EnSight/blow3_bin.case") 18reader.SetTimeValue(1) 19geom = vtk.vtkGeometryFilter() 20geom.SetInputConnection(reader.GetOutputPort()) 21mapper = vtk.vtkHierarchicalPolyDataMapper() 22mapper.SetInputConnection(geom.GetOutputPort()) 23mapper.SetColorModeToMapScalars() 24mapper.SetScalarModeToUsePointFieldData() 25mapper.ColorByArrayComponent("displacement",0) 26mapper.SetScalarRange(0,2.08) 27actor = vtk.vtkActor() 28actor.SetMapper(mapper) 29# assign our actor to the renderer 30ren1.AddActor(actor) 31# enable user interface interactor 32iren.Initialize() 33ren1.GetActiveCamera().SetPosition(99.3932,17.6571,-22.6071) 34ren1.GetActiveCamera().SetFocalPoint(3.5,12,1.5) 35ren1.GetActiveCamera().SetViewAngle(30) 36ren1.GetActiveCamera().SetViewUp(0.239617,-0.01054,0.97081) 37ren1.ResetCameraClippingRange() 38renWin.Render() 39# prevent the tk window from showing up then start the event loop 40reader.SetDefaultExecutivePrototype(None) 41# --- end of script -- 42