package require vtk

# This example demonstrates how to permute an image with vtkImageReslice.
# The advantage of using vtkImageReslice for this (rather than
# vtkImagePermute) is that you can specify negative or even oblique axes,
# and you can resample the image at the same time.

vtkImageReader reader
reader ReleaseDataFlagOff
reader SetDataByteOrderToLittleEndian
reader SetDataExtent 0 63 0 63 1 93
reader SetDataSpacing 3.2 3.2 1.5
reader SetFilePrefix "$VTK_DATA_ROOT/Data/headsq/quarter"
reader SetDataMask 0x7fff

vtkImageReslice reslice
reslice SetInputConnection [reader GetOutputPort]
# specify the new axes in terms of the original axes
reslice SetResliceAxesDirectionCosines  0 +1 0   0 0 -1   -1 0 0    
# resample the image to enlarge & to get the aspect right
reslice SetOutputSpacing 1.0 1.0 1.0

vtkImageViewer viewer
viewer SetInputConnection [reslice GetOutputPort]
viewer SetZSlice 100
viewer SetColorWindow 2000
viewer SetColorLevel 1000
viewer Render









