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