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()
8ren1.SetBackground(0,0,0)
9renWin = vtk.vtkRenderWindow()
10renWin.AddRenderer(ren1)
11renWin.SetSize(300,300)
12iren = vtk.vtkRenderWindowInteractor()
13iren.SetRenderWindow(renWin)
14# camera parameters
15camera = ren1.GetActiveCamera()
16camera.SetPosition(-54.8012,109.471,231.412)
17camera.SetFocalPoint(33,33,33)
18camera.SetViewUp(0.157687,0.942832,-0.293604)
19camera.SetViewAngle(30)
20camera.SetClippingRange(124.221,363.827)
21reader = vtk.vtkGenericEnSightReader()
22# Make sure all algorithms use the composite data pipeline
23cdp = vtk.vtkCompositeDataPipeline()
24reader.SetDefaultExecutivePrototype(cdp)
25reader.SetCaseFileName("" + str(VTK_DATA_ROOT) + "/Data/EnSight/ironProt_bin.case")
26Contour0 = vtk.vtkContourFilter()
27Contour0.SetInputConnection(reader.GetOutputPort())
28Contour0.SetValue(0,200)
29Contour0.SetComputeScalars(1)
30mapper = vtk.vtkHierarchicalPolyDataMapper()
31mapper.SetInputConnection(Contour0.GetOutputPort())
32mapper.SetScalarRange(0,1)
33mapper.SetScalarVisibility(1)
34actor = vtk.vtkActor()
35actor.SetMapper(mapper)
36actor.GetProperty().SetRepresentationToSurface()
37actor.GetProperty().SetInterpolationToGouraud()
38ren1.AddActor(actor)
39# enable user interface interactor
40iren.Initialize()
41# prevent the tk window from showing up then start the event loop
42reader.SetDefaultExecutivePrototype(None)
43# --- end of script --
44