Create an orthogonal slice view through a 3D vtkImageData object #python #vtk #imagedata #dicom #visualization

vtkimagedata_view_orthogonal.py

# Assuming that we have an object of type 'vtkImageData' under the name of
# 'image'

# Create a new vtkImageSliceMapper
mapper = vtk.vtkImageSliceMapper()
# Set 'image' as the input dataset
mapper.SetInput(image)
# Set the orthogonal slice orientation to be perpendicular to the 'X' axis. Use
# 'SetOrientationToY' and 'SetOrientationToZ' for the Y and Z axes. 
# Alternatively you can use the 'SetOrientation' method with a 'int' parameter
# of '1', '2', or '3' for the X, Y, or Z axis respectively
mapper.SetOrientationToX()
# Set the slice index to be visualized. Here we set it to the middle slice along
# the X axis 
mapper.SetSliceNumber(image.GetDimensions()[0]//2)

# Create a new 'vtkImageActor'
actor = vtk.vtkImageActor()
# Set the mapper
actor.SetMapper(mapper)
# Set the image opacity to 50%
actor.GetProperty().SetOpacity(0.5)