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