1#!/usr/bin/env python 2import vtk 3from vtk.util.misc import vtkGetDataRoot 4VTK_DATA_ROOT = vtkGetDataRoot() 5 6# Image pipeline 7reader = vtk.vtkImageReader() 8reader.GetExecutive().SetReleaseDataFlag(0,0) 9reader.SetDataByteOrderToLittleEndian() 10reader.SetDataExtent(0,63,0,63,1,93) 11reader.SetFilePrefix("" + str(VTK_DATA_ROOT) + "/Data/headsq/quarter") 12reader.SetDataMask(0x7fff) 13imageFloat = vtk.vtkImageCast() 14imageFloat.SetInputConnection(reader.GetOutputPort()) 15imageFloat.SetOutputScalarTypeToFloat() 16flipX = vtk.vtkImageFlip() 17flipX.SetInputConnection(imageFloat.GetOutputPort()) 18flipX.SetFilteredAxis(0) 19flipY = vtk.vtkImageFlip() 20flipY.SetInputConnection(imageFloat.GetOutputPort()) 21flipY.SetFilteredAxis(1) 22flipY.FlipAboutOriginOn() 23imageAppend = vtk.vtkImageAppend() 24imageAppend.AddInputConnection(imageFloat.GetOutputPort()) 25imageAppend.AddInputConnection(flipX.GetOutputPort()) 26imageAppend.AddInputConnection(flipY.GetOutputPort()) 27imageAppend.SetAppendAxis(0) 28viewer = vtk.vtkImageViewer() 29viewer.SetInputConnection(imageAppend.GetOutputPort()) 30viewer.SetZSlice(22) 31viewer.SetColorWindow(2000) 32viewer.SetColorLevel(1000) 33#make interface 34viewer.Render() 35# --- end of script -- 36