gdcmPython/demo/vtkGdcmReader.py.in

   1 import sys
   2 import os
   3
   4 # All paths must be added without system tests, because of ctest...
   5 sys.path.append('${GDCM_BINARY_DIR}')
   6 if os.name == 'posix':
   7    sys.path.append('${GDCM_BINARY_DIR}/bin')
   8 else:
   9    sys.path.append('${GDCM_BINARY_DIR}/bin/Release')
  10    sys.path.append('${GDCM_BINARY_DIR}/bin/Debug')
  11
  12 from gdcmPython.vtk import *
  13 from vtkpython import *
  14
  15 # Get the file names
  16 try:
  17    fileName = sys.argv[1]
  18 except IndexError:
  19    fileName = os.path.join(GDCM_DATA_PATH, "test.acr")
  20
  21 # read the image
  22 reader = vtkGdcmReader()
  23 reader.SetFileName(fileName)
  24 reader.Update()
  25
  26 # show the image
  27 viewer = vtkImageViewer()
  28 viewer.SetInput(reader.GetOutput())
  29 rng = reader.GetOutput().GetScalarRange()
  30 viewer.SetColorWindow(rng[1] - rng[0])
  31 viewer.SetColorLevel(0.5 * (rng[1] + rng[0]))
  32
  33 dim = reader.GetOutput().GetDimensions()
  34 viewer.SetSize(dim[0],dim[1])
  35 if(dim[2] != 1):
  36    # For multifame dicom, take a snapshot of the center slice (+/- 1)
  37    viewer.SetZSlice(dim[2] / 2)
  38 else:
  39    viewer.SetZSlice(0)
  40 viewer.Render()
  41
  42 # sys.stdin.read(1)
  43