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/blow2_ascii.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