]> Creatis software - gdcm.git/blob - gdcmPython/demo/vtkGdcmReader.py.in
aea5e5c879fda132c0fb83fa2faa3a4b2e78ef93
[gdcm.git] / gdcmPython / demo / vtkGdcmReader.py.in
1 from gdcmConfigDemo import *
2 from gdcmPython.vtk import *
3 from vtkpython import *
4
5 def PrintUse():
6    print ""
7    print "Use :"
8    print "-----"
9    print "%s <fileName>" % sys.argv[0]
10    print "   fileName : path to the DICOM file to see"
11    print ""
12    print ""
13
14 # Get the file names
15 try:
16    fileName = sys.argv[1]
17 except IndexError:
18    fileName = os.path.join(GDCM_DATA_ROOT, "test.acr")
19    
20 if( not os.path.isfile(fileName) ):
21    PrintUse()
22    sys.exit(1)
23
24 # read the image
25 reader = vtkGdcmReader()
26 reader.SetFileName(fileName)
27 reader.Update()
28
29 # show the image
30 viewer = vtkImageViewer()
31 viewer.SetInput(reader.GetOutput())
32 rng = reader.GetOutput().GetScalarRange()
33 viewer.SetColorWindow(rng[1] - rng[0])
34 viewer.SetColorLevel(0.5 * (rng[1] + rng[0]))
35
36 dim = reader.GetOutput().GetDimensions()
37 viewer.SetSize(dim[0],dim[1])
38 if(dim[2] != 1):
39    # For multifame dicom, take a snapshot of the center slice (+/- 1)
40    viewer.SetZSlice(dim[2] / 2)
41 else:
42    viewer.SetZSlice(0)
43 viewer.Render()
44
45 # sys.stdin.read(1)
46